X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FAliPHOSSDigitizer.cxx;h=9f9d44a229efbc5fa983f499a515914f1d7185cd;hb=9225a7c4941479e2c8a11fb9908e19a869acec61;hp=6cd2cea2a933bf5513f454c2239ea851b2d6fd49;hpb=7acf600847d7c3706c34409560cce379580d954e;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/AliPHOSSDigitizer.cxx b/PHOS/AliPHOSSDigitizer.cxx index 6cd2cea2a93..9f9d44a229e 100644 --- a/PHOS/AliPHOSSDigitizer.cxx +++ b/PHOS/AliPHOSSDigitizer.cxx @@ -7,343 +7,445 @@ * Permission to use, copy, modify and distribute this software and its * * documentation strictly for non-commercial purposes is hereby granted * * without fee, provided that the above copyright notice appears in all * - * copies and that both the copyright notice and this permission notice * - * appear in the supporting documentation. The authors make no claims * + * copies and that both the copyright notice and this permission notice *ÄQual * appear in the supporting documentation. The authors make no claims * * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ + /* $Id$ */ +/* History of cvs commits: + * + * $Log$ + * 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) + * + * Revision 1.50 2006/08/28 10:01:56 kharlov + * Effective C++ warnings fixed (Timur Pocheptsov) + * + * Revision 1.49 2006/05/10 06:42:53 kharlov + * Remove redundant loop over primaries + * + * Revision 1.48 2006/04/22 10:30:17 hristov + * Add fEnergy to AliPHOSDigit and operate with EMC amplitude in energy units (Yu.Kharlov) + * + * Revision 1.47 2005/05/28 14:19:04 schutz + * Compilation warnings fixed by T.P. + * + */ + //_________________________________________________________________________ // 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 +// the created SDigits +// The title of the TTAsk is the name of the file that contains the hits from +// which the SDigits are created // A Summable Digits is the sum of all hits originating // from one primary in one active cell // A threshold for assignment of the primary to SDigit is applied // SDigits are written to TreeS, branch "PHOS" // AliPHOSSDigitizer with all current parameters is written // to TreeS branch "AliPHOSSDigitizer". -// Both branches, "PHOS" and "AliPHOSSDigitizer", are written to the same -// file, and therefore, changing branch file name one can produce several -// versions of SDigitization from the same hits. -// +// Both branches have the same title. If necessary one can produce +// another set of SDigits with different parameters. Two versions +// can be distunguished using titles of the branches. +// User case: +// root [0] AliPHOSSDigitizer * s = new AliPHOSSDigitizer("galice.root") +// Warning in : object already instantiated +// root [1] s->ExecuteTask() +// // Makes SDigitis for all events stored in galice.root +// root [2] s->SetPedestalParameter(0.001) +// // One can change parameters of digitization +// root [3] s->SetSDigitsBranch("Pedestal 0.001") +// // and write them into the new branch +// root [4] s->ExecuteTask("deb all tim") +// // available parameters: +// deb - print # of produced SDigitis +// deb all - print # and list of produced SDigits +// tim - print benchmarking information // //*-- 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" + //#include "TObjectTable.h" + // --- Standard library --- -#include // --- AliRoot header files --- -#include "AliRun.h" +#include "AliLog.h" +#include "AliPHOSGeometry.h" #include "AliPHOSDigit.h" +#include "AliPHOSGetter.h" #include "AliPHOSHit.h" #include "AliPHOSSDigitizer.h" - +#include "AliPHOSQADataMaker.h" ClassImp(AliPHOSSDigitizer) //____________________________________________________________________________ - AliPHOSSDigitizer::AliPHOSSDigitizer():TTask("AliPHOSSDigitizer","") +AliPHOSSDigitizer::AliPHOSSDigitizer() : + TTask("",""), + fA(0.f), fB(0.f), + fPrimThreshold(0.f), + fDefaultInit(kTRUE), + fEventFolderName(""), + fInit(kFALSE), + fSDigitsInRun(0), + fFirstEvent(0), + fLastEvent(0), + fQADM (0x0) { // ctor - fA = 0; - fB = 10000000. ; - fPrimThreshold = 0.01 ; - fNevents = 0 ; - fSDigits = 0 ; - fHits = 0 ; - fIsInitialized = kFALSE ; - + // Intialize the quality assurance data maker } //____________________________________________________________________________ -AliPHOSSDigitizer::AliPHOSSDigitizer(const char* headerFile, const char *sDigitsTitle):TTask("AliPHOSSDigitizer","") +AliPHOSSDigitizer::AliPHOSSDigitizer(const char * alirunFileName, + const char * eventFolderName): + TTask("PHOS"+AliConfig::Instance()->GetSDigitizerTaskName(), alirunFileName), + fA(0.f), fB(0.f), + fPrimThreshold(0.f), + fDefaultInit(kFALSE), + fEventFolderName(eventFolderName), + fInit(kFALSE), + fSDigitsInRun(0), + fFirstEvent(0), + fLastEvent(0), + fQADM (0x0) { // ctor - fA = 0; - fB = 10000000.; - fPrimThreshold = 0.01 ; - fNevents = 0 ; - fSDigitsTitle = sDigitsTitle ; - fHeadersFile = headerFile ; - fSDigits = new TClonesArray("AliPHOSDigit",1000); - fHits = new TClonesArray("AliPHOSHit",1000); + InitParameters() ; + Init(); + fDefaultInit = kFALSE ; + // Intialize the quality assurance data maker + //FIXME: get the run number + Int_t run = 0 ; + //EMXIF + GetQADataMaker()->Init(AliQA::kHITS, run, fgkCycles) ; + GetQADataMaker()->StartOfCycle(AliQA::kHITS) ; + GetQADataMaker()->Init(AliQA::kSDIGITS, run) ; + GetQADataMaker()->StartOfCycle(AliQA::kSDIGITS, "same") ; +} - TFile * file = (TFile*) gROOT->GetFile(fHeadersFile.Data() ) ; - - //File was not opened yet - if(file == 0){ - file = new TFile(fHeadersFile.Data(),"update") ; - gAlice = (AliRun *) file->Get("gAlice") ; - } - - //add Task to //root/Tasks folder - TTask * roottasks = (TTask*)gROOT->GetRootFolder()->FindObject("Tasks") ; - roottasks->Add(this) ; - - fIsInitialized = kTRUE ; +//____________________________________________________________________________ +AliPHOSSDigitizer::AliPHOSSDigitizer(const AliPHOSSDigitizer& sd) : + TTask(sd.GetName(), sd.GetTitle()), + fA(sd.fA), fB(sd.fB), + fPrimThreshold(sd.fPrimThreshold), + fDefaultInit(kFALSE), + fEventFolderName(sd.fEventFolderName), + fInit(kFALSE), + fSDigitsInRun(sd.fSDigitsInRun), + fFirstEvent(sd.fFirstEvent), + fLastEvent(sd.fLastEvent), + fQADM (sd.fQADM) +{ + // cpy ctor + // Intialize the quality assurance data maker + //FIXME: get the run number + Int_t run = 0 ; + //EMXIF + GetQADataMaker()->Init(AliQA::kHITS, run, fgkCycles) ; + GetQADataMaker()->StartOfCycle(AliQA::kHITS) ; + GetQADataMaker()->Init(AliQA::kSDIGITS, run) ; + GetQADataMaker()->StartOfCycle(AliQA::kSDIGITS, "same") ; } -//____________________________________________________________________________ -AliPHOSSDigitizer::~AliPHOSSDigitizer() +//_____________________________________________________________________________ +AliPHOSSDigitizer& AliPHOSSDigitizer::operator = (const AliPHOSSDigitizer& qa) { - // dtor - if(fSDigits) - delete fSDigits ; - if(fHits) - delete fHits ; +// assignment operator + + this->~AliPHOSSDigitizer(); + new(this) AliPHOSSDigitizer(qa); + return *this; } + //____________________________________________________________________________ -void AliPHOSSDigitizer::Init(){ - //Initialization can not be done in the default constructor +AliPHOSSDigitizer::~AliPHOSSDigitizer() { + //dtor + // AliPHOSGetter * gime = + // AliPHOSGetter::Instance(GetTitle(),fEventFolderName.Data()); + AliPHOSGetter * gime = + AliPHOSGetter::Instance(); + gime->PhosLoader()->CleanSDigitizer(); + delete fQADM ; +} - if(!fIsInitialized){ +//____________________________________________________________________________ +void AliPHOSSDigitizer::Init() +{ + // Uses the getter 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 ; + } - if(fHeadersFile.IsNull()) - fHeadersFile="galice.root" ; + gime->PostSDigitizer(this); + gime->PhosLoader()->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE); + + fQADM = new AliPHOSQADataMaker() ; - TFile * file = (TFile*) gROOT->GetFile(fHeadersFile.Data() ) ; - - //if file was not opened yet, read gAlice - if(file == 0){ - file = new TFile(fHeadersFile.Data(),"update") ; - gAlice = (AliRun *) file->Get("gAlice") ; - } - - fHits = new TClonesArray("AliPHOSHit",1000); - fSDigits = new TClonesArray("AliPHOSDigit",1000); - - // add Task to //root/Tasks folder - TTask * roottasks = (TTask*)gROOT->GetRootFolder()->FindObject("Tasks") ; - roottasks->Add(this) ; - - fIsInitialized = kTRUE ; - } } + +//____________________________________________________________________________ +void AliPHOSSDigitizer::InitParameters() +{ + // initializes the parameters for digitization + fA = 0; + fB = 10000000.; + fPrimThreshold = 0.01 ; + fSDigitsInRun = 0 ; +} + //____________________________________________________________________________ -void AliPHOSSDigitizer::Exec(Option_t *option) { - //Collects all hits in the same active volume into digit +void AliPHOSSDigitizer::Exec(Option_t *option) +{ + // 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(!fIsInitialized) - Init() ; + if (strstr(option, "print") ) { + Print() ; + return ; + } if(strstr(option,"tim")) gBenchmark->Start("PHOSSDigitizer"); - fNevents = (Int_t) gAlice->TreeE()->GetEntries() ; + // check the QA result for RAWS + AliQA * qa = AliQA::Instance(AliQA::kPHOS) ; + if ( qa->IsSet(AliQA::kPHOS, AliQA::kRAW, AliQA::kFATAL)) { + AliFatal("QA status in RAW was Fatal") ; + } else if ( qa->IsSet(AliQA::kPHOS, AliQA::kRAW, AliQA::kERROR)) { + AliError("QA status in RAW was Error") ; + } else if ( qa->IsSet(AliQA::kPHOS, AliQA::kRAW, AliQA::kWARNING) ) { + AliWarning("QA status in RAW was Warning") ; + } else if ( qa->IsSet(AliQA::kPHOS, AliQA::kRAW, AliQA::kINFO) ) { + AliInfo("QA status in RAW was Info") ; + } + + AliPHOSGetter * gime = AliPHOSGetter::Instance() ; + + //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 ; + } + + if (fLastEvent == -1) + fLastEvent = gime->MaxEvent() - 1 ; + else + fLastEvent = TMath::Min(fFirstEvent, gime->MaxEvent()); // only ine event at the time + Int_t nEvents = fLastEvent - fFirstEvent + 1; - Int_t ievent ; - for(ievent = 0; ievent < fNevents; ievent++){ - gAlice->GetEvent(ievent) ; - gAlice->SetEvent(ievent) ; - - if(gAlice->TreeH()==0){ - cout << "AliPHOSSDigitizer: There is no Hit Tree" << endl; - return ; - } - - //set address of the hits - TBranch * branch = gAlice->TreeH()->GetBranch("PHOS"); - if (branch) - branch->SetAddress(&fHits); - else{ - cout << "ERROR in AliPHOSSDigitizer: "<< endl ; - cout << " no branch PHOS in TreeH"<< endl ; - cout << " do nothing " << endl ; - return ; - } - - fSDigits->Clear(); + Int_t ievent, i; + + //AliMemoryWatcher memwatcher; + + for (ievent = fFirstEvent; ievent <= fLastEvent; ievent++) { + gime->Event(ievent,"H") ; + TTree * treeS = gime->TreeS(); + TClonesArray * hits = gime->Hits() ; + TClonesArray * sdigits = gime->SDigits() ; + sdigits->Clear(); Int_t nSdigits = 0 ; - - - //Now made SDigits from hits, for PHOS it is the same, so just copy - Int_t itrack ; - for (itrack=0; itrackGetNtrack(); itrack++){ - - //=========== Get the Hits Tree for the Primary track itrack - gAlice->ResetHits(); - gAlice->TreeH()->GetEvent(itrack); + //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 - Int_t i; - for ( i = 0 ; i < fHits->GetEntries() ; i++ ) { - AliPHOSHit * hit = (AliPHOSHit*)fHits->At(i) ; - AliPHOSDigit * newdigit ; - - // Assign primary number only if contribution is significant - if( hit->GetEnergy() > fPrimThreshold) - newdigit = new AliPHOSDigit( hit->GetPrimary(), hit->GetId(), Digitize( hit->GetEnergy() ) ) ; - else - newdigit = new AliPHOSDigit( -1 , hit->GetId(), Digitize( hit->GetEnergy() ) ) ; - - new((*fSDigits)[nSdigits]) AliPHOSDigit(* newdigit) ; - nSdigits++ ; - - delete newdigit ; - } + 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 tracks - - fSDigits->Sort() ; - - nSdigits = fSDigits->GetEntriesFast() ; - fSDigits->Expand(nSdigits) ; - - Int_t i ; + } + + sdigits->Sort() ; + + nSdigits = sdigits->GetEntriesFast() ; + + fSDigitsInRun += nSdigits ; + sdigits->Expand(nSdigits) ; + for (i = 0 ; i < nSdigits ; i++) { - AliPHOSDigit * digit = (AliPHOSDigit *) fSDigits->At(i) ; + AliPHOSDigit * digit = dynamic_cast(sdigits->At(i)) ; digit->SetIndexInList(i) ; } - if(gAlice->TreeS() == 0) - gAlice->MakeTree("S") ; - - //check, if this branch already exits? - TBranch * sdigitsBranch = 0; - TBranch * sdigitizerBranch = 0; - - TObjArray * branches = gAlice->TreeS()->GetListOfBranches() ; - Int_t ibranch; - Bool_t phosNotFound = kTRUE ; - Bool_t sdigitizerNotFound = kTRUE ; - - for(ibranch = 0;ibranch GetEntries();ibranch++){ - - if(phosNotFound){ - sdigitsBranch=(TBranch *) branches->At(ibranch) ; - if( (strcmp("PHOS",sdigitsBranch->GetName())==0 ) && - (fSDigitsTitle.CompareTo(sdigitsBranch->GetTitle()) == 0) ) - phosNotFound = kFALSE ; - } - if(sdigitizerNotFound){ - sdigitizerBranch = (TBranch *) branches->At(ibranch) ; - if( (strcmp(sdigitizerBranch->GetName(),"AliPHOSSDigitizer") == 0)&& - (fSDigitsTitle.CompareTo(sdigitizerBranch->GetTitle()) == 0) ) - sdigitizerNotFound = kFALSE ; - } - } + // make Quality Assurance data + + if (GetQADataMaker()->IsCycleDone() ) { + GetQADataMaker()->EndOfCycle(AliQA::kHITS) ; + GetQADataMaker()->EndOfCycle(AliQA::kSDIGITS) ; + GetQADataMaker()->StartOfCycle(AliQA::kHITS) ; + GetQADataMaker()->StartOfCycle(AliQA::kSDIGITS, "same") ; + } + GetQADataMaker()->Exec(AliQA::kHITS, hits) ; + GetQADataMaker()->Exec(AliQA::kSDIGITS, sdigits) ; + GetQADataMaker()->Increment() ; + + //Now write SDigits - if(!(sdigitizerNotFound && phosNotFound)){ - cout << "AliPHOSSdigitizer error:" << endl ; - cout << "Can not overwrite existing branches: do not write" << endl ; - return ; - } - - //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; //First list of sdigits - Int_t bufferSize = 32000 ; - sdigitsBranch = gAlice->TreeS()->Branch("PHOS",&fSDigits,bufferSize); - sdigitsBranch->SetTitle(fSDigitsTitle.Data()); - if (file) { - sdigitsBranch->SetFile(file); - TIter next( sdigitsBranch->GetListOfBranches()); - while ((sdigitsBranch=(TBranch*)next())) { - sdigitsBranch->SetFile(file); - } - cwd->cd(); - } - - //second - SDigitizer - Int_t splitlevel = 0 ; - AliPHOSSDigitizer * sd = this ; - sdigitizerBranch = gAlice->TreeS()->Branch("AliPHOSSDigitizer","AliPHOSSDigitizer", - &sd,bufferSize,splitlevel); - sdigitizerBranch->SetTitle(fSDigitsTitle.Data()); - if (file) { - sdigitizerBranch->SetFile(file); - TIter next( sdigitizerBranch->GetListOfBranches()); - while ((sdigitizerBranch=(TBranch*)next())) { - sdigitizerBranch->SetFile(file); - } - cwd->cd(); - delete file; - } - gAlice->TreeS()->Fill() ; - gAlice->TreeS()->Write(0,TObject::kOverwrite) ; - + Int_t bufferSize = 32000 ; + TBranch * sdigitsBranch = treeS->Branch("PHOS",&sdigits,bufferSize); + + sdigitsBranch->Fill() ; + + gime->WriteSDigits("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(AliQA::kHITS) ; + GetQADataMaker()->EndOfCycle(AliQA::kSDIGITS) ; + GetQADataMaker()->Finish(AliQA::kHITS) ; + GetQADataMaker()->Finish(AliQA::kSDIGITS) ; + + Unload(); + + // gime->PhosLoader()->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE); if(strstr(option,"tim")){ gBenchmark->Stop("PHOSSDigitizer"); - cout << "AliPHOSSDigitizer:" << endl ; - cout << " took " << gBenchmark->GetCpuTime("PHOSSDigitizer") << " seconds for SDigitizing " - << gBenchmark->GetCpuTime("PHOSSDigitizer")/fNevents << " seconds per event " << endl ; - cout << endl ; + Info("Exec"," took %f seconds for SDigitizing %f seconds per event", + gBenchmark->GetCpuTime("PHOSSDigitizer"), gBenchmark->GetCpuTime("PHOSSDigitizer")/nEvents) ; } - -} -//__________________________________________________________________ -void AliPHOSSDigitizer::SetSDigitsBranch(const char * title ){ - //Seting title to branch SDigits - if(!fSDigitsTitle.IsNull()) - cout << "AliPHOSSdigitizer: changing SDigits file from " <GetEntriesFast() << endl ; - cout << endl ; +void AliPHOSSDigitizer::PrintSDigits(Option_t * option) +{ + // Prints list of digits produced in the current pass of AliPHOSDigitizer + + + AliPHOSGetter * gime = AliPHOSGetter::Instance() ; + const TClonesArray * sdigits = gime->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 = gime->PHOSGeometry()->GetNModules()*gime->PHOSGeometry()->GetNCristalsInModule() ; Int_t index ; - for (index = 0 ; index < fSDigits->GetEntries() ; index++) { - digit = (AliPHOSDigit * ) fSDigits->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 = gime->PHOSGeometry()->GetNModules()*gime->PHOSGeometry()->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::Unload() const +{ + // Unloads the objects from the folder + AliPHOSGetter * gime = AliPHOSGetter::Instance() ; + AliPHOSLoader * loader = gime->PhosLoader() ; + loader->UnloadHits() ; + loader->UnloadSDigits() ; +}