2 /**************************************************************************
3 * Copyright(c) 1998-2000, ALICE Experiment at CERN, All rights reserved. *
5 * Author: The ALICE Off-line Project. *
6 * Contributors are mentioned in the code where appropriate. *
8 * Permission to use, copy, modify and distribute this software and its *
9 * documentation strictly for non-commercial purposes is hereby granted *
10 * without fee, provided that the above copyright notice appears in all *
11 * copies and that both the copyright notice and this permission notice *
12 * appear in the supporting documentation. The authors make no claims *
13 * about the suitability of this software for any purpose. It is *
14 * provided "as is" without express or implied warranty. *
15 **************************************************************************/
20 #include <TDirectory.h>
28 #include "AliSTARTDigitizer.h"
30 #include "AliSTARThit.h"
31 #include "AliSTARTdigit.h"
32 #include "AliRunDigitizer.h"
33 #include <AliDetector.h>
35 #include <AliLoader.h>
36 #include <AliRunLoader.h>
38 #include <Riostream.h>
39 #include <Riostream.h>
40 #include "AliSTARTParameters.h"
41 #include "AliCDBLocal.h"
42 #include "AliCDBStorage.h"
43 #include "AliCDBManager.h"
44 #include "AliCDBEntry.h"
46 ClassImp(AliSTARTDigitizer)
48 //___________________________________________
49 AliSTARTDigitizer::AliSTARTDigitizer() :AliDigitizer()
51 // Default ctor - don't use it
55 //___________________________________________
56 AliSTARTDigitizer::AliSTARTDigitizer(AliRunDigitizer* manager)
57 :AliDigitizer(manager),
66 // ctor which should be used
68 AliDebug(1,"processed");
74 ftimeCFD = new TArrayI(24);
75 fADC = new TArrayI(24);
76 ftimeLED = new TArrayI(24);
77 fADC0 = new TArrayI(24);
82 //------------------------------------------------------------------------
83 AliSTARTDigitizer::~AliSTARTDigitizer()
95 //------------------------------------------------------------------------
96 Bool_t AliSTARTDigitizer::Init()
104 //---------------------------------------------------------------------
106 void AliSTARTDigitizer::Exec(Option_t* /*option*/)
110 Produde digits from hits
111 digits is TObject and includes
112 We are writing array if left & right TDC
113 left & right ADC (will need for slow simulation)
114 TOF first particle left & right
115 mean time and time difference (vertex position)
120 AliRunLoader *outRL = AliRunLoader::GetRunLoader(fManager->GetOutputFolderName());
121 AliLoader * pOutStartLoader = outRL->GetLoader("STARTLoader");
123 AliDebug(1,"start...");
126 // From hits to digits
130 Int_t volume, pmt, trCFD, trLED;
132 Int_t bestRightTDC, bestLeftTDC;
133 Float_t time[24], besttime[24], timeGaus[24] ;
134 // Float_t channelWidth=25.; //ps
135 //Q->T-> coefficients !!!! should be asked!!!
136 Float_t gain[24],timeDelayCFD[24], timeDelayLED[24];
137 Int_t threshold =50; //photoelectrons
138 Int_t qtThreshold=13, qtCh;
139 Float_t zdetA, zdetC;
140 TObjArray slewingLED;
141 AliSTARTParameters* param = AliSTARTParameters::Instance();
142 Int_t ph2Mip = param->GetPh2Mip();
143 Int_t channelWidth = param->GetChannelWidth() ;
146 for (Int_t i=0; i<24; i++){
147 timeDelayCFD[i] = param->GetTimeDelayCFD(i);
148 timeDelayLED[i] = param->GetTimeDelayLED(i);
149 gain[i] = param->GetGain(i);
150 TGraph* gr = param ->GetSlew(i);
151 slewingLED.AddAtAndExpand(gr,i);
152 TGraph* grEff = param ->GetPMTeff(i);
153 fEffPMT.AddAtAndExpand(grEff,i);
155 zdetC = param->GetZposition(0);
156 zdetA = param->GetZposition(1);
158 AliSTARThit *startHit;
161 Int_t nFiles=fManager->GetNinputs();
162 for (Int_t inputFile=0; inputFile<nFiles; inputFile++) {
163 if (inputFile < nFiles-1) {
164 AliWarning(Form("ignoring input stream %d", inputFile));
169 Float_t besttimeright=99999.;
170 Float_t besttimeleft=99999.;
171 Int_t pmtBestRight=9999;
172 Int_t pmtBestLeft=9999;
173 Int_t timeDiff=99999, meanTime=0;
179 cout<<" before reading event "<<endl;
180 for (Int_t i0=0; i0<24; i0++)
182 time[i0]=besttime[i0]=timeGaus[i0]=99999; countE[i0]=0;
184 AliRunLoader * inRL = AliRunLoader::GetRunLoader(fManager->GetInputFolderName(inputFile));
185 AliLoader * pInStartLoader = inRL->GetLoader("STARTLoader");
186 if (!inRL->GetAliRun()) inRL->LoadgAlice();
187 AliSTART *fSTART = (AliSTART*)inRL ->GetAliRun()->GetDetector("START");
190 pInStartLoader->LoadHits("READ");//probably it is necessary to load them before
191 TClonesArray *fHits = fSTART->Hits ();
192 TTree *th = pInStartLoader->TreeH();
193 brHits = th->GetBranch("START");
195 fSTART->SetHitsAddressBranch(brHits);
197 AliError("Branch START hit not found");
200 Int_t ntracks = (Int_t) th->GetEntries();
202 if (ntracks<=0) return;
203 // Start loop on tracks in the hits containers
204 for (Int_t track=0; track<ntracks;track++) {
205 brHits->GetEntry(track);
206 nhits = fHits->GetEntriesFast();
207 for (hit=0;hit<nhits;hit++)
209 startHit = (AliSTARThit*) fHits->UncheckedAt(hit);
211 AliError("The unchecked hit doesn't exist");
216 Double_t e=startHit->Etot();
217 volume = startHit->Volume();
219 if(e>0 && RegisterPhotoE(numpmt,e)) {
221 besttime[numpmt] = startHit->Time();
222 if(besttime[numpmt]<time[numpmt])
224 time[numpmt]=besttime[numpmt];
226 } //photoelectron accept
229 cout<<" hits was read "<<endl;
231 //spread time right&left by 25ps && besttime
232 Float_t c = 0.0299792; // cm/ps
234 Float_t koef=(zdetA-zdetC)/c; //correction position difference by cable
235 for (Int_t ipmt=0; ipmt<12; ipmt++){
236 if(countE[ipmt] > threshold) {
237 timeGaus[ipmt]=gRandom->Gaus(time[ipmt],25)+koef;
238 if(timeGaus[ipmt]<besttimeleft){
239 besttimeleft=timeGaus[ipmt]; //timeleft
243 cout<<"left "<<besttimeleft<<" "<<pmtBestLeft<<" koef "<<koef<<endl;
244 for ( Int_t ipmt=12; ipmt<24; ipmt++){
245 if(countE[ipmt] > threshold) {
246 timeGaus[ipmt]=gRandom->Gaus(time[ipmt],25);
247 if(timeGaus[ipmt]<besttimeright) {
248 besttimeright=timeGaus[ipmt]; //timeright
252 cout<<"right"<<besttimeright<<" "<<pmtBestRight<<endl;
254 //folding with alignmentz position distribution
256 // Float_t koef=(zdetA-zdetC)/c;
257 // Float_t besttimeleftR= besttimeleft;
258 // besttimeleft=koef+besttimeleft;
259 cout<<" besttimeleft "<< besttimeleft<<" delay "<< timeDelayCFD[pmtBestLeft]<<" channelWidth "<<channelWidth<<endl;
260 bestLeftTDC=Int_t ((besttimeleft+1000*timeDelayCFD[pmtBestLeft])
262 cout<<" bestLeftTDC "<<bestLeftTDC<<" delay "<<timeDelayCFD[pmtBestLeft]<<endl;
263 bestRightTDC=Int_t ((besttimeright+1000*timeDelayCFD[pmtBestLeft])
265 timeDiff=Int_t (((besttimeleft+1000*timeDelayCFD[pmtBestLeft])-
266 (besttimeright+1000*timeDelayCFD[pmtBestLeft])
268 meanTime=Int_t (((besttimeright+1000*timeDelayCFD[pmtBestLeft]+
269 besttimeleft+1000*timeDelayCFD[pmtBestLeft])/2.)
271 cout<<"bestLeftTDC"<<bestLeftTDC<<" bestRightTDC "<< bestRightTDC<<
272 " timeDiff "<<timeDiff<<" meanTime= "<<meanTime<<endl;
273 AliDebug(2,Form(" time in ns %f ", besttimeleft));
274 for (Int_t i=0; i<24; i++)
276 Float_t al = countE[i];
277 if (al>threshold && timeGaus[i]<50000 ) {
280 // phe -> mV 0.3; 1MIP ->500phe -> ln (amp (mV)) = 5;
281 // max 200ns, HIJING mean 50000phe -> 15000mv -> ln = 15 (s zapasom)
283 qt= 50.*al*gain[i]/ph2Mip; // 50mv/Mip amp in mV
285 trCFD = Int_t (timeGaus[i] + 1000.*timeDelayCFD[i])/channelWidth;
286 trLED= Int_t (timeGaus[i] + 1000.*timeDelayLED[i]);
287 sl = ((TGraph*)slewingLED.At(i))->Eval(qt);
288 trLED = Int_t(( trLED + 1000*sl )/channelWidth);
289 qtCh=Int_t (1000.*TMath::Log(qt)) / channelWidth;
292 ftimeCFD->AddAt(Int_t (trCFD),i);
293 ftimeLED->AddAt(trLED,i);
294 sumMult += Int_t ((al*gain[i]/ph2Mip)*50) ;
295 cout<<al<<" "<<qt<<" trCFD "<<trCFD<<" sl "<<sl<<" trLED "<<trLED<<
296 " qtThreshold "<<qtThreshold<<" qtCh "<< qtCh<<endl;
300 if (sumMult > threshold){
301 fSumMult = Int_t (TMath::Log(sumMult));
303 fSTART->AddDigit(bestRightTDC,bestLeftTDC,meanTime,timeDiff,fSumMult,
304 ftimeCFD,fADC,ftimeLED,fADC0);
305 pOutStartLoader->UnloadHits();
306 } //input streams loop
309 pOutStartLoader->LoadDigits("UPDATE");
310 TTree *treeD = pOutStartLoader->TreeD();
312 pOutStartLoader->MakeTree("D");
313 treeD = pOutStartLoader->TreeD();
316 fSTART = (AliSTART*)outRL ->GetAliRun()->GetDetector("START");
317 // Make a branch in the tree
318 fSTART->MakeBranch("D");
321 cout<<" treeD->Fill(); "<<endl;
323 pOutStartLoader->WriteDigits("OVERWRITE");
325 cout<<" pOutStartLoader->WriteDigits "<<endl;
326 fSTART->ResetDigits();
327 pOutStartLoader->UnloadDigits();
332 //------------------------------------------------------------------------
333 Bool_t AliSTARTDigitizer::RegisterPhotoE(Int_t ipmt,Double_t energy)
337 // Float_t hc=197.326960*1.e6; //mev*nm
338 Double_t hc=1.973*1.e-6; //gev*nm
339 Float_t lambda=hc/energy;
340 Float_t eff = ((TGraph*) fEffPMT.At(ipmt))->Eval(lambda);
341 Double_t p = gRandom->Rndm();
349 //----------------------------------------------------------------------------