// --- ROOT system ---
-#include "TSystem.h"
-#include "TROOT.h"
-
+#include <TFile.h>
+#include <TROOT.h>
+#include <TSystem.h>
+#include <TParticle.h>
+#include <TF1.h>
+#include <TGraph.h>
+//#include <TCanvas.h>
+//#include <TFrame.h>
// --- Standard library ---
// --- AliRoot header files ---
-#include "AliPHOSGetter.h"
+#include "AliESD.h"
+#include "AliHeader.h"
+#include "AliMC.h"
#include "AliPHOS.h"
+#include "AliPHOSBeamTestEvent.h"
+#include "AliPHOSGetter.h"
+#include "AliPHOSLoader.h"
#include "AliRunLoader.h"
#include "AliStack.h"
-#include "AliPHOSLoader.h"
-#include "AliPHOSBeamTestEvent.h"
-#include "AliMC.h"
-#include "AliESD.h"
+#include "AliPHOSRawStream.h"
+#include "AliRawReaderFile.h"
+#include "AliLog.h"
ClassImp(AliPHOSGetter)
if (rl->GetAliRun() == 0x0) {
rl->LoadgAlice();
gAlice = rl->GetAliRun(); // should be removed
+ rl->LoadHeader();
}
}
fgPhosLoader = dynamic_cast<AliPHOSLoader*>(rl->GetLoader("PHOSLoader"));
else
fgPhosLoader->SetTitle(version);
-
// initialize data members
SetDebug(0) ;
fBTE = 0 ;
fPrimaries = 0 ;
fLoadingStatus = "" ;
- fESDFileName = "AliESDs.root" ;
+
+ fESDFileName = rl->GetFileName() ; // this should be the galice.root file
+ fESDFileName.ReplaceAll("galice.root", "AliESDs.root") ;
fESDFile = 0 ;
+ fESD = 0 ;
+ fESDTree = 0 ;
+ fRawDigits = kFALSE ;
}
//____________________________________________________________________________
AliPHOSGetter::~AliPHOSGetter()
{
// dtor
- delete fgPhosLoader ;
- fgPhosLoader = 0 ;
- delete fBTE ;
- fBTE = 0 ;
- fPrimaries->Delete() ;
- delete fPrimaries ;
+ if(fgPhosLoader){
+ delete fgPhosLoader ;
+ fgPhosLoader = 0 ;
+ }
+ if(fBTE){
+ delete fBTE ;
+ fBTE = 0 ;
+ }
+ if(fPrimaries){
+ fPrimaries->Delete() ;
+ delete fPrimaries ;
+ }
+ if (fESD)
+ delete fESD ;
+ if (fESDTree)
+ delete fESDTree ;
+
+ fgObjGetter = 0;
+}
+
+//____________________________________________________________________________
+void AliPHOSGetter::Reset()
+{
+ // resets things in case the getter is called consecutively with different files
+ // the PHOS Loader is already deleted by the Run Loader
+
+ if (fPrimaries) {
+ fPrimaries->Delete() ;
+ delete fPrimaries ;
+ }
+ fgPhosLoader = 0;
+ fgObjGetter = 0;
}
//____________________________________________________________________________
{
// Reads the content of all Tree's S, D and R
- if ( event >= MaxEvent() ) {
- Error("Event", "%d not found in TreeE !", event) ;
- return ;
- }
+// if ( event >= MaxEvent() ) {
+// Error("Event", "%d not found in TreeE !", event) ;
+// return ;
+// }
AliRunLoader * rl = AliRunLoader::GetRunLoader(PhosLoader()->GetTitle());
- // checks if we are dealing with test-beam data
- TBranch * btb = rl->TreeE()->GetBranch("AliPHOSBeamTestEvent") ;
- if(btb){
- if(!fBTE)
- fBTE = new AliPHOSBeamTestEvent() ;
- btb->SetAddress(&fBTE) ;
- btb->GetEntry(event) ;
- }
- else{
- if(fBTE){
- delete fBTE ;
- fBTE = 0 ;
- }
- }
+// // checks if we are dealing with test-beam data
+// TBranch * btb = rl->TreeE()->GetBranch("AliPHOSBeamTestEvent") ;
+// if(btb){
+// if(!fBTE)
+// fBTE = new AliPHOSBeamTestEvent() ;
+// btb->SetAddress(&fBTE) ;
+// btb->GetEntry(event) ;
+// }
+// else{
+// if(fBTE){
+// delete fBTE ;
+// fBTE = 0 ;
+// }
+// }
// Loads the type of object(s) requested
rl->GetEvent(event) ;
- if( strstr(opt,"X") || (strcmp(opt,"")==0) )
+ if(strstr(opt,"X") || (strcmp(opt,"")==0)){
ReadPrimaries() ;
-
- if(strstr(opt,"H") )
+ }
+
+ if(strstr(opt,"H") ){
ReadTreeH();
-
- if(strstr(opt,"S") )
+ }
+
+ if(strstr(opt,"S") ){
ReadTreeS() ;
-
- if( strstr(opt,"D") )
+ }
+
+ if(strstr(opt,"D") ){
ReadTreeD() ;
-
- if( strstr(opt,"R") )
+ }
+
+ if(strstr(opt,"R") ){
ReadTreeR() ;
+ }
- if( strstr(opt,"T") )
+ if( strstr(opt,"T") ){
ReadTreeT() ;
+ }
- if( strstr(opt,"P") )
+ if( strstr(opt,"P") ){
ReadTreeP() ;
-
-// if( strstr(opt,"Q") )
-// ReadTreeQA() ;
+ }
+
+ if( strstr(opt,"E") ){
+ ReadTreeE(event) ;
+ }
+
+}
+
+
+//____________________________________________________________________________
+void AliPHOSGetter::Event(AliRawReader *rawReader, const char* opt)
+{
+ // Reads the raw event from rawReader
+
+ if( strstr(opt,"W") ){
+ ReadRaw(rawReader) ;
+ }
}
// Creates and returns the pointer of the unique instance
// Must be called only when the environment has changed
- //::Info("Instance","alirunFileName=%s version=%s openingOption=%s",alirunFileName,version,openingOption);
-
if(!fgObjGetter){ // first time the getter is called
fgObjGetter = new AliPHOSGetter(alirunFileName, version, openingOption) ;
}
}
}
}
- else
+ else {
+ AliRunLoader * rl = AliRunLoader::GetRunLoader(fgPhosLoader->GetTitle()) ;
+ if ( strstr(version, AliConfig::GetDefaultEventFolderName()) ) // false in case of merging
+ delete rl ;
fgObjGetter = new AliPHOSGetter(alirunFileName, version, openingOption) ;
+ }
}
if (!fgObjGetter)
::Error("AliPHOSGetter::Instance", "Failed to create the PHOS Getter object") ;
AliPHOS * AliPHOSGetter:: PHOS() const
{
// returns the PHOS object
- AliPHOS * phos = dynamic_cast<AliPHOS*>(PhosLoader()->GetModulesFolder()->FindObject("PHOS")) ;
+ AliPHOSLoader * loader = 0;
+ static AliPHOSLoader * oldloader = 0;
+ static AliPHOS * phos = 0;
+
+ loader = PhosLoader();
+
+ if ( loader != oldloader) {
+ phos = dynamic_cast<AliPHOS*>(loader->GetModulesFolder()->FindObject("PHOS")) ;
+ oldloader = loader;
+ }
if (!phos)
if (fgDebug)
Warning("PHOS", "PHOS module not found in module folders: %s", PhosLoader()->GetModulesFolder()->GetName() ) ;
AliRunLoader * rl = AliRunLoader::GetRunLoader(PhosLoader()->GetTitle());
// gets kine tree from the root file (Kinematics.root)
- if ( ! rl->TreeK() ) // load treeK the first time
+ if ( ! rl->TreeK() ) { // load treeK the first time
rl->LoadKinematics() ;
+ }
- fNPrimaries = rl->Stack()->GetNtrack() ;
-
+ fNPrimaries = (rl->GetHeader())->GetNtrack();
if (fgDebug)
Info( "ReadTreeK", "Found %d particles in event # %d", fNPrimaries, EventNumber() ) ;
}
//____________________________________________________________________________
-AliESD * AliPHOSGetter::ESD(Int_t event)
-{
- //Read the ESD
- if (!fESDFile)
- OpenESDFile() ;
-
- TString esdEvent("ESD") ;
- esdEvent+= event ;
- AliESD * esd = dynamic_cast<AliESD *>(fESDFile->Get(esdEvent)) ;
- return esd ;
-}
-
-//____________________________________________________________________________
-Bool_t AliPHOSGetter::OpenESDFile(TString name)
+Bool_t AliPHOSGetter::OpenESDFile()
{
+ //Open the ESD file
Bool_t rv = kTRUE ;
- fESDFileName = name ;
- if (!fESDFile)
- fESDFile = new TFile(fESDFileName) ;
+ if (!fESDFile) {
+ fESDFile = TFile::Open(fESDFileName) ;
+ if (!fESDFile )
+ return kFALSE ;
+ }
else if (fESDFile->IsOpen()) {
fESDFile->Close() ;
fESDFile = TFile::Open(fESDFileName) ;
return rv ;
}
+//____________________________________________________________________________
+void AliPHOSGetter::FitRaw(Bool_t lowGainFlag, TGraph * gLowGain, TGraph * gHighGain, TF1* signalF, Int_t & amp, Double_t & time)
+{
+ // Fits the raw signal time distribution
+
+ const Int_t kNoiseThreshold = 0 ;
+ Double_t timezero1 = 0., timezero2 = 0., timemax = 0. ;
+ Double_t signal = 0., signalmax = 0. ;
+ Double_t energy = time = 0. ;
+
+ if (lowGainFlag) {
+ timezero1 = timezero2 = signalmax = timemax = 0. ;
+ signalF->FixParameter(0, PHOS()->GetRawFormatLowCharge()) ;
+ signalF->FixParameter(1, PHOS()->GetRawFormatLowGain()) ;
+ Int_t index ;
+ for (index = 0; index < PHOS()->GetRawFormatTimeBins(); index++) {
+ gLowGain->GetPoint(index, time, signal) ;
+ if (signal > kNoiseThreshold && timezero1 == 0.)
+ timezero1 = time ;
+ if (signal <= kNoiseThreshold && timezero1 > 0. && timezero2 == 0.)
+ timezero2 = time ;
+ if (signal > signalmax) {
+ signalmax = signal ;
+ timemax = time ;
+ }
+ }
+ signalmax /= PHOS()->RawResponseFunctionMax(PHOS()->GetRawFormatLowCharge(),
+ PHOS()->GetRawFormatLowGain()) ;
+ if ( timezero1 + PHOS()->GetRawFormatTimePeak() < PHOS()->GetRawFormatTimeMax() * 0.4 ) { // else its noise
+ signalF->SetParameter(2, signalmax) ;
+ signalF->SetParameter(3, timezero1) ;
+ gLowGain->Fit(signalF, "QRON", "", 0., timezero2); //, "QRON") ;
+ energy = signalF->GetParameter(2) ;
+ time = signalF->GetMaximumX() - PHOS()->GetRawFormatTimePeak() - PHOS()->GetRawFormatTimeTrigger() ;
+ }
+ } else {
+ timezero1 = timezero2 = signalmax = timemax = 0. ;
+ signalF->FixParameter(0, PHOS()->GetRawFormatHighCharge()) ;
+ signalF->FixParameter(1, PHOS()->GetRawFormatHighGain()) ;
+ Int_t index ;
+ for (index = 0; index < PHOS()->GetRawFormatTimeBins(); index++) {
+ gHighGain->GetPoint(index, time, signal) ;
+ if (signal > kNoiseThreshold && timezero1 == 0.)
+ timezero1 = time ;
+ if (signal <= kNoiseThreshold && timezero1 > 0. && timezero2 == 0.)
+ timezero2 = time ;
+ if (signal > signalmax) {
+ signalmax = signal ;
+ timemax = time ;
+ }
+ }
+ signalmax /= PHOS()->RawResponseFunctionMax(PHOS()->GetRawFormatHighCharge(),
+ PHOS()->GetRawFormatHighGain()) ;;
+ if ( timezero1 + PHOS()->GetRawFormatTimePeak() < PHOS()->GetRawFormatTimeMax() * 0.4 ) { // else its noise
+ signalF->SetParameter(2, signalmax) ;
+ signalF->SetParameter(3, timezero1) ;
+ gHighGain->Fit(signalF, "QRON", "", 0., timezero2) ;
+ energy = signalF->GetParameter(2) ;
+ time = signalF->GetMaximumX() - PHOS()->GetRawFormatTimePeak() - PHOS()->GetRawFormatTimeTrigger() ;
+ }
+ }
+
+ if (time == 0. && energy == 0.)
+ amp = 0 ;
+ else {
+ AliPHOSDigitizer * digitizer = Digitizer() ;
+ amp = static_cast<Int_t>( (energy - digitizer->GetEMCpedestal()) / digitizer->GetEMCchannel() + 0.5 ) ;
+ }
+}
+
+//____________________________________________________________________________
+Int_t AliPHOSGetter::ReadRaw(AliRawReader *rawReader)
+{
+ // reads the raw format data, converts it into digits format and store digits in Digits()
+ // container.
+
+ AliPHOSRawStream in(rawReader);
+
+ Bool_t first = kTRUE ;
+
+ TF1 * signalF = new TF1("signal", AliPHOS::RawResponseFunction, 0, PHOS()->GetRawFormatTimeMax(), 4);
+ signalF->SetParNames("Charge", "Gain", "Amplitude", "TimeZero") ;
+
+ Int_t relId[4], id ;
+ Bool_t lowGainFlag = kFALSE ;
+
+ TClonesArray * digits = Digits() ;
+ digits->Clear() ;
+ Int_t idigit = 0 ;
+ Int_t amp = 0 ;
+ Double_t time = 0. ;
+
+ TGraph * gLowGain = new TGraph(PHOS()->GetRawFormatTimeBins()) ;
+ TGraph * gHighGain= new TGraph(PHOS()->GetRawFormatTimeBins()) ;
+
+ while ( in.Next() ) { // PHOS entries loop
+ if ( (in.IsNewRow() || in.IsNewColumn() || in.IsNewModule()) ) {
+ if (!first) {
+ FitRaw(lowGainFlag, gLowGain, gHighGain, signalF, amp, time) ;
+ if (amp > 0) {
+ new((*digits)[idigit]) AliPHOSDigit( -1, id, amp, time) ;
+ idigit++ ;
+ }
+ Int_t index ;
+ for (index = 0; index < PHOS()->GetRawFormatTimeBins(); index++) {
+ gLowGain->SetPoint(index, index * PHOS()->GetRawFormatTimeMax() / PHOS()->GetRawFormatTimeBins(), 0) ;
+ gHighGain->SetPoint(index, index * PHOS()->GetRawFormatTimeMax() / PHOS()->GetRawFormatTimeBins(), 0) ;
+ }
+ }
+ first = kFALSE ;
+ relId[0] = in.GetModule() ;
+ if ( relId[0] >= PHOS()->GetRawFormatLowGainOffset() ) {
+ relId[0] -= PHOS()->GetRawFormatLowGainOffset() ;
+ lowGainFlag = kTRUE ;
+ } else
+ lowGainFlag = kFALSE ;
+ relId[1] = 0 ;
+ relId[2] = in.GetRow() ;
+ relId[3] = in.GetColumn() ;
+ PHOSGeometry()->RelToAbsNumbering(relId, id) ;
+ }
+ if (lowGainFlag)
+ gLowGain->SetPoint(in.GetTime(),
+ in.GetTime()* PHOS()->GetRawFormatTimeMax() / PHOS()->GetRawFormatTimeBins(),
+ in.GetSignal()) ;
+ else
+ gHighGain->SetPoint(in.GetTime(),
+ in.GetTime() * PHOS()->GetRawFormatTimeMax() / PHOS()->GetRawFormatTimeBins(),
+ in.GetSignal() ) ;
+
+ } // PHOS entries loop
+
+ FitRaw(lowGainFlag, gLowGain, gHighGain, signalF, amp, time) ;
+ if (amp > 0 ) {
+ new((*digits)[idigit]) AliPHOSDigit( -1, id, amp, time) ;
+ idigit++ ;
+ }
+ digits->Sort() ;
+
+ delete signalF ;
+ delete gLowGain, gHighGain ;
+
+ return digits->GetEntriesFast() ;
+}
+
+// TClonesArray * digits = Digits() ;
+// digits->Clear() ;
+// Int_t idigit = 0 ;
+
+// while ( in.Next() ) { // PHOS entries loop
+
+// Int_t amp = in.GetSignal() ;
+// Double_t time = in.GetTime() ;
+// Int_t relId[4], id ;
+
+// relId[0] = in.GetModule() ;
+// if ( relId[0] >= PHOS()->GetRawFormatLowGainOffset() ) {
+// relId[0] -= PHOS()->GetRawFormatLowGainOffset() ;
+// }
+// relId[1] = 0 ;
+// relId[2] = in.GetRow() ;
+// relId[3] = in.GetColumn() ;
+// PHOSGeometry()->RelToAbsNumbering(relId, id) ;
+
+// if (amp > 0 && id >=0 ) {
+// new((*digits)[idigit]) AliPHOSDigit( -1, id, amp, time) ;
+// idigit++ ;
+// }
+
+// } // PHOS entries loop
+
+// digits->Sort() ;
+
+// return digits->GetEntriesFast() ;
+// }
//____________________________________________________________________________
Int_t AliPHOSGetter::ReadTreeD()
{
// Read the Digits
-
- // gets TreeD from the root file (PHOS.SDigits.root)
- if ( !IsLoaded("D") ) {
- PhosLoader()->LoadDigits("UPDATE") ;
- PhosLoader()->LoadDigitizer("UPDATE") ;
- SetLoaded("D") ;
- }
+ PhosLoader()->CleanDigits() ;
+ PhosLoader()->LoadDigits("UPDATE") ;
+ PhosLoader()->LoadDigitizer("UPDATE") ;
return Digits()->GetEntries() ;
}
Int_t AliPHOSGetter::ReadTreeH()
{
// Read the Hits
-
+ PhosLoader()->CleanHits() ;
// gets TreeH from the root file (PHOS.Hit.root)
- if ( !IsLoaded("H") ) {
+ //if ( !IsLoaded("H") ) {
PhosLoader()->LoadHits("UPDATE") ;
- SetLoaded("H") ;
- }
+ // SetLoaded("H") ;
+ //}
return Hits()->GetEntries() ;
}
{
// Read the RecPoints
-
+ PhosLoader()->CleanRecPoints() ;
// gets TreeR from the root file (PHOS.RecPoints.root)
- if ( !IsLoaded("R") ) {
+ //if ( !IsLoaded("R") ) {
PhosLoader()->LoadRecPoints("UPDATE") ;
PhosLoader()->LoadClusterizer("UPDATE") ;
- SetLoaded("R") ;
- }
+ // SetLoaded("R") ;
+ //}
return EmcRecPoints()->GetEntries() ;
}
{
// Read the TrackSegments
-
+ PhosLoader()->CleanTracks() ;
// gets TreeT from the root file (PHOS.TrackSegments.root)
- if ( !IsLoaded("T") ) {
+ //if ( !IsLoaded("T") ) {
PhosLoader()->LoadTracks("UPDATE") ;
PhosLoader()->LoadTrackSegmentMaker("UPDATE") ;
- SetLoaded("T") ;
- }
+ // SetLoaded("T") ;
+ //}
return TrackSegments()->GetEntries() ;
}
//____________________________________________________________________________
Int_t AliPHOSGetter::ReadTreeP()
{
- // Read the TrackSegments
-
+ // Read the RecParticles
+ PhosLoader()->CleanRecParticles() ;
+
// gets TreeT from the root file (PHOS.TrackSegments.root)
- if ( !IsLoaded("P") ) {
+ // if ( !IsLoaded("P") ) {
PhosLoader()->LoadRecParticles("UPDATE") ;
PhosLoader()->LoadPID("UPDATE") ;
- SetLoaded("P") ;
- }
+ // SetLoaded("P") ;
+ //}
return RecParticles()->GetEntries() ;
}
{
// Read the SDigits
-
+ PhosLoader()->CleanSDigits() ;
// gets TreeS from the root file (PHOS.SDigits.root)
- if ( !IsLoaded("S") ) {
- PhosLoader()->LoadSDigits("UPDATE") ;
- PhosLoader()->LoadSDigitizer("UPDATE") ;
- SetLoaded("S") ;
- }
+ //if ( !IsLoaded("S") ) {
+ PhosLoader()->LoadSDigits("READ") ;
+ PhosLoader()->LoadSDigitizer("READ") ;
+ // SetLoaded("S") ;
+ //}
return SDigits()->GetEntries() ;
}
+//____________________________________________________________________________
+Int_t AliPHOSGetter::ReadTreeE(Int_t event)
+{
+ // Read the ESD
+
+ // gets esdTree from the root file (AliESDs.root)
+ if (!fESDFile)
+ if ( !OpenESDFile() )
+ return -1 ;
+
+ fESDTree = static_cast<TTree*>(fESDFile->Get("esdTree")) ;
+ fESD = new AliESD;
+ if (!fESDTree) {
+
+ Error("ReadTreeE", "no ESD tree found");
+ return -1;
+ }
+ fESDTree->SetBranchAddress("ESD", &fESD);
+ fESDTree->GetEvent(event);
+
+ return event ;
+}
+
//____________________________________________________________________________
TClonesArray * AliPHOSGetter::SDigits()
{
if ( opt == "sdigits") {
// add the version name to the root file name
TString fileName( PhosLoader()->GetSDigitsFileName() ) ;
- if (version != AliConfig::fgkDefaultEventFolderName) // only if not the default folder name
+ if (version != AliConfig::GetDefaultEventFolderName()) // only if not the default folder name
fileName = fileName.ReplaceAll(".root", "") + "_" + version + ".root" ;
if ( !(gSystem->AccessPathName(fileName)) ) {
Warning("VersionExists", "The file %s already exists", fileName.Data()) ;
if ( opt == "digits") {
// add the version name to the root file name
TString fileName( PhosLoader()->GetDigitsFileName() ) ;
- if (version != AliConfig::fgkDefaultEventFolderName) // only if not the default folder name
+ if (version != AliConfig::GetDefaultEventFolderName()) // only if not the default folder name
fileName = fileName.ReplaceAll(".root", "") + "_" + version + ".root" ;
if ( !(gSystem->AccessPathName(fileName)) ) {
Warning("VersionExists", "The file %s already exists", fileName.Data()) ;