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