1 /* $Id: AliPhysicsSelection.h 35782 2009-10-22 11:54:31Z jgrosseo $ */
3 #ifndef ALIPHYSICSSELECTION_H
4 #define ALIPHYSICSSELECTION_H
6 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
7 * See cxx source for full Copyright notice */
9 //-------------------------------------------------------------------------
10 // Implementation of Class AliPhysicsSelection
12 // This class selects collision candidates from data runs, applying selection cuts on triggers
13 // and background rejection based on the content of the ESD
15 // Origin: Jan Fiete Grosse-Oetringhaus, CERN
16 // Michele Floris, CERN
17 //-------------------------------------------------------------------------
19 #include <AliAnalysisCuts.h>
21 #include "TObjString.h"
22 #include "AliVEvent.h"
24 //#define VERBOSE_STAT
30 class AliTriggerAnalysis;
31 class AliAnalysisTaskSE;
34 class AliPhysicsSelection : public AliAnalysisCuts
39 enum {kStatTriggerClass=1,kStatHWTrig,kStatV0ABG,kStatV0CBG,kStatMB1,kStatMB1Prime,kStatFMD,kStatFO1,kStatFO2,kStatFO2L1,kStatV0A,kStatV0C,kStatZDCA,kStatZDCC,kStatZDCAC,kStatZDCTime,kStatV0,kStatOffline,kStatBG,kStatAccepted};
42 enum {kStatRowAllB=0, kStatRowAllAC, kStatRowAllE, kStatRowBG,kStatRowAcc,kStatRowBGFrac,kStatRowAccFrac,kStatRowErrGoodFrac,kStatRowGoodFrac,kStatRowErrGood,kStatRowGood}; // offset wrt fBGStatOffset
44 enum {kStatRowAllB=0, kStatRowAllAC, kStatRowAllE, kStatRowBG, kStatRowAcc,kStatRowGood}; // offset wrt fBGStatOffset
47 enum {kStatIdxAll=0,kStatIdxBin0=1};
48 enum ETriggerLogic { kCINT1 = 0, kCMBS2A, kCMBS2C, kCMBAC, kCMBACS2, kHighMultL1 };
50 typedef Bool_t (*Bin0Callback_t)(const AliESDEvent *);
52 AliPhysicsSelection();
53 virtual ~AliPhysicsSelection();
55 // AliAnalysisCuts interface
56 virtual UInt_t GetSelectionMask(const TObject* obj) { return IsCollisionCandidate((const AliESDEvent*) obj); }
57 virtual Bool_t IsSelected(TList*) { return kFALSE; }
58 virtual Bool_t IsSelected(TObject*) {return kFALSE;}
60 Int_t GetCurrentRun() const {return fCurrentRun;}
61 UInt_t IsCollisionCandidate(const AliESDEvent* aEsd);
62 Bool_t Initialize(const AliESDEvent* aEsd);
63 Bool_t Initialize(Int_t runNumber, Bool_t pp);
65 void SetAnalyzeMC(Bool_t flag = kTRUE) { fMC = flag; }
66 void SetSkipTriggerClassSelection(Bool_t flag = kTRUE) { fSkipTriggerClassSelection = flag; }
67 void SetSkipV0(Bool_t flag=kTRUE) { fSkipV0 = flag;}
68 void SetSkipZDCTime(Bool_t flag=kTRUE) { fSkipZDCTime = flag;}
70 void AddBackgroundIdentification(AliAnalysisCuts* background) { fBackgroundIdentification = background; }
72 virtual void Print(Option_t* option = "") const;
73 virtual Long64_t Merge(TCollection* list);
74 void SaveHistograms(const char* folder = 0);
76 const TList* GetCollisionTriggerClasses() const { return &fCollTrigClasses; }
77 const TList* GetBGTriggerClasses() const { return &fBGTrigClasses; }
78 void AddCollisionTriggerClass(const char* className){ fCollTrigClasses.Add(new TObjString(className)); fUsingCustomClasses = kTRUE; }
79 void AddBGTriggerClass(const char* className) { fBGTrigClasses.Add(new TObjString(className)); fUsingCustomClasses = kTRUE; }
81 AliTriggerAnalysis* GetTriggerAnalysis(Int_t i = 0) { return (fTriggerAnalysis.GetEntries() > 0) ? (AliTriggerAnalysis*) fTriggerAnalysis.At(i) : 0; }
83 const TH2F* GetStatisticsHistogram(Int_t idx=kStatIdxAll) const { return fHistStatistics[idx]; }
84 const TH2F* GetBunchCrossingHistogram() const { return fHistBunchCrossing; }
85 virtual TObject *GetStatistics(Option_t *option) const;
87 void SetBIFactors(const AliESDEvent * aESD);
89 void SetUseBXNumbers(Bool_t flag = kTRUE) {fUseBXNumbers = flag;}
90 void SetComputeBG (UInt_t flag = AliVEvent::kMB) {fComputeBG = flag; if(flag) fUseBXNumbers = flag;}
91 void SetUseMuonTriggers(Bool_t) { Printf("ERROR: Use of AliPhysicsSelection::SetUseMuonTriggers is deprecated. Use AliVEvent::kMUON selection instead."); }
92 void SetBin0Callback( const char * cb) {fBin0CallBack = cb;}
93 void SetBin0CallbackViaPointer( Bin0Callback_t cb) {fBin0CallBackPointer = cb;}// WARNING: THIS SHOULD NOT BE USED, WILL BE REMOVED SOON
96 UInt_t CheckTriggerClass(const AliESDEvent* aEsd, const char* trigger, Int_t& triggerLogic) const;
97 Int_t GetTriggerScheme(UInt_t runNumber) const;
98 const char * GetBXIDs(UInt_t runNumber, const char * trigger ) ;
99 const char * GetFillingScheme(UInt_t runNumber) ;
100 TH2F * BookHistStatistics(const char * tag) ;
101 Int_t GetStatRow(const char * triggerBXClass, UInt_t offlineTriggerType, UInt_t ** rowIDs) const;
103 Int_t fCurrentRun; // run number for which the object is initialized
104 Bool_t fMC; // flag if MC is analyzed
105 TList fCollTrigClasses; // trigger class identifying collision candidates
106 TList fBGTrigClasses; // trigger classes identifying background events
108 TList fTriggerAnalysis; // list of offline trigger objects (several are needed to keep the control histograms separate per trigger class)
110 AliAnalysisCuts* fBackgroundIdentification; // class that performs additional background identification
112 TH2F* fHistStatistics[2]; // how many events are cut away why {all,bin 0}
113 TH2F* fHistBunchCrossing; // histograms of accepted bunch crossing numbers
114 TH1F* fHistTriggerPattern; // Pattern of the individual detectors in the MB1 trigger. Can reveal inconsistencies/inefficiencies in the trigger
116 Bool_t fSkipTriggerClassSelection; // flag that determines if the trigger classs selection is skipped
117 Bool_t fUsingCustomClasses; // flag that is set if costum trigger classes are defined
118 Bool_t fSkipV0; // ignore information from v0
119 Bool_t fSkipZDCTime; // ignore ZDC timing cut (used in HI)
121 Float_t fBIFactorA; // ratio of interacting over non interacting bunch intensities for beam 1
122 Float_t fBIFactorC; // ratio of interacting over non interacting bunch intensities for beam 2
123 Float_t fBIFactorAC; // ratio of interacting over non interacting bunch intensities for beam 1&2:
124 // for some trigger classes we do not distinguish betweem 1 and 2
127 UInt_t fComputeBG; // Switch on computation of background and filling of relevant stat table entries. If you enable this you can only process one run at a time (the relative bunch intensity used to compute this chages from run to run); This is a mask of the trigger types to be used in the selection, as defined in AliVEvent
128 Int_t fBGStatOffset; // Offset of the BG statistics computed at the end of the processing. It depends on how many trigger classes you have selected
129 Bool_t fUseBXNumbers; // Explicitely select "good" bunch crossing numbers (exclude pilot, afterpulses and fakes). If you anable this you can only process runs within the same filling scheme.
130 Bool_t fUseMuonTriggers; // if true, also use the muon triggers
131 TString fFillingScheme; // stores the filling scheme of the current run.
133 TString fBin0CallBack; // callback used to determine if an event is in the bin0 (name of the task where the callback is implemented);
134 Bin0Callback_t fBin0CallBackPointer; //! don't stream this. TO BE REMOVED SOON
136 Bool_t fIsPP; // True if processing pp run, false if heavy ion
138 ClassDef(AliPhysicsSelection, 12)
141 AliPhysicsSelection(const AliPhysicsSelection&);
142 AliPhysicsSelection& operator=(const AliPhysicsSelection&);