1 #ifndef ALIANALYSISTASKCHECKCASCADE_H
2 #define ALIANALYSISTASKCHECKCASCADE_H
4 /* See cxx source for full Copyright notice */
6 //-----------------------------------------------------------------
7 // AliAnalysisTaskCheckCascade class
8 // (AliAnalysisTaskCheckCascade)
9 // This task has three roles :
10 // 1. QAing the Cascades from ESD and AOD
11 // Origin: AliAnalysisTaskESDCheckV0 by B.H. Nov2007, hippolyt@in2p3.fr
12 // 2. Prepare the plots which stand as raw material for yield extraction
13 // 3. Rough azimuthal correlation study (Eta, Phi)
14 // Adapted to Cascade : A.Maire Mar2008, antonin.maire@ires.in2p3.fr
15 // Modified : A.Maire Aug2009, antonin.maire@ires.in2p3.fr
16 //-----------------------------------------------------------------
28 #include "AliAnalysisTaskSE.h"
30 class AliAnalysisTaskCheckCascade : public AliAnalysisTaskSE {
32 AliAnalysisTaskCheckCascade();
33 AliAnalysisTaskCheckCascade(const char *name);
34 virtual ~AliAnalysisTaskCheckCascade() {}
36 virtual void UserCreateOutputObjects();
37 virtual void UserExec(Option_t *option);
38 void DoAngularCorrelation(const Char_t *lCascType,
39 Double_t lInvMassCascade,
40 const Int_t *lArrTrackID,
41 TVector3 &lTVect3MomXi,
43 virtual void Terminate(Option_t *);
45 void SetCollidingSystems(Short_t collidingSystems = 0) {fCollidingSystems = collidingSystems;}
46 void SetAnalysisType (const char* analysisType = "ESD") {fAnalysisType = analysisType;}
49 TString fAnalysisType; // "ESD" or "AOD" analysis type
50 Short_t fCollidingSystems; // 0 = pp collisions or 1 = AA collisions
53 TList *fListHistCascade; //! List of Cascade histograms
55 // - General histos (filled for any event)
56 TH1F *fHistTrackMultiplicity; //! Track multiplicity distribution
57 TH1F *fHistCascadeMultiplicity; //! Cascade multiplicity distribution
61 TH1F *fHistVtxStatus; //! Is there a tracking vertex in the event ?
63 // Vtx coming from the full tracking
64 TH1F *fHistPosTrkgPrimaryVtxX; //! primary vertex position distribution in x
65 TH1F *fHistPosTrkgPrimaryVtxY; //! primary vertex position distribution in y
66 TH1F *fHistPosTrkgPrimaryVtxZ; //! primary vertex position distribution in z
67 TH1F *fHistTrkgPrimaryVtxRadius; //! primary vertex (3D) radius distribution
69 // Best primary Vtx available for the event
70 TH1F *fHistPosBestPrimaryVtxX; //! (best) primary vertex position distribution in x
71 TH1F *fHistPosBestPrimaryVtxY; //! (best) primary vertex position distribution in y
72 TH1F *fHistPosBestPrimaryVtxZ; //! (best) primary vertex position distribution in z
73 TH1F *fHistBestPrimaryVtxRadius; //! (best) primary vertex radius distribution
75 // Correlation Best Vtx / Full Tracking Vtx
76 TH2F *f2dHistTrkgPrimVtxVsBestPrimVtx; //! Radius of prim. Vtx from tracks Vs Radius of best Prim. Vtx
79 // PART 1 : Adavanced QA
80 // - Typical histos on the variables used for the selection of cascades
81 TH1F *fHistEffMassXi; //! reconstructed cascade effective mass
82 TH1F *fHistChi2Xi; //! chi2 value
83 TH1F *fHistDcaXiDaughters; //! dca between Xi's daughters
84 TH1F *fHistDcaBachToPrimVertex; //! dca of the bachelor track to primary vertex
85 TH1F *fHistXiCosineOfPointingAngle; //! cosine of Xi pointing angle in a cascade
86 TH1F *fHistXiRadius; //! (transverse) radius of the cascade vertex
88 // - Histos about ~ the "V0 selection part" of the cascade, coming by inheritance from AliESDv0
89 TH1F *fHistMassLambdaAsCascDghter; //! Test Invariant Mass of Lambda coming from Cascade
90 TH1F *fHistV0Chi2Xi; //! V0 chi2 distribution, for the V0 associated to a cascade
91 TH1F *fHistDcaV0DaughtersXi; //! Dca between V0 daughters, for the V0 associated to a cascade
92 TH1F *fHistDcaV0ToPrimVertexXi; //! Dca of V0 to primary vertex, for the V0 associated to a cascade
93 TH1F *fHistV0CosineOfPointingAngleXi; //! Cosine of V0 pointing angle, for the V0 associated to a cascade
94 TH1F *fHistV0RadiusXi; //! V0 (transverse) distance distribution, for the V0 associated to a cascade
96 TH1F *fHistDcaPosToPrimVertexXi; //! Dca of V0 positive daughter to primary vertex, for the V0 associated to a cascade
97 TH1F *fHistDcaNegToPrimVertexXi; //! Dca of V0 negative daughter to primary vertex, for the V0 associated to a cascade
100 // - Effective mass histos for cascades.
101 TH1F *fHistMassXiMinus; //! reconstructed cascade effective mass, under Xi- hyp.
102 TH1F *fHistMassXiPlus; //! reconstructed cascade effective mass, under Xi+ hyp.
103 TH1F *fHistMassOmegaMinus; //! reconstructed cascade effective mass, under Omega- hyp.
104 TH1F *fHistMassOmegaPlus; //! reconstructed cascade effective mass, under Omega+ hyp.
106 TH1F *fHistMassWithCombPIDXiMinus; //! reconstructed Xi- effective mass, with bach. comb PID
107 TH1F *fHistMassWithCombPIDXiPlus; //! reconstructed Xi+ effective mass, with bach. comb PID
108 TH1F *fHistMassWithCombPIDOmegaMinus; //! reconstructed Omega- effective mass, with bach. comb PID
109 TH1F *fHistMassWithCombPIDOmegaPlus; //! reconstructed Omega+ effective mass, with bach. comb PID
111 // - Complements for QA
112 TH1F *fHistXiTransvMom; //! Xi transverse momentum
113 TH1F *fHistXiTotMom; //! Xi momentum norm
115 TH1F *fHistBachTransvMom; //! bachelor transverse momentum
116 TH1F *fHistBachTotMom; //! bachelor momentum norm
118 TH1F *fHistChargeXi; //! Charge sign of the cascade candidate
119 TH1F *fHistV0toXiCosineOfPointingAngle; //! Cos. of Pointing angle between the V0 mom and the Xi-V0 vtx line
121 TH1F *fHistRapXi; //! rapidity of Xi candidates
122 TH1F *fHistRapOmega; //! rapidity of Omega candidates
123 TH1F *fHistEta; //! eta distrib. of all the cascade candidates
124 TH1F *fHistTheta; //! theta distrib. of all the cascade candidates
125 TH1F *fHistPhi; //! phi distrib. of all the cascade candidates
127 TH2F *f2dHistArmenteros; //! alpha(casc. cand.) Vs PtArm(casc. cand.)
129 TH2F *f2dHistEffMassLambdaVsEffMassXiMinus; //! Xi- Eff mass Vs V0 Eff mass, under Xi- hyp.
130 TH2F *f2dHistEffMassXiVsEffMassOmegaMinus; //! Xi- Eff mass Vs Omega- Eff mass, for negative cascades
131 TH2F *f2dHistEffMassLambdaVsEffMassXiPlus; //! Xi+ Eff mass Vs V0 Eff mass, under Xi+ hyp.
132 TH2F *f2dHistEffMassXiVsEffMassOmegaPlus; //! Xi+ Eff mass Vs Omega+ Eff mass, for positive cascades
134 TH2F *f2dHistXiRadiusVsEffMassXiMinus; //! transv. casc. decay radius Vs Xi- Eff mass, under Xi- hyp.
135 TH2F *f2dHistXiRadiusVsEffMassXiPlus; //! transv. casc. decay radius Vs Xi+ Eff mass, under Xi+ hyp.
136 TH2F *f2dHistXiRadiusVsEffMassOmegaMinus; //! transv. casc. decay radius Vs Omega- Eff mass, under Omega- hyp.
137 TH2F *f2dHistXiRadiusVsEffMassOmegaPlus; //! transv. casc. decay radius Vs Omega+ Eff mass, under Omega+ hyp.
140 // PART 2 : TH3F needed for pt spectrum and yield extraction
142 TH3F *f3dHistXiPtVsEffMassVsYXiMinus; //! casc. transv. momemtum Vs Xi- Eff mass Vs Y
143 TH3F *f3dHistXiPtVsEffMassVsYXiPlus; //! casc. transv. momemtum Vs Xi+ Eff mass Vs Y
144 TH3F *f3dHistXiPtVsEffMassVsYOmegaMinus; //! casc. transv. momemtum Vs Omega- Eff mass Vs Y
145 TH3F *f3dHistXiPtVsEffMassVsYOmegaPlus; //! casc. transv. momemtum Vs Omega+ Eff mass Vs Y
147 // With single PID : proton PID for Xi (pion = useless) / bachelor PID for Omega
148 // = a priori, the most important PID info for each species
149 TH3F *f3dHistXiPtVsEffMassVsYWithCombPIDXiMinus; //! casc. transv. momemtum Vs Xi- Eff mass Vs Y
150 TH3F *f3dHistXiPtVsEffMassVsYWithCombPIDXiPlus; //! casc. transv. momemtum Vs Xi+ Eff mass Vs Y
151 TH3F *f3dHistXiPtVsEffMassVsYWithCombPIDOmegaMinus; //! casc. transv. momemtum Vs Omega- Eff mass Vs Y
152 TH3F *f3dHistXiPtVsEffMassVsYWithCombPIDOmegaPlus; //! casc. transv. momemtum Vs Omega+ Eff mass Vs Y
154 // With double PID : proton PID + bachelor PID for Omega
155 // = "second order" refinement for omegas
156 TH3F *f3dHistXiPtVsEffMassVsYWith2CombPIDOmegaMinus; //! casc. transv. momemtum Vs Omega- Eff mass Vs Y
157 TH3F *f3dHistXiPtVsEffMassVsYWith2CombPIDOmegaPlus; //! casc. transv. momemtum Vs Omega+ Eff mass Vs Y
161 // PART 3 : Azimuthal correlation study
162 THnSparseF *fHnSpAngularCorrXiMinus; //! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
163 THnSparseF *fHnSpAngularCorrXiPlus; //! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
164 THnSparseF *fHnSpAngularCorrOmegaMinus; //! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
165 THnSparseF *fHnSpAngularCorrOmegaPlus; //! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
168 AliAnalysisTaskCheckCascade(const AliAnalysisTaskCheckCascade&); // not implemented
169 AliAnalysisTaskCheckCascade& operator=(const AliAnalysisTaskCheckCascade&); // not implemented
171 ClassDef(AliAnalysisTaskCheckCascade, 6);