X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOSSDigitizer.cxx;h=78fe69dc90c7ed7654088bacf4220ab6ab49e692;hb=7e21be4955fcc878edc0e41895ffce541636d2bf;hp=113a6367794e6d9ccec532befa3b1621a970e1f4;hpb=7b7c15335990bc27772c67f9ee7bc90cc48bc9a3;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOSSDigitizer.cxx b/PHOS/AliPHOSSDigitizer.cxx index 113a6367794..78fe69dc90c 100644 --- a/PHOS/AliPHOSSDigitizer.cxx +++ b/PHOS/AliPHOSSDigitizer.cxx @@ -13,8 +13,31 @@ * provided "as is" without express or implied warranty. * **************************************************************************/ + /* $Id$ */ +/* 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) + */ + //_________________________________________________________________________ // This is a TTask that makes SDigits out of Hits // The name of the TTask is also the title of the branch that will contain @@ -45,328 +68,350 @@ // 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 "TFile.h" -#include "TTask.h" -#include "TTree.h" -#include "TSystem.h" -#include "TROOT.h" -#include "TFolder.h" #include "TBenchmark.h" + // --- Standard library --- -#include // --- AliRoot header files --- -#include "AliRun.h" +#include "AliLog.h" +#include "AliPHOSGeometry.h" #include "AliPHOSDigit.h" -#include "AliPHOSGetter.h" +#include "AliPHOSLoader.h" #include "AliPHOSHit.h" #include "AliPHOSSDigitizer.h" - ClassImp(AliPHOSSDigitizer) //____________________________________________________________________________ - AliPHOSSDigitizer::AliPHOSSDigitizer():TTask("","") +AliPHOSSDigitizer::AliPHOSSDigitizer() : + TTask("",""), + fPrimThreshold(0.f), + fDefaultInit(kTRUE), + fEventFolderName(""), + fInit(kFALSE), + fSDigitsInRun(0), + fFirstEvent(0), + fLastEvent(0) { // ctor - fA = 0; - fB = 10000000.; - fPrimThreshold = 0.01 ; + // Intialize the quality assurance data maker } //____________________________________________________________________________ -AliPHOSSDigitizer::AliPHOSSDigitizer(const char * headerFile, const char * sDigitsTitle):TTask(sDigitsTitle, headerFile) +AliPHOSSDigitizer::AliPHOSSDigitizer(const char * alirunFileName, + const char * eventFolderName): + TTask("PHOS"+AliConfig::Instance()->GetSDigitizerTaskName(), alirunFileName), + fPrimThreshold(0.f), + fDefaultInit(kFALSE), + fEventFolderName(eventFolderName), + fInit(kFALSE), + fSDigitsInRun(0), + fFirstEvent(0), + fLastEvent(0) { // ctor - fA = 0; - fB = 10000000.; - fPrimThreshold = 0.01 ; + InitParameters() ; Init(); + fDefaultInit = kFALSE ; } +//____________________________________________________________________________ +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 + AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; + if(rl){ + AliPHOSLoader * phosLoader = + dynamic_cast(rl->GetLoader("PHOSLoader")); + phosLoader->CleanSDigitizer() ; + } +} //____________________________________________________________________________ void AliPHOSSDigitizer::Init() { - // Initialization: open root-file, allocate arrays for hits and sdigits, - // attach task SDigitizer to the list of PHOS tasks - // - // Initialization can not be done in the default constructor - //============================================================= YS - // The initialisation is now done by AliPHOSGetter + // Uses the Loader to access the required files - if( strcmp(GetTitle(), "") == 0 ) - SetTitle("galice.root") ; + fInit = kTRUE ; - AliPHOSGetter * gime = AliPHOSGetter::GetInstance(GetTitle(), GetName()) ; - if ( gime == 0 ) { - cerr << "ERROR: AliPHOSSDigitizer::Init -> Could not obtain the Getter object !" << endl ; - return ; - } - - //add Task to //YSAlice/tasks/SDiditizer/PHOS - TTask * aliceSD = (TTask*)gROOT->FindObjectAny("YSAlice/tasks/SDigitizer") ; - TTask * phosSD = (TTask*)aliceSD->GetListOfTasks()->FindObject("PHOS") ; - if ( ! phosSD->GetListOfTasks()->FindObject(GetName()) ) - phosSD->Add(this) ; - // create a folder on the white board //YSAlice/WhiteBoard/SDigits/PHOS/headerFile/sdigitsTitle - gime->Post(GetTitle(), "S", GetName() ) ; - + //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 + fPrimThreshold = 0.01 ; + fSDigitsInRun = 0 ; } //____________________________________________________________________________ void AliPHOSSDigitizer::Exec(Option_t *option) { - // Collects all hits in the same active volume into digit - if( strcmp(GetName(), "") == 0 ) - Init() ; + // Steering method to produce summable digits for events + // in the range from fFirstEvent to fLastEvent. + // This range is optionally set by SetEventRange(). + // if fLastEvent=-1 (by default), then process events until the end. + // + // Summable digit is a sum of all hits in the same active + // volume into digit if (strstr(option, "print") ) { - Print("") ; + Print() ; return ; } if(strstr(option,"tim")) gBenchmark->Start("PHOSSDigitizer"); - - //Check, if this branch already exits - TObjArray * lob = (TObjArray*)gAlice->TreeS()->GetListOfBranches() ; - TIter next(lob) ; - TBranch * branch = 0 ; - Bool_t phosfound = kFALSE, sdigitizerfound = kFALSE ; - while ( (branch = (TBranch*)next()) && (!phosfound || !sdigitizerfound) ) { - if ( (strcmp(branch->GetName(), "PHOS")==0) && (strcmp(branch->GetTitle(), GetName())==0) ) - phosfound = kTRUE ; - - else if ( (strcmp(branch->GetName(), "AliPHOSSDigitizer")==0) && (strcmp(branch->GetTitle(), GetName())==0) ) - sdigitizerfound = kTRUE ; +/* + // 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")); - if ( phosfound || sdigitizerfound ) { - cerr << "WARNING: AliPHOSSDigitizer::Exec -> SDigits and/or SDigitizer branch with name " << GetName() - << " already exits" << endl ; - return ; - } - - - AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ; - - TClonesArray * hits = gime->Hits() ; + //switch off reloading of this task while getting event + if (!fInit) { // to prevent overwrite existing file + AliError( Form("Give a version name different from %s", fEventFolderName.Data()) ) ; + return ; + } - Int_t nevents = (Int_t) gAlice->TreeE()->GetEntries() ; + if (fLastEvent == -1) + fLastEvent = rl->GetNumberOfEvents() - 1 ; + else + fLastEvent = TMath::Min(fFirstEvent, rl->GetNumberOfEvents()); // only one event at the time + Int_t nEvents = fLastEvent - fFirstEvent + 1; - Int_t ievent ; - for(ievent = 0; ievent < nevents; ievent++){ - gAlice->GetEvent(ievent) ; - gAlice->SetEvent(ievent) ; - - if(gAlice->TreeH()==0){ - cerr << "ERROR: AliPHOSSDigitizer::Exec There is no Hit Tree" << endl; - return ; - } - - //set address of the hits - TBranch * branch = gAlice->TreeH()->GetBranch("PHOS"); - if (branch) - branch->SetAddress(&hits); - else{ - cerr << "ERROR: AliPHOSSDigitizer::Exec -> No branch PHOS in TreeH " << endl ; - cerr << " Do nothing " << endl ; - return ; + Int_t ievent, i; + + for (ievent = fFirstEvent; ievent <= fLastEvent; ievent++) { + rl->GetEvent(ievent) ; + TTree * treeS = phosLoader->TreeS(); + if(!treeS){ + phosLoader->MakeTree("S"); + treeS = phosLoader->TreeS(); } - TClonesArray * sdigits = gime->SDigits(GetName()) ; + 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 itrack ; - for (itrack=0; itrack < gAlice->GetNtrack(); itrack++){ - - //=========== Get the PHOS branch from Hits Tree for the Primary track itrack - branch->GetEntry(itrack,0); - Int_t i; - for ( i = 0 ; i < hits->GetEntries() ; i++ ) { - const AliPHOSHit * hit = gime->Hit(i) ; - // Assign primary number only if contribution is significant - - if( hit->GetEnergy() > fPrimThreshold) - new((*sdigits)[nSdigits]) AliPHOSDigit( hit->GetPrimary(), hit->GetId(), Digitize( hit->GetEnergy() ) ) ; - else - new((*sdigits)[nSdigits]) AliPHOSDigit( -1 , hit->GetId(), Digitize( hit->GetEnergy() ) ) ; - - nSdigits++ ; - - } - } // loop over tracks - + //Now make SDigits from hits, for PHOS it is the same, so just copy + 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++ ; + + } + sdigits->Sort() ; - + nSdigits = sdigits->GetEntriesFast() ; + + fSDigitsInRun += nSdigits ; sdigits->Expand(nSdigits) ; - Int_t i ; + for (i = 0 ; i < nSdigits ; i++) { - AliPHOSDigit * digit = (AliPHOSDigit *) sdigits->At(i) ; + AliPHOSDigit * digit = dynamic_cast(sdigits->At(i)) ; digit->SetIndexInList(i) ; } - if(gAlice->TreeS() == 0) - gAlice->MakeTree("S") ; - - //Make (if necessary) branches - char * file =0; - if(gSystem->Getenv("CONFIG_SPLIT_FILE")){ //generating file name - file = new char[strlen(gAlice->GetBaseFile())+20] ; - sprintf(file,"%s/PHOS.SDigits.root",gAlice->GetBaseFile()) ; - } - - TDirectory *cwd = gDirectory; +// // 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 + //First list of sdigits - Int_t bufferSize = 32000 ; - TBranch * sdigitsBranch = gAlice->TreeS()->Branch("PHOS",&sdigits,bufferSize); - sdigitsBranch->SetTitle(GetName()); - if (file) { - sdigitsBranch->SetFile(file); - TIter next( sdigitsBranch->GetListOfBranches()); - TBranch * subbr; - while ((subbr=(TBranch*)next())) { - subbr->SetFile(file); - } - cwd->cd(); - } - - //Next - SDigitizer - Int_t splitlevel = 0 ; - AliPHOSSDigitizer * sd = this ; - TBranch * sdigitizerBranch = gAlice->TreeS()->Branch("AliPHOSSDigitizer","AliPHOSSDigitizer", - &sd,bufferSize,splitlevel); - sdigitizerBranch->SetTitle(GetName()); - if (file) { - sdigitizerBranch->SetFile(file); - TIter next( sdigitizerBranch->GetListOfBranches()); - TBranch * subbr ; - while ((subbr=(TBranch*)next())) { - subbr->SetFile(file); - } - cwd->cd(); - delete file; - } + Int_t bufferSize = 32000 ; + TBranch * sdigitsBranch = treeS->Branch("PHOS",&sdigits,bufferSize); sdigitsBranch->Fill() ; - sdigitizerBranch->Fill() ; - gAlice->TreeS()->Write(0,TObject::kOverwrite) ; - + + phosLoader->WriteSDigits("OVERWRITE"); + phosLoader->WriteSDigitizer("OVERWRITE"); + if(strstr(option,"deb")) PrintSDigits(option) ; - - } - - if(strstr(option,"tim")){ - gBenchmark->Stop("PHOSSDigitizer"); - cout << "AliPHOSSDigitizer:" << endl ; - cout << " took " << gBenchmark->GetCpuTime("PHOSSDigitizer") << " seconds for SDigitizing " - << gBenchmark->GetCpuTime("PHOSSDigitizer")/nevents << " seconds per event " << endl ; - cout << endl ; - } - + + }// event loop -} -//__________________________________________________________________ -void AliPHOSSDigitizer::SetSDigitsBranch(const char * title ) -{ - // Setting title to branch SDigits +// //Write the quality assurance data +// GetQADataMaker()->EndOfCycle(AliQAv1::kHITS) ; +// GetQADataMaker()->EndOfCycle(AliQAv1::kSDIGITS) ; +// GetQADataMaker()->Finish() ; - TString stitle(title) ; + Unload(); - // check if branch with title already exists - TBranch * sdigitsBranch = (TBranch*)gAlice->TreeS()->GetListOfBranches()->FindObject("PHOS") ; - TBranch * sdigitizerBranch = (TBranch*)gAlice->TreeS()->GetListOfBranches()->FindObject("AliPHOSSDigitizer") ; - const char * sdigitsTitle = sdigitsBranch ->GetTitle() ; - const char * sdigitizerTitle = sdigitizerBranch ->GetTitle() ; - if ( stitle.CompareTo(sdigitsTitle)==0 || stitle.CompareTo(sdigitizerTitle)==0 ){ - cerr << "ERROR: AliPHOSSdigitizer::SetSDigitsBranch -> Cannot overwrite existing branch with title " << title << endl ; - return ; + 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) ; } - cout << "AliPHOSSdigitizer::SetSDigitsBranch -> Changing SDigits file from " << GetName() << " to " << title << endl ; - - SetName(title) ; - - // Post to the WhiteBoard - AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ; - gime->Post(GetTitle(), "S", GetName()) ; } //__________________________________________________________________ -void AliPHOSSDigitizer::Print(Option_t* option)const +void AliPHOSSDigitizer::Print(const Option_t *)const { // Prints parameters of SDigitizer - cout << "------------------- "<< GetName() << " -------------" << endl ; - cout << " Writing SDigitis to branch with title " << GetName() << endl ; - cout << " with digitization parameters A = " << fA << endl ; - cout << " B = " << fB << endl ; - cout << " Threshold for Primary assignment= " << fPrimThreshold << endl ; - cout << "---------------------------------------------------"<SDigits(GetName()) ; + 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",""); - cout << "AliPHOSSDigitizer: " << endl ; - cout << " Number of entries in SDigits list " << sdigits->GetEntriesFast() << endl ; - cout << endl ; + const TClonesArray * sdigits = phosLoader->SDigits() ; - if(strstr(option,"all")){// print all digits + Info( "\nPrintSDigits", "event # %d %d sdigits", + gAlice->GetEvNumber(), sdigits->GetEntriesFast() ) ; + + if(strstr(option,"all")||strstr(option,"EMC")){ //loop over digits AliPHOSDigit * digit; - cout << "SDigit Id " << " Amplitude " << " Index " << " Nprim " << " Primaries list " << endl; + printf("\nEMC sdigits\n") ; + Int_t maxEmc = geom->GetNModules() * geom->GetNCristalsInModule() ; Int_t index ; - for (index = 0 ; index < sdigits->GetEntries() ; index++) { - digit = (AliPHOSDigit * ) sdigits->At(index) ; - cout << setw(8) << digit->GetId() << " " << setw(3) << digit->GetAmp() << " " - << setw(6) << digit->GetIndexInList() << " " - << setw(5) << digit->GetNprimary() <<" "; - + 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", // 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++) - cout << setw(5) << digit->GetPrimary(iprimary+1) << " "; - cout << endl; - } + for (iprimary=0; iprimaryGetNprimary(); iprimary++) { + printf("%d ",digit->GetPrimary(iprimary+1) ) ; + } + } + } + + if(strstr(option,"all")||strstr(option,"CPV")){ + //loop over CPV digits + AliPHOSDigit * digit; + printf("\nCPV sdigits\n") ; + Int_t maxEmc = geom->GetNModules() * geom->GetNCristalsInModule() ; + Int_t index ; + for (index = 0 ; index < sdigits->GetEntriesFast(); index++) { + digit = dynamic_cast( sdigits->At(index) ) ; + if(digit->GetId() > maxEmc){ + printf("\n%6d %8d %4d %2d :", + digit->GetId(), digit->GetAmp(), digit->GetIndexInList(), digit->GetNprimary()) ; + Int_t iprimary; + for (iprimary=0; iprimaryGetNprimary(); iprimary++) { + printf("%d ",digit->GetPrimary(iprimary+1) ) ; + } + } + } } } //____________________________________________________________________________ -void AliPHOSSDigitizer::UseHitsFrom(const char * filename) +void AliPHOSSDigitizer::Unload() const { - SetTitle(filename) ; - Init() ; + // Unloads the objects from the folder + AliRunLoader* rl = AliRunLoader::GetRunLoader(fEventFolderName) ; + AliPHOSLoader * phosLoader = dynamic_cast(rl->GetLoader("PHOSLoader")); + phosLoader->UnloadHits() ; + phosLoader->UnloadSDigits() ; }