]>
Commit | Line | Data |
---|---|---|
990119d6 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
c65c502a | 10 | * copies and that both the copyright notice and this permission notice * |
11 | * appear in the supporting documentation. The authors make no claims * | |
990119d6 | 12 | * about the suitability of this software for any purpose. It is * |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
c093f031 | 16 | |
990119d6 | 17 | /* $Id$ */ |
18 | ||
702ab87e | 19 | /* History of cvs commits: |
20 | * | |
21 | * $Log$ | |
04236e67 | 22 | * Revision 1.56 2007/10/19 18:04:29 schutz |
23 | * The standalone QA data maker is called from AliSimulation and AliReconstruction outside the event loop; i.e. re-reading the data. The QA data making in the event loop has been commented out. | |
24 | * | |
c65c502a | 25 | * Revision 1.55 2007/10/14 21:08:10 schutz |
26 | * Introduced the checking of QA results from previous step before entering the event loop | |
27 | * | |
8661738e | 28 | * Revision 1.54 2007/10/10 09:05:10 schutz |
29 | * Changing name QualAss to QA | |
30 | * | |
b8274834 | 31 | * Revision 1.53 2007/09/30 17:08:20 schutz |
32 | * Introducing the notion of QA data acquisition cycle (needed by online) | |
33 | * | |
5b188f2f | 34 | * Revision 1.52 2007/09/26 14:22:18 cvetan |
35 | * Important changes to the reconstructor classes. Complete elimination of the run-loaders, which are now steered only from AliReconstruction. Removal of the corresponding Reconstruct() and FillESD() methods. | |
36 | * | |
d76c31f4 | 37 | * Revision 1.51 2007/08/07 14:12:03 kharlov |
38 | * Quality assurance added (Yves Schutz) | |
702ab87e | 39 | */ |
40 | ||
990119d6 | 41 | //_________________________________________________________________________ |
7acf6008 | 42 | // This is a TTask that makes SDigits out of Hits |
7b7c1533 | 43 | // The name of the TTask is also the title of the branch that will contain |
44 | // the created SDigits | |
45 | // The title of the TTAsk is the name of the file that contains the hits from | |
46 | // which the SDigits are created | |
7acf6008 | 47 | // A Summable Digits is the sum of all hits originating |
48 | // from one primary in one active cell | |
49 | // A threshold for assignment of the primary to SDigit is applied | |
50 | // SDigits are written to TreeS, branch "PHOS" | |
51 | // AliPHOSSDigitizer with all current parameters is written | |
52 | // to TreeS branch "AliPHOSSDigitizer". | |
f035f6ce | 53 | // Both branches have the same title. If necessary one can produce |
54 | // another set of SDigits with different parameters. Two versions | |
55 | // can be distunguished using titles of the branches. | |
56 | // User case: | |
a4e98857 | 57 | // root [0] AliPHOSSDigitizer * s = new AliPHOSSDigitizer("galice.root") |
58 | // Warning in <TDatabasePDG::TDatabasePDG>: object already instantiated | |
f21fc003 | 59 | // root [1] s->Digitize() |
f035f6ce | 60 | // // Makes SDigitis for all events stored in galice.root |
a4e98857 | 61 | // root [2] s->SetPedestalParameter(0.001) |
f035f6ce | 62 | // // One can change parameters of digitization |
3de072dc | 63 | // root [3] s->SetSDigitsBranch("Pedestal 0.001") |
f035f6ce | 64 | // // and write them into the new branch |
3de072dc | 65 | // root [4] s->ExecuteTask("deb all tim") |
f035f6ce | 66 | // // available parameters: |
67 | // deb - print # of produced SDigitis | |
68 | // deb all - print # and list of produced SDigits | |
69 | // tim - print benchmarking information | |
990119d6 | 70 | // |
6c8cd883 | 71 | //-- Author : Dmitri Peressounko (SUBATECH & KI) |
990119d6 | 72 | ////////////////////////////////////////////////////////////////////////////// |
73 | ||
f035f6ce | 74 | |
990119d6 | 75 | // --- ROOT system --- |
7acf6008 | 76 | #include "TBenchmark.h" |
106fc2fa | 77 | |
990119d6 | 78 | // --- Standard library --- |
79 | ||
80 | // --- AliRoot header files --- | |
9d6866ec | 81 | #include "AliLog.h" |
e957fea8 | 82 | #include "AliPHOSGeometry.h" |
990119d6 | 83 | #include "AliPHOSDigit.h" |
bfae5a5d | 84 | #include "AliPHOSLoader.h" |
990119d6 | 85 | #include "AliPHOSHit.h" |
990119d6 | 86 | #include "AliPHOSSDigitizer.h" |
87 | ||
990119d6 | 88 | ClassImp(AliPHOSSDigitizer) |
89 | ||
90 | ||
91 | //____________________________________________________________________________ | |
3663622c | 92 | AliPHOSSDigitizer::AliPHOSSDigitizer() : |
f21fc003 | 93 | TNamed("",""), |
3663622c | 94 | fPrimThreshold(0.f), |
95 | fDefaultInit(kTRUE), | |
96 | fEventFolderName(""), | |
97 | fInit(kFALSE), | |
98 | fSDigitsInRun(0), | |
99 | fFirstEvent(0), | |
c65c502a | 100 | fLastEvent(0) |
548f0134 | 101 | { |
990119d6 | 102 | // ctor |
ddd1a39c | 103 | // Intialize the quality assurance data maker |
990119d6 | 104 | } |
7acf6008 | 105 | |
990119d6 | 106 | //____________________________________________________________________________ |
e191bb57 | 107 | AliPHOSSDigitizer::AliPHOSSDigitizer(const char * alirunFileName, |
108 | const char * eventFolderName): | |
f21fc003 | 109 | TNamed("PHOSSDigitizer", alirunFileName), |
3663622c | 110 | fPrimThreshold(0.f), |
111 | fDefaultInit(kFALSE), | |
112 | fEventFolderName(eventFolderName), | |
113 | fInit(kFALSE), | |
114 | fSDigitsInRun(0), | |
115 | fFirstEvent(0), | |
c65c502a | 116 | fLastEvent(0) |
990119d6 | 117 | { |
118 | // ctor | |
8d0f3f77 | 119 | InitParameters() ; |
2bd5457f | 120 | Init(); |
92f521a9 | 121 | fDefaultInit = kFALSE ; |
ddd1a39c | 122 | } |
123 | ||
124 | //____________________________________________________________________________ | |
125 | AliPHOSSDigitizer::AliPHOSSDigitizer(const AliPHOSSDigitizer& sd) : | |
f21fc003 | 126 | TNamed(sd.GetName(), sd.GetTitle()), |
ddd1a39c | 127 | fPrimThreshold(sd.fPrimThreshold), |
128 | fDefaultInit(kFALSE), | |
129 | fEventFolderName(sd.fEventFolderName), | |
130 | fInit(kFALSE), | |
131 | fSDigitsInRun(sd.fSDigitsInRun), | |
132 | fFirstEvent(sd.fFirstEvent), | |
c65c502a | 133 | fLastEvent(sd.fLastEvent) |
ddd1a39c | 134 | { |
135 | // cpy ctor | |
ddd1a39c | 136 | } |
137 | ||
ddd1a39c | 138 | //_____________________________________________________________________________ |
139 | AliPHOSSDigitizer& AliPHOSSDigitizer::operator = (const AliPHOSSDigitizer& qa) | |
140 | { | |
141 | // assignment operator | |
142 | ||
143 | this->~AliPHOSSDigitizer(); | |
144 | new(this) AliPHOSSDigitizer(qa); | |
145 | return *this; | |
990119d6 | 146 | } |
147 | ||
7acf6008 | 148 | //____________________________________________________________________________ |
f21fc003 | 149 | AliPHOSSDigitizer::~AliPHOSSDigitizer() |
150 | { | |
797e311f | 151 | //dtor |
797e311f | 152 | } |
ddd1a39c | 153 | |
797e311f | 154 | //____________________________________________________________________________ |
a4e98857 | 155 | void AliPHOSSDigitizer::Init() |
156 | { | |
bfae5a5d | 157 | // Uses the Loader to access the required files |
7b7c1533 | 158 | |
88cb7938 | 159 | fInit = kTRUE ; |
160 | ||
bfae5a5d | 161 | //to prevent cleaning of this object while GetEvent is called |
162 | AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; | |
95a38b7d | 163 | if (!rl) AliRunLoader::Open(GetTitle(), fEventFolderName) ; |
7acf6008 | 164 | } |
7b7c1533 | 165 | |
8d0f3f77 | 166 | //____________________________________________________________________________ |
167 | void AliPHOSSDigitizer::InitParameters() | |
168 | { | |
88cb7938 | 169 | // initializes the parameters for digitization |
8d0f3f77 | 170 | fPrimThreshold = 0.01 ; |
171 | fSDigitsInRun = 0 ; | |
8d0f3f77 | 172 | } |
173 | ||
990119d6 | 174 | //____________________________________________________________________________ |
f21fc003 | 175 | void AliPHOSSDigitizer::Digitize(Option_t *option) |
a4e98857 | 176 | { |
212d1c0f | 177 | // Steering method to produce summable digits for events |
178 | // in the range from fFirstEvent to fLastEvent. | |
179 | // This range is optionally set by SetEventRange(). | |
180 | // if fLastEvent=-1 (by default), then process events until the end. | |
181 | // | |
182 | // Summable digit is a sum of all hits in the same active | |
183 | // volume into digit | |
7b7c1533 | 184 | |
185 | if (strstr(option, "print") ) { | |
88cb7938 | 186 | Print() ; |
7b7c1533 | 187 | return ; |
188 | } | |
990119d6 | 189 | |
7acf6008 | 190 | if(strstr(option,"tim")) |
191 | gBenchmark->Start("PHOSSDigitizer"); | |
fbf811ec | 192 | |
04236e67 | 193 | /* |
194 | // check the QA result for RAWS | |
4e25ac79 | 195 | AliQAv1 * qa = AliQAv1::Instance(AliQAv1::kPHOS) ; |
196 | if ( qa->IsSet(AliQAv1::kPHOS, AliQAv1::kRAW, AliQAv1::kFATAL)) { | |
8661738e | 197 | AliFatal("QA status in RAW was Fatal") ; |
4e25ac79 | 198 | } else if ( qa->IsSet(AliQAv1::kPHOS, AliQAv1::kRAW, AliQAv1::kERROR)) { |
8661738e | 199 | AliError("QA status in RAW was Error") ; |
4e25ac79 | 200 | } else if ( qa->IsSet(AliQAv1::kPHOS, AliQAv1::kRAW, AliQAv1::kWARNING) ) { |
8661738e | 201 | AliWarning("QA status in RAW was Warning") ; |
4e25ac79 | 202 | } else if ( qa->IsSet(AliQAv1::kPHOS, AliQAv1::kRAW, AliQAv1::kINFO) ) { |
8661738e | 203 | AliInfo("QA status in RAW was Info") ; |
204 | } | |
04236e67 | 205 | */ |
bfae5a5d | 206 | AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; |
daa897f5 | 207 | AliPHOSLoader * phosLoader = static_cast<AliPHOSLoader*>(rl->GetLoader("PHOSLoader")); |
825dc01a | 208 | |
88cb7938 | 209 | //switch off reloading of this task while getting event |
210 | if (!fInit) { // to prevent overwrite existing file | |
9d6866ec | 211 | AliError( Form("Give a version name different from %s", fEventFolderName.Data()) ) ; |
88cb7938 | 212 | return ; |
213 | } | |
daa897f5 | 214 | |
212d1c0f | 215 | if (fLastEvent == -1) |
bfae5a5d | 216 | fLastEvent = rl->GetNumberOfEvents() - 1 ; |
212d1c0f | 217 | else |
bfae5a5d | 218 | fLastEvent = TMath::Min(fFirstEvent, rl->GetNumberOfEvents()); // only one event at the time |
212d1c0f | 219 | Int_t nEvents = fLastEvent - fFirstEvent + 1; |
220 | ||
e3310962 | 221 | Int_t ievent, i; |
9d6866ec | 222 | |
9d6866ec | 223 | for (ievent = fFirstEvent; ievent <= fLastEvent; ievent++) { |
bfae5a5d | 224 | rl->GetEvent(ievent) ; |
225 | TTree * treeS = phosLoader->TreeS(); | |
226 | if(!treeS){ | |
227 | phosLoader->MakeTree("S"); | |
228 | treeS = phosLoader->TreeS(); | |
229 | } | |
230 | ||
231 | phosLoader->CleanHits() ; | |
232 | phosLoader->LoadHits("READ") ; | |
233 | ||
234 | TClonesArray * hits = phosLoader->Hits() ; | |
235 | TClonesArray * sdigits = phosLoader->SDigits() ; | |
236 | if( !sdigits ) { | |
237 | phosLoader->MakeSDigitsArray() ; | |
238 | sdigits = phosLoader->SDigits() ; | |
239 | } | |
7b7c1533 | 240 | sdigits->Clear(); |
990119d6 | 241 | Int_t nSdigits = 0 ; |
bfae5a5d | 242 | |
7b7c1533 | 243 | //Now make SDigits from hits, for PHOS it is the same, so just copy |
e3310962 | 244 | for ( i = 0 ; i < hits->GetEntries() ; i++ ) { |
c78764fb | 245 | |
daa897f5 | 246 | AliPHOSHit * hit = static_cast<AliPHOSHit *>(hits->At(i)) ; |
e3310962 | 247 | // Assign primary number only if contribution is significant |
248 | ||
249 | if( hit->GetEnergy() > fPrimThreshold) | |
250 | new((*sdigits)[nSdigits]) AliPHOSDigit(hit->GetPrimary(),hit->GetId(), | |
251 | hit->GetEnergy() ,hit->GetTime()) ; | |
252 | else | |
253 | new((*sdigits)[nSdigits]) AliPHOSDigit(-1 ,hit->GetId(), | |
254 | hit->GetEnergy() ,hit->GetTime()) ; | |
255 | nSdigits++ ; | |
256 | ||
257 | } | |
88cb7938 | 258 | |
7b7c1533 | 259 | sdigits->Sort() ; |
88cb7938 | 260 | |
7b7c1533 | 261 | nSdigits = sdigits->GetEntriesFast() ; |
88cb7938 | 262 | |
69866bef | 263 | fSDigitsInRun += nSdigits ; |
7b7c1533 | 264 | sdigits->Expand(nSdigits) ; |
daa897f5 | 265 | |
266 | for (i = 0 ; i < nSdigits ; i++) { | |
267 | AliPHOSDigit * digit = static_cast<AliPHOSDigit *>(sdigits->At(i)) ; | |
990119d6 | 268 | digit->SetIndexInList(i) ; |
269 | } | |
daa897f5 | 270 | |
c65c502a | 271 | // // make Quality Assurance data |
272 | // | |
273 | // if (GetQADataMaker()->IsCycleDone() ) { | |
4e25ac79 | 274 | // GetQADataMaker()->EndOfCycle(AliQAv1::kHITS) ; |
275 | // GetQADataMaker()->EndOfCycle(AliQAv1::kSDIGITS) ; | |
276 | // GetQADataMaker()->StartOfCycle(AliQAv1::kHITS) ; | |
277 | // GetQADataMaker()->StartOfCycle(AliQAv1::kSDIGITS, kTRUE) ; | |
c65c502a | 278 | // } |
4e25ac79 | 279 | // GetQADataMaker()->Exec(AliQAv1::kHITS, hits) ; |
280 | // GetQADataMaker()->Exec(AliQAv1::kSDIGITS, sdigits) ; | |
c65c502a | 281 | // GetQADataMaker()->Increment() ; |
5b188f2f | 282 | |
fbf811ec | 283 | //Now write SDigits |
fbf811ec | 284 | |
88cb7938 | 285 | |
7acf6008 | 286 | //First list of sdigits |
88cb7938 | 287 | |
fbf811ec | 288 | Int_t bufferSize = 32000 ; |
88cb7938 | 289 | TBranch * sdigitsBranch = treeS->Branch("PHOS",&sdigits,bufferSize); |
fbf811ec | 290 | sdigitsBranch->Fill() ; |
fbf811ec | 291 | |
bfae5a5d | 292 | phosLoader->WriteSDigits("OVERWRITE"); |
88cb7938 | 293 | |
7acf6008 | 294 | if(strstr(option,"deb")) |
295 | PrintSDigits(option) ; | |
9d6866ec | 296 | |
9d6866ec | 297 | }// event loop |
fbf811ec | 298 | |
c65c502a | 299 | // //Write the quality assurance data |
4e25ac79 | 300 | // GetQADataMaker()->EndOfCycle(AliQAv1::kHITS) ; |
301 | // GetQADataMaker()->EndOfCycle(AliQAv1::kSDIGITS) ; | |
c65c502a | 302 | // GetQADataMaker()->Finish() ; |
ddd1a39c | 303 | |
88cb7938 | 304 | Unload(); |
305 | ||
7acf6008 | 306 | if(strstr(option,"tim")){ |
307 | gBenchmark->Stop("PHOSSDigitizer"); | |
21cd0c07 | 308 | Info("Exec"," took %f seconds for SDigitizing %f seconds per event", |
bfae5a5d | 309 | gBenchmark->GetCpuTime("PHOSSDigitizer"), |
310 | gBenchmark->GetCpuTime("PHOSSDigitizer")/nEvents) ; | |
7acf6008 | 311 | } |
9d6866ec | 312 | |
990119d6 | 313 | } |
106fc2fa | 314 | |
990119d6 | 315 | //__________________________________________________________________ |
702ab87e | 316 | void AliPHOSSDigitizer::Print(const Option_t *)const |
a4e98857 | 317 | { |
318 | // Prints parameters of SDigitizer | |
88cb7938 | 319 | Info("Print", "\n------------------- %s -------------", GetName() ) ; |
320 | printf(" Writing SDigits to branch with title %s\n", fEventFolderName.Data()) ; | |
88cb7938 | 321 | printf(" Threshold for Primary assignment= %f\n", fPrimThreshold) ; |
322 | printf("---------------------------------------------------\n") ; | |
7acf6008 | 323 | |
990119d6 | 324 | } |
548f0134 | 325 | |
990119d6 | 326 | //__________________________________________________________________ |
a4e98857 | 327 | Bool_t AliPHOSSDigitizer::operator==( AliPHOSSDigitizer const &sd )const |
328 | { | |
3de072dc | 329 | // Equal operator. |
330 | // SDititizers are equal if their pedestal, slope and threshold are equal | |
331 | ||
f898e0f3 | 332 | if(fPrimThreshold==sd.fPrimThreshold) |
990119d6 | 333 | return kTRUE ; |
334 | else | |
335 | return kFALSE ; | |
336 | } | |
548f0134 | 337 | |
a6eedfad | 338 | //__________________________________________________________________ |
a4e98857 | 339 | void AliPHOSSDigitizer::PrintSDigits(Option_t * option) |
340 | { | |
341 | // Prints list of digits produced in the current pass of AliPHOSDigitizer | |
7b7c1533 | 342 | |
bfae5a5d | 343 | AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; |
c78764fb | 344 | AliPHOSLoader * phosLoader = static_cast<AliPHOSLoader*>(rl->GetLoader("PHOSLoader")); |
6c8cd883 | 345 | |
346 | // Get PHOS Geometry object | |
347 | AliPHOSGeometry *geom; | |
348 | if (!(geom = AliPHOSGeometry::GetInstance())) | |
349 | geom = AliPHOSGeometry::GetInstance("IHEP",""); | |
350 | ||
bfae5a5d | 351 | const TClonesArray * sdigits = phosLoader->SDigits() ; |
21cd0c07 | 352 | |
bfae5a5d | 353 | Info( "\nPrintSDigits", "event # %d %d sdigits", |
354 | gAlice->GetEvNumber(), sdigits->GetEntriesFast() ) ; | |
71dfb0d8 | 355 | |
a6eedfad | 356 | if(strstr(option,"all")||strstr(option,"EMC")){ |
7acf6008 | 357 | |
358 | //loop over digits | |
359 | AliPHOSDigit * digit; | |
71dfb0d8 | 360 | printf("\nEMC sdigits\n") ; |
6c8cd883 | 361 | Int_t maxEmc = geom->GetNModules() * geom->GetNCristalsInModule() ; |
7acf6008 | 362 | Int_t index ; |
a6eedfad | 363 | for (index = 0 ; (index < sdigits->GetEntriesFast()) && |
c78764fb | 364 | ((static_cast<AliPHOSDigit *> (sdigits->At(index)))->GetId() <= maxEmc) ; index++) { |
365 | digit = static_cast<AliPHOSDigit *>( sdigits->At(index) ) ; | |
71dfb0d8 | 366 | // if(digit->GetNprimary() == 0) |
367 | // continue; | |
27a73a5d | 368 | // printf("%6d %8d %6.5e %4d %2d :\n", // YVK |
369 | printf("%6d %.4f %6.5e %4d %2d :\n", | |
370 | digit->GetId(), digit->GetEnergy(), digit->GetTime(), digit->GetIndexInList(), digit->GetNprimary()) ; | |
7acf6008 | 371 | Int_t iprimary; |
11f9c5ff | 372 | for (iprimary=0; iprimary<digit->GetNprimary(); iprimary++) { |
71dfb0d8 | 373 | printf("%d ",digit->GetPrimary(iprimary+1) ) ; |
11f9c5ff | 374 | } |
a6eedfad | 375 | } |
a6eedfad | 376 | } |
377 | ||
378 | if(strstr(option,"all")||strstr(option,"CPV")){ | |
7acf6008 | 379 | |
a6eedfad | 380 | //loop over CPV digits |
381 | AliPHOSDigit * digit; | |
71dfb0d8 | 382 | printf("\nCPV sdigits\n") ; |
6c8cd883 | 383 | Int_t maxEmc = geom->GetNModules() * geom->GetNCristalsInModule() ; |
a6eedfad | 384 | Int_t index ; |
385 | for (index = 0 ; index < sdigits->GetEntriesFast(); index++) { | |
daa897f5 | 386 | digit = static_cast<AliPHOSDigit *>( sdigits->At(index) ) ; |
a6eedfad | 387 | if(digit->GetId() > maxEmc){ |
71dfb0d8 | 388 | printf("\n%6d %8d %4d %2d :", |
11f9c5ff | 389 | digit->GetId(), digit->GetAmp(), digit->GetIndexInList(), digit->GetNprimary()) ; |
a6eedfad | 390 | Int_t iprimary; |
21cd0c07 | 391 | for (iprimary=0; iprimary<digit->GetNprimary(); iprimary++) { |
71dfb0d8 | 392 | printf("%d ",digit->GetPrimary(iprimary+1) ) ; |
21cd0c07 | 393 | } |
a6eedfad | 394 | } |
395 | } | |
7acf6008 | 396 | } |
397 | } | |
7b7c1533 | 398 | |
399 | //____________________________________________________________________________ | |
88cb7938 | 400 | void AliPHOSSDigitizer::Unload() const |
7b7c1533 | 401 | { |
e957fea8 | 402 | // Unloads the objects from the folder |
bfae5a5d | 403 | AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; |
daa897f5 | 404 | AliPHOSLoader * phosLoader = static_cast<AliPHOSLoader*>(rl->GetLoader("PHOSLoader")); |
c78764fb | 405 | |
bfae5a5d | 406 | phosLoader->UnloadHits() ; |
407 | phosLoader->UnloadSDigits() ; | |
7b7c1533 | 408 | } |