Lines Matching refs:TS
136 Source TS;
519 extern Source TS;
548 extern Source TS;
556 timeOffsetFlag=findTimeOffset(ss,&TS.timeOffset);
572 TS.ZAxis+=(RPts[3+j]-RPts[j])*(RPts[3+j]-RPts[j]);
573 TS.XAxis+=(RPts[6+j]-RPts[3+j])*(RPts[6+j]-RPts[3+j]);
577 TS.ZAxis=sqrt(TS.ZAxis)/100.0; // convert to metres from cm
578 TS.XAxis=sqrt(TS.XAxis)/100.0;
581 TS.ZAxis=yh;
582 TS.XAxis=xw;
584 fprintf(stderr,"Time off sec == %g \n",TS.timeOffset);
585 fprintf(stderr,"mod size z == %g \n",TS.ZAxis);
586 // TS.rdumMid=calcRDum(RPts);
587 TS.rdumMid=TS.timeOffset; // Goran
615 extern Source TS;
631 TS.nTime=0;
632 TS.nEnergy=0;
647 TS.nTime=tIndex;
663 if (TS.nTime==0)
676 TS.nEnergy=eIndex+2;
677 if (!TS.nTime && tIndex) // SMALL TEST if reading first energy bin
678 TS.nTime=tIndex;
679 // printf("tIndex %d %d %d %d \n",tIndex,eIndex,TS.nEnergy,TS.nTime);
684 TS.Flux=matrix(TS.nEnergy,TS.nTime);
685 TS.EInt=(double*) malloc(TS.nEnergy*sizeof(double)); // Runing integral of energy
686 TS.TimeBin=(double*) malloc(TS.nTime*sizeof(double));
687 TS.EnergyBin=(double*) malloc(TS.nEnergy*sizeof(double));
688 if (verbose) fprintf(stderr,"NUMBER %d %d \n",TS.nEnergy,TS.nTime);
696 TS.EInt[0]=0.0;
708 TS.EnergyBin[eIndex]=(Einit>Ea) ? Einit : Ea;
726 TS.EInt[eIndex+1]=Tsum;
728 …if (verbose) fprintf(stderr,"Energy %g %g \n",EtoLambda(Ea),(TS.EInt[eIndex+1]-TS.EInt[eIndex])*3.…
736 TS.TimeBin[tIndex]=T/1e8; // convert Time into second (from shakes)
738 TS.Flux[eIndex][tIndex]=Tsum;
745 TS.EnergyBin[eIndex+1]=Eend;
746 TS.Total=Tsum;
748 // printf("tIndex %d %d %d \n",tIndex,eIndex,TS.nTime);
750 // fprintf(stderr,"ebin1 ebinN %g %g\n",TS.EnergyBin[0],TS.EnergyBin[TS.nEnergy-1]);
760 Uses TS table to find the point
769 extern Source TS;
786 Rend=R=TS.Total*R0;
788 Epnt=binSearch(TS.nEnergy-1,TS.EInt,R);
789 Tpnt=binSearch(TS.nTime-1,TS.Flux[Epnt-1],R);
791 if(Epnt<0 || R<TS.Flux[Epnt-1][Tpnt-1] || R >TS.Flux[Epnt-1][Tpnt] )
794 TS.Flux[Epnt-1][Tpnt-1],R,TS.Flux[Epnt-1][Tpnt]);
800 Espread=TS.EInt[0];
801 *EV=TS.EnergyBin[1];
805 Estart=TS.EInt[Epnt-1];
806 Espread=TS.EInt[Epnt]-TS.EInt[Epnt-1];
807 *EV=(Epnt>TS.nEnergy-1) ? TS.EnergyBin[Epnt+1] : TS.EnergyBin[Epnt];
810 if (Tpnt==0 || Epnt==0 || Tpnt==TS.nTime)
815 *TV=TS.TimeBin[Tpnt-1];
816 TRange=TS.TimeBin[Tpnt]-TS.TimeBin[Tpnt-1];
817 Tspread=TS.Flux[Epnt-1][Tpnt]-TS.Flux[Epnt-1][Tpnt-1];
818 R-=TS.Flux[Epnt-1][Tpnt-1];
824 R1=TS.EInt[Epnt-1]+Espread*rand01();
826 iEnd=TS.nEnergy>Epnt+highPoly ? Epnt+highPoly : TS.nEnergy-1; // min(nEnergy-1,Epnt+highPoly
828 *EV=polInterp(TS.EInt+iStart,TS.EnergyBin+iStart,1+iEnd-iStart,R1);
830 if(*TV < TS.TimeBin[Tpnt-1] || *TV > TS.TimeBin[Tpnt])
832 fprintf(stderr,"%d Tpnt %d Tval %g Epnt %d \n",TS.nTime,Tpnt,*TV,Epnt);
833 … fprintf(stderr,"TBoundary == %12.6e,%g , %12.6e \n\n",TS.TimeBin[Tpnt-1],*TV,TS.TimeBin[Tpnt]);
840 Epnt,Tpnt,TS.EnergyBin[Epnt-1],*EV,TS.EnergyBin[Epnt]);
842 fprintf(stderr,"TS == %g %g :: %d %d \n",TS.EInt[Epnt-1],TS.EInt[Epnt],iStart,iEnd);
845 fprintf(stderr,"Time start %g %g",*(TS.EInt+i+iStart),*(TS.EnergyBin+iStart+i));
989 Control Values TS.XAxis TS.ZAxis focus_xw focus_yh
1003 fprintf(stderr,"TS axis == %g %g\n",TS.XAxis,TS.ZAxis);
1007 Mx=(i+0.5)*TS.XAxis/(NStep*2.0);
1010 My=(j+0.5)*TS.ZAxis/(NStep*2.0);
1022 A*= (TS.XAxis*TS.ZAxis)/(NStep*NStep*(2.0*NStep+1.0)*(2.0*NStep+1.0));
1030 // A*=TS.XAxis*TS.ZAxis/(TS.XAxis*TS.ZAxis-projArea);
1033 …stderr,"Viewport == %g %g Moderator size == (%g * %g) m^2 \n",focus_xw,focus_yh,TS.XAxis,TS.ZAxis);
1064 fprintf(stderr,"Number of Energy Points == %d\n",TS.nEnergy);
1067 dist=TS.rdumMid;
1074 fprintf(stderr,"Totals:: %g %d %d \n",TS.Total,TS.nEnergy,TS.nTime);
1086 x += TS.XAxis*(0.5-rand01()); /* Get point on shutter enterance */
1087 y += TS.ZAxis*(0.5-rand01()); /* Get point on shutter enterance */
1089 // if (modPosition) z=TS.rdumMid; // Goran
1096 // Ival=TS.Total*6.2415093e+12; // * 1.1879451; /* ( of proton in 1uAmp) * (1-cos(30))*2*Pi */
1097 // Ival=TS.Total*5.9294338e+12; // Number of proton in 1uAmp * 0.95 (0.95 because of muon target l…
1098 Ival=TS.Total*6.2415093e+12; // Number of proton in 1uAmp
1109 t=Tval-(TS.rdumMid-TS.timeOffset)/vz; ///vz;
1113 t+=TS.rdumMid/vz;
1119 // if (!modPosition) p*=TS.timeOffset;
1126 fprintf(stderr,"PlaneAxis %g %g \n",TS.XAxis,fullAngle);
1127 fprintf(stderr,"RMID %g \n",TS.rdumMid);
1128 fprintf(stderr,"TimeMid[%d]=> %g\n",Ncount,TS.rdumMid);
1129 fprintf(stderr,"TimeOffset[%d]=> %g\n",Ncount,TS.timeOffset);
1144 multiline(5,-0.5*TS.XAxis,-0.5*TS.ZAxis,0.0,
1145 0.5*TS.XAxis,-0.5*TS.ZAxis,0.0,
1146 0.5*TS.XAxis,0.5*TS.ZAxis,0.0,
1147 -0.5*TS.XAxis,0.5*TS.ZAxis,0.0,
1148 -0.5*TS.XAxis,-0.5*TS.ZAxis,0.0);