X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOSSDigitizer.cxx;h=78fe69dc90c7ed7654088bacf4220ab6ab49e692;hb=a8c512112de63fcf3e88eaa15e4d5f3d78427e9d;hp=65a73193f25143a1065f057b5046d4ab8cd1f7c9;hpb=702ab87e845509d9379dca12e47c7e8a37eba7f6;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOSSDigitizer.cxx b/PHOS/AliPHOSSDigitizer.cxx index 65a73193f25..78fe69dc90c 100644 --- a/PHOS/AliPHOSSDigitizer.cxx +++ b/PHOS/AliPHOSSDigitizer.cxx @@ -19,6 +19,23 @@ /* History of cvs commits: * * $Log$ + * Revision 1.56 2007/10/19 18:04:29 schutz + * 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. + * + * Revision 1.55 2007/10/14 21:08:10 schutz + * Introduced the checking of QA results from previous step before entering the event loop + * + * Revision 1.54 2007/10/10 09:05:10 schutz + * Changing name QualAss to QA + * + * Revision 1.53 2007/09/30 17:08:20 schutz + * Introducing the notion of QA data acquisition cycle (needed by online) + * + * Revision 1.52 2007/09/26 14:22:18 cvetan + * 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. + * + * Revision 1.51 2007/08/07 14:12:03 kharlov + * Quality assurance added (Yves Schutz) */ //_________________________________________________________________________ @@ -51,13 +68,12 @@ // deb all - print # and list of produced SDigits // tim - print benchmarking information // -//*-- Author : Dmitri Peressounko (SUBATECH & KI) +//-- Author : Dmitri Peressounko (SUBATECH & KI) ////////////////////////////////////////////////////////////////////////////// // --- ROOT system --- #include "TBenchmark.h" - //#include "TObjectTable.h" // --- Standard library --- @@ -65,93 +81,102 @@ #include "AliLog.h" #include "AliPHOSGeometry.h" #include "AliPHOSDigit.h" -#include "AliPHOSGetter.h" +#include "AliPHOSLoader.h" #include "AliPHOSHit.h" #include "AliPHOSSDigitizer.h" - //#include "AliMemoryWatcher.h" ClassImp(AliPHOSSDigitizer) //____________________________________________________________________________ - AliPHOSSDigitizer::AliPHOSSDigitizer():TTask("","") +AliPHOSSDigitizer::AliPHOSSDigitizer() : + TTask("",""), + fPrimThreshold(0.f), + fDefaultInit(kTRUE), + fEventFolderName(""), + fInit(kFALSE), + fSDigitsInRun(0), + fFirstEvent(0), + fLastEvent(0) { // ctor - fFirstEvent = fLastEvent = 0 ; - fDefaultInit = kTRUE ; + // Intialize the quality assurance data maker } //____________________________________________________________________________ AliPHOSSDigitizer::AliPHOSSDigitizer(const char * alirunFileName, const char * eventFolderName): TTask("PHOS"+AliConfig::Instance()->GetSDigitizerTaskName(), alirunFileName), - fEventFolderName(eventFolderName) + fPrimThreshold(0.f), + fDefaultInit(kFALSE), + fEventFolderName(eventFolderName), + fInit(kFALSE), + fSDigitsInRun(0), + fFirstEvent(0), + fLastEvent(0) { - // ctor - fFirstEvent = fLastEvent = 0 ; // runs one event by defaut InitParameters() ; Init(); fDefaultInit = kFALSE ; } -//____________________________________________________________________________ -AliPHOSSDigitizer::AliPHOSSDigitizer(const AliPHOSSDigitizer & sd) - : TTask(sd) -{ - //cpy ctor - - fFirstEvent = sd.fFirstEvent ; - fLastEvent = sd.fLastEvent ; - fA = sd.fA ; - fB = sd.fB ; - fPrimThreshold = sd.fPrimThreshold ; - fSDigitsInRun = sd.fSDigitsInRun ; - SetName(sd.GetName()) ; - SetTitle(sd.GetTitle()) ; - fEventFolderName = sd.fEventFolderName; +//____________________________________________________________________________ +AliPHOSSDigitizer::AliPHOSSDigitizer(const AliPHOSSDigitizer& sd) : + TTask(sd.GetName(), sd.GetTitle()), + fPrimThreshold(sd.fPrimThreshold), + fDefaultInit(kFALSE), + fEventFolderName(sd.fEventFolderName), + fInit(kFALSE), + fSDigitsInRun(sd.fSDigitsInRun), + fFirstEvent(sd.fFirstEvent), + fLastEvent(sd.fLastEvent) +{ + // cpy ctor } +//_____________________________________________________________________________ +AliPHOSSDigitizer& AliPHOSSDigitizer::operator = (const AliPHOSSDigitizer& qa) +{ +// assignment operator + + this->~AliPHOSSDigitizer(); + new(this) AliPHOSSDigitizer(qa); + return *this; +} //____________________________________________________________________________ AliPHOSSDigitizer::~AliPHOSSDigitizer() { //dtor - // AliPHOSGetter * gime = - // AliPHOSGetter::Instance(GetTitle(),fEventFolderName.Data()); - AliPHOSGetter * gime = - AliPHOSGetter::Instance(); - gime->PhosLoader()->CleanSDigitizer(); + AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; + if(rl){ + AliPHOSLoader * phosLoader = + dynamic_cast(rl->GetLoader("PHOSLoader")); + phosLoader->CleanSDigitizer() ; + } } + //____________________________________________________________________________ void AliPHOSSDigitizer::Init() { - // Uses the getter to access the required files + // Uses the Loader to access the required files fInit = kTRUE ; - AliPHOSGetter * gime = AliPHOSGetter::Instance(GetTitle(), fEventFolderName.Data()); - if ( gime == 0 ) { - Fatal("Init" ,"Could not obtain the Getter object for file %s and event %s !", GetTitle(), fEventFolderName.Data()) ; - return ; - } - - TString opt("SDigits") ; - if(gime->VersionExists(opt) ) { - Error( "Init", "Give a version name different from %s", fEventFolderName.Data() ) ; - fInit = kFALSE ; - } - - gime->PostSDigitizer(this); - gime->PhosLoader()->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE); - + //to prevent cleaning of this object while GetEvent is called + AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; + if (!rl) + rl = AliRunLoader::Open(GetTitle(), fEventFolderName) ; + + AliPHOSLoader * phosLoader = dynamic_cast(rl->GetLoader("PHOSLoader")); + phosLoader->PostSDigitizer(this); + phosLoader->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE); } //____________________________________________________________________________ void AliPHOSSDigitizer::InitParameters() { // initializes the parameters for digitization - fA = 0; - fB = 10000000.; fPrimThreshold = 0.01 ; fSDigitsInRun = 0 ; } @@ -175,7 +200,21 @@ void AliPHOSSDigitizer::Exec(Option_t *option) if(strstr(option,"tim")) gBenchmark->Start("PHOSSDigitizer"); - AliPHOSGetter * gime = AliPHOSGetter::Instance() ; +/* + // check the QA result for RAWS + AliQAv1 * qa = AliQAv1::Instance(AliQAv1::kPHOS) ; + if ( qa->IsSet(AliQAv1::kPHOS, AliQAv1::kRAW, AliQAv1::kFATAL)) { + AliFatal("QA status in RAW was Fatal") ; + } else if ( qa->IsSet(AliQAv1::kPHOS, AliQAv1::kRAW, AliQAv1::kERROR)) { + AliError("QA status in RAW was Error") ; + } else if ( qa->IsSet(AliQAv1::kPHOS, AliQAv1::kRAW, AliQAv1::kWARNING) ) { + AliWarning("QA status in RAW was Warning") ; + } else if ( qa->IsSet(AliQAv1::kPHOS, AliQAv1::kRAW, AliQAv1::kINFO) ) { + AliInfo("QA status in RAW was Info") ; + } +*/ + AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; + AliPHOSLoader * phosLoader = dynamic_cast(rl->GetLoader("PHOSLoader")); //switch off reloading of this task while getting event if (!fInit) { // to prevent overwrite existing file @@ -184,48 +223,48 @@ void AliPHOSSDigitizer::Exec(Option_t *option) } if (fLastEvent == -1) - fLastEvent = gime->MaxEvent() - 1 ; + fLastEvent = rl->GetNumberOfEvents() - 1 ; else - fLastEvent = TMath::Min(fFirstEvent, gime->MaxEvent()); // only ine event at the time + fLastEvent = TMath::Min(fFirstEvent, rl->GetNumberOfEvents()); // only one event at the time Int_t nEvents = fLastEvent - fFirstEvent + 1; - Int_t ievent ; - - //AliMemoryWatcher memwatcher; + Int_t ievent, i; for (ievent = fFirstEvent; ievent <= fLastEvent; ievent++) { - gime->Event(ievent,"H") ; - TTree * treeS = gime->TreeS(); - TClonesArray * hits = gime->Hits() ; - TClonesArray * sdigits = gime->SDigits() ; + rl->GetEvent(ievent) ; + TTree * treeS = phosLoader->TreeS(); + if(!treeS){ + phosLoader->MakeTree("S"); + treeS = phosLoader->TreeS(); + } + + phosLoader->CleanHits() ; + phosLoader->LoadHits("READ") ; + + TClonesArray * hits = phosLoader->Hits() ; + TClonesArray * sdigits = phosLoader->SDigits() ; + if( !sdigits ) { + phosLoader->MakeSDigitsArray() ; + sdigits = phosLoader->SDigits() ; + } sdigits->Clear(); Int_t nSdigits = 0 ; + //Now make SDigits from hits, for PHOS it is the same, so just copy - Int_t nPrim = static_cast((gime->TreeH())->GetEntries()) ; - // Attention nPrim is the number of primaries tracked by Geant - // and this number could be different to the number of Primaries in TreeK; - Int_t iprim ; - - for (iprim = 0 ; iprim < nPrim ; iprim ++) { - //=========== Get the PHOS branch from Hits Tree for the Primary iprim - gime->Track(iprim) ; - Int_t i; - for ( i = 0 ; i < hits->GetEntries() ; i++ ) { - AliPHOSHit * hit = dynamic_cast(hits->At(i)) ; - // Assign primary number only if contribution is significant - - if( hit->GetEnergy() > fPrimThreshold) - new((*sdigits)[nSdigits]) AliPHOSDigit(hit->GetPrimary(),hit->GetId(), - Digitize(hit->GetEnergy()), hit->GetTime()) ; - else - new((*sdigits)[nSdigits]) AliPHOSDigit( -1 , hit->GetId(), - Digitize(hit->GetEnergy()), hit->GetTime()) ; - nSdigits++ ; - - } + for ( i = 0 ; i < hits->GetEntries() ; i++ ) { + AliPHOSHit * hit = dynamic_cast(hits->At(i)) ; + // Assign primary number only if contribution is significant + + if( hit->GetEnergy() > fPrimThreshold) + new((*sdigits)[nSdigits]) AliPHOSDigit(hit->GetPrimary(),hit->GetId(), + hit->GetEnergy() ,hit->GetTime()) ; + else + new((*sdigits)[nSdigits]) AliPHOSDigit(-1 ,hit->GetId(), + hit->GetEnergy() ,hit->GetTime()) ; + nSdigits++ ; + + } - } // loop over iprim - sdigits->Sort() ; nSdigits = sdigits->GetEntriesFast() ; @@ -233,12 +272,23 @@ void AliPHOSSDigitizer::Exec(Option_t *option) fSDigitsInRun += nSdigits ; sdigits->Expand(nSdigits) ; - Int_t i ; for (i = 0 ; i < nSdigits ; i++) { AliPHOSDigit * digit = dynamic_cast(sdigits->At(i)) ; digit->SetIndexInList(i) ; } +// // make Quality Assurance data +// +// if (GetQADataMaker()->IsCycleDone() ) { +// GetQADataMaker()->EndOfCycle(AliQAv1::kHITS) ; +// GetQADataMaker()->EndOfCycle(AliQAv1::kSDIGITS) ; +// GetQADataMaker()->StartOfCycle(AliQAv1::kHITS) ; +// GetQADataMaker()->StartOfCycle(AliQAv1::kSDIGITS, kTRUE) ; +// } +// GetQADataMaker()->Exec(AliQAv1::kHITS, hits) ; +// GetQADataMaker()->Exec(AliQAv1::kSDIGITS, sdigits) ; +// GetQADataMaker()->Increment() ; + //Now write SDigits @@ -246,35 +296,30 @@ void AliPHOSSDigitizer::Exec(Option_t *option) Int_t bufferSize = 32000 ; TBranch * sdigitsBranch = treeS->Branch("PHOS",&sdigits,bufferSize); - sdigitsBranch->Fill() ; - gime->WriteSDigits("OVERWRITE"); + phosLoader->WriteSDigits("OVERWRITE"); + phosLoader->WriteSDigitizer("OVERWRITE"); - //Next - SDigitizer - - gime->WriteSDigitizer("OVERWRITE"); - //gObjectTable->Print() ; - if(strstr(option,"deb")) PrintSDigits(option) ; - //memwatcher.Watch(ievent); }// event loop +// //Write the quality assurance data +// GetQADataMaker()->EndOfCycle(AliQAv1::kHITS) ; +// GetQADataMaker()->EndOfCycle(AliQAv1::kSDIGITS) ; +// GetQADataMaker()->Finish() ; + Unload(); - // gime->PhosLoader()->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE); - if(strstr(option,"tim")){ gBenchmark->Stop("PHOSSDigitizer"); Info("Exec"," took %f seconds for SDigitizing %f seconds per event", - gBenchmark->GetCpuTime("PHOSSDigitizer"), gBenchmark->GetCpuTime("PHOSSDigitizer")/nEvents) ; + gBenchmark->GetCpuTime("PHOSSDigitizer"), + gBenchmark->GetCpuTime("PHOSSDigitizer")/nEvents) ; } - //TFile f("out.root","RECREATE"); - //memwatcher.WriteToFile(); - //f.Close(); } //__________________________________________________________________ @@ -283,8 +328,6 @@ void AliPHOSSDigitizer::Print(const Option_t *)const // Prints parameters of SDigitizer Info("Print", "\n------------------- %s -------------", GetName() ) ; printf(" Writing SDigits to branch with title %s\n", fEventFolderName.Data()) ; - printf(" with digitization parameters A = %f\n", fA) ; - printf(" B = %f\n", fB) ; printf(" Threshold for Primary assignment= %f\n", fPrimThreshold) ; printf("---------------------------------------------------\n") ; @@ -296,7 +339,7 @@ Bool_t AliPHOSSDigitizer::operator==( AliPHOSSDigitizer const &sd )const // Equal operator. // SDititizers are equal if their pedestal, slope and threshold are equal - if( (fA==sd.fA)&&(fB==sd.fB)&&(fPrimThreshold==sd.fPrimThreshold)) + if(fPrimThreshold==sd.fPrimThreshold) return kTRUE ; else return kFALSE ; @@ -307,26 +350,34 @@ void AliPHOSSDigitizer::PrintSDigits(Option_t * option) { // Prints list of digits produced in the current pass of AliPHOSDigitizer + AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; + AliPHOSLoader * phosLoader = dynamic_cast(rl->GetLoader("PHOSLoader")); + + // Get PHOS Geometry object + AliPHOSGeometry *geom; + if (!(geom = AliPHOSGeometry::GetInstance())) + geom = AliPHOSGeometry::GetInstance("IHEP",""); - AliPHOSGetter * gime = AliPHOSGetter::Instance() ; - const TClonesArray * sdigits = gime->SDigits() ; + const TClonesArray * sdigits = phosLoader->SDigits() ; - Info( "\nPrintSDigits", "event # %d %d sdigits", gAlice->GetEvNumber(), sdigits->GetEntriesFast() ) ; + Info( "\nPrintSDigits", "event # %d %d sdigits", + gAlice->GetEvNumber(), sdigits->GetEntriesFast() ) ; if(strstr(option,"all")||strstr(option,"EMC")){ //loop over digits AliPHOSDigit * digit; printf("\nEMC sdigits\n") ; - Int_t maxEmc = gime->PHOSGeometry()->GetNModules()*gime->PHOSGeometry()->GetNCristalsInModule() ; + Int_t maxEmc = geom->GetNModules() * geom->GetNCristalsInModule() ; Int_t index ; for (index = 0 ; (index < sdigits->GetEntriesFast()) && ((dynamic_cast (sdigits->At(index)))->GetId() <= maxEmc) ; index++) { digit = dynamic_cast( sdigits->At(index) ) ; // if(digit->GetNprimary() == 0) // continue; - printf("%6d %8d %6.5e %4d %2d :\n", - digit->GetId(), digit->GetAmp(), digit->GetTime(), digit->GetIndexInList(), digit->GetNprimary()) ; +// printf("%6d %8d %6.5e %4d %2d :\n", // YVK + printf("%6d %.4f %6.5e %4d %2d :\n", + digit->GetId(), digit->GetEnergy(), digit->GetTime(), digit->GetIndexInList(), digit->GetNprimary()) ; Int_t iprimary; for (iprimary=0; iprimaryGetNprimary(); iprimary++) { printf("%d ",digit->GetPrimary(iprimary+1) ) ; @@ -339,7 +390,7 @@ void AliPHOSSDigitizer::PrintSDigits(Option_t * option) //loop over CPV digits AliPHOSDigit * digit; printf("\nCPV sdigits\n") ; - Int_t maxEmc = gime->PHOSGeometry()->GetNModules()*gime->PHOSGeometry()->GetNCristalsInModule() ; + Int_t maxEmc = geom->GetNModules() * geom->GetNCristalsInModule() ; Int_t index ; for (index = 0 ; index < sdigits->GetEntriesFast(); index++) { digit = dynamic_cast( sdigits->At(index) ) ; @@ -359,8 +410,8 @@ void AliPHOSSDigitizer::PrintSDigits(Option_t * option) void AliPHOSSDigitizer::Unload() const { // Unloads the objects from the folder - AliPHOSGetter * gime = AliPHOSGetter::Instance() ; - AliPHOSLoader * loader = gime->PhosLoader() ; - loader->UnloadHits() ; - loader->UnloadSDigits() ; + AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; + AliPHOSLoader * phosLoader = dynamic_cast(rl->GetLoader("PHOSLoader")); + phosLoader->UnloadHits() ; + phosLoader->UnloadSDigits() ; }