]>
Commit | Line | Data |
---|---|---|
61899827 | 1 | /* $Id: AliPhysicsSelection.h 35782 2009-10-22 11:54:31Z jgrosseo $ */ |
2 | ||
3 | #ifndef ALIPHYSICSSELECTION_H | |
4 | #define ALIPHYSICSSELECTION_H | |
5 | ||
6 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
7 | * See cxx source for full Copyright notice */ | |
8 | ||
9 | //------------------------------------------------------------------------- | |
10 | // Implementation of Class AliPhysicsSelection | |
11 | // | |
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 | |
14 | // | |
15 | // Origin: Jan Fiete Grosse-Oetringhaus, CERN | |
73cc8654 | 16 | // Michele Floris, CERN |
61899827 | 17 | //------------------------------------------------------------------------- |
18 | ||
296dd262 | 19 | #include <AliAnalysisCuts.h> |
61899827 | 20 | #include <TList.h> |
91bea6e7 | 21 | #include "TObjString.h" |
141265a2 | 22 | #include "AliVEvent.h" |
61899827 | 23 | |
8dec6e35 | 24 | //#define VERBOSE_STAT |
73cc8654 | 25 | |
61899827 | 26 | class AliESDEvent; |
17ba346c | 27 | class TH2F; |
73cc8654 | 28 | class TH1F; |
61899827 | 29 | class TCollection; |
30 | class AliTriggerAnalysis; | |
85c71ba7 | 31 | class AliAnalysisTaskSE; |
61899827 | 32 | |
141265a2 | 33 | |
296dd262 | 34 | class AliPhysicsSelection : public AliAnalysisCuts |
61899827 | 35 | { |
9b0cb3c3 | 36 | public: |
85c71ba7 | 37 | |
5f337f3d | 38 | public: |
869f9767 | 39 | enum {kStatTriggerClass=1,kStatHWTrig,kStatV0ABG,kStatV0CBG,kStatMB1,kStatMB1Prime,kStatFMD,kStatFO1,kStatFO2,kStatFO2L1,kStatV0A,kStatV0C,kStatZDCA,kStatZDCC,kStatZDCAC,kStatZDCTime,kStatV0,kStatOffline,kStatBG,kStatAccepted}; |
85c71ba7 | 40 | |
41 | #ifdef VERBOSE_STAT | |
8dec6e35 | 42 | enum {kStatRowAllB=0, kStatRowAllAC, kStatRowAllE, kStatRowBG,kStatRowAcc,kStatRowBGFrac,kStatRowAccFrac,kStatRowErrGoodFrac,kStatRowGoodFrac,kStatRowErrGood,kStatRowGood}; // offset wrt fBGStatOffset |
85c71ba7 | 43 | #else |
8dec6e35 | 44 | enum {kStatRowAllB=0, kStatRowAllAC, kStatRowAllE, kStatRowBG, kStatRowAcc,kStatRowGood}; // offset wrt fBGStatOffset |
85c71ba7 | 45 | #endif |
46 | ||
47 | enum {kStatIdxAll=0,kStatIdxBin0=1}; | |
05e1da5d | 48 | enum ETriggerLogic { kCINT1 = 0, kCMBS2A, kCMBS2C, kCMBAC, kCMBACS2, kHighMultL1 }; |
85c71ba7 | 49 | |
ff097e3f | 50 | typedef Bool_t (*Bin0Callback_t)(const AliESDEvent *); |
85c71ba7 | 51 | |
52 | AliPhysicsSelection(); | |
53 | virtual ~AliPhysicsSelection(); | |
61899827 | 54 | |
85c71ba7 | 55 | // AliAnalysisCuts interface |
0c6c629b | 56 | virtual UInt_t GetSelectionMask(const TObject* obj) { return IsCollisionCandidate((const AliESDEvent*) obj); } |
85c71ba7 | 57 | virtual Bool_t IsSelected(TList*) { return kFALSE; } |
0c6c629b | 58 | virtual Bool_t IsSelected(TObject*) {return kFALSE;} |
296dd262 | 59 | |
518b6d59 | 60 | Int_t GetCurrentRun() const {return fCurrentRun;} |
0c6c629b | 61 | UInt_t IsCollisionCandidate(const AliESDEvent* aEsd); |
1107e2ce | 62 | Bool_t Initialize(const AliESDEvent* aEsd); |
63 | Bool_t Initialize(Int_t runNumber, Bool_t pp); | |
decf6fd4 | 64 | |
85c71ba7 | 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;} | |
8dec6e35 | 68 | void SetSkipZDCTime(Bool_t flag=kTRUE) { fSkipZDCTime = flag;} |
91bea6e7 | 69 | |
85c71ba7 | 70 | void AddBackgroundIdentification(AliAnalysisCuts* background) { fBackgroundIdentification = background; } |
296dd262 | 71 | |
85c71ba7 | 72 | virtual void Print(Option_t* option = "") const; |
73 | virtual Long64_t Merge(TCollection* list); | |
8dec6e35 | 74 | void SaveHistograms(const char* folder = 0); |
61899827 | 75 | |
85c71ba7 | 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; } | |
91bea6e7 | 80 | |
c3a12f21 | 81 | AliTriggerAnalysis* GetTriggerAnalysis(Int_t i = 0) { return (fTriggerAnalysis.GetEntries() > 0) ? (AliTriggerAnalysis*) fTriggerAnalysis.At(i) : 0; } |
61899827 | 82 | |
85c71ba7 | 83 | const TH2F* GetStatisticsHistogram(Int_t idx=kStatIdxAll) const { return fHistStatistics[idx]; } |
84 | const TH2F* GetBunchCrossingHistogram() const { return fHistBunchCrossing; } | |
2d45a1a7 | 85 | virtual TObject *GetStatistics(Option_t *option) const; |
17ba346c | 86 | |
141265a2 | 87 | void SetBIFactors(const AliESDEvent * aESD); |
85c71ba7 | 88 | |
89 | void SetUseBXNumbers(Bool_t flag = kTRUE) {fUseBXNumbers = flag;} | |
141265a2 | 90 | void SetComputeBG (UInt_t flag = AliVEvent::kMB) {fComputeBG = flag; if(flag) fUseBXNumbers = flag;} |
85fc0208 | 91 | void SetUseMuonTriggers(Bool_t) { Printf("ERROR: Use of AliPhysicsSelection::SetUseMuonTriggers is deprecated. Use AliVEvent::kMUON selection instead."); } |
85c71ba7 | 92 | void SetBin0Callback( const char * cb) {fBin0CallBack = cb;} |
ff097e3f | 93 | void SetBin0CallbackViaPointer( Bin0Callback_t cb) {fBin0CallBackPointer = cb;}// WARNING: THIS SHOULD NOT BE USED, WILL BE REMOVED SOON |
73cc8654 | 94 | |
85c71ba7 | 95 | protected: |
c2ba5a61 | 96 | UInt_t CheckTriggerClass(const AliESDEvent* aEsd, const char* trigger, Int_t& triggerLogic) const; |
85c71ba7 | 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) ; | |
85c71ba7 | 100 | TH2F * BookHistStatistics(const char * tag) ; |
141265a2 | 101 | Int_t GetStatRow(const char * triggerBXClass, UInt_t offlineTriggerType, UInt_t ** rowIDs) const; |
85c71ba7 | 102 | |
85c71ba7 | 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 | |
296dd262 | 107 | |
85c71ba7 | 108 | TList fTriggerAnalysis; // list of offline trigger objects (several are needed to keep the control histograms separate per trigger class) |
61899827 | 109 | |
85c71ba7 | 110 | AliAnalysisCuts* fBackgroundIdentification; // class that performs additional background identification |
61899827 | 111 | |
85c71ba7 | 112 | TH2F* fHistStatistics[2]; // how many events are cut away why {all,bin 0} |
113 | TH2F* fHistBunchCrossing; // histograms of accepted bunch crossing numbers | |
73cc8654 | 114 | TH1F* fHistTriggerPattern; // Pattern of the individual detectors in the MB1 trigger. Can reveal inconsistencies/inefficiencies in the trigger |
61899827 | 115 | |
85c71ba7 | 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 | |
8dec6e35 | 119 | Bool_t fSkipZDCTime; // ignore ZDC timing cut (used in HI) |
85c71ba7 | 120 | |
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 | |
141265a2 | 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 | |
85c71ba7 | 125 | |
85c71ba7 | 126 | |
141265a2 | 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 | |
73cc8654 | 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. | |
85c71ba7 | 132 | |
133 | TString fBin0CallBack; // callback used to determine if an event is in the bin0 (name of the task where the callback is implemented); | |
ff097e3f | 134 | Bin0Callback_t fBin0CallBackPointer; //! don't stream this. TO BE REMOVED SOON |
91bea6e7 | 135 | |
869f9767 | 136 | Bool_t fIsPP; // True if processing pp run, false if heavy ion |
137 | ||
8dec6e35 | 138 | ClassDef(AliPhysicsSelection, 12) |
61899827 | 139 | |
85c71ba7 | 140 | private: |
141 | AliPhysicsSelection(const AliPhysicsSelection&); | |
142 | AliPhysicsSelection& operator=(const AliPhysicsSelection&); | |
61899827 | 143 | }; |
144 | ||
145 | #endif |