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, qtCh;
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 Float_t zdetA, zdetC;
139 TObjArray slewingLED;
140 AliSTARTParameters* param = AliSTARTParameters::Instance();
141 Int_t ph2Mip = param->GetPh2Mip();
142 Int_t channelWidth = param->GetChannelWidth() ;
145 for (Int_t i=0; i<24; i++){
146 timeDelayCFD[i] = param->GetTimeDelayCFD(i);
147 timeDelayLED[i] = param->GetTimeDelayLED(i);
148 gain[i] = param->GetGain(i);
149 TGraph* gr = param ->GetSlew(i);
150 slewingLED.AddAtAndExpand(gr,i);
151 TGraph* grEff = param ->GetPMTeff(i);
152 fEffPMT.AddAtAndExpand(grEff,i);
154 zdetC = param->GetZposition(0);
155 zdetA = param->GetZposition(1);
157 AliSTARThit *startHit;
160 Int_t nFiles=fManager->GetNinputs();
161 for (Int_t inputFile=0; inputFile<nFiles; inputFile++) {
162 if (inputFile < nFiles-1) {
163 AliWarning(Form("ignoring input stream %d", inputFile));
168 Float_t besttimeright=99999.;
169 Float_t besttimeleft=99999.;
170 Int_t pmtBestRight=9999;
171 Int_t pmtBestLeft=9999;
172 Int_t timeDiff=99999, meanTime=0;
178 for (Int_t i0=0; i0<24; i0++)
180 time[i0]=besttime[i0]=timeGaus[i0]=99999; countE[i0]=0;
182 AliRunLoader * inRL = AliRunLoader::GetRunLoader(fManager->GetInputFolderName(inputFile));
183 AliLoader * pInStartLoader = inRL->GetLoader("STARTLoader");
184 if (!inRL->GetAliRun()) inRL->LoadgAlice();
185 AliSTART *fSTART = (AliSTART*)inRL ->GetAliRun()->GetDetector("START");
188 pInStartLoader->LoadHits("READ");//probably it is necessary to load them before
189 TClonesArray *fHits = fSTART->Hits ();
190 TTree *th = pInStartLoader->TreeH();
191 brHits = th->GetBranch("START");
193 fSTART->SetHitsAddressBranch(brHits);
195 AliError("Branch START hit not found");
198 Int_t ntracks = (Int_t) th->GetEntries();
200 if (ntracks<=0) return;
201 // Start loop on tracks in the hits containers
202 for (Int_t track=0; track<ntracks;track++) {
203 brHits->GetEntry(track);
204 nhits = fHits->GetEntriesFast();
205 for (hit=0;hit<nhits;hit++)
207 startHit = (AliSTARThit*) fHits->UncheckedAt(hit);
209 AliError("The unchecked hit doesn't exist");
214 Double_t e=startHit->Etot();
215 volume = startHit->Volume();
217 if(e>0 && RegisterPhotoE(numpmt,e)) {
219 besttime[numpmt] = startHit->Time();
220 if(besttime[numpmt]<time[numpmt])
222 time[numpmt]=besttime[numpmt];
224 } //photoelectron accept
228 //spread time right&left by 25ps && besttime
229 Float_t c = 0.0299792; // cm/ps
231 Float_t koef=(zdetA-zdetC)/c; //correction position difference by cable
232 for (Int_t ipmt=0; ipmt<12; ipmt++){
233 if(countE[ipmt] > threshold) {
234 timeGaus[ipmt]=gRandom->Gaus(time[ipmt],25)+koef;
235 if(timeGaus[ipmt]<besttimeleft){
236 besttimeleft=timeGaus[ipmt]; //timeleft
240 for ( Int_t ipmt=12; ipmt<24; ipmt++){
241 if(countE[ipmt] > threshold) {
242 timeGaus[ipmt]=gRandom->Gaus(time[ipmt],25);
243 if(timeGaus[ipmt]<besttimeright) {
244 besttimeright=timeGaus[ipmt]; //timeright
249 //folding with alignmentz position distribution
251 bestLeftTDC=Int_t ((besttimeleft+1000*timeDelayCFD[pmtBestLeft])
254 bestRightTDC=Int_t ((besttimeright+1000*timeDelayCFD[pmtBestLeft])
256 timeDiff=Int_t (((besttimeleft-besttimeright)+1000*timeDelayCFD[pmtBestRight])
258 meanTime=Int_t (((besttimeright+1000*timeDelayCFD[pmtBestLeft]+
259 besttimeleft+1000*timeDelayCFD[pmtBestLeft])/2.)
262 AliDebug(2,Form(" time in ns %f %f ",besttimeleft,besttimeright ));
263 for (Int_t i=0; i<24; i++)
265 Float_t al = countE[i];
266 if (al>threshold && timeGaus[i]<50000 ) {
269 // phe -> mV 0.3; 1MIP ->500phe -> ln (amp (mV)) = 5;
270 // max 200ns, HIJING mean 50000phe -> 15000mv -> ln = 15 (s zapasom)
272 qt= 50.*al*gain[i]/ph2Mip; // 50mv/Mip amp in mV
274 trCFD = Int_t (timeGaus[i] + 1000.*timeDelayCFD[i])/channelWidth;
275 trLED= Int_t (timeGaus[i] + 1000.*timeDelayLED[i]);
276 sl = ((TGraph*)slewingLED.At(i))->Eval(qt);
277 trLED = Int_t(( trLED + 1000*sl )/channelWidth);
278 qtCh=Int_t (1000.*TMath::Log(qt)) / channelWidth;
281 ftimeCFD->AddAt(Int_t (trCFD),i);
282 ftimeLED->AddAt(trLED,i);
283 sumMult += Int_t ((al*gain[i]/ph2Mip)*50) ;
287 if (sumMult > threshold){
288 fSumMult = Int_t (TMath::Log(sumMult));
290 fSTART->AddDigit(bestRightTDC,bestLeftTDC,meanTime,timeDiff,fSumMult,
291 ftimeCFD,fADC,ftimeLED,fADC0);
292 pOutStartLoader->UnloadHits();
293 } //input streams loop
296 pOutStartLoader->LoadDigits("UPDATE");
297 TTree *treeD = pOutStartLoader->TreeD();
299 pOutStartLoader->MakeTree("D");
300 treeD = pOutStartLoader->TreeD();
303 fSTART = (AliSTART*)outRL ->GetAliRun()->GetDetector("START");
304 // Make a branch in the tree
305 fSTART->MakeBranch("D");
308 pOutStartLoader->WriteDigits("OVERWRITE");
310 fSTART->ResetDigits();
311 pOutStartLoader->UnloadDigits();
316 //------------------------------------------------------------------------
317 Bool_t AliSTARTDigitizer::RegisterPhotoE(Int_t ipmt,Double_t energy)
321 // Float_t hc=197.326960*1.e6; //mev*nm
322 Double_t hc=1.973*1.e-6; //gev*nm
323 Float_t lambda=hc/energy;
324 Float_t eff = ((TGraph*) fEffPMT.At(ipmt))->Eval(lambda);
325 Double_t p = gRandom->Rndm();
333 //----------------------------------------------------------------------------