/* History of cvs commits:
*
* $Log$
+ * Revision 1.99 2006/11/07 16:49:51 kharlov
+ * Corrections for next event switch in case of raw data (B.Polichtchouk)
+ *
* Revision 1.98 2006/10/27 17:14:27 kharlov
* Introduce AliDebug and AliLog (B.Polichtchouk)
*
digit->SetIndexInList(i) ;
}
+ //Overwrite digits tree
+ AliPHOSGetter* gime = AliPHOSGetter::Instance();
+ TTree * treeD = gime->TreeD();
+ treeD->Branch("PHOS", &digits);
+ treeD->Fill() ;
+ gime->WriteDigits("OVERWRITE");
+ gime->PhosLoader()->UnloadDigits() ;
}
//____________________________________________________________________________
Bool_t AliPHOSClusterizerv1::IsInEmc(AliPHOSDigit * digit) const
Int_t eventNumber = runLoader->GetEventNumber() ;
- if(eventNumber==0) {
- rawReader->RewindEvents();
- rawReader->NextEvent();
- }
-
AliPHOSGetter *gime = AliPHOSGetter::Instance();
-
- Bool_t isOldRCUFormat = kFALSE;
- TString opt = GetOption();
- if(opt.Contains("OldRCUFormat"))
- isOldRCUFormat = kTRUE;
-
- gime->ReadRaw(rawReader,isOldRCUFormat) ;
+ gime->Event(eventNumber, "DRTP") ;
TClonesArray *recParticles = gime->RecParticles();
Int_t nOfRecParticles = recParticles->GetEntries();
AliFatal(Form("Digit not found at the expected position %d!",iDigit));
}
else {
- amplList[iDigit] = (UShort_t)(digit->GetEnergy()*500); // Energy in units of GeV/500
+ amplList[iDigit] = (UShort_t)digit->GetEnergy();
digiList[iDigit] = (UShort_t)(digit->GetId());
}
}
ec->SetM20(emcRP->GetM2z()) ; //second moment M2z
ec->SetNExMax(emcRP->GetNExMax()); //number of local maxima
ec->SetNumberOfDigits(digitMult); //digit multiplicity
- ec->SetDigitAmplitude(amplList); //energies in 1/500 of GeV
+ ec->SetDigitAmplitude(amplList); //digit energies
ec->SetDigitIndex(digiList); //abs id of the cell
ec->SetEmcCpvDistance(-1); //not yet implemented
ec->SetClusterChi2(-1); //not yet implemented