]>
Commit | Line | Data |
---|---|---|
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 | |
16 | // Michele Floris, CERN | |
17 | //------------------------------------------------------------------------- | |
18 | ||
19 | #include <AliAnalysisCuts.h> | |
20 | #include <TList.h> | |
21 | #include "TObjString.h" | |
22 | #include "AliVEvent.h" | |
23 | ||
24 | //#define VERBOSE_STAT | |
25 | ||
26 | class AliESDEvent; | |
27 | class TH2F; | |
28 | class TH1F; | |
29 | class TCollection; | |
30 | class AliTriggerAnalysis; | |
31 | class AliAnalysisTaskSE; | |
32 | ||
33 | ||
34 | class AliPhysicsSelection : public AliAnalysisCuts | |
35 | { | |
36 | public: | |
37 | ||
38 | public: | |
39 | enum {kStatTriggerClass=1,kStatHWTrig,kStatV0ABG,kStatV0CBG,kStatMB1,kStatMB1Prime,kStatFMD,kStatFO1,kStatFO2,kStatFO2L1,kStatV0A,kStatV0C,kStatZDCA,kStatZDCC,kStatZDCAC,kStatZDCTime,kStatV0,kStatOffline,kStatBG,kStatAccepted}; | |
40 | ||
41 | #ifdef VERBOSE_STAT | |
42 | enum {kStatRowAllB=0, kStatRowAllAC, kStatRowAllE, kStatRowBG,kStatRowAcc,kStatRowBGFrac,kStatRowAccFrac,kStatRowErrGoodFrac,kStatRowGoodFrac,kStatRowErrGood,kStatRowGood}; // offset wrt fBGStatOffset | |
43 | #else | |
44 | enum {kStatRowAllB=0, kStatRowAllAC, kStatRowAllE, kStatRowBG, kStatRowAcc,kStatRowGood}; // offset wrt fBGStatOffset | |
45 | #endif | |
46 | ||
47 | enum {kStatIdxAll=0,kStatIdxBin0=1}; | |
48 | enum ETriggerLogic { kCINT1 = 0, kCMBS2A, kCMBS2C, kCMBAC, kCMBACS2, kHighMultL1 }; | |
49 | ||
50 | typedef Bool_t (*Bin0Callback_t)(const AliESDEvent *); | |
51 | ||
52 | AliPhysicsSelection(); | |
53 | virtual ~AliPhysicsSelection(); | |
54 | ||
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;} | |
59 | ||
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); | |
64 | ||
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;} | |
69 | ||
70 | void AddBackgroundIdentification(AliAnalysisCuts* background) { fBackgroundIdentification = background; } | |
71 | ||
72 | virtual void Print(Option_t* option = "") const; | |
73 | virtual Long64_t Merge(TCollection* list); | |
74 | void SaveHistograms(const char* folder = 0); | |
75 | ||
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; } | |
80 | ||
81 | AliTriggerAnalysis* GetTriggerAnalysis(Int_t i = 0) { return (fTriggerAnalysis.GetEntries() > 0) ? (AliTriggerAnalysis*) fTriggerAnalysis.At(i) : 0; } | |
82 | ||
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; | |
86 | ||
87 | void SetBIFactors(const AliESDEvent * aESD); | |
88 | ||
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 | |
94 | ||
95 | protected: | |
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; | |
102 | ||
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 | |
107 | ||
108 | TList fTriggerAnalysis; // list of offline trigger objects (several are needed to keep the control histograms separate per trigger class) | |
109 | ||
110 | AliAnalysisCuts* fBackgroundIdentification; // class that performs additional background identification | |
111 | ||
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 | |
115 | ||
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) | |
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 | |
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 | |
125 | ||
126 | ||
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. | |
132 | ||
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 | |
135 | ||
136 | Bool_t fIsPP; // True if processing pp run, false if heavy ion | |
137 | ||
138 | ClassDef(AliPhysicsSelection, 12) | |
139 | ||
140 | private: | |
141 | AliPhysicsSelection(const AliPhysicsSelection&); | |
142 | AliPhysicsSelection& operator=(const AliPhysicsSelection&); | |
143 | }; | |
144 | ||
145 | #endif |