]>
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" |
61899827 | 22 | |
73cc8654 | 23 | #define VERBOSE_STAT |
24 | ||
61899827 | 25 | class AliESDEvent; |
17ba346c | 26 | class TH2F; |
73cc8654 | 27 | class TH1F; |
61899827 | 28 | class TCollection; |
29 | class AliTriggerAnalysis; | |
85c71ba7 | 30 | class AliAnalysisTaskSE; |
61899827 | 31 | |
296dd262 | 32 | class AliPhysicsSelection : public AliAnalysisCuts |
61899827 | 33 | { |
85c71ba7 | 34 | |
35 | enum {kStatTriggerClass=1,kStatHWTrig,kStatV0ABG,kStatV0CBG,kStatMB1,kStatMB1Prime,kStatFMD,kStatFO1,kStatFO2,kStatV0A,kStatV0C,kStatSSD1,kStatFO1AndV0,kStatV0,kStatAny2Hits,kStatOffline,kStatBG,kStatAccepted}; | |
36 | ||
37 | #ifdef VERBOSE_STAT | |
ca234d7a | 38 | enum {kStatRowBG=5,kStatRowAcc,kStatRowBGFrac,kStatRowAccFrac,kStatRowErrGoodFrac,kStatRowGoodFrac,kStatRowErrGood,kStatRowGood}; |
85c71ba7 | 39 | #else |
40 | enum {kStatRowBG=5,kStatRowAcc,kStatRowGood}; | |
41 | #endif | |
42 | ||
43 | enum {kStatIdxAll=0,kStatIdxBin0=1}; | |
44 | ||
45 | ||
46 | public: | |
ff097e3f | 47 | typedef Bool_t (*Bin0Callback_t)(const AliESDEvent *); |
85c71ba7 | 48 | |
49 | AliPhysicsSelection(); | |
50 | virtual ~AliPhysicsSelection(); | |
61899827 | 51 | |
85c71ba7 | 52 | // AliAnalysisCuts interface |
53 | virtual Bool_t IsSelected(TObject* obj) { return IsCollisionCandidate((const AliESDEvent*) obj); } | |
54 | virtual Bool_t IsSelected(TList*) { return kFALSE; } | |
296dd262 | 55 | |
85c71ba7 | 56 | Bool_t IsCollisionCandidate(const AliESDEvent* aEsd); |
57 | Bool_t Initialize(Int_t runNumber); | |
decf6fd4 | 58 | |
85c71ba7 | 59 | void SetAnalyzeMC(Bool_t flag = kTRUE) { fMC = flag; } |
60 | void SetSkipTriggerClassSelection(Bool_t flag = kTRUE) { fSkipTriggerClassSelection = flag; } | |
61 | void SetSkipV0(Bool_t flag=kTRUE) { fSkipV0 = flag;} | |
91bea6e7 | 62 | |
85c71ba7 | 63 | void AddBackgroundIdentification(AliAnalysisCuts* background) { fBackgroundIdentification = background; } |
296dd262 | 64 | |
85c71ba7 | 65 | virtual void Print(Option_t* option = "") const; |
66 | virtual Long64_t Merge(TCollection* list); | |
67 | void SaveHistograms(const char* folder = 0) const; | |
61899827 | 68 | |
85c71ba7 | 69 | const TList* GetCollisionTriggerClasses() const { return &fCollTrigClasses; } |
70 | const TList* GetBGTriggerClasses() const { return &fBGTrigClasses; } | |
71 | void AddCollisionTriggerClass(const char* className){ fCollTrigClasses.Add(new TObjString(className)); fUsingCustomClasses = kTRUE; } | |
72 | void AddBGTriggerClass(const char* className) { fBGTrigClasses.Add(new TObjString(className)); fUsingCustomClasses = kTRUE; } | |
91bea6e7 | 73 | |
c3a12f21 | 74 | AliTriggerAnalysis* GetTriggerAnalysis(Int_t i = 0) { return (fTriggerAnalysis.GetEntries() > 0) ? (AliTriggerAnalysis*) fTriggerAnalysis.At(i) : 0; } |
61899827 | 75 | |
85c71ba7 | 76 | const TH2F* GetStatisticsHistogram(Int_t idx=kStatIdxAll) const { return fHistStatistics[idx]; } |
77 | const TH2F* GetBunchCrossingHistogram() const { return fHistBunchCrossing; } | |
17ba346c | 78 | |
85c71ba7 | 79 | void SetBIFactors(Int_t run); |
80 | ||
81 | void SetUseBXNumbers(Bool_t flag = kTRUE) {fUseBXNumbers = flag;} | |
82 | void SetComputeBG (Bool_t flag = kTRUE) {fComputeBG = flag; if(flag) fUseBXNumbers = flag;} | |
73cc8654 | 83 | void SetUseMuonTriggers(Bool_t flag = kTRUE) {fUseMuonTriggers = flag;} |
85c71ba7 | 84 | void SetBin0Callback( const char * cb) {fBin0CallBack = cb;} |
ff097e3f | 85 | void SetBin0CallbackViaPointer( Bin0Callback_t cb) {fBin0CallBackPointer = cb;}// WARNING: THIS SHOULD NOT BE USED, WILL BE REMOVED SOON |
73cc8654 | 86 | |
85c71ba7 | 87 | |
88 | protected: | |
89 | Bool_t CheckTriggerClass(const AliESDEvent* aEsd, const char* trigger) const; | |
90 | Int_t GetTriggerScheme(UInt_t runNumber) const; | |
91 | const char * GetBXIDs(UInt_t runNumber, const char * trigger ) ; | |
92 | const char * GetFillingScheme(UInt_t runNumber) ; | |
93 | Int_t GetRatioBBBE(Int_t runNumber); | |
94 | TH2F * BookHistStatistics(const char * tag) ; | |
95 | ||
296dd262 | 96 | |
85c71ba7 | 97 | Int_t fCurrentRun; // run number for which the object is initialized |
98 | Bool_t fMC; // flag if MC is analyzed | |
99 | TList fCollTrigClasses; // trigger class identifying collision candidates | |
100 | TList fBGTrigClasses; // trigger classes identifying background events | |
296dd262 | 101 | |
85c71ba7 | 102 | TList fTriggerAnalysis; // list of offline trigger objects (several are needed to keep the control histograms separate per trigger class) |
61899827 | 103 | |
85c71ba7 | 104 | AliAnalysisCuts* fBackgroundIdentification; // class that performs additional background identification |
61899827 | 105 | |
85c71ba7 | 106 | TH2F* fHistStatistics[2]; // how many events are cut away why {all,bin 0} |
107 | TH2F* fHistBunchCrossing; // histograms of accepted bunch crossing numbers | |
73cc8654 | 108 | TH1F* fHistTriggerPattern; // Pattern of the individual detectors in the MB1 trigger. Can reveal inconsistencies/inefficiencies in the trigger |
61899827 | 109 | |
85c71ba7 | 110 | Bool_t fSkipTriggerClassSelection; // flag that determines if the trigger classs selection is skipped |
111 | Bool_t fUsingCustomClasses; // flag that is set if costum trigger classes are defined | |
112 | Bool_t fSkipV0; // ignore information from v0 | |
113 | ||
114 | Float_t fBIFactorA; // ratio of interacting over non interacting bunch intensities for beam 1 | |
115 | Float_t fBIFactorC; // ratio of interacting over non interacting bunch intensities for beam 2 | |
116 | ||
73cc8654 | 117 | Int_t fRatioBEEE; // ratio between the number of BX in the Beam-Empty and the Empty-Empty. Depends on the filling and on the trigger scheme |
85c71ba7 | 118 | |
73cc8654 | 119 | Bool_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) |
120 | 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. | |
121 | Bool_t fUseMuonTriggers; // if true, also use the muon triggers | |
122 | TString fFillingScheme; // stores the filling scheme of the current run. | |
85c71ba7 | 123 | |
124 | TString fBin0CallBack; // callback used to determine if an event is in the bin0 (name of the task where the callback is implemented); | |
ff097e3f | 125 | Bin0Callback_t fBin0CallBackPointer; //! don't stream this. TO BE REMOVED SOON |
91bea6e7 | 126 | |
73cc8654 | 127 | ClassDef(AliPhysicsSelection, 8) |
61899827 | 128 | |
85c71ba7 | 129 | private: |
130 | AliPhysicsSelection(const AliPhysicsSelection&); | |
131 | AliPhysicsSelection& operator=(const AliPhysicsSelection&); | |
61899827 | 132 | }; |
133 | ||
134 | #endif |