1 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. */
2 /* See cxx source for full Copyright notice */
7 // origin: Mikolaj Krzewicki (mikolaj.krzewicki@cern.ch)
9 #ifndef ALIFLOWEVENTCUTS_H
10 #define ALIFLOWEVENTCUTS_H
21 #include "AliTriggerAnalysis.h"
22 #include "AliFlowTrackCuts.h"
23 #include "AliFlowEventSimpleCuts.h"
25 class AliFlowEventCuts : public AliFlowEventSimpleCuts {
28 enum refMultMethod { kTPConly, kSPDtracklets, kVZERO, kV0=kVZERO, kSPD1clusters };
31 AliFlowEventCuts(const char* name, const char* title = "AliFlowEventCuts");
32 AliFlowEventCuts(const AliFlowEventCuts& someCuts);
33 AliFlowEventCuts& operator=(const AliFlowEventCuts& someCuts);
34 virtual ~AliFlowEventCuts();
36 virtual Bool_t IsSelected(TObject* obj, TObject *objmc);
38 Bool_t PassesCuts(AliVEvent* event, AliMCEvent *mcevent);
40 static AliFlowEventCuts* StandardCuts();
42 void SetNumberOfTracksMax(Int_t value) {fNumberOfTracksMax=value;fCutNumberOfTracks=kTRUE;}
43 void SetNumberOfTracksMin(Int_t value) {fNumberOfTracksMin=value;fCutNumberOfTracks=kTRUE;}
44 void SetNumberOfTracksRange(Int_t min, Int_t max) {fNumberOfTracksMin=min;fNumberOfTracksMax=max;fCutNumberOfTracks=kTRUE;}
45 void SetRefMultMax(Int_t value) {fRefMultMax=value;fCutRefMult=kTRUE;}
46 void SetRefMultMin(Int_t value) {fRefMultMin=value;fCutRefMult=kTRUE;}
47 void SetRefMultRange(Int_t min, Int_t max) {fRefMultMin=min;fRefMultMax=max;fCutRefMult=kTRUE;}
48 void SetImpactParameterMax(Double_t value) {fImpactParameterMax=value;fCutImpactParameter=kTRUE;}
49 void SetImpactParameterMin(Double_t value) {fImpactParameterMin=value;fCutImpactParameter=kTRUE;}
50 void SetImpactParameterRange(Double_t min, Double_t max) {fImpactParameterMin=min;fImpactParameterMax=max;fCutImpactParameter=kTRUE;}
51 void SetPrimaryVertexXrange(Double_t min, Double_t max)
52 {fCutPrimaryVertexX=kTRUE; fPrimaryVertexXmin=min; fPrimaryVertexXmax=max;}
53 void SetPrimaryVertexYrange(Double_t min, Double_t max)
54 {fCutPrimaryVertexY=kTRUE; fPrimaryVertexYmin=min; fPrimaryVertexYmax=max;}
55 void SetPrimaryVertexZrange(Double_t min, Double_t max)
56 {fCutPrimaryVertexZ=kTRUE; fPrimaryVertexZmin=min; fPrimaryVertexZmax=max;}
57 void SetNContributorsRange(Int_t min, Int_t max=INT_MAX)
58 {fCutNContributors=kTRUE; fNContributorsMin=min; fNContributorsMax=max;}
59 void SetMeanPtRange(Double_t min, Double_t max) {fCutMeanPt=kTRUE; fMeanPtMax=max; fMeanPtMin=min;}
60 void SetCutSPDvertexerAnomaly(Bool_t b=kTRUE) {fCutSPDvertexerAnomaly=b;}
61 void SetCutZDCtiming(Bool_t c=kTRUE) {fCutZDCtiming=c;}
62 void SetCutSPDTRKVtxZ(Bool_t b=kTRUE) {fCutSPDTRKVtxZ=b;}
63 void SetCutTPCmultiplicityOutliers(Bool_t b=kTRUE) {fCutTPCmultiplicityOutliers=b;}
65 Int_t GetNumberOfTracksMax() const {return fNumberOfTracksMax;}
66 Int_t GetNumberOfTracksMin() const {return fNumberOfTracksMin;}
67 Int_t GetRefMultMax() const {return fRefMultMax;}
68 Int_t GetRefMultMin() const {return fRefMultMin;}
69 void SetRefMultMethod(refMultMethod m) {fRefMultMethod=m;}
70 void SetRefMultMethod(AliESDtrackCuts::MultEstTrackType m) { fRefMultMethodAliESDtrackCuts=m;
71 fUseAliESDtrackCutsRefMult=kTRUE; }
72 refMultMethod GetRefMultMethod() const {return fRefMultMethod;}
73 void SetRefMultCuts( AliFlowTrackCuts* cuts ) {fRefMultCuts=static_cast<AliFlowTrackCuts*>(cuts->Clone());}
74 void SetMeanPtCuts( AliFlowTrackCuts* cuts ) {fMeanPtCuts=static_cast<AliFlowTrackCuts*>(cuts->Clone());}
75 AliFlowTrackCuts* GetRefMultCuts() const {return fRefMultCuts;}
76 void DefineHistograms();
77 void SetQA(Bool_t b=kTRUE) {if (b) DefineHistograms();}
78 TList* GetQA() const {return fQA;}
79 TH1* QAbefore(Int_t i) {return static_cast<TH1*>(static_cast<TList*>(fQA->At(0))->At(i));}
80 TH1* QAafter(Int_t i) {return static_cast<TH1*>(static_cast<TList*>(fQA->At(1))->At(i));}
82 Int_t RefMult(AliVEvent* event, AliMCEvent *mcEvent = 0x0);
83 //Int_t GetRefMult() {return fRefMult;}
84 Int_t GetReferenceMultiplicity(AliVEvent* event, AliMCEvent *mcEvent) {return RefMult(event,mcEvent);}
85 const char* CentrMethName(refMultMethod method) const;
86 void SetCentralityPercentileMethod( refMultMethod m) {fCentralityPercentileMethod=m;}
87 void SetUseCentralityUnchecked(Bool_t b=kTRUE) {fUseCentralityUnchecked=b;}
89 Float_t GetCentrality(AliVEvent* event, AliMCEvent* mcEvent);
90 void SetUsedDataset(Bool_t b=kTRUE) {fData2011=b;}
92 void Browse(TBrowser* b);
93 Long64_t Merge(TCollection* list);
94 TH2F *GetCorrelationTPCvsGlobalMultiplicity() {return fhistTPCvsGlobalMult;}
99 Bool_t fCutNumberOfTracks;//cut on # of tracks
100 Int_t fNumberOfTracksMax; //limits
101 Int_t fNumberOfTracksMin; //limits
102 Bool_t fCutRefMult; //cut on refmult
103 refMultMethod fRefMultMethod; //how do we calculate refmult?
104 Bool_t fUseAliESDtrackCutsRefMult; //use AliESDtrackCuts for refmult calculation
105 AliESDtrackCuts::MultEstTrackType fRefMultMethodAliESDtrackCuts;
106 Int_t fRefMultMax; //max refmult
107 Int_t fRefMultMin; //min refmult
108 AliFlowTrackCuts* fRefMultCuts; //cuts
109 AliFlowTrackCuts* fMeanPtCuts; //mean pt cuts
110 AliFlowTrackCuts* fStandardTPCcuts; //Standard TPC cuts
111 AliFlowTrackCuts* fStandardGlobalCuts; //StandardGlobalCuts
112 Bool_t fCutPrimaryVertexX; //cut on x of prim vtx
113 Double_t fPrimaryVertexXmax; //max x prim vtx
114 Double_t fPrimaryVertexXmin; //min x prim vtx
115 Bool_t fCutPrimaryVertexY; //cut on y of prim vtx
116 Double_t fPrimaryVertexYmax; //max y prim vtx
117 Double_t fPrimaryVertexYmin; //min y prim vtx
118 Bool_t fCutPrimaryVertexZ; //cut on z of prim vtx
119 Double_t fPrimaryVertexZmax; //max z prim vtx
120 Double_t fPrimaryVertexZmin; //min z prim vtx
121 Bool_t fCutNContributors; //cut on number of contributors
122 Int_t fNContributorsMax; //maximal number of contrib
123 Int_t fNContributorsMin; //minimal number of contrib
124 Bool_t fCutMeanPt; //cut on mean pt
125 Double_t fMeanPtMax; //max mean pt
126 Double_t fMeanPtMin; //min mean pt
127 Bool_t fCutSPDvertexerAnomaly; //cut on the spd vertexer anomaly
128 Bool_t fCutSPDTRKVtxZ; //require compatibility between SPDvertexz TRKvertexz
129 Bool_t fCutTPCmultiplicityOutliers; //cut TPC multiplicity outliers
130 Bool_t fUseCentralityUnchecked; //use the unchecked method
131 refMultMethod fCentralityPercentileMethod; //where to get the percentile from
132 Bool_t fCutZDCtiming; //cut on ZDC timing
133 AliTriggerAnalysis fTrigAna; //trigger analysis object
134 Bool_t fCutImpactParameter; //cut on impact parameter (MC header)
135 Double_t fImpactParameterMin; // min impact parameter
136 Double_t fImpactParameterMax; // max impact parameter
137 TH2F *fhistTPCvsGlobalMult; //!correlation between TPCMult and GlobalMult
138 Bool_t fData2011; //2011 data is used
140 ClassDef(AliFlowEventCuts,5)