**************************************************************************/
//_____________________________________________________________________________
-// Steering class for particle (gamma, hadron) identification and correlation
-// analysis. It is called by the task class AliAnalysisTaskCaloTrackCorrelation
-// and it connects the input (ESD/AOD/MonteCarlo) got with AliCaloTrackReader
-// (produces TClonesArrays of AODs (TParticles in MC case if requested)), with
+// Steering class for particle (gamma, hadron) identification and correlation
+// analysis. It is called by the task class AliAnalysisTaskCaloTrackCorrelation
+// and it connects the input (ESD/AOD/MonteCarlo) got with AliCaloTrackReader
+// (produces TClonesArrays of AODs (TParticles in MC case if requested)), with
// the analysis classes that derive from AliAnaCaloTrackCorrBaseClass
//
// -- Author: Gustavo Conesa (INFN-LNF, LPSC-Grenoble)
#include <cstdlib>
// --- ROOT system ---
-#include "TClonesArray.h"
-#include "TList.h"
-#include "TH1F.h"
+#include <TClonesArray.h>
+#include <TList.h>
+#include <TH1F.h>
//#include <TObjectTable.h>
+#include <TGeoGlobalMagField.h>
//---- AliRoot system ----
#include "AliAnalysisManager.h"
//__________________________________________________
-AliAnaCaloTrackCorrMaker::AliAnaCaloTrackCorrMaker() :
+AliAnaCaloTrackCorrMaker::AliAnaCaloTrackCorrMaker() :
TObject(),
fReader(0), fCaloUtils(0),
fOutputContainer(new TList ), fAnalysisContainer(new TList ),
-fMakeHisto(kFALSE), fMakeAOD(kFALSE),
-fAnaDebug(0), fCuts(new TList),
+fMakeHisto(kFALSE), fMakeAOD(kFALSE),
+fAnaDebug(0), fCuts(new TList),
fScaleFactor(-1),
+fFillDataControlHisto(kTRUE),
+// Control histograms
fhNEvents(0), fhNExoticEvents(0),
fhNEventsNoTriggerFound(0),
fhNPileUpEvents(0), fhNPileUpEventsTriggerBC0(0),
//Default Ctor
if(fAnaDebug > 1 ) printf("*** Analysis Maker Constructor *** \n");
+ for(Int_t i = 0; i < 3; i++)
+ {
+ fhClusterTriggerBCUnMatchReMatch [0] = 0;
+ fhClusterTriggerBCExoticUnMatchReMatch[0] = 0;
+ }
+
//Initialize parameters, pointers and histograms
InitParameters();
}
//________________________________________________________________________________________
-AliAnaCaloTrackCorrMaker::AliAnaCaloTrackCorrMaker(const AliAnaCaloTrackCorrMaker & maker) :
+AliAnaCaloTrackCorrMaker::AliAnaCaloTrackCorrMaker(const AliAnaCaloTrackCorrMaker & maker) :
TObject(),
fReader(), //(new AliCaloTrackReader(*maker.fReader)),
fCaloUtils(),//(new AliCalorimeterUtils(*maker.fCaloUtils)),
-fOutputContainer(new TList()), fAnalysisContainer(new TList()),
+fOutputContainer(new TList()), fAnalysisContainer(new TList()),
fMakeHisto(maker.fMakeHisto), fMakeAOD(maker.fMakeAOD),
fAnaDebug(maker.fAnaDebug), fCuts(new TList()),
fScaleFactor(maker.fScaleFactor),
+fFillDataControlHisto(maker.fFillDataControlHisto),
fhNEvents(maker.fhNEvents),
fhNExoticEvents(maker.fhNExoticEvents),
fhNEventsNoTriggerFound(maker.fhNEventsNoTriggerFound),
fhTrackMult(maker.fhTrackMult),
fhCentrality(maker.fhCentrality),
fhEventPlaneAngle(maker.fhEventPlaneAngle),
-fhNMergedFiles(maker.fhNMergedFiles),
+fhNMergedFiles(maker.fhNMergedFiles),
fhScaleFactor(maker.fhScaleFactor),
fhEMCalBCEvent(maker.fhEMCalBCEvent),
fhEMCalBCEventCut(maker.fhEMCalBCEventCut),
fhClusterTriggerBCExoticEventBCUnMatch(maker.fhClusterTriggerBCExoticEventBCUnMatch)
{
+ for(Int_t i = 0; i < 3; i++)
+ {
+ fhClusterTriggerBCUnMatchReMatch [i] = maker.fhClusterTriggerBCUnMatchReMatch [i];
+ fhClusterTriggerBCExoticUnMatchReMatch[i] = maker.fhClusterTriggerBCExoticUnMatchReMatch[i];
+ }
// cpy ctor
}
//___________________________________________________
-AliAnaCaloTrackCorrMaker::~AliAnaCaloTrackCorrMaker()
+AliAnaCaloTrackCorrMaker::~AliAnaCaloTrackCorrMaker()
{
// Remove all owned pointers.
- // Do not delete it here, already done somewhere else, need to understand where.
+ // Do not delete it here, already done somewhere else, need to understand where.
// if (fOutputContainer) {
// fOutputContainer->Clear();
// delete fOutputContainer ;
- // }
+ // }
if (fAnalysisContainer)
{
fAnalysisContainer->Delete();
delete fAnalysisContainer ;
- }
+ }
if (fReader) delete fReader ;
if (fCaloUtils) delete fCaloUtils ;
}
//__________________________________________________________________
-void AliAnaCaloTrackCorrMaker::AddAnalysis(TObject* ana, Int_t n)
+void AliAnaCaloTrackCorrMaker::AddAnalysis(TObject* ana, Int_t n)
{
// Add analysis depending on AliAnaCaloTrackCorrBaseClass to list
if ( fAnalysisContainer)
- {
- fAnalysisContainer->AddAt(ana,n);
+ {
+ fAnalysisContainer->AddAt(ana,n);
}
else
- {
+ {
printf("AliAnaCaloTrackCorrMaker::AddAnalysis() - AnalysisContainer not initialized\n");
abort();
}
-}
+}
//_________________________________________________________
TList * AliAnaCaloTrackCorrMaker::FillAndGetAODBranchList()
-{
+{
// Get any new output AOD branches from analysis and put them in a list
// The list is filled in the maker, and new branch passed to the analysis frame
TList *aodBranchList = fReader->GetAODBranchList() ;
for(Int_t iana = 0; iana < fAnalysisContainer->GetEntries(); iana++)
- {
- AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
- if(ana->NewOutputAOD()) aodBranchList->Add(ana->GetCreateOutputAODBranch());
- }
+ {
+ AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
+ if(ana->NewOutputAOD()) aodBranchList->Add(ana->GetCreateOutputAODBranch());
+ }
return aodBranchList ;
}
-//_________________________________________________________
+//____________________________________________________
void AliAnaCaloTrackCorrMaker::FillControlHistograms()
{
// Event control histograms
fhNEvents ->Fill(0); // Number of events analyzed
- if( fReader->IsPileUpFromSPD())
- fhNPileUpEvents->Fill(0.5);
- //if( event->IsPileupFromSPDInMultBins())
- // fhNPileUpEvents->Fill(1.5);
- if( fReader->IsPileUpFromEMCal())
- fhNPileUpEvents->Fill(2.5);
- if( fReader->IsPileUpFromSPDOrEMCal() )
- fhNPileUpEvents->Fill(3.5);
- if( fReader->IsPileUpFromSPDAndEMCal() )
- fhNPileUpEvents->Fill(4.5);
- if( fReader->IsPileUpFromSPDAndNotEMCal() )
- fhNPileUpEvents->Fill(5.5);
- if( fReader->IsPileUpFromEMCalAndNotSPD() )
- fhNPileUpEvents->Fill(6.5);
- if( fReader->IsPileUpFromNotSPDAndNotEMCal() )
- fhNPileUpEvents->Fill(7.5);
-
- Int_t triggerBC = fReader->GetTriggerClusterBC() ;
- if( triggerBC == 0 &&
- !fReader->IsExoticEvent() &&
- !fReader->IsBadCellTriggerEvent())
+ Double_t v[3];
+ event->GetPrimaryVertex()->GetXYZ(v) ;
+ fhXVertex->Fill(v[0]);
+ fhYVertex->Fill(v[1]);
+ fhZVertex->Fill(v[2]);
+
+ fhTrackMult ->Fill(fReader->GetTrackMultiplicity());
+ fhCentrality ->Fill(fReader->GetEventCentrality ());
+ fhEventPlaneAngle ->Fill(fReader->GetEventPlaneAngle ());
+
+ if(fFillDataControlHisto)
{
if( fReader->IsPileUpFromSPD())
- fhNPileUpEventsTriggerBC0->Fill(0.5);
+ fhNPileUpEvents->Fill(0.5);
//if( event->IsPileupFromSPDInMultBins())
- // fhNPileUpEventsTriggerBC0->Fill(1.5);
+ // fhNPileUpEvents->Fill(1.5);
if( fReader->IsPileUpFromEMCal())
- fhNPileUpEventsTriggerBC0->Fill(2.5);
+ fhNPileUpEvents->Fill(2.5);
if( fReader->IsPileUpFromSPDOrEMCal() )
- fhNPileUpEventsTriggerBC0->Fill(3.5);
+ fhNPileUpEvents->Fill(3.5);
if( fReader->IsPileUpFromSPDAndEMCal() )
- fhNPileUpEventsTriggerBC0->Fill(4.5);
+ fhNPileUpEvents->Fill(4.5);
if( fReader->IsPileUpFromSPDAndNotEMCal() )
- fhNPileUpEventsTriggerBC0->Fill(5.5);
+ fhNPileUpEvents->Fill(5.5);
if( fReader->IsPileUpFromEMCalAndNotSPD() )
- fhNPileUpEventsTriggerBC0->Fill(6.5);
+ fhNPileUpEvents->Fill(6.5);
if( fReader->IsPileUpFromNotSPDAndNotEMCal() )
- fhNPileUpEventsTriggerBC0->Fill(7.5);
+ fhNPileUpEvents->Fill(7.5);
+
+ Int_t triggerBC = fReader->GetTriggerClusterBC() ;
+ if( triggerBC == 0 &&
+ !fReader->IsExoticEvent() &&
+ !fReader->IsBadCellTriggerEvent())
+ {
+ if( fReader->IsPileUpFromSPD())
+ fhNPileUpEventsTriggerBC0->Fill(0.5);
+ //if( event->IsPileupFromSPDInMultBins())
+ // fhNPileUpEventsTriggerBC0->Fill(1.5);
+ if( fReader->IsPileUpFromEMCal())
+ fhNPileUpEventsTriggerBC0->Fill(2.5);
+ if( fReader->IsPileUpFromSPDOrEMCal() )
+ fhNPileUpEventsTriggerBC0->Fill(3.5);
+ if( fReader->IsPileUpFromSPDAndEMCal() )
+ fhNPileUpEventsTriggerBC0->Fill(4.5);
+ if( fReader->IsPileUpFromSPDAndNotEMCal() )
+ fhNPileUpEventsTriggerBC0->Fill(5.5);
+ if( fReader->IsPileUpFromEMCalAndNotSPD() )
+ fhNPileUpEventsTriggerBC0->Fill(6.5);
+ if( fReader->IsPileUpFromNotSPDAndNotEMCal() )
+ fhNPileUpEventsTriggerBC0->Fill(7.5);
+ }
+
+ if(fReader->IsPileUpFromSPD())
+ fhPileUpClusterMultAndSPDPileUp ->Fill(fReader->GetNPileUpClusters());
+
+ fhPileUpClusterMult ->Fill(fReader->GetNPileUpClusters ());
+
+ for(Int_t i = 0; i < 19; i++)
+ {
+ if(fReader->GetTrackEventBC(i)) fhTrackBCEvent ->Fill(i);
+ if(fReader->GetTrackEventBCcut(i))fhTrackBCEventCut->Fill(i);
+ if(fReader->GetEMCalEventBC(i)) fhEMCalBCEvent ->Fill(i);
+ if(fReader->GetEMCalEventBCcut(i))fhEMCalBCEventCut->Fill(i);
+ }
+
+ Int_t bc = fReader->GetVertexBC();
+ if(bc!=AliVTrack::kTOFBCNA)fhPrimaryVertexBC->Fill(bc);
+
+
+ // N pile up vertices
+ Int_t nVerticesSPD = -1;
+ Int_t nVerticesTracks = -1;
+
+ if (esdevent)
+ {
+ nVerticesSPD = esdevent->GetNumberOfPileupVerticesSPD();
+ nVerticesTracks = esdevent->GetNumberOfPileupVerticesTracks();
+
+ }//ESD
+ else if (aodevent)
+ {
+ nVerticesSPD = aodevent->GetNumberOfPileupVerticesSPD();
+ nVerticesTracks = aodevent->GetNumberOfPileupVerticesTracks();
+ }//AOD
+
+ fhNPileUpVertSPD ->Fill(nVerticesSPD);
+ fhNPileUpVertTracks->Fill(nVerticesTracks);
+
+ // Time stamp
+ if(fReader->IsSelectEventTimeStampOn() && esdevent)
+ {
+ Int_t timeStamp = esdevent->GetTimeStamp();
+ Float_t timeStampFrac = 1.*(timeStamp-fReader->GetRunTimeStampMin()) /
+ (fReader->GetRunTimeStampMax()-fReader->GetRunTimeStampMin());
+
+ //printf("stamp %d, min %d, max %d, frac %f\n", timeStamp, fReader->GetRunTimeStampMin(), fReader->GetRunTimeStampMax(), timeStampFrac);
+
+ fhTimeStampFraction->Fill(timeStampFrac);
+ }
}
+}
+
+//___________________________________________________________
+void AliAnaCaloTrackCorrMaker::FillTriggerControlHistograms()
+{
+ if(!fFillDataControlHisto) return;
- if(fReader->IsPileUpFromSPD())
- fhPileUpClusterMultAndSPDPileUp ->Fill(fReader->GetNPileUpClusters());
-
- fhPileUpClusterMult ->Fill(fReader->GetNPileUpClusters ());
- fhTrackMult ->Fill(fReader->GetTrackMultiplicity());
- fhCentrality ->Fill(fReader->GetEventCentrality ());
- fhEventPlaneAngle ->Fill(fReader->GetEventPlaneAngle ());
+ Int_t triggerBC = fReader->GetTriggerClusterBC() ;
+ Bool_t exotic = fReader->IsExoticEvent();
+ Bool_t badCluster = fReader->IsBadCellTriggerEvent();
+ Bool_t badCell = fReader->IsBadMaxCellTriggerEvent();
+ Bool_t triggerMatch= fReader->IsTriggerMatched();
+ Bool_t triggerBCOK = kTRUE;
+ Int_t triggerId = fReader->GetTriggerClusterId() ;
+ Bool_t reMatchOpenTime = fReader->IsTriggerMatchedOpenCuts(0);
+ Bool_t reMatchNeigbour = fReader->IsTriggerMatchedOpenCuts(1);
+ Bool_t reMatchBoth = fReader->IsTriggerMatchedOpenCuts(2);
- for(Int_t i = 0; i < 19; i++)
+ if(triggerId < 0)
{
- if(fReader->GetTrackEventBC(i)) fhTrackBCEvent ->Fill(i);
- if(fReader->GetTrackEventBCcut(i))fhTrackBCEventCut->Fill(i);
- if(fReader->GetEMCalEventBC(i)) fhEMCalBCEvent ->Fill(i);
- if(fReader->GetEMCalEventBCcut(i))fhEMCalBCEventCut->Fill(i);
+ //printf("Trigger id %d\n",triggerId);
+ if(triggerId == -2)fhNEventsNoTriggerFound->Fill(0);
+ triggerBCOK = kFALSE;
}
- Double_t v[3];
- event->GetPrimaryVertex()->GetXYZ(v) ;
- fhXVertex->Fill(v[0]);
- fhYVertex->Fill(v[1]);
- fhZVertex->Fill(v[2]);
-
- Int_t bc = fReader->GetVertexBC();
- if(bc!=AliVTrack::kTOFBCNA)fhPrimaryVertexBC->Fill(bc);
+ if(exotic)
+ {
+ fhNExoticEvents->Fill(0) ;
+ Double_t v[3];
+ fReader->GetInputEvent()->GetPrimaryVertex()->GetXYZ(v) ;
+ fhXVertexExotic->Fill(v[0]);
+ fhYVertexExotic->Fill(v[1]);
+ fhZVertexExotic->Fill(v[2]);
+ }
+ //if(fReader->IsExoticEvent()) printf("Maker: EXOTIC Cluster trigger\n");
+ if(!triggerBCOK) return;
- // N pile up vertices
- Int_t nVerticesSPD = -1;
- Int_t nVerticesTracks = -1;
+ Int_t eventBC = fReader->GetInputEvent()->GetBunchCrossNumber();
+ if(eventBC%4 < 0 || eventBC%4 > 3 )
+ printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - STRANGE: Trigger BC %d - Event BC %d, modulo4 %d \n",triggerBC,eventBC,eventBC%4);
- if (esdevent)
+ if(triggerMatch)
{
- nVerticesSPD = esdevent->GetNumberOfPileupVerticesSPD();
- nVerticesTracks = esdevent->GetNumberOfPileupVerticesTracks();
+ if (!exotic && !badCluster) fhClusterTriggerBC->Fill(triggerBC);
+ else if( exotic && badCluster)
+ {
+ fhClusterTriggerBCBadClusterExotic->Fill(triggerBC);
+ if(badCell) fhClusterTriggerBCBadCellExotic->Fill(triggerBC);
+ }
+ else if( exotic && !badCluster) fhClusterTriggerBCExotic->Fill(triggerBC);
+ else if( badCluster && !exotic )
+ {
+ fhClusterTriggerBCBadCluster ->Fill(triggerBC);
+ if(badCell) fhClusterTriggerBCBadCell->Fill(triggerBC);
+ }
- }//ESD
- else if (aodevent)
- {
- nVerticesSPD = aodevent->GetNumberOfPileupVerticesSPD();
- nVerticesTracks = aodevent->GetNumberOfPileupVerticesTracks();
- }//AOD
-
- fhNPileUpVertSPD ->Fill(nVerticesSPD);
- fhNPileUpVertTracks->Fill(nVerticesTracks);
-
- // Time stamp
- if(fReader->IsSelectEventTimeStampOn() && esdevent)
+ if(!exotic) fhClusterTriggerBCEventBC ->Fill(triggerBC,eventBC%4);
+ else fhClusterTriggerBCExoticEventBC->Fill(triggerBC,eventBC%4);
+ }
+ else
{
- Int_t timeStamp = esdevent->GetTimeStamp();
- Float_t timeStampFrac = 1.*(timeStamp-fReader->GetRunTimeStampMin()) /
- (fReader->GetRunTimeStampMax()-fReader->GetRunTimeStampMin());
+ if (!exotic && !badCluster)
+ {
+ fhClusterTriggerBCUnMatch->Fill(triggerBC);
+ if(reMatchOpenTime) fhClusterTriggerBCUnMatchReMatch[0]->Fill(triggerBC);
+ if(reMatchNeigbour) fhClusterTriggerBCUnMatchReMatch[1]->Fill(triggerBC);
+ if(reMatchBoth) fhClusterTriggerBCUnMatchReMatch[2]->Fill(triggerBC);
+ }
+ else if( exotic && badCluster)
+ {
+ fhClusterTriggerBCBadClusterExoticUnMatch->Fill(triggerBC);
+ if(badCell) fhClusterTriggerBCBadCellExoticUnMatch ->Fill(triggerBC);
+ }
+ else if( exotic && !badCluster)
+ {
+ fhClusterTriggerBCExoticUnMatch->Fill(triggerBC);
+ if(reMatchOpenTime) fhClusterTriggerBCExoticUnMatchReMatch[0]->Fill(triggerBC);
+ if(reMatchNeigbour) fhClusterTriggerBCExoticUnMatchReMatch[1]->Fill(triggerBC);
+ if(reMatchBoth) fhClusterTriggerBCExoticUnMatchReMatch[2]->Fill(triggerBC);
+ }
+ else if( badCluster && !exotic )
+ {
+ fhClusterTriggerBCBadClusterUnMatch->Fill(triggerBC);
+ if(badCell)fhClusterTriggerBCBadCellUnMatch->Fill(triggerBC);
+ }
- //printf("stamp %d, min %d, max %d, frac %f\n", timeStamp, fReader->GetRunTimeStampMin(), fReader->GetRunTimeStampMax(), timeStampFrac);
-
- fhTimeStampFraction->Fill(timeStampFrac);
+ if(!exotic) fhClusterTriggerBCEventBCUnMatch ->Fill(triggerBC,eventBC%4);
+ else fhClusterTriggerBCExoticEventBCUnMatch->Fill(triggerBC,eventBC%4);
}
+
}
+
//_______________________________________________________
TList * AliAnaCaloTrackCorrMaker::GetListOfAnalysisCuts()
-{
+{
// Get the list of the cuts used for the analysis
// The list is filled in the maker, called by the task in LocalInit() and posted there
for(Int_t iana = 0; iana < fAnalysisContainer->GetEntries(); iana++)
- {
- AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
- TObjString * objstring = ana->GetAnalysisCuts();
-
- if(objstring)fCuts->Add(objstring);
- }
+ {
+ AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
+ TObjString * objstring = ana->GetAnalysisCuts();
+
+ if(objstring)fCuts->Add(objstring);
+ }
return fCuts ;
fhNEvents->SetYTitle("# events");
fOutputContainer->Add(fhNEvents);
- fhNExoticEvents = new TH1F("hNExoticEvents", "Number of analyzed events triggered by exotic cluster" , 1 , 0 , 1 ) ;
- fhNExoticEvents->SetYTitle("# exotic events");
- fOutputContainer->Add(fhNExoticEvents);
-
- fhNEventsNoTriggerFound = new TH1F("hNEventsNoTriggerFound", "Number of analyzed events triggered but no trigger found" , 1 , 0 , 1 ) ;
- fhNEventsNoTriggerFound->SetYTitle("# exotic events");
- fOutputContainer->Add(fhNEventsNoTriggerFound);
-
-
- Int_t nbin = 11;
- Float_t minbin =-5.5;
- Float_t maxbin = 5.5;
- Int_t labelshift = 6;
-
- fhClusterTriggerBCEventBC = new TH2F("hClusterTriggerBCEventBC", "Found trigger BC and Event BC",
- nbin , minbin ,maxbin,4,0, 4) ;
- fhClusterTriggerBCEventBC->SetXTitle("cluster trigger BC");
- for(Int_t i = 0; i < 4; i++)
- fhClusterTriggerBCEventBC->GetYaxis()->SetBinLabel(i+1 ,Form("BC/4=%d",i));
- fhClusterTriggerBCEventBC->SetXTitle("cluster trigger BC");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCEventBC->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fhClusterTriggerBCEventBC->SetYTitle("Event BC%4");
- fOutputContainer->Add(fhClusterTriggerBCEventBC);
-
- fhClusterTriggerBCExoticEventBC = new TH2F("hClusterTriggerBCExoticEventBC", "Found exotic trigger BC and Event BC",
- nbin , minbin ,maxbin,4,1, 4) ;
- for(Int_t i = 0; i < 4; i++)
- fhClusterTriggerBCExoticEventBC->GetYaxis()->SetBinLabel(i+1 ,Form("BC/4=%d",i));
- fhClusterTriggerBCExoticEventBC->SetXTitle("cluster trigger BC");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCExoticEventBC->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fhClusterTriggerBCExoticEventBC->SetYTitle("Event BC%4");
- fOutputContainer->Add(fhClusterTriggerBCExoticEventBC);
-
- fhClusterTriggerBCEventBCUnMatch = new TH2F("hClusterTriggerBCEventBCUnMatch", "Found unmatched trigger BC and Event BC",
- nbin , minbin ,maxbin,4,1, 4) ;
- for(Int_t i = 0; i < 4; i++)
- fhClusterTriggerBCEventBCUnMatch->GetYaxis()->SetBinLabel(i+1 ,Form("BC/4=%d",i));
- fhClusterTriggerBCEventBCUnMatch->SetXTitle("cluster trigger BC");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCEventBCUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fhClusterTriggerBCEventBCUnMatch->SetYTitle("Event BC%4");
- fOutputContainer->Add(fhClusterTriggerBCEventBCUnMatch);
-
- fhClusterTriggerBCExoticEventBCUnMatch = new TH2F("hClusterTriggerExoticBCEventBCUnMatch", "Found unmatched trigger BC and Event BC",
- nbin , minbin ,maxbin,4,1, 4) ;
- for(Int_t i = 0; i < 4; i++)
- fhClusterTriggerBCExoticEventBCUnMatch->GetYaxis()->SetBinLabel(i+1 ,Form("BC/4=%d",i));
- fhClusterTriggerBCExoticEventBCUnMatch->SetXTitle("cluster trigger BC");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCExoticEventBCUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fhClusterTriggerBCExoticEventBCUnMatch->SetYTitle("Event BC%4");
- fOutputContainer->Add(fhClusterTriggerBCExoticEventBCUnMatch);
-
- fhClusterTriggerBC = new TH1F("hClusterTriggerBC",
- "Number of analyzed events triggered by a cluster in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBC->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBC->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBC);
-
- fhClusterTriggerBCExotic = new TH1F("hClusterTriggerBCExotic",
- "Number of analyzed events triggered by a exotic cluster in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBCExotic->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCExotic->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCExotic);
-
-
- fhClusterTriggerBCBadCell = new TH1F("hClusterTriggerBCBadCell",
- "Number of analyzed events triggered by a bad cell in a given BC",
- nbin , minbin ,maxbin) ;
-
- fhClusterTriggerBCBadCell->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCBadCell->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCBadCell);
-
- fhClusterTriggerBCBadCellExotic = new TH1F("hClusterTriggerBCBadCellExotic",
- "Number of analyzed events triggered by a bad cell & exotic cluster in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBCBadCellExotic->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCBadCellExotic->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCBadCellExotic);
-
- fhClusterTriggerBCBadCluster = new TH1F("hClusterTriggerBCBadCluster",
- "Number of analyzed events triggered by a bad cluster in a given BC",
- nbin , minbin ,maxbin) ;
-
- fhClusterTriggerBCBadCluster->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCBadCluster->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCBadCluster);
-
-
- fhClusterTriggerBCBadClusterExotic = new TH1F("hClusterTriggerBCBadClusterExotic",
- "Number of analyzed events triggered by a bad cluster & exotic cluster in a given BC",
- nbin , minbin ,maxbin) ;
-
- fhClusterTriggerBCBadClusterExotic->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCBadClusterExotic->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCBadClusterExotic);
-
- fhClusterTriggerBCUnMatch = new TH1F("hClusterTriggerBCUnMatch",
- "Number of analyzed events triggered by a cluster (no trigger patch match) in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBCUnMatch->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCUnMatch);
-
- fhClusterTriggerBCExoticUnMatch = new TH1F("hClusterTriggerBCExoticUnMatch",
- "Number of analyzed events triggered by a exotic cluster (no trigger patch match) in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBCExoticUnMatch->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCExoticUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCExoticUnMatch);
-
-
- fhClusterTriggerBCBadCellUnMatch = new TH1F("hClusterTriggerBCBadCellUnMatch",
- "Number of analyzed events triggered by a bad cluster (no trigger patch match) in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBCBadCellUnMatch->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCBadCellUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCBadCellUnMatch);
-
-
- fhClusterTriggerBCBadCellExoticUnMatch = new TH1F("hClusterTriggerBCBadCellExoticUnMatch",
- "Number of analyzed events triggered by a bad&exotic cluster (no trigger patch match) in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBCBadCellExoticUnMatch->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCBadCellExoticUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCBadCellExoticUnMatch);
-
-
- fhClusterTriggerBCBadClusterUnMatch = new TH1F("hClusterTriggerBCBadClusterUnMatch",
- "Number of analyzed events triggered by a bad cluster (no trigger patch match) in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBCBadClusterUnMatch->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCBadClusterUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCBadClusterUnMatch);
-
-
- fhClusterTriggerBCBadClusterExoticUnMatch = new TH1F("hClusterTriggerBCBadClusterExoticUnMatch",
- "Number of analyzed events triggered by a bad&exotic cluster (no trigger patch match) in a given BC",
- nbin , minbin ,maxbin) ;
- fhClusterTriggerBCBadClusterExoticUnMatch->SetYTitle("# events");
- for(Int_t i = 1; i < 12; i++)
- fhClusterTriggerBCBadClusterExoticUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
- fOutputContainer->Add(fhClusterTriggerBCBadClusterExoticUnMatch);
-
-
- fhNPileUpEvents = new TH1F("hNPileUpEvents", "Number of events considered as pile-up", 8 , 0 , 8 ) ;
- fhNPileUpEvents->SetYTitle("# events");
- fhNPileUpEvents->GetXaxis()->SetBinLabel(1 ,"SPD");
- fhNPileUpEvents->GetXaxis()->SetBinLabel(2 ,"Multi SPD");
- fhNPileUpEvents->GetXaxis()->SetBinLabel(3 ,"EMCal");
- fhNPileUpEvents->GetXaxis()->SetBinLabel(4 ,"EMCal || SPD");
- fhNPileUpEvents->GetXaxis()->SetBinLabel(5 ,"EMCal && SPD");
- fhNPileUpEvents->GetXaxis()->SetBinLabel(6 ,"!EMCal && SPD");
- fhNPileUpEvents->GetXaxis()->SetBinLabel(7 ,"EMCal && !SPD");
- fhNPileUpEvents->GetXaxis()->SetBinLabel(8 ,"!EMCal && !SPD");
- fOutputContainer->Add(fhNPileUpEvents);
-
- fhNPileUpEventsTriggerBC0 = new TH1F("hNPileUpEventsTriggerBC0","Number of events considered as pile-up, trigger cluster in BC=0", 8 , 0 , 8 ) ;
- fhNPileUpEventsTriggerBC0->SetYTitle("# events");
- fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(1 ,"SPD");
- fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(2 ,"Multi SPD");
- fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(3 ,"EMCal");
- fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(4 ,"EMCal || SPD");
- fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(5 ,"EMCal && SPD");
- fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(6 ,"!EMCal && SPD");
- fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(7 ,"EMCal && !SPD");
- fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(8 ,"!EMCal && !SPD");
- fOutputContainer->Add(fhNPileUpEventsTriggerBC0);
-
-
- fhTrackBCEvent = new TH1F("hTrackBCEvent", "Number of events with at least 1 track in a bunch crossing ", 19 , 0 , 19 ) ;
- fhTrackBCEvent->SetYTitle("# events");
- fhTrackBCEvent->SetXTitle("Bunch crossing");
- for(Int_t i = 1; i < 20; i++)
- fhTrackBCEvent->GetXaxis()->SetBinLabel(i ,Form("%d",i-10));
- fOutputContainer->Add(fhTrackBCEvent);
-
- fhTrackBCEventCut = new TH1F("hTrackBCEventCut", "Number of events with at least 1 track in a bunch crossing ", 19 , 0 , 19 ) ;
- fhTrackBCEventCut->SetYTitle("# events");
- fhTrackBCEventCut->SetXTitle("Bunch crossing");
- for(Int_t i = 1; i < 20; i++)
- fhTrackBCEventCut->GetXaxis()->SetBinLabel(i ,Form("%d",i-10));
- fOutputContainer->Add(fhTrackBCEventCut);
-
- fhPrimaryVertexBC = new TH1F("hPrimaryVertexBC", "Number of primary vertex per bunch crossing ", 41 , -20 , 20 ) ;
- fhPrimaryVertexBC->SetYTitle("# events");
- fhPrimaryVertexBC->SetXTitle("Bunch crossing");
- fOutputContainer->Add(fhPrimaryVertexBC);
-
- fhEMCalBCEvent = new TH1F("hEMCalBCEvent", "Number of events with at least 1 cluster in a bunch crossing ", 19 , 0 , 19 ) ;
- fhEMCalBCEvent->SetYTitle("# events");
- fhEMCalBCEvent->SetXTitle("Bunch crossing");
- for(Int_t i = 1; i < 20; i++)
- fhEMCalBCEvent->GetXaxis()->SetBinLabel(i ,Form("%d",i-10));
- fOutputContainer->Add(fhEMCalBCEvent);
-
- fhEMCalBCEventCut = new TH1F("hEMCalBCEventCut", "Number of events with at least 1 cluster in a bunch crossing", 19 , 0 , 19 ) ;
- fhEMCalBCEventCut->SetYTitle("# events");
- fhEMCalBCEventCut->SetXTitle("Bunch crossing");
- for(Int_t i = 1; i < 20; i++)
- fhEMCalBCEventCut->GetXaxis()->SetBinLabel(i ,Form("%d",i-10));
- fOutputContainer->Add(fhEMCalBCEventCut);
-
fhXVertex = new TH1F("hXVertex", " X vertex distribution" , 200 , -4 , 4 ) ;
fhXVertex->SetXTitle("v_{x} (cm)");
fOutputContainer->Add(fhXVertex);
fhZVertex = new TH1F("hZVertex", " Z vertex distribution" , 200 , -50 , 50 ) ;
fhZVertex->SetXTitle("v_{z} (cm)");
fOutputContainer->Add(fhZVertex);
-
- fhXVertexExotic = new TH1F("hXVertexExotic", " X vertex distribution in exotic events" , 200 , -4 , 4 ) ;
- fhXVertexExotic->SetXTitle("v_{x} (cm)");
- fOutputContainer->Add(fhXVertexExotic);
-
- fhYVertexExotic = new TH1F("hYVertexExotic", " Y vertex distribution in exotic events" , 200 , -4 , 4 ) ;
- fhYVertexExotic->SetXTitle("v_{y} (cm)");
- fOutputContainer->Add(fhYVertexExotic);
-
- fhZVertexExotic = new TH1F("hZVertexExotic", " Z vertex distribution in exotic events" , 200 , -50 , 50 ) ;
- fhZVertexExotic->SetXTitle("v_{z} (cm)");
- fOutputContainer->Add(fhZVertexExotic);
-
- fhTrackMult = new TH1F("hTrackMult", "Number of tracks per events" , 2000 , 0 , 2000 ) ;
- fhTrackMult->SetXTitle("# tracks");
- fOutputContainer->Add(fhTrackMult);
-
- fhPileUpClusterMult = new TH1F("hPileUpClusterMult", "Number of clusters per event with large time (|t| > 20 ns)" , 100 , 0 , 100 ) ;
- fhPileUpClusterMult->SetXTitle("# clusters");
- fOutputContainer->Add(fhPileUpClusterMult);
-
- fhPileUpClusterMultAndSPDPileUp = new TH1F("hPileUpClusterMultAndSPDPileUp", "Number of clusters per event with large time (|t| > 20 ns, events tagged as pile-up by SPD)" , 100 , 0 , 100 ) ;
- fhPileUpClusterMultAndSPDPileUp->SetXTitle("# clusters");
- fOutputContainer->Add(fhPileUpClusterMultAndSPDPileUp);
-
- fhNPileUpVertSPD = new TH1F ("hNPileUpVertSPD","N pile-up SPD vertex", 50,0,50);
- fhNPileUpVertSPD->SetYTitle("# vertex ");
- fOutputContainer->Add(fhNPileUpVertSPD);
-
- fhNPileUpVertTracks = new TH1F ("hNPileUpVertTracks","N pile-up Tracks vertex", 50,0,50);
- fhNPileUpVertTracks->SetYTitle("# vertex ");
- fOutputContainer->Add(fhNPileUpVertTracks);
-
+
fhCentrality = new TH1F("hCentrality","Number of events in centrality bin",100,0.,100) ;
fhCentrality->SetXTitle("Centrality bin");
- fOutputContainer->Add(fhCentrality) ;
+ fOutputContainer->Add(fhCentrality) ;
fhEventPlaneAngle=new TH1F("hEventPlaneAngle","Number of events in event plane",100,0.,TMath::Pi()) ;
fhEventPlaneAngle->SetXTitle("EP angle (rad)");
fOutputContainer->Add(fhEventPlaneAngle) ;
-
- if(fReader->IsSelectEventTimeStampOn())
+
+ fhTrackMult = new TH1F("hTrackMult", "Number of tracks per events" , 2000 , 0 , 2000 ) ;
+ fhTrackMult->SetXTitle("# tracks");
+ fOutputContainer->Add(fhTrackMult);
+
+ if(fFillDataControlHisto)
{
- fhTimeStampFraction = new TH1F("hTimeStampFraction","Fraction of events within a given time stamp range",150, -1, 2) ;
- fhTimeStampFraction->SetXTitle("fraction");
- fOutputContainer->Add(fhTimeStampFraction) ;
+ fhNExoticEvents = new TH1F("hNExoticEvents", "Number of analyzed events triggered by exotic cluster" , 1 , 0 , 1 ) ;
+ fhNExoticEvents->SetYTitle("# exotic events");
+ fOutputContainer->Add(fhNExoticEvents);
+
+ fhNEventsNoTriggerFound = new TH1F("hNEventsNoTriggerFound", "Number of analyzed events triggered but no trigger found" , 1 , 0 , 1 ) ;
+ fhNEventsNoTriggerFound->SetYTitle("# exotic events");
+ fOutputContainer->Add(fhNEventsNoTriggerFound);
+
+
+ Int_t nbin = 11;
+ Float_t minbin =-5.5;
+ Float_t maxbin = 5.5;
+ Int_t labelshift = 6;
+
+ fhClusterTriggerBCEventBC = new TH2F("hClusterTriggerBCEventBC", "Found trigger BC and Event BC",
+ nbin , minbin ,maxbin,4,0, 4) ;
+ fhClusterTriggerBCEventBC->SetXTitle("cluster trigger BC");
+ for(Int_t i = 0; i < 4; i++)
+ fhClusterTriggerBCEventBC->GetYaxis()->SetBinLabel(i+1 ,Form("BC/4=%d",i));
+ fhClusterTriggerBCEventBC->SetXTitle("cluster trigger BC");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCEventBC->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fhClusterTriggerBCEventBC->SetYTitle("Event BC%4");
+ fOutputContainer->Add(fhClusterTriggerBCEventBC);
+
+ fhClusterTriggerBCExoticEventBC = new TH2F("hClusterTriggerBCExoticEventBC", "Found exotic trigger BC and Event BC",
+ nbin , minbin ,maxbin,4,1, 4) ;
+ for(Int_t i = 0; i < 4; i++)
+ fhClusterTriggerBCExoticEventBC->GetYaxis()->SetBinLabel(i+1 ,Form("BC/4=%d",i));
+ fhClusterTriggerBCExoticEventBC->SetXTitle("cluster trigger BC");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCExoticEventBC->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fhClusterTriggerBCExoticEventBC->SetYTitle("Event BC%4");
+ fOutputContainer->Add(fhClusterTriggerBCExoticEventBC);
+
+ fhClusterTriggerBCEventBCUnMatch = new TH2F("hClusterTriggerBCEventBCUnMatch", "Found unmatched trigger BC and Event BC",
+ nbin , minbin ,maxbin,4,1, 4) ;
+ for(Int_t i = 0; i < 4; i++)
+ fhClusterTriggerBCEventBCUnMatch->GetYaxis()->SetBinLabel(i+1 ,Form("BC/4=%d",i));
+ fhClusterTriggerBCEventBCUnMatch->SetXTitle("cluster trigger BC");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCEventBCUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fhClusterTriggerBCEventBCUnMatch->SetYTitle("Event BC%4");
+ fOutputContainer->Add(fhClusterTriggerBCEventBCUnMatch);
+
+ fhClusterTriggerBCExoticEventBCUnMatch = new TH2F("hClusterTriggerExoticBCEventBCUnMatch", "Found unmatched trigger BC and Event BC",
+ nbin , minbin ,maxbin,4,1, 4) ;
+ for(Int_t i = 0; i < 4; i++)
+ fhClusterTriggerBCExoticEventBCUnMatch->GetYaxis()->SetBinLabel(i+1 ,Form("BC/4=%d",i));
+ fhClusterTriggerBCExoticEventBCUnMatch->SetXTitle("cluster trigger BC");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCExoticEventBCUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fhClusterTriggerBCExoticEventBCUnMatch->SetYTitle("Event BC%4");
+ fOutputContainer->Add(fhClusterTriggerBCExoticEventBCUnMatch);
+
+ fhClusterTriggerBC = new TH1F("hClusterTriggerBC",
+ "Number of analyzed events triggered by a cluster in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBC->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBC->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBC);
+
+ fhClusterTriggerBCExotic = new TH1F("hClusterTriggerBCExotic",
+ "Number of analyzed events triggered by a exotic cluster in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCExotic->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCExotic->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCExotic);
+
+
+ fhClusterTriggerBCBadCell = new TH1F("hClusterTriggerBCBadCell",
+ "Number of analyzed events triggered by a bad cell in a given BC",
+ nbin , minbin ,maxbin) ;
+
+ fhClusterTriggerBCBadCell->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCBadCell->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCBadCell);
+
+ fhClusterTriggerBCBadCellExotic = new TH1F("hClusterTriggerBCBadCellExotic",
+ "Number of analyzed events triggered by a bad cell & exotic cluster in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCBadCellExotic->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCBadCellExotic->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCBadCellExotic);
+
+ fhClusterTriggerBCBadCluster = new TH1F("hClusterTriggerBCBadCluster",
+ "Number of analyzed events triggered by a bad cluster in a given BC",
+ nbin , minbin ,maxbin) ;
+
+ fhClusterTriggerBCBadCluster->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCBadCluster->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCBadCluster);
+
+
+ fhClusterTriggerBCBadClusterExotic = new TH1F("hClusterTriggerBCBadClusterExotic",
+ "Number of analyzed events triggered by a bad cluster & exotic cluster in a given BC",
+ nbin , minbin ,maxbin) ;
+
+ fhClusterTriggerBCBadClusterExotic->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCBadClusterExotic->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCBadClusterExotic);
+
+ fhClusterTriggerBCUnMatch = new TH1F("hClusterTriggerBCUnMatch",
+ "Number of analyzed events triggered by a cluster (no trigger patch match) in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCUnMatch->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCUnMatch);
+
+ fhClusterTriggerBCExoticUnMatch = new TH1F("hClusterTriggerBCExoticUnMatch",
+ "Number of analyzed events triggered by a exotic cluster (no trigger patch match) in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCExoticUnMatch->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCExoticUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCExoticUnMatch);
+
+
+ fhClusterTriggerBCBadCellUnMatch = new TH1F("hClusterTriggerBCBadCellUnMatch",
+ "Number of analyzed events triggered by a bad cluster (no trigger patch match) in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCBadCellUnMatch->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCBadCellUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCBadCellUnMatch);
+
+
+ fhClusterTriggerBCBadCellExoticUnMatch = new TH1F("hClusterTriggerBCBadCellExoticUnMatch",
+ "Number of analyzed events triggered by a bad&exotic cluster (no trigger patch match) in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCBadCellExoticUnMatch->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCBadCellExoticUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCBadCellExoticUnMatch);
+
+
+ fhClusterTriggerBCBadClusterUnMatch = new TH1F("hClusterTriggerBCBadClusterUnMatch",
+ "Number of analyzed events triggered by a bad cluster (no trigger patch match) in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCBadClusterUnMatch->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCBadClusterUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCBadClusterUnMatch);
+
+
+ fhClusterTriggerBCBadClusterExoticUnMatch = new TH1F("hClusterTriggerBCBadClusterExoticUnMatch",
+ "Number of analyzed events triggered by a bad&exotic cluster (no trigger patch match) in a given BC",
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCBadClusterExoticUnMatch->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCBadClusterExoticUnMatch->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCBadClusterExoticUnMatch);
+
+ TString rematch[] = {"OpenTime","CheckNeighbours","Both"};
+ for(Int_t j = 0; j < 3; j++)
+ {
+ fhClusterTriggerBCUnMatchReMatch[j] = new TH1F(Form("hClusterTriggerBCUnMatch_ReMatch_%s",rematch[j].Data()),
+ Form("Number of analyzed events triggered by a cluster (no trigger patch match) in a given BC, re-match %s",rematch[j].Data()),
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCUnMatchReMatch[j]->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCUnMatchReMatch[j]->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCUnMatchReMatch[j]);
+
+ fhClusterTriggerBCExoticUnMatchReMatch[j] = new TH1F(Form("hClusterTriggerBCExoticUnMatch_ReMatch_%s",rematch[j].Data()),
+ Form("Number of analyzed events triggered by a exotic cluster (no trigger patch match) in a given BC, re-match %s",rematch[j].Data()),
+ nbin , minbin ,maxbin) ;
+ fhClusterTriggerBCExoticUnMatchReMatch[j]->SetYTitle("# events");
+ for(Int_t i = 1; i < 12; i++)
+ fhClusterTriggerBCExoticUnMatchReMatch[j]->GetXaxis()->SetBinLabel(i ,Form("BC%d",i-labelshift));
+ fOutputContainer->Add(fhClusterTriggerBCExoticUnMatchReMatch[j]);
+ }
+
+ fhNPileUpEvents = new TH1F("hNPileUpEvents", "Number of events considered as pile-up", 8 , 0 , 8 ) ;
+ fhNPileUpEvents->SetYTitle("# events");
+ fhNPileUpEvents->GetXaxis()->SetBinLabel(1 ,"SPD");
+ fhNPileUpEvents->GetXaxis()->SetBinLabel(2 ,"Multi SPD");
+ fhNPileUpEvents->GetXaxis()->SetBinLabel(3 ,"EMCal");
+ fhNPileUpEvents->GetXaxis()->SetBinLabel(4 ,"EMCal || SPD");
+ fhNPileUpEvents->GetXaxis()->SetBinLabel(5 ,"EMCal && SPD");
+ fhNPileUpEvents->GetXaxis()->SetBinLabel(6 ,"!EMCal && SPD");
+ fhNPileUpEvents->GetXaxis()->SetBinLabel(7 ,"EMCal && !SPD");
+ fhNPileUpEvents->GetXaxis()->SetBinLabel(8 ,"!EMCal && !SPD");
+ fOutputContainer->Add(fhNPileUpEvents);
+
+ fhNPileUpEventsTriggerBC0 = new TH1F("hNPileUpEventsTriggerBC0","Number of events considered as pile-up, trigger cluster in BC=0", 8 , 0 , 8 ) ;
+ fhNPileUpEventsTriggerBC0->SetYTitle("# events");
+ fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(1 ,"SPD");
+ fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(2 ,"Multi SPD");
+ fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(3 ,"EMCal");
+ fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(4 ,"EMCal || SPD");
+ fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(5 ,"EMCal && SPD");
+ fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(6 ,"!EMCal && SPD");
+ fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(7 ,"EMCal && !SPD");
+ fhNPileUpEventsTriggerBC0->GetXaxis()->SetBinLabel(8 ,"!EMCal && !SPD");
+ fOutputContainer->Add(fhNPileUpEventsTriggerBC0);
+
+
+ fhTrackBCEvent = new TH1F("hTrackBCEvent", "Number of events with at least 1 track in a bunch crossing ", 19 , 0 , 19 ) ;
+ fhTrackBCEvent->SetYTitle("# events");
+ fhTrackBCEvent->SetXTitle("Bunch crossing");
+ for(Int_t i = 1; i < 20; i++)
+ fhTrackBCEvent->GetXaxis()->SetBinLabel(i ,Form("%d",i-10));
+ fOutputContainer->Add(fhTrackBCEvent);
+
+ fhTrackBCEventCut = new TH1F("hTrackBCEventCut", "Number of events with at least 1 track in a bunch crossing ", 19 , 0 , 19 ) ;
+ fhTrackBCEventCut->SetYTitle("# events");
+ fhTrackBCEventCut->SetXTitle("Bunch crossing");
+ for(Int_t i = 1; i < 20; i++)
+ fhTrackBCEventCut->GetXaxis()->SetBinLabel(i ,Form("%d",i-10));
+ fOutputContainer->Add(fhTrackBCEventCut);
+
+ fhPrimaryVertexBC = new TH1F("hPrimaryVertexBC", "Number of primary vertex per bunch crossing ", 41 , -20 , 20 ) ;
+ fhPrimaryVertexBC->SetYTitle("# events");
+ fhPrimaryVertexBC->SetXTitle("Bunch crossing");
+ fOutputContainer->Add(fhPrimaryVertexBC);
+
+ fhEMCalBCEvent = new TH1F("hEMCalBCEvent", "Number of events with at least 1 cluster in a bunch crossing ", 19 , 0 , 19 ) ;
+ fhEMCalBCEvent->SetYTitle("# events");
+ fhEMCalBCEvent->SetXTitle("Bunch crossing");
+ for(Int_t i = 1; i < 20; i++)
+ fhEMCalBCEvent->GetXaxis()->SetBinLabel(i ,Form("%d",i-10));
+ fOutputContainer->Add(fhEMCalBCEvent);
+
+ fhEMCalBCEventCut = new TH1F("hEMCalBCEventCut", "Number of events with at least 1 cluster in a bunch crossing", 19 , 0 , 19 ) ;
+ fhEMCalBCEventCut->SetYTitle("# events");
+ fhEMCalBCEventCut->SetXTitle("Bunch crossing");
+ for(Int_t i = 1; i < 20; i++)
+ fhEMCalBCEventCut->GetXaxis()->SetBinLabel(i ,Form("%d",i-10));
+ fOutputContainer->Add(fhEMCalBCEventCut);
+
+ fhXVertexExotic = new TH1F("hXVertexExotic", " X vertex distribution in exotic events" , 200 , -4 , 4 ) ;
+ fhXVertexExotic->SetXTitle("v_{x} (cm)");
+ fOutputContainer->Add(fhXVertexExotic);
+
+ fhYVertexExotic = new TH1F("hYVertexExotic", " Y vertex distribution in exotic events" , 200 , -4 , 4 ) ;
+ fhYVertexExotic->SetXTitle("v_{y} (cm)");
+ fOutputContainer->Add(fhYVertexExotic);
+
+ fhZVertexExotic = new TH1F("hZVertexExotic", " Z vertex distribution in exotic events" , 200 , -50 , 50 ) ;
+ fhZVertexExotic->SetXTitle("v_{z} (cm)");
+ fOutputContainer->Add(fhZVertexExotic);
+
+ fhPileUpClusterMult = new TH1F("hPileUpClusterMult", "Number of clusters per event with large time (|t| > 20 ns)" , 100 , 0 , 100 ) ;
+ fhPileUpClusterMult->SetXTitle("# clusters");
+ fOutputContainer->Add(fhPileUpClusterMult);
+
+ fhPileUpClusterMultAndSPDPileUp = new TH1F("hPileUpClusterMultAndSPDPileUp", "Number of clusters per event with large time (|t| > 20 ns, events tagged as pile-up by SPD)" , 100 , 0 , 100 ) ;
+ fhPileUpClusterMultAndSPDPileUp->SetXTitle("# clusters");
+ fOutputContainer->Add(fhPileUpClusterMultAndSPDPileUp);
+
+ fhNPileUpVertSPD = new TH1F ("hNPileUpVertSPD","N pile-up SPD vertex", 50,0,50);
+ fhNPileUpVertSPD->SetYTitle("# vertex ");
+ fOutputContainer->Add(fhNPileUpVertSPD);
+
+ fhNPileUpVertTracks = new TH1F ("hNPileUpVertTracks","N pile-up Tracks vertex", 50,0,50);
+ fhNPileUpVertTracks->SetYTitle("# vertex ");
+ fOutputContainer->Add(fhNPileUpVertTracks);
+
+ if(fReader->IsSelectEventTimeStampOn())
+ {
+ fhTimeStampFraction = new TH1F("hTimeStampFraction","Fraction of events within a given time stamp range",150, -1, 2) ;
+ fhTimeStampFraction->SetXTitle("fraction");
+ fOutputContainer->Add(fhTimeStampFraction) ;
+ }
}
if(fScaleFactor > 0)
fhScaleFactor = new TH1F("hScaleFactor", "Number of merged output files" , 1 , 0 , 1 ) ;
fhScaleFactor->SetYTitle("scale factor");
- fhScaleFactor->SetBinContent(1,fScaleFactor); // Fill here
- fOutputContainer->Add(fhScaleFactor);
+ fhScaleFactor->SetBinContent(1,fScaleFactor); // Fill here
+ fOutputContainer->Add(fhScaleFactor);
}
if(!fAnalysisContainer || fAnalysisContainer->GetEntries()==0)
if(fMakeHisto) // Analysis with histograms as output on
{
- //Fill container with appropriate histograms
- TList * templist = ana ->GetCreateOutputObjects();
+ //Fill container with appropriate histograms
+ TList * templist = ana ->GetCreateOutputObjects();
templist->SetOwner(kFALSE); //Owner is fOutputContainer.
for(Int_t i = 0; i < templist->GetEntries(); i++)
{
//Add only to the histogram name the name of the task
- if( strcmp((templist->At(i))->ClassName(),"TObjString") )
+ if( strcmp((templist->At(i))->ClassName(),"TObjString") )
{
snprintf(newname,buffersize, "%s%s", (ana->GetAddedHistogramsStringToName()).Data(), (templist->At(i))->GetName());
//printf("name %s, new name %s\n",(templist->At(i))->GetName(),newname);
//___________________________________
void AliAnaCaloTrackCorrMaker::Init()
-{
+{
//Init container histograms and other common variables
// Fill the output list of histograms during the CreateOutputObjects stage.
//_____________________________________________
void AliAnaCaloTrackCorrMaker::InitParameters()
-{
+{
//Init data members
fMakeHisto = kTRUE;
- fMakeAOD = kTRUE;
+ fMakeAOD = kTRUE;
fAnaDebug = 0; // No debugging info displayed by default
}
//______________________________________________________________
void AliAnaCaloTrackCorrMaker::Print(const Option_t * opt) const
-{
+{
//Print some relevant parameters set for the analysis
if(! opt)
printf("Number of analysis tasks = %d\n", fAnalysisContainer->GetEntries()) ;
if(!strcmp("all",opt))
+ {
+ printf("Print analysis Tasks settings :\n") ;
+ for(Int_t iana = 0; iana<fAnalysisContainer->GetEntries(); iana++)
{
- printf("Print analysis Tasks settings :\n") ;
- for(Int_t iana = 0; iana<fAnalysisContainer->GetEntries(); iana++)
- {
- ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana))->Print("");
- }
-
- printf("Print analysis Reader settings :\n") ;
- fReader->Print("");
- printf("Print analysis Calorimeter Utils settings :\n") ;
- fCaloUtils->Print("");
-
+ ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana))->Print("");
+ }
+
+ printf("Print analysis Reader settings :\n") ;
+ fReader->Print("");
+ printf("Print analysis Calorimeter Utils settings :\n") ;
+ fCaloUtils->Print("");
+
}
-}
+}
-//_______________________________________________________________________
-void AliAnaCaloTrackCorrMaker::ProcessEvent(const Int_t iEntry,
- const char * currentFileName)
+//_____________________________________________________________________________________
+void AliAnaCaloTrackCorrMaker::ProcessEvent(Int_t iEntry, const char * currentFileName)
{
//Process analysis for this event
if(fMakeHisto && !fOutputContainer)
- {
- printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Histograms not initialized\n");
- abort();
- }
+ {
+ printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Histograms not initialized\n");
+ abort();
+ }
if(fAnaDebug >= 0 )
+ {
+ printf("*** AliAnaCaloTrackCorrMaker::ProcessEvent() Event %d *** \n",iEntry);
+ if(fAnaDebug > 1 )
{
- printf("*** AliAnaCaloTrackCorrMaker::ProcessEvent() Event %d *** \n",iEntry);
- if(fAnaDebug > 1 )
- {
- printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Current File Name : %s\n", currentFileName);
- //printf("fAODBranchList %p, entries %d\n",fAODBranchList,fAODBranchList->GetEntries());
- }
+ printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - Current File Name : %s\n", currentFileName);
+ //printf("fAODBranchList %p, entries %d\n",fAODBranchList,fAODBranchList->GetEntries());
}
+ }
//Each event needs an empty branch
TList * aodList = fReader->GetAODBranchList();
}
//Set geometry matrices before filling arrays, in case recalibration/position calculation etc is needed
- fCaloUtils->AccessGeometry(fReader->GetInputEvent());
+ fCaloUtils->AccessGeometry(fReader->GetInputEvent());
//Set the AODB calibration, bad channels etc. parameters at least once
- fCaloUtils->AccessOADB(fReader->GetInputEvent());
+ fCaloUtils->AccessOADB(fReader->GetInputEvent());
//Tell the reader to fill the data in the 3 detector lists
Bool_t ok = fReader->FillInputEvent(iEntry, currentFileName);
- Int_t triggerBC = fReader->GetTriggerClusterBC() ;
- Bool_t exotic = fReader->IsExoticEvent();
- Bool_t badCluster = fReader->IsBadCellTriggerEvent();
- Bool_t badCell = fReader->IsBadMaxCellTriggerEvent();
- Bool_t triggerMatch= fReader->IsTriggerMatched();
- Bool_t triggerBCOK = kTRUE;
- Int_t triggerId = fReader->GetTriggerClusterId() ;
-
- if(triggerId < 0)
- {
- //printf("Trigger id %d\n",triggerId);
- if(triggerId == -2)fhNEventsNoTriggerFound->Fill(0);
- triggerBCOK = kFALSE;
- }
-
- if(exotic)
- {
- fhNExoticEvents->Fill(0) ;
- Double_t v[3];
- fReader->GetInputEvent()->GetPrimaryVertex()->GetXYZ(v) ;
- fhXVertexExotic->Fill(v[0]);
- fhYVertexExotic->Fill(v[1]);
- fhZVertexExotic->Fill(v[2]);
- }
- //if(fReader->IsExoticEvent()) printf("Maker: EXOTIC Cluster trigger\n");
-
- if(triggerBCOK)
- {
- Int_t eventBC = fReader->GetInputEvent()->GetBunchCrossNumber();
- if(eventBC%4 < 0 || eventBC%4 > 3 )
- printf("AliAnaCaloTrackCorrMaker::ProcessEvent() - STRANGE: Trigger BC %d - Event BC %d, modulo4 %d \n",triggerBC,eventBC,eventBC%4);
-
- if(triggerMatch)
- {
- if (!exotic && !badCluster) fhClusterTriggerBC->Fill(triggerBC);
- else if( exotic && badCluster)
- {
- fhClusterTriggerBCBadClusterExotic->Fill(triggerBC);
- if(badCell) fhClusterTriggerBCBadCellExotic->Fill(triggerBC);
- }
- else if( exotic && !badCluster) fhClusterTriggerBCExotic->Fill(triggerBC);
- else if( badCluster && !exotic )
- {
- fhClusterTriggerBCBadCluster ->Fill(triggerBC);
- if(badCell) fhClusterTriggerBCBadCell->Fill(triggerBC);
- }
-
- if(!exotic) fhClusterTriggerBCEventBC ->Fill(triggerBC,eventBC%4);
- else fhClusterTriggerBCExoticEventBC->Fill(triggerBC,eventBC%4);
- }
- else
- {
- if (!exotic && !badCluster) fhClusterTriggerBCUnMatch->Fill(triggerBC);
- else if( exotic && badCluster)
- {
- fhClusterTriggerBCBadClusterExoticUnMatch->Fill(triggerBC);
- if(badCell) fhClusterTriggerBCBadCellExoticUnMatch ->Fill(triggerBC);
- }
- else if( exotic && !badCluster) fhClusterTriggerBCExoticUnMatch->Fill(triggerBC);
- else if( badCluster && !exotic )
- {
- fhClusterTriggerBCBadClusterUnMatch->Fill(triggerBC);
- if(badCell)fhClusterTriggerBCBadCellUnMatch->Fill(triggerBC);
- }
-
- if(!exotic) fhClusterTriggerBCEventBCUnMatch ->Fill(triggerBC,eventBC%4);
- else fhClusterTriggerBCExoticEventBCUnMatch->Fill(triggerBC,eventBC%4);
- }
- }
+ FillTriggerControlHistograms();
if(!ok)
- {
+ {
if(fAnaDebug >= 1 )printf("*** Skip event *** %d \n",iEntry);
fReader->ResetLists();
return ;
isTrigger = inputHandler->IsEventSelected() & fReader->GetEventTriggerMask();
}
+ // Init mag field for tracks in case of ESDs, not really necessary
+ if (!TGeoGlobalMagField::Instance()->GetField() && ((AliESDEvent*) fReader->GetInputEvent()))
+ ((AliESDEvent*)fReader->GetInputEvent())->InitMagneticField();
+
//Loop on analysis algorithms
if(fAnaDebug > 0 ) printf("*** Begin analysis *** \n");
Int_t nana = fAnalysisContainer->GetEntries() ;
for(Int_t iana = 0; iana < nana; iana++)
{
- AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
+ AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
ana->ConnectInputOutputAODBranches(); //Sets branches for each analysis
}
fReader->ResetLists();
-
+
// In case of mixing analysis, non triggered events are used,
// do not fill control histograms for a non requested triggered event
if(!fReader->IsEventTriggerAtSEOn() && !isTrigger)
- {
+ {
if(fAnaDebug > 0 ) printf("AliAnaCaloTrackMaker::ProcessEvent() - *** End analysis, MB for mixing *** \n");
return;
}
//__________________________________________________________
void AliAnaCaloTrackCorrMaker::Terminate(TList * outputList)
-{
+{
//Execute Terminate of analysis
//Do some final plots.
- if (!outputList)
+ if (!outputList)
{
Error("Terminate", "No output list");
return;
}
for(Int_t iana = 0; iana < fAnalysisContainer->GetEntries(); iana++)
- {
-
- AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
- if(ana->MakePlotsOn())ana->Terminate(outputList);
-
- }//Loop on analysis defined
+ {
+
+ AliAnaCaloTrackCorrBaseClass * ana = ((AliAnaCaloTrackCorrBaseClass *) fAnalysisContainer->At(iana)) ;
+ if(ana->MakePlotsOn())ana->Terminate(outputList);
+
+ }//Loop on analysis defined
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff