- // Time to TDC signal
- // 256 channels for timediff, range 1ns
- iTimeAv=(Int_t)timeAv;
- timeDiff= 512+1000*timeDiff/channelWidth; // time channel numbres
- iTimeDiff=(Int_t)timeDiff;
- // fill digits
- fdigits->SetTimeBestLeft(bestLeftADC);
- fdigits->SetTimeBestRight(bestRightADC);
- fdigits->SetMeanTime(iTimeAv);
- fdigits->SetTimeDiff(iTimeDiff);
- for (Int_t i=0; i<12; i++)
- {
- // fill TDC
- timeright[i+1]=gRandom->Gaus(timeright[i+1],0.05);
- timeleft[i+1]=gRandom->Gaus(timeleft[i+1],0.05);
- tr= Int_t (timeright[i+1]*1000/channelWidth);
- if(tr<200) tr=0;
- tl= Int_t (timeleft[i+1]*1000/channelWidth);
- if(tl<1000) tl=0;
-
- ftimeRightTDC->AddAt(tr,i);
- ftimeLeftTDC->AddAt(tl,i);
- //fill ADC
- Int_t al=( Int_t ) countEl[i+1]/ channelWidthADC;
- Int_t ar=( Int_t ) countEr[i+1]/ channelWidthADC;
- fRightADC->AddAt(ar,i);
- fLeftADC ->AddAt(al,i);
- sumRight+=countEr[i+1];
- }
- fdigits->SetTimeRight(*ftimeRightTDC);
- fdigits->SetTimeLeft(*ftimeLeftTDC);
- fdigits->SetADCRight(*fRightADC);
- fdigits->SetADCLeft(*fLeftADC);
- fdigits->SetSumADCRight(sumRight);
- }
- else
- {timeAv=999999; timeDiff=99999;}
-
-// trick to find out output dir:
+ //fill ADC
+ Int_t al= countE[i];
+ // QTC procedure:
+ // phe -> mV 0.3; 1MIP ->500phe -> ln (amp (mV)) = 5;
+ // max 200ns, HIJING mean 50000phe -> 15000mv -> ln = 15 (s zapasom)
+ // channel 25ps
+ if (al>threshold) {
+ qt=Int_t (TMath::Log(al*ph2mV) * mV2channel);
+ qtAmp=Int_t (TMath::Log(al*10*ph2mV) * mV2channel);
+ fADC->AddAt(qt,i);
+ ftimeTDC->AddAt(tr,i);
+ fADCAmp->AddAt(qtAmp,i);
+ ftimeTDCAmp->AddAt(tr,i); //now is the same as non-amplified but will be change
+ }
+ } //pmt loop
+ for (Int_t im=0; im<3; im++)
+ {
+ if (sumMult[im]>threshold){
+ Int_t sum=Int_t (TMath::Log(sumMult[im]*ph2mV) * mV2channel);
+ Int_t sumAmp=Int_t (TMath::Log(10*sumMult[im]*ph2mV) * mV2channel);
+ fSumMult->AddAt(sum,im);
+ fSumMult->AddAt(sumAmp,im+3);
+ }
+ }