X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=EMCAL%2FAliEMCALSDigitizer.cxx;h=1ba3a8fc279ce7d05ef1e8f5bb061bcfff105ef5;hb=6b4892386b6faf4675cc086e591bd800c1e8c5c4;hp=f1d88b262d9b7e2fd3deb60b7fe1de04ec1a5d69;hpb=70479d0eaf6647221605280c1c435802bb47bfb8;p=u%2Fmrichter%2FAliRoot.git diff --git a/EMCAL/AliEMCALSDigitizer.cxx b/EMCAL/AliEMCALSDigitizer.cxx index f1d88b262d9..1ba3a8fc279 100644 --- a/EMCAL/AliEMCALSDigitizer.cxx +++ b/EMCAL/AliEMCALSDigitizer.cxx @@ -20,6 +20,13 @@ // A Summable Digits is the sum of all hits originating // from one in one tower of the EMCAL // A threshold for assignment of the primary to SDigit is applied +// +// JLK 26-Jun-2008 Added explanation: +// SDigits need to hold the energy sum of the hits, but AliEMCALDigit +// can (should) only store amplitude. Therefore, the SDigit energy is +// "digitized" before being stored and must be "calibrated" back to an +// energy before SDigits are summed to form true Digits +// // SDigits are written to TreeS, branch "EMCAL" // AliEMCALSDigitizer with all current parameters is written // to TreeS branch "AliEMCALSDigitizer". @@ -43,59 +50,91 @@ // //*-- Author : Sahal Yacoob (LBL) // based on : AliPHOSSDigitzer -// Modif: -// August 2002 Yves Schutz: clone PHOS as closely as possible and intoduction -// of new IO (à la PHOS) ////////////////////////////////////////////////////////////////////////////// // --- ROOT system --- -#include "TFile.h" -#include "TTask.h" -#include "TTree.h" -#include "TSystem.h" -#include "TROOT.h" -#include "TFolder.h" -#include "TBenchmark.h" -#include "TGeometry.h" +#include +#include +//#include +#include +#include // --- Standard library --- -#include +#include "stdlib.h" // --- AliRoot header files --- -#include "AliRun.h" -#include "AliHeader.h" +#include "AliLog.h" +#include "AliRunLoader.h" +#include "AliStack.h" #include "AliEMCALDigit.h" -#include "AliEMCALGeometry.h" -#include "AliEMCALGetter.h" +#include "AliEMCALLoader.h" #include "AliEMCALHit.h" #include "AliEMCALSDigitizer.h" +#include "AliEMCALGeometry.h" +#include "AliEMCALSimParam.h" ClassImp(AliEMCALSDigitizer) //____________________________________________________________________________ - AliEMCALSDigitizer::AliEMCALSDigitizer():TTask("AliEMCALSDigitizer","") +AliEMCALSDigitizer::AliEMCALSDigitizer() + : TTask("",""), + fA(0.),fB(0.),fECPrimThreshold(0.), + fDefaultInit(kTRUE), + fEventFolderName(0), + fInit(0), + fSDigitsInRun(0), + fFirstEvent(0), + fLastEvent(0), + fSampling(0.) { // ctor - InitParameters() ; - fDefaultInit = kTRUE ; + InitParameters(); } //____________________________________________________________________________ -AliEMCALSDigitizer::AliEMCALSDigitizer(const char* headerFile, const char *sDigitsTitle, const Bool_t toSplit): -TTask(sDigitsTitle, headerFile) +AliEMCALSDigitizer::AliEMCALSDigitizer(const char * alirunFileName, + const char * eventFolderName) + : TTask("EMCAL"+AliConfig::Instance()->GetSDigitizerTaskName(), alirunFileName), + fA(0.),fB(0.),fECPrimThreshold(0.), + fDefaultInit(kFALSE), + fEventFolderName(eventFolderName), + fInit(0), + fSDigitsInRun(0), + fFirstEvent(0), + fLastEvent(0), + fSampling(0.) { // ctor - InitParameters() ; - fToSplit = toSplit ; Init(); - fDefaultInit = kFALSE ; + InitParameters() ; + } + //____________________________________________________________________________ -AliEMCALSDigitizer::~AliEMCALSDigitizer() +AliEMCALSDigitizer::AliEMCALSDigitizer(const AliEMCALSDigitizer & sd) + : TTask(sd.GetName(),sd.GetTitle()), + fA(sd.fA), + fB(sd.fB), + fECPrimThreshold(sd.fECPrimThreshold), + fDefaultInit(sd.fDefaultInit), + fEventFolderName(sd.fEventFolderName), + fInit(sd.fInit), + fSDigitsInRun(sd.fSDigitsInRun), + fFirstEvent(sd.fFirstEvent), + fLastEvent(sd.fLastEvent), + fSampling(sd.fSampling) { - // dtor - fSplitFile = 0 ; + //cpy ctor +} + + +//____________________________________________________________________________ +AliEMCALSDigitizer::~AliEMCALSDigitizer() { + //dtor + AliLoader *emcalLoader=0; + if ((emcalLoader = AliRunLoader::Instance()->GetDetectorLoader("EMCAL"))) + emcalLoader->CleanSDigitizer(); } //____________________________________________________________________________ @@ -107,290 +146,265 @@ void AliEMCALSDigitizer::Init(){ //============================================================= YS // The initialisation is now done by the getter - if( strcmp(GetTitle(), "") == 0 ) - SetTitle("galice.root") ; + fInit = kTRUE; - AliEMCALGetter * gime = AliEMCALGetter::GetInstance(GetTitle(), GetName(), fToSplit) ; - if ( gime == 0 ) { - cerr << "ERROR: AliEMCALSDigitizer::Init -> Could not obtain the Getter object !" - << endl ; + AliEMCALLoader *emcalLoader = dynamic_cast(AliRunLoader::Instance()->GetDetectorLoader("EMCAL")); + + if ( emcalLoader == 0 ) { + Fatal("Init", "Could not obtain the AliEMCALLoader"); return ; } - gime->PostSDigits( GetName(), GetTitle() ) ; + emcalLoader->PostSDigitizer(this); + emcalLoader->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE); - fSplitFile = 0 ; - if(fToSplit){ - // construct the name of the file as /path/EMCAL.SDigits.root - // First - extract full path if necessary - TString sDigitsFileName(GetTitle()) ; - Ssiz_t islash = sDigitsFileName.Last('/') ; - if(islash(gROOT->GetFile(sDigitsFileName.Data())); - if(!fSplitFile) - fSplitFile = TFile::Open(sDigitsFileName.Data(),"update") ; - } - - TString sdname(GetName() ); - sdname.Append(":") ; - sdname.Append(GetTitle() ) ; - SetName(sdname) ; - gime->PostSDigitizer(this) ; } //____________________________________________________________________________ void AliEMCALSDigitizer::InitParameters() { - fA = 0; - fB = 10000000.; - fTowerPrimThreshold = 0.01 ; - fPreShowerPrimThreshold = 0.0001 ; - fPhotonElectronFactor = 5000. ; // photoelectrons per GeV - fSplitFile = 0 ; - fToSplit = kFALSE ; + //initialize parameters for sdigitization + + const AliEMCALGeometry * geom = AliEMCALGeometry::GetInstance(); + if (geom->GetSampling() == 0.) { + Fatal("InitParameters", "Sampling factor not set !") ; + } + + // + //JLK 26-Jun-2008 THIS SHOULD HAVE BEEN EXPLAINED AGES AGO: + // + //In order to be able to store SDigit Energy info into + //AliEMCALDigit, we need to convert it temporarily to an ADC amplitude + //and later when summing SDigits to form digits, convert it back to + //energy. These fA and fB parameters accomplish this through the + //Digitize() and Calibrate() methods + // + // Initializes parameters + fA = AliEMCALSimParam::GetInstance()->GetA(); //0; + fB = AliEMCALSimParam::GetInstance()->GetB(); //1.e+6; // Changed 24 Apr 2007. Dynamic range now 2 TeV + fSampling = geom->GetSampling(); + + // threshold for deposit energy of hit + fECPrimThreshold = AliEMCALSimParam::GetInstance()->GetECPrimaryThreshold();//0.05;// GeV // 22-may-07 was 0// 24-nov-04 - was 1.e-6; + + AliDebug(2,Form("Print: \n------------------- %s -------------\n",GetName())); + AliDebug(2,Form(" fInit %i\n", int(fInit))); + AliDebug(2,Form(" fFirstEvent %i\n", fFirstEvent)); + AliDebug(2,Form(" fLastEvent %i\n", fLastEvent)); + AliDebug(2,Form(" Writing SDigits to branch with title %s\n", fEventFolderName.Data())); + AliDebug(2,Form(" with digitization parameters A = %f\n", fA)); + AliDebug(2,Form(" B = %f\n", fB)); + AliDebug(2,Form(" Threshold for EC Primary assignment = %f\n", fECPrimThreshold)); + AliDebug(2,Form(" Sampling = %f\n", fSampling)); + AliDebug(2,Form("---------------------------------------------------\n")); + } //____________________________________________________________________________ void AliEMCALSDigitizer::Exec(Option_t *option) { - // Collects all hits in the same active volume into digit - - if( strcmp(GetName(), "") == 0 ) - Init() ; - + // Collects all hit of the same tower into digits + TString o(option); o.ToUpper(); if (strstr(option, "print") ) { - Print("") ; + + AliDebug(2,Form("Print: \n------------------- %s -------------\n",GetName())); + AliDebug(2,Form(" fInit %i\n", int(fInit))); + AliDebug(2,Form(" fFirstEvent %i\n", fFirstEvent)); + AliDebug(2,Form(" fLastEvent %i\n", fLastEvent)); + AliDebug(2,Form(" Writing SDigits to branch with title %s\n", fEventFolderName.Data())); + AliDebug(2,Form(" with digitization parameters A = %f\n", fA)); + AliDebug(2,Form(" B = %f\n", fB)); + AliDebug(2,Form(" Threshold for EC Primary assignment = %f\n", fECPrimThreshold)); + AliDebug(2,Form(" Sampling = %f\n", fSampling)); + AliDebug(2,Form("---------------------------------------------------\n")); + return ; } if(strstr(option,"tim")) gBenchmark->Start("EMCALSDigitizer"); - //Check, if this branch already exits - AliEMCALGetter * gime = AliEMCALGetter::GetInstance() ; - if(gime->BranchExists("SDigits") ) - return; + AliRunLoader *rl = AliRunLoader::Instance(); + AliEMCALLoader *emcalLoader = dynamic_cast(rl->GetDetectorLoader("EMCAL")); - TString sdname(GetName()) ; - sdname.Remove(sdname.Index(GetTitle())-1) ; + //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 = gime->MaxEvent() ; - Int_t ievent ; - - for(ievent = 0; ievent < nevents; ievent++){ - gime->Event(ievent,"H") ; - - const TClonesArray * hits = gime->Hits() ; - - TClonesArray * sdigits = gime->SDigits(sdname.Data()) ; + if (fLastEvent == -1) + fLastEvent = rl->GetNumberOfEvents() - 1 ; + else { + fLastEvent = TMath::Min(fLastEvent, rl->GetNumberOfEvents()-1); + } + Int_t nEvents = fLastEvent - fFirstEvent + 1; + + Int_t ievent; + Float_t energy=0.; // de * fSampling - 23-nov-04 + rl->LoadKinematics(); + rl->LoadHits("EMCAL"); + + for (ievent = fFirstEvent; ievent <= fLastEvent; ievent++) { + rl->GetEvent(ievent); + TTree * treeS = emcalLoader->TreeS(); + if ( !treeS ) { + emcalLoader->MakeSDigitsContainer(); + treeS = emcalLoader->TreeS(); + } + TClonesArray * hits = emcalLoader->Hits() ; + TClonesArray * sdigits = emcalLoader->SDigits() ; sdigits->Clear(); - Int_t nSdigits = 0 ; - //Collects all hits in the same active volume into digit - - //Now make SDigits from hits, for EMCAL it is the same, so just copy - Int_t nPrim = static_cast((gAlice->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 EMCAL branch from Hits Tree for the Primary iprim - gime->Track(iprim) ; - Int_t i; - for ( i = 0 ; i < hits->GetEntries() ; i++ ) { - AliEMCALHit * hit = dynamic_cast(hits->At(i)) ; - AliEMCALDigit * curSDigit = 0 ; - AliEMCALDigit * sdigit = 0 ; - Bool_t newsdigit = kTRUE; + Int_t nSdigits = 0 ; + Int_t i; + AliEMCALGeometry *geom = AliEMCALGeometry::GetInstance(); + for ( i = 0 ; i < hits->GetEntries() ; i++ ) { + AliEMCALHit * hit = dynamic_cast(hits->At(i)) ; + AliEMCALDigit * curSDigit = 0 ; + AliEMCALDigit * sdigit = 0 ; + Bool_t newsdigit = kTRUE; + + // hit->GetId() - Absolute Id number EMCAL segment + if(geom->CheckAbsCellId(hit->GetId())) { // was IsInECA(hit->GetId()) + energy = hit->GetEnergy() * fSampling; // 23-nov-04 + if(energy > fECPrimThreshold ) // Assign primary number only if deposited energy is significant - - if( (!hit->IsInPreShower() && hit->GetEnergy() > fTowerPrimThreshold) || - (hit->IsInPreShower() && hit->GetEnergy() > fPreShowerPrimThreshold)) - curSDigit = new AliEMCALDigit( hit->GetPrimary(), - hit->GetIparent(),Layer2TowerID(hit->GetId(),hit->IsInPreShower()), - Digitize(hit->GetEnergy()), hit->GetTime() ) ; - else + curSDigit = new AliEMCALDigit( hit->GetPrimary(), + hit->GetIparent(), hit->GetId(), + Digitize(energy), hit->GetTime(), + -1, energy ) ; + else curSDigit = new AliEMCALDigit( -1 , -1 , - Layer2TowerID(hit->GetId(),hit->IsInPreShower()), - Digitize(hit->GetEnergy()), hit->GetTime() ) ; - Int_t check = 0 ; - for(check= 0; check < nSdigits ; check++) { - sdigit = dynamic_cast(sdigits->At(check)) ; - if( sdigit->GetId() == curSDigit->GetId()) { // Are we in the same tower or the same preshower ? - *sdigit = *sdigit + *curSDigit; - newsdigit = kFALSE; - } - } - if (newsdigit) { - new((*sdigits)[nSdigits]) AliEMCALDigit(*curSDigit); - nSdigits++ ; - } - delete curSDigit ; - } // loop over all hits (hit = deposited energy/layer/entering particle) - } // loop over iprim - - sdigits->Sort() ; - - nSdigits = sdigits->GetEntriesFast() ; - fSDigitsInRun += nSdigits ; - sdigits->Expand(nSdigits) ; - - const AliEMCALGeometry * geom = gime->EMCALGeometry() ; + hit->GetId(), + Digitize(energy), hit->GetTime(), + -1, energy ) ; + } else { + Warning("Exec"," abs id %i is bad \n", hit->GetId()); + newsdigit = kFALSE; + curSDigit = 0; + } - if (nSdigits != 0 ) { - Int_t lastPreShowerIndex = nSdigits - 1 ; - - - if (!(dynamic_cast(sdigits->At(lastPreShowerIndex))->IsInPreShower())) + if(curSDigit != 0){ + for(Int_t check= 0; check < nSdigits ; check++) { + sdigit = dynamic_cast(sdigits->At(check)) ; - lastPreShowerIndex = -2; - - Int_t firstPreShowerIndex = 100000 ; - Int_t index ; - AliEMCALDigit * sdigit = 0 ; - - for ( index = 0; index < nSdigits ; index++) { - sdigit = dynamic_cast(sdigits->At(index) ) ; - if (sdigit->IsInPreShower() ){ - firstPreShowerIndex = index ; - break ; - } - } - - AliEMCALDigit * preshower ; - AliEMCALDigit * tower ; - Int_t lastIndex = lastPreShowerIndex +1 ; - - for (index = firstPreShowerIndex ; index <= lastPreShowerIndex; index++) { - preshower = dynamic_cast(sdigits->At(index) ); - Bool_t towerFound = kFALSE ; - Int_t jndex ; - for (jndex = 0; jndex < firstPreShowerIndex; jndex++) { - tower = dynamic_cast(sdigits->At(jndex) ); - if ( (preshower->GetId() - (geom->GetNZ() * geom->GetNPhi()) ) == tower->GetId() ) { - Float_t towerEnergy = static_cast(tower->GetAmp()) ; - Float_t preshoEnergy = static_cast(preshower->GetAmp()) ; - towerEnergy +=preshoEnergy ; - *tower = *tower + *preshower ; // and add preshower multiplied by layer ratio to tower - tower->SetAmp(static_cast(TMath::Ceil(towerEnergy))) ; - towerFound = kTRUE ; - } - } - if ( !towerFound ) { - new((*sdigits)[lastIndex]) AliEMCALDigit(*preshower); - AliEMCALDigit * temp = dynamic_cast(sdigits->At(lastIndex)) ; - temp->SetId(temp->GetId() - (geom->GetNZ() * geom->GetNPhi()) ) ; - lastIndex++ ; + if( sdigit->GetId() == curSDigit->GetId()) { // Are we in the same ECAL tower ? + *sdigit = *sdigit + *curSDigit; + newsdigit = kFALSE; } } - sdigits->Sort() ; - Int_t NPrimarymax = -1 ; - Int_t i ; - for (i = 0 ; i < sdigits->GetEntriesFast() ; i++) { - sdigit = dynamic_cast(sdigits->At(i)) ; - sdigit->SetIndexInList(i) ; - } - - for (i = 0 ; i < sdigits->GetEntriesFast() ; i++) { - if (((dynamic_cast(sdigits->At(i)))->GetNprimary()) > NPrimarymax) - NPrimarymax = ((dynamic_cast(sdigits->At(i)))->GetNprimary()) ; - } } - - //Now write SDigits - - if(gAlice->TreeS() == 0 || (fSplitFile)) //<--- To be checked: we should not create TreeS if it is already here - gAlice->MakeTree("S",fSplitFile); - - if(fSplitFile) - fSplitFile->cd() ; - - //First list of sdigits - Int_t bufferSize = 32000 ; - TBranch * sdigitsBranch = gAlice->TreeS()->Branch("EMCAL",&sdigits,bufferSize); - sdigitsBranch->SetTitle(sdname); - - //NEXT - SDigitizer - Int_t splitlevel = 0 ; - AliEMCALSDigitizer * sd = this ; - TBranch * sdigitizerBranch = gAlice->TreeS()->Branch("AliEMCALSDigitizer","AliEMCALSDigitizer", - &sd,bufferSize,splitlevel); - sdigitizerBranch->SetTitle(sdname); - - sdigitsBranch->Fill() ; - sdigitizerBranch->Fill() ; - gAlice->TreeS()->AutoSave() ; + if (newsdigit) { + new((*sdigits)[nSdigits]) AliEMCALDigit(*curSDigit); + nSdigits++ ; + } + delete curSDigit ; + } // loop over all hits (hit = deposited energy/entering particle) + sdigits->Sort() ; + + nSdigits = sdigits->GetEntriesFast() ; + fSDigitsInRun += nSdigits ; - if(strstr(option,"deb")) - PrintSDigits(option) ; - + for (i = 0 ; i < sdigits->GetEntriesFast() ; i++) { + AliEMCALDigit * sdigit = dynamic_cast(sdigits->At(i)) ; + sdigit->SetIndexInList(i) ; + } + + // Now write SDigits + + Int_t bufferSize = 32000 ; + TBranch * sdigitsBranch = treeS->GetBranch("EMCAL"); + if (sdigitsBranch) + sdigitsBranch->SetAddress(&sdigits); + else + treeS->Branch("EMCAL",&sdigits,bufferSize); + + treeS->Fill(); + + emcalLoader->WriteSDigits("OVERWRITE"); + + //NEXT - SDigitizer + emcalLoader->WriteSDigitizer("OVERWRITE"); // why in event cycle ? + + if(strstr(option,"deb")) + PrintSDigits(option) ; } - + + Unload(); + + emcalLoader->GetSDigitsDataLoader()->GetBaseTaskLoader()->SetDoNotReload(kTRUE); + if(strstr(option,"tim")){ - gBenchmark->Stop("EMCALSDigitizer"); - cout << "AliEMCALSDigitizer:" << endl ; - cout << " took " << gBenchmark->GetCpuTime("EMCALSDigitizer") << " seconds for SDigitizing " - << gBenchmark->GetCpuTime("EMCALSDigitizer") << " seconds per event " << endl ; - cout << endl ; - } + gBenchmark->Stop("EMCALSDigitizer"); + printf("\n Exec: took %f seconds for SDigitizing %f seconds per event\n", + gBenchmark->GetCpuTime("EMCALSDigitizer"), gBenchmark->GetCpuTime("EMCALSDigitizer")/nEvents ) ; + } } //__________________________________________________________________ -void AliEMCALSDigitizer::SetSDigitsBranch(const char * title ){ - - // Setting title to branch SDigits - - TString stitle(title) ; - - // check if branch with title already exists - TBranch * sdigitsBranch = - static_cast(gAlice->TreeS()->GetListOfBranches()->FindObject("EMCAL")) ; - TBranch * sdigitizerBranch = - static_cast(gAlice->TreeS()->GetListOfBranches()->FindObject("AliEMCALSDigitizer")) ; - const char * sdigitsTitle = sdigitsBranch ->GetTitle() ; - const char * sdigitizerTitle = sdigitizerBranch ->GetTitle() ; - if ( stitle.CompareTo(sdigitsTitle)==0 || stitle.CompareTo(sdigitizerTitle)==0 ){ - cerr << "ERROR: AliEMCALSdigitizer::SetSDigitsBranch -> Cannot overwrite existing branch with title " << title << endl ; - return ; +Int_t AliEMCALSDigitizer::Digitize(Float_t energy)const { + // Digitize the energy + // + //JLK 26-Jun-2008 EXPLANATION LONG OVERDUE: + // + //We have to digitize the SDigit energy so that it can be stored in + //AliEMCALDigit, which has only an ADC amplitude member and + //(rightly) no energy member. This method converts the energy to an + //integer which can be re-converted back to an energy with the + //Calibrate(energy) method when it is time to create Digits from SDigits + // + Double_t aSignal = fA + energy*fB; + if (TMath::Abs(aSignal)>2147483647.0) { + //PH 2147483647 is the max. integer + //PH This apparently is a problem which needs investigation + AliWarning(Form("Too big or too small energy %f",aSignal)); + aSignal = TMath::Sign((Double_t)2147483647,aSignal); } - - cout << "AliEMCALSdigitizer::SetSDigitsBranch -> Changing SDigits file from " << GetName() << " to " << title << endl ; + return (Int_t ) aSignal; +} + +//__________________________________________________________________ +Float_t AliEMCALSDigitizer::Calibrate(Int_t amp)const { + // + // Convert the amplitude back to energy in GeV + // + //JLK 26-Jun-2008 EXPLANATION LONG OVERDUE: + // + //We have to digitize the SDigit energy with the method Digitize() + //so that it can be stored in AliEMCALDigit, which has only an ADC + //amplitude member and (rightly) no energy member. This method is + //just the reverse of Digitize(): it converts the stored amplitude + //back to an energy value in GeV so that the SDigit energies can be + //summed before adding noise and creating digits out of them + // + return (Float_t)(amp - fA)/fB; - SetName(title) ; - - // Post to the WhiteBoard - AliEMCALGetter * gime = AliEMCALGetter::GetInstance() ; - gime->PostSDigits( title, GetTitle()) ; } + +//__________________________________________________________________ +void AliEMCALSDigitizer::Print1(Option_t * option) +{ + Print(); + PrintSDigits(option); +} //__________________________________________________________________ -void AliEMCALSDigitizer::Print(Option_t* option)const +void AliEMCALSDigitizer::Print(Option_t *option) 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 in Tower = " << fTowerPrimThreshold << endl ; - cout << " Threshold for Primary assignment in PreShower = " << fPreShowerPrimThreshold << endl ; - cout << "---------------------------------------------------"<(AliRunLoader::Instance()->GetDetectorLoader("EMCAL")); + const TClonesArray * sdigits = rl->SDigits() ; - AliEMCALGetter * gime = AliEMCALGetter::GetInstance() ; - TString sdname(GetName()) ; - sdname.Remove(sdname.Index(GetTitle())-1) ; - const TClonesArray * sdigits = gime->SDigits(sdname.Data()) ; - - cout << "AliEMCALSDigitiser: event " << gAlice->GetEvNumber() << endl ; - cout << " Number of entries in SDigits list " << sdigits->GetEntriesFast() << endl ; - cout << endl ; + printf("\n") ; + printf("event %i", rl->GetRunLoader()->GetEventNumber()); + printf(" Number of entries in SDigits list %i", sdigits->GetEntriesFast()); if(strstr(option,"all")||strstr(option,"EMC")){ - + //loop over digits AliEMCALDigit * digit; - cout << "SDigit Id " << " Amplitude " << " Time " << " Index " << " Nprim " << " Primaries list " << endl; - Int_t index ; + printf("\n Id Amplitude Time Index Nprim: Primaries list \n") ; + Int_t index, isum=0; + char * tempo = new char[8192]; for (index = 0 ; index < sdigits->GetEntries() ; index++) { digit = dynamic_cast( sdigits->At(index) ) ; - cout << setw(6) << digit->GetId() << " " << setw(10) << digit->GetAmp() << " " << digit->GetTime() - << setw(6) << digit->GetIndexInList() << " " - << setw(5) << digit->GetNprimary() <<" "; + sprintf(tempo, "\n%6d %8d %6.5e %4d %2d :", + digit->GetId(), digit->GetAmp(), digit->GetTime(), digit->GetIndexInList(), digit->GetNprimary()) ; + printf("%s",tempo); + isum += digit->GetAmp(); Int_t iprimary; - for (iprimary=0; iprimaryGetNprimary(); iprimary++) - cout << " " << digit->GetPrimary(iprimary+1) << " "; - cout << endl; + for (iprimary=0; iprimaryGetNprimary(); iprimary++) { + sprintf(tempo, "%d ",digit->GetPrimary(iprimary+1) ) ; + printf("%s",tempo); + } } - cout <EMCALGeometry(); - Int_t ieta = ((ihit-1)/geom->GetNPhi())%geom->GetNZ(); // eta Tower Index - Int_t iphi = (ihit-1)%(geom->GetNPhi())+1; //phi Tower Index - Int_t it = -10; - Int_t ipre = 0; - - if (preshower)ipre = 1; - if (iphi > 0 && ieta >= 0){ - it = iphi+ieta*geom->GetNPhi() + ipre*geom->GetNPhi()*geom->GetNZ(); - return it; - }else{ - cerr << " AliEMCALSDigitizer::Layer2TowerID() -- there is an error "<< endl << "Eta number = " - << ieta << "Phi number = " << iphi << endl ; - return it; - } // end if iphi>0 && ieta>0 + // Unload Hits and SDigits from the folder + AliEMCALLoader *rl = dynamic_cast(AliRunLoader::Instance()->GetDetectorLoader("EMCAL")); + rl->UnloadHits() ; + rl->UnloadSDigits() ; } -//_______________________________________________________________________________________ -// void AliEMCALSDigitizer::TestTowerID(void) -// { -// Int_t j; - -// Bool_t preshower = kFALSE; -// for (j = 0 ; j < 10 ; j++){ // loop over hit id -// Int_t i; -// for (i = 0 ; i <= 2 ; i++){ // loop over -// Int_t k = i*96*144+j*144+1; -// cout << " Hit Index = " << k << " " << j*10 << " TOWERID = " << Layer2TowerID(k, preshower) << endl ; -// } -// } -// } //____________________________________________________________________________ -void AliEMCALSDigitizer::UseHitsFrom(const char * filename) +void AliEMCALSDigitizer::Browse(TBrowser* b) { - SetTitle(filename) ; - Init() ; + TTask::Browse(b); }