X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ANALYSIS%2FAliESDtrackCuts.h;h=c2ad374bf44fc8a1b20b4c7cbea3a640ec376bbb;hb=624cda230152b9ee2efcc0b12669a92290025ec0;hp=1e282d5ee8f6b8df6cd6406ad4786b2ea371db8b;hpb=07b1b131167ee0e0f656c363e83394f4d0189618;p=u%2Fmrichter%2FAliRoot.git diff --git a/ANALYSIS/AliESDtrackCuts.h b/ANALYSIS/AliESDtrackCuts.h index 1e282d5ee8f..c2ad374bf44 100644 --- a/ANALYSIS/AliESDtrackCuts.h +++ b/ANALYSIS/AliESDtrackCuts.h @@ -41,6 +41,9 @@ class AliESDtrackCuts : public AliAnalysisCuts public: enum ITSClusterRequirement { kOff = 0, kNone, kAny, kFirst, kOnlyFirst, kSecond, kOnlySecond, kBoth }; enum Detector { kSPD = 0, kSDD, kSSD }; + enum MultEstTrackCuts { kMultEstTrackCutGlobal = 0, kMultEstTrackCutITSSA, kMultEstTrackCutDCAwSPD, kMultEstTrackCutDCAwoSPD, kNMultEstTrackCuts /* this must always be the last */}; + enum MultEstTrackType { kTrackletsITSTPC = 0, kTrackletsITSSA, kTracklets }; + enum VertexType { kVertexTracks = 0x1, kVertexSPD = 0x2, kVertexTPC = 0x4 }; AliESDtrackCuts(const Char_t* name = "AliESDtrackCuts", const Char_t* title = ""); virtual ~AliESDtrackCuts(); @@ -49,13 +52,15 @@ public: {return AcceptTrack((AliESDtrack*)obj);} virtual Bool_t IsSelected(TList* /*list*/) {return kTRUE;} - Bool_t AcceptTrack(AliESDtrack* esdTrack); - TObjArray* GetAcceptedTracks(AliESDEvent* esd, Bool_t bTPC = kFALSE); - Int_t CountAcceptedTracks(AliESDEvent* const esd); + Bool_t AcceptTrack(const AliESDtrack* esdTrack); + TObjArray* GetAcceptedTracks(const AliESDEvent* esd, Bool_t bTPC = kFALSE); + Int_t CountAcceptedTracks(const AliESDEvent* const esd); - static Int_t GetReferenceMultiplicity(AliESDEvent* esd, Bool_t tpcOnly); + static Int_t GetReferenceMultiplicity(const AliESDEvent* esd, Bool_t tpcOnly); + static Int_t GetReferenceMultiplicity(const AliESDEvent* esd, MultEstTrackType trackType = kTrackletsITSTPC, Float_t etaRange = 0.5); + static AliESDtrackCuts* GetMultEstTrackCuts(MultEstTrackCuts cut); - static AliESDtrack* GetTPCOnlyTrack(AliESDEvent* esd, Int_t iTrack); + static AliESDtrack* GetTPCOnlyTrack(const AliESDEvent* esd, Int_t iTrack); // Standard cut definitions static AliESDtrackCuts* GetStandardTPCOnlyTrackCuts(); @@ -63,8 +68,11 @@ public: static AliESDtrackCuts* GetStandardITSTPCTrackCuts2010(Bool_t selPrimaries=kTRUE, Int_t clusterCut=0); static AliESDtrackCuts* GetStandardITSSATrackCuts2009(Bool_t selPrimaries=kTRUE, Bool_t useForPid=kTRUE); static AliESDtrackCuts* GetStandardITSSATrackCuts2010(Bool_t selPrimaries=kTRUE, Bool_t useForPid=kTRUE); + static AliESDtrackCuts* GetStandardITSSATrackCutsPbPb2010(Bool_t selPrimaries=kTRUE, Bool_t useForPid=kTRUE); static AliESDtrackCuts* GetStandardITSPureSATrackCuts2009(Bool_t selPrimaries=kTRUE, Bool_t useForPid=kTRUE); static AliESDtrackCuts* GetStandardITSPureSATrackCuts2010(Bool_t selPrimaries=kTRUE, Bool_t useForPid=kTRUE); + // Standard cuts for daughter tracks + static AliESDtrackCuts* GetStandardV0DaughterCuts(); virtual Long64_t Merge(TCollection* list); virtual void Copy(TObject &c) const; @@ -74,12 +82,16 @@ public: //###################################################### // track quality cut setters void SetMinNClustersTPC(Int_t min=-1) {fCutMinNClusterTPC=min;} + void SetMinNClustersTPCPtDep(TFormula *f1=0x0, Float_t ptmax=0.); void SetMinNClustersITS(Int_t min=-1) {fCutMinNClusterITS=min;} void SetMinNCrossedRowsTPC(Float_t min=-1) { fCutMinNCrossedRowsTPC=min;} void SetMinRatioCrossedRowsOverFindableClustersTPC(Float_t min = -1) { fCutMinRatioCrossedRowsOverFindableClustersTPC=min;} void SetClusterRequirementITS(Detector det, ITSClusterRequirement req = kOff) { fCutClusterRequirementITS[det] = req; } void SetMaxChi2PerClusterTPC(Float_t max=1e10) {fCutMaxChi2PerClusterTPC=max;} void SetMaxChi2PerClusterITS(Float_t max=1e10) {fCutMaxChi2PerClusterITS=max;} + void SetMaxChi2TPCConstrainedGlobal(Float_t max=1e10) {fCutMaxChi2TPCConstrainedVsGlobal = max; } + void SetMaxChi2TPCConstrainedGlobalVertexType(Int_t vertexType = kVertexTracks | kVertexSPD) { fCutMaxChi2TPCConstrainedVsGlobalVertexType = vertexType; } + void SetMaxNOfMissingITSPoints(Int_t max=6) {fCutMaxMissingITSPoints=max;} void SetRequireTPCRefit(Bool_t b=kFALSE) {fCutRequireTPCRefit=b;} void SetRequireTPCStandAlone(Bool_t b=kFALSE) {fCutRequireTPCStandAlone=b;} void SetRequireITSRefit(Bool_t b=kFALSE) {fCutRequireITSRefit=b;} @@ -114,9 +126,13 @@ public: Int_t GetMinNClusterTPC() const { return fCutMinNClusterTPC;} Int_t GetMinNClustersITS() const { return fCutMinNClusterITS;} + TFormula *GetMinNClustersTPCPtDep() const { return f1CutMinNClustersTPCPtDep;} ITSClusterRequirement GetClusterRequirementITS(Detector det) const { return fCutClusterRequirementITS[det]; } Float_t GetMaxChi2PerClusterTPC() const { return fCutMaxChi2PerClusterTPC;} Float_t GetMaxChi2PerClusterITS() const { return fCutMaxChi2PerClusterITS;} + Float_t GetMaxChi2TPCConstrainedGlobal() const { return fCutMaxChi2TPCConstrainedVsGlobal; } + Int_t GetMaxChi2TPCConstrainedGlobalVertexType() const { return fCutMaxChi2TPCConstrainedVsGlobalVertexType; } + Int_t GetMaxNOfMissingITSPoints() const { return fCutMaxMissingITSPoints;} Bool_t GetRequireTPCRefit() const { return fCutRequireTPCRefit;} Bool_t GetRequireTPCStandAlone() const { return fCutRequireTPCStandAlone;} Bool_t GetRequireITSRefit() const { return fCutRequireITSRefit;} @@ -163,14 +179,16 @@ public: void SaveHistograms(const Char_t* dir = 0); void DrawHistograms(); - static Float_t GetSigmaToVertex(AliESDtrack* const esdTrack); + static Float_t GetSigmaToVertex(const AliESDtrack* const esdTrack); static void EnableNeededBranches(TTree* tree); // void SaveQualityCuts(Char_t* file) // void LoadQualityCuts(Char_t* file) - TH1F* GetDZNormalized(Int_t i) const { return fhDZNormalized[i]; } + TH1F* GetDZNormalized(Int_t i) const { return fhDZNormalized[i]; } + TH1F* GetNClustersTPC(Int_t i) const { return fhNClustersTPC[i]; } + TH1F* GetPtHist(Int_t i) const { return fhPt[i]; } protected: void Init(); // sets everything to 0 @@ -178,21 +196,27 @@ protected: Bool_t CheckPtDepDCA(TString dist,Bool_t print=kFALSE) const; void SetPtDepDCACuts(Double_t pt); - enum { kNCuts = 38 }; + enum { kNCuts = 40 }; //###################################################### // esd track quality cuts static const Char_t* fgkCutNames[kNCuts]; //! names of cuts (for internal use) + static AliESDtrackCuts* fgMultEstTrackCuts[kNMultEstTrackCuts]; //! track cuts used for the multiplicity estimate Int_t fCutMinNClusterTPC; // min number of tpc clusters Int_t fCutMinNClusterITS; // min number of its clusters Float_t fCutMinNCrossedRowsTPC; // min number of tpc crossed rows Float_t fCutMinRatioCrossedRowsOverFindableClustersTPC; // min ratio crossed rows / findable clusters + TFormula *f1CutMinNClustersTPCPtDep; // pt dependent tpc clusters cut + Float_t fCutMaxPtDepNClustersTPC; // maximum pt for pt dependend TPC cluster cut. For pt=>ptmax NClusterMin = f1CutMinNClustersTPCPtDep->Eval(fCutMaxPtDepNClustersTPC). ITSClusterRequirement fCutClusterRequirementITS[3]; // detailed ITS cluster requirements for (SPD, SDD, SSD) Float_t fCutMaxChi2PerClusterTPC; // max tpc fit chi2 per tpc cluster Float_t fCutMaxChi2PerClusterITS; // max its fit chi2 per its cluster + Float_t fCutMaxChi2TPCConstrainedVsGlobal; // max chi2 TPC track constrained with vtx vs. global track + Int_t fCutMaxChi2TPCConstrainedVsGlobalVertexType; // vertex type for max chi2 TPC track constrained with vtx vs. global track (can be configured to accept several vertex types) + Int_t fCutMaxMissingITSPoints; // max n. of missing ITS points Float_t fCutMaxC11; // max cov. matrix diag. elements (res. y^2) Float_t fCutMaxC22; // max cov. matrix diag. elements (res. z^2) @@ -249,11 +273,15 @@ protected: TH1F* fhNClustersITS[2]; //-> TH1F* fhNClustersTPC[2]; //-> + TH1F* fhNSharedClustersTPC[2]; //-> TH1F* fhNCrossedRowsTPC[2]; //-> TH1F* fhRatioCrossedRowsOverFindableClustersTPC[2]; // -> TH1F* fhChi2PerClusterITS[2]; //-> TH1F* fhChi2PerClusterTPC[2]; //-> + TH1F* fhChi2TPCConstrainedVsGlobal[2]; //-> + TH1F* fhNClustersForITSPID[2]; //-> number of points in SDD+SSD (ITS PID selection) + TH1F* fhNMissingITSPoints[2]; //-> number of missing ITS points TH1F* fhC11[2]; //-> TH1F* fhC22[2]; //-> @@ -281,7 +309,7 @@ protected: TH1F* fhCutStatistics; //-> statistics of what cuts the tracks did not survive TH2F* fhCutCorrelation; //-> 2d statistics plot - ClassDef(AliESDtrackCuts, 13) + ClassDef(AliESDtrackCuts, 18) };