[u/mrichter/AliRoot.git] / PWGLF / STRANGENESS / Cascades / AliAnalysisTaskCheckCascadePbPb.h

#ifndef ALIANALYSISTASKCHECKCASCADEPBPB_H
#define ALIANALYSISTASKCHECKCASCADEPBPB_H

/*  See cxx source for full Copyright notice */

//-----------------------------------------------------------------
//            AliAnalysisTaskCheckCascadePbPb class
//              Origin AliAnalysisTaskCheckCascade
//              This task has four roles :
//                1. QAing the Cascades from ESD and AOD
//                   Origin:  AliAnalysisTaskESDCheckV0 by Boris Hippolyte Nov2007, hippolyt@in2p3.fr
//                2. Prepare the plots which stand as raw material for yield extraction (wi/wo PID)
//                3. Supply an AliCFContainer meant to define the optimised topological selections
//                4. Rough azimuthal correlation study (Eta, Phi)
//                Adapted to Cascade : A.Maire Mar2008, antonin.maire@ires.in2p3.fr
//                Modified :           A.Maire Mar2010, antonin.maire@ires.in2p3.fr
//                Modified for PbPb analysis: M. Nicassio Feb 2011, maria.nicassio@ba.infn.it
//-----------------------------------------------------------------

class TList;
class TH1F;
class TH2F;
class TH3F;
class TVector3;
class THnSparse;
 
class AliESDEvent;
class AliPhysicsSelection;
class AliCFContainer;
class AliPIDResponse;

#include "TString.h"

#include "AliAnalysisTaskSE.h"

class AliAnalysisTaskCheckCascadePbPb : public AliAnalysisTaskSE {
 public:
  AliAnalysisTaskCheckCascadePbPb();
  AliAnalysisTaskCheckCascadePbPb(const char *name);
  virtual ~AliAnalysisTaskCheckCascadePbPb();
  
  virtual void   UserCreateOutputObjects();
  virtual void   UserExec(Option_t *option);
/*	  void   DoAngularCorrelation(const Char_t   *lCascType, 
				            Double_t  lInvMassCascade, 
				      const Int_t    *lArrTrackID,
				            TVector3 &lTVect3MomXi, 
				            Double_t  lEtaXi);*/
  virtual Int_t  DoESDTrackWithTPCrefitMultiplicity(const AliESDEvent *lESDevent);
  virtual void   Terminate(Option_t *);
  
  void SetAnalysisType               (const char* analysisType          = "ESD") { fAnalysisType                = analysisType;               }
  void SetRelaunchV0CascVertexers    (Bool_t rerunV0CascVertexers       = 0    ) { fkRerunV0CascVertexers       = rerunV0CascVertexers;       }
  void SetQualityCutZprimVtxPos      (Bool_t qualityCutZprimVtxPos      = kTRUE) { fkQualityCutZprimVtxPos      = qualityCutZprimVtxPos;      }
  void SetQualityCutNoTPConlyPrimVtx (Bool_t qualityCutNoTPConlyPrimVtx = kTRUE) { fkQualityCutNoTPConlyPrimVtx = qualityCutNoTPConlyPrimVtx; }
  void SetQualityCutTPCrefit         (Bool_t qualityCutTPCrefit         = kTRUE) { fkQualityCutTPCrefit         = qualityCutTPCrefit;         }
  void SetQualityCutnTPCcls          (Bool_t qualityCutnTPCcls          = kTRUE) { fkQualityCutnTPCcls          = qualityCutnTPCcls;          }
  void SetQualityCutMinnTPCcls       (Int_t minnTPCcls                  = 70   ) { fMinnTPCcls                  = minnTPCcls;                 }
  void SetExtraSelections            (Bool_t extraSelections            = 0    ) { fkExtraSelections            = extraSelections;            }
  void SetCentralityLowLim           (Float_t centrlowlim               = 0.   ) { fCentrLowLim                 = centrlowlim;                }  
  void SetCentralityUpLim            (Float_t centruplim                = 100. ) { fCentrUpLim                  = centruplim;                 }
  void SetCentralityEst              (TString   centrest                = "V0M") { fCentrEstimator              = centrest;                   }
  void SetUseCleaning                (Bool_t   usecleaning              = kTRUE) { fkUseCleaning                = usecleaning;                }
  void SetVertexRange                (Float_t vtxrange                  = 0.   ) { fVtxRange                    = vtxrange;                   }
  void SetMinptCutOnDaughterTracks   (Float_t minptdaughtrks            = 0.   ) { fMinPtCutOnDaughterTracks    = minptdaughtrks;             }
  void SetEtaCutOnDaughterTracks     (Float_t etadaughtrks              = 0.   ) { fEtaCutOnDaughterTracks      = etadaughtrks;               }
  //Setters for the V0 and cascade Vertexer Parameters
  void SetV0VertexerMaxChisquare           (Double_t lParameter){ fV0VertexerSels[0] = lParameter; }
  void SetV0VertexerDCAFirstToPV           (Double_t lParameter){ fV0VertexerSels[1] = lParameter; }
  void SetV0VertexerDCASecondtoPV          (Double_t lParameter){ fV0VertexerSels[2] = lParameter; }
  void SetV0VertexerDCAV0Daughters         (Double_t lParameter){ fV0VertexerSels[3] = lParameter; }
  void SetV0VertexerCosinePA               (Double_t lParameter){ fV0VertexerSels[4] = lParameter; }
  void SetV0VertexerMinRadius              (Double_t lParameter){ fV0VertexerSels[5] = lParameter; }
  void SetV0VertexerMaxRadius              (Double_t lParameter){ fV0VertexerSels[6] = lParameter; }
  void SetCascVertexerMaxChisquare         (Double_t lParameter){ fCascadeVertexerSels[0] = lParameter; }
  void SetCascVertexerMinV0ImpactParameter (Double_t lParameter){ fCascadeVertexerSels[1] = lParameter; }
  void SetCascVertexerV0MassWindow         (Double_t lParameter){ fCascadeVertexerSels[2] = lParameter; }
  void SetCascVertexerDCABachToPV          (Double_t lParameter){ fCascadeVertexerSels[3] = lParameter; }
  void SetCascVertexerDCACascadeDaughters  (Double_t lParameter){ fCascadeVertexerSels[4] = lParameter; }
  void SetCascVertexerCascadeCosinePA      (Double_t lParameter){ fCascadeVertexerSels[5] = lParameter; }
  void SetCascVertexerCascadeMinRadius     (Double_t lParameter){ fCascadeVertexerSels[6] = lParameter; }
  void SetCascVertexerCascadeMaxRadius     (Double_t lParameter){ fCascadeVertexerSels[7] = lParameter; } 


 private:
        // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
        // your data member object is created on the worker nodes and streaming is not needed.
        // http://root.cern.ch/download/doc/11InputOutput.pdf, page 14


        TString         fAnalysisType;                  // "ESD" or "AOD" analysis type	
        AliESDtrackCuts *fESDtrackCuts;                 // ESD track cuts used for primary track definition
        //TPaveText       *fPaveTextBookKeeping;          // TString to store all the relevant info necessary for book keeping (v0 cuts, cascade cuts, quality cuts, ...)
        AliPIDResponse *fPIDResponse;                   //! PID response object

        Bool_t          fkRerunV0CascVertexers;         // Boolean : kTRUE = relaunch both V0 + Cascade vertexers
        Bool_t          fkQualityCutZprimVtxPos;        // Boolean : kTRUE = cut on the prim.vtx  z-position
        Bool_t          fkQualityCutNoTPConlyPrimVtx;   // Boolean : kTRUE = prim vtx should be SPD or Tracking vertex
        Bool_t          fkQualityCutTPCrefit;           // Boolean : kTRUE = ask for TPCrefit for the 3 daughter tracks
        Bool_t          fkQualityCutnTPCcls;            // Boolean : kTRUE = ask for at least n TPC clusters for each daughter track
        Int_t           fMinnTPCcls;                    // minimum number of TPC cluster for daughter tracks
        Bool_t          fkExtraSelections;              // Boolean : kTRUE = apply tighter selections, before starting the analysis
        Float_t         fCentrLowLim;                   // Lower limit for centrality percentile selection
        Float_t         fCentrUpLim;                    // Upper limit for centrality percentile selection
        TString         fCentrEstimator;                // string for the centrality estimator
        Bool_t          fkUseCleaning;                  // Boolean : kTRUE = uses all the cleaning criteria of centrality selections (vertex cut + outliers) otherwise only outliers
        Float_t         fVtxRange;                      // to select events with |zvtx|<fVtxRange cm
        Float_t         fMinPtCutOnDaughterTracks;      // minimum pt cut on daughter tracks
        Float_t         fEtaCutOnDaughterTracks;        // pseudorapidity cut on daughter tracks
       
        Double_t        fV0VertexerSels[7];             // Array to store the 7 values for the different selections V0 related (if fkRerunV0CascVertexers)
        Double_t        fCascadeVertexerSels[8];        // Array to store the 8 values for the different selections Casc. related (if fkRerunV0CascVertexers)

        TList      *fListHistCascade;                   //! List of Cascade histograms
        
        // - General histos (filled before the trigger selection)
        TH2F    *fHistEvtsInCentralityBinsvsNtracks;    //! Events in centrality bins vs N ESDtracks
        TH1F    *fHistCascadeMultiplicityBeforeEvSel;   //! Cascade multiplicity distribution
         
        // - General histos (filled for any triggered event)
        TH1F    *fHistCascadeMultiplicityForCentrEvt;              //! Cascade multiplicity distribution
        TH1F    *fHistTrackMultiplicityForCentrEvt;                //! Track multiplicity distribution (without any cut = include ITS stand-alone + global tracks)
        TH1F    *fHistTPCrefitTrackMultiplicityForCentrEvt;        //! Track multiplicity distribution for tracks with TPCrefit
        
        // - General histos (filled for events selected in this analysis (|z(prim. vtx)| < 10 cm + prim vtx = SPD or Tracking) )
        TH1F    *fHistCascadeMultiplicityForSelEvt;     //! Cascade multiplicity distribution
        TH1F    *fHistPosBestPrimaryVtxXForSelEvt;      //! (best) primary vertex position distribution in x 
        TH1F    *fHistPosBestPrimaryVtxYForSelEvt;      //! (best) primary vertex position distribution in y
        TH1F    *fHistPosBestPrimaryVtxZForSelEvt;      //! (best) primary vertex position distribution in z
        
        
        // - Characteristics for event with >1 cascade : Track Multiplicity, TPC clusters + Prim. vertex positions
        TH1F	*fHistTPCrefitTrackMultiplicityForCascadeEvt;   //! TPCrefit Track multiplicity distribution for event with >1 cascade candidate (NB: after quality sel. within the task)

        TH1F    *fHistPosV0TPCClusters;                 //! TPC clusters distribution for Positive V0 daughter track
        TH1F    *fHistNegV0TPCClusters;                 //! TPC clusters distribution for Negative V0 daughter track
        TH1F    *fHistBachTPCClusters;                  //! TPC clusters distribution for Bachelor             track

        TH1F    *fHistVtxStatus;                        //! Is there a tracking vertex in the cascade event ?

                // Vtx coming from the full tracking, for events containing at least a cascade
        TH1F    *fHistPosTrkgPrimaryVtxXForCascadeEvt;  //! primary vertex position distribution in x 
        TH1F    *fHistPosTrkgPrimaryVtxYForCascadeEvt;  //! primary vertex position distribution in y
        TH1F    *fHistPosTrkgPrimaryVtxZForCascadeEvt;  //! primary vertex position distribution in z
        TH1F    *fHistTrkgPrimaryVtxRadius;             //! primary vertex (3D) radius distribution 

                // Best primary Vtx available, for events containing at least a cascade
        TH1F    *fHistPosBestPrimaryVtxXForCascadeEvt;  //! (best) primary vertex position distribution in x 
        TH1F    *fHistPosBestPrimaryVtxYForCascadeEvt;  //! (best) primary vertex position distribution in y
        TH1F    *fHistPosBestPrimaryVtxZForCascadeEvt;  //! (best) primary vertex position distribution in z
        TH1F    *fHistBestPrimaryVtxRadius;             //! (best) primary vertex radius distribution 

                // Correlation Best Vtx / Full Tracking Vtx
        TH2F    *f2dHistTrkgPrimVtxVsBestPrimVtx;       //!  Radius of prim. Vtx from tracks Vs Radius of best Prim. Vtx


// PART 1 : Adavanced QA
// - Typical histos on the variables used for the selection of cascades
        TH1F    *fHistEffMassXi;      			//! reconstructed cascade effective mass
        TH1F    *fHistChi2Xi;         			//! chi2 value
        TH1F    *fHistDcaXiDaughters; 			//! dca between Xi's daughters
        TH1F    *fHistDcaBachToPrimVertex;		//! dca of the bachelor track to primary vertex
        TH1F    *fHistXiCosineOfPointingAngle;		//! cosine of Xi pointing angle in a cascade
        TH1F    *fHistXiRadius;                         //! (transverse) radius of the cascade vertex 
		
	// - Histos about ~ the "V0 selection part" of the cascade,  coming by inheritance from AliESDv0
	TH1F	*fHistMassLambdaAsCascDghter;		//! Test Invariant Mass of Lambda coming from Cascade
	TH1F	*fHistV0Chi2Xi;				//! V0 chi2 distribution, for the V0 associated to a cascade
	TH1F	*fHistDcaV0DaughtersXi;			//! Dca between V0 daughters, for the V0 associated to a cascade
	TH1F	*fHistDcaV0ToPrimVertexXi;		//! Dca of V0 to primary vertex, for the V0 associated to a cascade	
	TH1F	*fHistV0CosineOfPointingAngle;	        //! Cosine of V0 pointing angle, for the V0 associated to a cascade
	TH1F	*fHistV0RadiusXi;			//! V0 (transverse) distance distribution, for the V0 associated to a cascade

	TH1F	*fHistDcaPosToPrimVertexXi;		//! Dca of V0 positive daughter to primary vertex, for the V0 associated to a cascade
	TH1F	*fHistDcaNegToPrimVertexXi;		//! Dca of V0 negative daughter to primary vertex, for the V0 associated to a cascade
	

	// - Effective mass histos for cascades.
	TH1F	*fHistMassXiMinus;			//! reconstructed cascade effective mass, under Xi- hyp.
	TH1F	*fHistMassXiPlus;			//! reconstructed cascade effective mass, under Xi+ hyp.
	TH1F	*fHistMassOmegaMinus;			//! reconstructed cascade effective mass, under Omega- hyp.
	TH1F	*fHistMassOmegaPlus;			//! reconstructed cascade effective mass, under Omega+ hyp.
	
	TH1F	*fHistMassWithCombPIDXiMinus;		//! reconstructed Xi- effective mass, with bach. comb PID
	TH1F	*fHistMassWithCombPIDXiPlus;		//! reconstructed Xi+ effective mass, with bach. comb PID
	TH1F	*fHistMassWithCombPIDOmegaMinus;	//! reconstructed Omega- effective mass, with bach. comb PID
	TH1F	*fHistMassWithCombPIDOmegaPlus;		//! reconstructed Omega+ effective mass, with bach. comb PID

	// - Complements for QA
        TH1F	*fHistXiTransvMom;                      //! Xi transverse momentum, around the mass peak of Xi-/+ 
	TH1F	*fHistXiTotMom;                         //! Xi momentum norm, around the mass peak of Xi-/+
	
	TH1F	*fHistBachTransvMomXi;                  //! bachelor transverse momentum, for cand. around the mass peak of Xi-/+
	TH1F	*fHistBachTotMomXi;                     //! bachelor momentum norm, for cand. around the mass peak of Xi-/+
        TH1F    *fHistPosTransvMomXi;                   //! positive daughter transverse momentum, for cand. around the mass peak of Xi-/+
        TH1F    *fHistNegTransvMomXi;                   //! negative daughter transverse momentum, for cand. around the mass peak of Xi-/+

	TH1F	*fHistChargeXi;				//! Charge sign of the cascade candidate
	TH1F	*fHistV0toXiCosineOfPointingAngle;	//! Cos. of Pointing angle between the V0 mom and the Xi-V0 vtx line
  
	TH1F	*fHistRapXi;                            //! rapidity of Xi candidates, around the mass peak of Xi-/+
	TH1F	*fHistRapOmega;                         //! rapidity of Omega candidates, around the mass peak of Omega-/+
	TH1F	*fHistEtaXi;                            //! eta distrib. of all the cascade candidates, around the mass peak of Xi-/+
	TH1F	*fHistThetaXi;                          //! theta distrib. of all the cascade candidates, around the mass peak of Xi-/+
	TH1F	*fHistPhiXi;                            //! phi distrib. of all the cascade candidates, around the mass peak of Xi-/+
	
	TH2F	*f2dHistArmenteros;			//! alpha(casc. cand.) Vs PtArm(casc. cand.)
	
	TH2F	*f2dHistEffMassLambdaVsEffMassXiMinus;	//! Xi- Eff mass Vs V0 Eff mass, under Xi- hyp.
	TH2F	*f2dHistEffMassXiVsEffMassOmegaMinus;	//! Xi- Eff mass Vs Omega- Eff mass, for negative cascades
	TH2F	*f2dHistEffMassLambdaVsEffMassXiPlus;	//! Xi+ Eff mass Vs V0 Eff mass, under Xi+ hyp.
	TH2F	*f2dHistEffMassXiVsEffMassOmegaPlus;	//! Xi+ Eff mass Vs Omega+ Eff mass, for positive cascades
	
	TH2F	*f2dHistXiRadiusVsEffMassXiMinus;	//! transv. casc. decay radius Vs Xi- Eff mass, under Xi- hyp.
	TH2F	*f2dHistXiRadiusVsEffMassXiPlus;	//! transv. casc. decay radius Vs Xi+ Eff mass, under Xi+ hyp.
	TH2F	*f2dHistXiRadiusVsEffMassOmegaMinus;	//! transv. casc. decay radius Vs Omega- Eff mass, under Omega- hyp.
	TH2F	*f2dHistXiRadiusVsEffMassOmegaPlus;	//! transv. casc. decay radius Vs Omega+ Eff mass, under Omega+ hyp.
        
        TH2F    *f2dHistTPCdEdxOfCascDghters;           //! TPC Bethe-Bloch curve, populated with the cascade daughters
	
	
	// PART 2 : TH3F needed for pt spectrum and yield extraction
	// Without any PID
/*	TH3F	*f3dHistXiPtVsEffMassVsYXiMinus;        //! casc. transv. momemtum Vs Xi- Eff mass Vs Y
	TH3F	*f3dHistXiPtVsEffMassVsYXiPlus;         //! casc. transv. momemtum Vs Xi+ Eff mass Vs Y
	TH3F	*f3dHistXiPtVsEffMassVsYOmegaMinus;     //! casc. transv. momemtum Vs Omega- Eff mass Vs Y
	TH3F	*f3dHistXiPtVsEffMassVsYOmegaPlus;      //! casc. transv. momemtum Vs Omega+ Eff mass Vs Y
*/	
	// Compilation of all PID plots (3D = casc. transv. momemtum Vs Casc Eff mass Vs Y), stored into an AliCFContainer
	AliCFContainer  *fCFContCascadePIDXiMinus;      //! for Xi-   : Container to store any 3D histos with the different PID flavours
	AliCFContainer  *fCFContCascadePIDXiPlus;       //! for Xi+   : Container to store any 3D histos with the different PID flavours
	AliCFContainer  *fCFContCascadePIDOmegaMinus;   //! for Omega-: Container to store any 3D histos with the different PID flavours
	AliCFContainer  *fCFContCascadePIDOmegaPlus;    //! for Omega+: Container to store any 3D histos with the different PID flavours
	
	
	// PART 3 : Towards the optimisation of topological selections / systematics
	AliCFContainer  *fCFContCascadeCuts;            //! Container meant to store all the relevant distributions corresponding to the cut variables
	
	
	// PART 4 :  Azimuthal correlation study
/*	THnSparseF	*fHnSpAngularCorrXiMinus;	//! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
	THnSparseF	*fHnSpAngularCorrXiPlus;	//! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
	THnSparseF	*fHnSpAngularCorrOmegaMinus;	//! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
	THnSparseF	*fHnSpAngularCorrOmegaPlus;	//! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
*/
        TH1F *fV0Ampl;                                  //! histo to check the V0 amplitude distribution  

        TH2F    *fHistDcaXiDaughtersvsInvMass;          //! cut variables vs inv. mass
        TH2F    *fHistDcaBachToPrimVertexvsInvMass;     //! cut variables vs inv. mass
        TH2F    *fHistXiCosineOfPointingAnglevsInvMass; //! cut variables vs inv. mass
        TH2F    *fHistMassLambdaAsCascDghtervsInvMass;  //! cut variables vs inv. mass
        TH2F    *fHistDcaV0DaughtersXivsInvMass;        //! cut variables vs inv. mass
        TH2F    *fHistDcaV0ToPrimVertexXivsInvMass;     //! cut variables vs inv. mass

        // Control plots for reco pseudorapidity of daughter tracks (Xi- only)

        TH1F    *fHistEtaBachXiM;                       //! bachelor pseudorapidity
        TH1F    *fHistEtaPosXiM;                        //! positive daughter pseudorapidity
        TH1F    *fHistEtaNegXiM;                        //! negative daughter pseudorapidity


  AliAnalysisTaskCheckCascadePbPb(const AliAnalysisTaskCheckCascadePbPb&);            // not implemented
  AliAnalysisTaskCheckCascadePbPb& operator=(const AliAnalysisTaskCheckCascadePbPb&); // not implemented
  
  ClassDef(AliAnalysisTaskCheckCascadePbPb, 7);
};

#endif
Commit	Line	Data
f9a6cab5	1	#ifndef ALIANALYSISTASKCHECKCASCADEPBPB_H
	2	#define ALIANALYSISTASKCHECKCASCADEPBPB_H
	3
	4	/* See cxx source for full Copyright notice */
	5
	6	//-----------------------------------------------------------------
	7	// AliAnalysisTaskCheckCascadePbPb class
	8	// Origin AliAnalysisTaskCheckCascade
	9	// This task has four roles :
	10	// 1. QAing the Cascades from ESD and AOD
	11	// Origin: AliAnalysisTaskESDCheckV0 by Boris Hippolyte Nov2007, hippolyt@in2p3.fr
	12	// 2. Prepare the plots which stand as raw material for yield extraction (wi/wo PID)
	13	// 3. Supply an AliCFContainer meant to define the optimised topological selections
	14	// 4. Rough azimuthal correlation study (Eta, Phi)
	15	// Adapted to Cascade : A.Maire Mar2008, antonin.maire@ires.in2p3.fr
	16	// Modified : A.Maire Mar2010, antonin.maire@ires.in2p3.fr
	17	// Modified for PbPb analysis: M. Nicassio Feb 2011, maria.nicassio@ba.infn.it
	18	//-----------------------------------------------------------------
	19
	20	class TList;
	21	class TH1F;
	22	class TH2F;
	23	class TH3F;
	24	class TVector3;
	25	class THnSparse;
	26
f9a6cab5	27	class AliESDEvent;
	28	class AliPhysicsSelection;
	29	class AliCFContainer;
11bcd1e4	30	class AliPIDResponse;
f9a6cab5	31
	32	#include "TString.h"
	33
	34	#include "AliAnalysisTaskSE.h"
	35
	36	class AliAnalysisTaskCheckCascadePbPb : public AliAnalysisTaskSE {
	37	public:
	38	AliAnalysisTaskCheckCascadePbPb();
	39	AliAnalysisTaskCheckCascadePbPb(const char *name);
	40	virtual ~AliAnalysisTaskCheckCascadePbPb();
	41
	42	virtual void UserCreateOutputObjects();
	43	virtual void UserExec(Option_t *option);
	44	/* void DoAngularCorrelation(const Char_t *lCascType,
	45	Double_t lInvMassCascade,
	46	const Int_t *lArrTrackID,
	47	TVector3 &lTVect3MomXi,
	48	Double_t lEtaXi);*/
	49	virtual Int_t DoESDTrackWithTPCrefitMultiplicity(const AliESDEvent *lESDevent);
f9a6cab5	50	virtual void Terminate(Option_t *);
f9a6cab5	51
f9a6cab5	52	void SetAnalysisType (const char* analysisType = "ESD") { fAnalysisType = analysisType; }
	53	void SetRelaunchV0CascVertexers (Bool_t rerunV0CascVertexers = 0 ) { fkRerunV0CascVertexers = rerunV0CascVertexers; }
	54	void SetQualityCutZprimVtxPos (Bool_t qualityCutZprimVtxPos = kTRUE) { fkQualityCutZprimVtxPos = qualityCutZprimVtxPos; }
	55	void SetQualityCutNoTPConlyPrimVtx (Bool_t qualityCutNoTPConlyPrimVtx = kTRUE) { fkQualityCutNoTPConlyPrimVtx = qualityCutNoTPConlyPrimVtx; }
	56	void SetQualityCutTPCrefit (Bool_t qualityCutTPCrefit = kTRUE) { fkQualityCutTPCrefit = qualityCutTPCrefit; }
79313da9	57	void SetQualityCutnTPCcls (Bool_t qualityCutnTPCcls = kTRUE) { fkQualityCutnTPCcls = qualityCutnTPCcls; }
79313da9	58	void SetQualityCutMinnTPCcls (Int_t minnTPCcls = 70 ) { fMinnTPCcls = minnTPCcls; }
f9a6cab5	59	void SetExtraSelections (Bool_t extraSelections = 0 ) { fkExtraSelections = extraSelections; }
	60	void SetCentralityLowLim (Float_t centrlowlim = 0. ) { fCentrLowLim = centrlowlim; }
	61	void SetCentralityUpLim (Float_t centruplim = 100. ) { fCentrUpLim = centruplim; }
	62	void SetCentralityEst (TString centrest = "V0M") { fCentrEstimator = centrest; }
314cd92c	63	void SetUseCleaning (Bool_t usecleaning = kTRUE) { fkUseCleaning = usecleaning; }
f9a6cab5	64	void SetVertexRange (Float_t vtxrange = 0. ) { fVtxRange = vtxrange; }
f27a407b	65	void SetMinptCutOnDaughterTracks (Float_t minptdaughtrks = 0. ) { fMinPtCutOnDaughterTracks = minptdaughtrks; }
433d53fd	66	void SetEtaCutOnDaughterTracks (Float_t etadaughtrks = 0. ) { fEtaCutOnDaughterTracks = etadaughtrks; }
	67	//Setters for the V0 and cascade Vertexer Parameters
	68	void SetV0VertexerMaxChisquare (Double_t lParameter){ fV0VertexerSels[0] = lParameter; }
	69	void SetV0VertexerDCAFirstToPV (Double_t lParameter){ fV0VertexerSels[1] = lParameter; }
	70	void SetV0VertexerDCASecondtoPV (Double_t lParameter){ fV0VertexerSels[2] = lParameter; }
	71	void SetV0VertexerDCAV0Daughters (Double_t lParameter){ fV0VertexerSels[3] = lParameter; }
	72	void SetV0VertexerCosinePA (Double_t lParameter){ fV0VertexerSels[4] = lParameter; }
	73	void SetV0VertexerMinRadius (Double_t lParameter){ fV0VertexerSels[5] = lParameter; }
	74	void SetV0VertexerMaxRadius (Double_t lParameter){ fV0VertexerSels[6] = lParameter; }
	75	void SetCascVertexerMaxChisquare (Double_t lParameter){ fCascadeVertexerSels[0] = lParameter; }
	76	void SetCascVertexerMinV0ImpactParameter (Double_t lParameter){ fCascadeVertexerSels[1] = lParameter; }
	77	void SetCascVertexerV0MassWindow (Double_t lParameter){ fCascadeVertexerSels[2] = lParameter; }
	78	void SetCascVertexerDCABachToPV (Double_t lParameter){ fCascadeVertexerSels[3] = lParameter; }
	79	void SetCascVertexerDCACascadeDaughters (Double_t lParameter){ fCascadeVertexerSels[4] = lParameter; }
	80	void SetCascVertexerCascadeCosinePA (Double_t lParameter){ fCascadeVertexerSels[5] = lParameter; }
	81	void SetCascVertexerCascadeMinRadius (Double_t lParameter){ fCascadeVertexerSels[6] = lParameter; }
	82	void SetCascVertexerCascadeMaxRadius (Double_t lParameter){ fCascadeVertexerSels[7] = lParameter; }
	83
	84
f9a6cab5	85
	86	private:
	87	// Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
	88	// your data member object is created on the worker nodes and streaming is not needed.
	89	// http://root.cern.ch/download/doc/11InputOutput.pdf, page 14
	90
	91
	92	TString fAnalysisType; // "ESD" or "AOD" analysis type
39b18694	93	AliESDtrackCuts *fESDtrackCuts; // ESD track cuts used for primary track definition
f9a6cab5	94	//TPaveText *fPaveTextBookKeeping; // TString to store all the relevant info necessary for book keeping (v0 cuts, cascade cuts, quality cuts, ...)
11bcd1e4	95	AliPIDResponse *fPIDResponse; //! PID response object
f9a6cab5	96
	97	Bool_t fkRerunV0CascVertexers; // Boolean : kTRUE = relaunch both V0 + Cascade vertexers
	98	Bool_t fkQualityCutZprimVtxPos; // Boolean : kTRUE = cut on the prim.vtx z-position
	99	Bool_t fkQualityCutNoTPConlyPrimVtx; // Boolean : kTRUE = prim vtx should be SPD or Tracking vertex
	100	Bool_t fkQualityCutTPCrefit; // Boolean : kTRUE = ask for TPCrefit for the 3 daughter tracks
79313da9	101	Bool_t fkQualityCutnTPCcls; // Boolean : kTRUE = ask for at least n TPC clusters for each daughter track
79313da9	102	Int_t fMinnTPCcls; // minimum number of TPC cluster for daughter tracks
f9a6cab5	103	Bool_t fkExtraSelections; // Boolean : kTRUE = apply tighter selections, before starting the analysis
	104	Float_t fCentrLowLim; // Lower limit for centrality percentile selection
	105	Float_t fCentrUpLim; // Upper limit for centrality percentile selection
	106	TString fCentrEstimator; // string for the centrality estimator
314cd92c	107	Bool_t fkUseCleaning; // Boolean : kTRUE = uses all the cleaning criteria of centrality selections (vertex cut + outliers) otherwise only outliers
f9a6cab5	108	Float_t fVtxRange; // to select events with \|zvtx\|<fVtxRange cm
f27a407b	109	Float_t fMinPtCutOnDaughterTracks; // minimum pt cut on daughter tracks
6fd25d0b	110	Float_t fEtaCutOnDaughterTracks; // pseudorapidity cut on daughter tracks
f9a6cab5	111
433d53fd	112	Double_t fV0VertexerSels[7]; // Array to store the 7 values for the different selections V0 related (if fkRerunV0CascVertexers)
433d53fd	113	Double_t fCascadeVertexerSels[8]; // Array to store the 8 values for the different selections Casc. related (if fkRerunV0CascVertexers)
f9a6cab5	114
39b18694	115	TList *fListHistCascade; //! List of Cascade histograms
f9a6cab5	116
f9a6cab5	117	// - General histos (filled before the trigger selection)
f27a407b	118	TH2F *fHistEvtsInCentralityBinsvsNtracks; //! Events in centrality bins vs N ESDtracks
f27a407b	119	TH1F *fHistCascadeMultiplicityBeforeEvSel; //! Cascade multiplicity distribution
f9a6cab5	120
	121	// - General histos (filled for any triggered event)
	122	TH1F *fHistCascadeMultiplicityForCentrEvt; //! Cascade multiplicity distribution
	123	TH1F *fHistTrackMultiplicityForCentrEvt; //! Track multiplicity distribution (without any cut = include ITS stand-alone + global tracks)
	124	TH1F *fHistTPCrefitTrackMultiplicityForCentrEvt; //! Track multiplicity distribution for tracks with TPCrefit
	125
f9a6cab5	126	// - General histos (filled for events selected in this analysis (\|z(prim. vtx)\| < 10 cm + prim vtx = SPD or Tracking) )
	127	TH1F *fHistCascadeMultiplicityForSelEvt; //! Cascade multiplicity distribution
	128	TH1F *fHistPosBestPrimaryVtxXForSelEvt; //! (best) primary vertex position distribution in x
	129	TH1F *fHistPosBestPrimaryVtxYForSelEvt; //! (best) primary vertex position distribution in y
	130	TH1F *fHistPosBestPrimaryVtxZForSelEvt; //! (best) primary vertex position distribution in z
	131
	132
	133
	134
	135	// - Characteristics for event with >1 cascade : Track Multiplicity, TPC clusters + Prim. vertex positions
	136	TH1F *fHistTPCrefitTrackMultiplicityForCascadeEvt; //! TPCrefit Track multiplicity distribution for event with >1 cascade candidate (NB: after quality sel. within the task)
	137
	138	TH1F *fHistPosV0TPCClusters; //! TPC clusters distribution for Positive V0 daughter track
	139	TH1F *fHistNegV0TPCClusters; //! TPC clusters distribution for Negative V0 daughter track
	140	TH1F *fHistBachTPCClusters; //! TPC clusters distribution for Bachelor track
	141
	142	TH1F *fHistVtxStatus; //! Is there a tracking vertex in the cascade event ?
	143
	144	// Vtx coming from the full tracking, for events containing at least a cascade
	145	TH1F *fHistPosTrkgPrimaryVtxXForCascadeEvt; //! primary vertex position distribution in x
	146	TH1F *fHistPosTrkgPrimaryVtxYForCascadeEvt; //! primary vertex position distribution in y
	147	TH1F *fHistPosTrkgPrimaryVtxZForCascadeEvt; //! primary vertex position distribution in z
	148	TH1F *fHistTrkgPrimaryVtxRadius; //! primary vertex (3D) radius distribution
	149
	150	// Best primary Vtx available, for events containing at least a cascade
	151	TH1F *fHistPosBestPrimaryVtxXForCascadeEvt; //! (best) primary vertex position distribution in x
	152	TH1F *fHistPosBestPrimaryVtxYForCascadeEvt; //! (best) primary vertex position distribution in y
	153	TH1F *fHistPosBestPrimaryVtxZForCascadeEvt; //! (best) primary vertex position distribution in z
	154	TH1F *fHistBestPrimaryVtxRadius; //! (best) primary vertex radius distribution
	155
	156	// Correlation Best Vtx / Full Tracking Vtx
	157	TH2F *f2dHistTrkgPrimVtxVsBestPrimVtx; //! Radius of prim. Vtx from tracks Vs Radius of best Prim. Vtx
	158
	159
	160	// PART 1 : Adavanced QA
	161	// - Typical histos on the variables used for the selection of cascades
	162	TH1F *fHistEffMassXi; //! reconstructed cascade effective mass
	163	TH1F *fHistChi2Xi; //! chi2 value
	164	TH1F *fHistDcaXiDaughters; //! dca between Xi's daughters
	165	TH1F *fHistDcaBachToPrimVertex; //! dca of the bachelor track to primary vertex
	166	TH1F *fHistXiCosineOfPointingAngle; //! cosine of Xi pointing angle in a cascade
	167	TH1F *fHistXiRadius; //! (transverse) radius of the cascade vertex
	168
	169	// - Histos about ~ the "V0 selection part" of the cascade, coming by inheritance from AliESDv0
	170	TH1F *fHistMassLambdaAsCascDghter; //! Test Invariant Mass of Lambda coming from Cascade
	171	TH1F *fHistV0Chi2Xi; //! V0 chi2 distribution, for the V0 associated to a cascade
	172	TH1F *fHistDcaV0DaughtersXi; //! Dca between V0 daughters, for the V0 associated to a cascade
	173	TH1F *fHistDcaV0ToPrimVertexXi; //! Dca of V0 to primary vertex, for the V0 associated to a cascade
988a1924	174	TH1F *fHistV0CosineOfPointingAngle; //! Cosine of V0 pointing angle, for the V0 associated to a cascade
f9a6cab5	175	TH1F *fHistV0RadiusXi; //! V0 (transverse) distance distribution, for the V0 associated to a cascade
	176
	177	TH1F *fHistDcaPosToPrimVertexXi; //! Dca of V0 positive daughter to primary vertex, for the V0 associated to a cascade
	178	TH1F *fHistDcaNegToPrimVertexXi; //! Dca of V0 negative daughter to primary vertex, for the V0 associated to a cascade
	179
	180
	181	// - Effective mass histos for cascades.
	182	TH1F *fHistMassXiMinus; //! reconstructed cascade effective mass, under Xi- hyp.
	183	TH1F *fHistMassXiPlus; //! reconstructed cascade effective mass, under Xi+ hyp.
	184	TH1F *fHistMassOmegaMinus; //! reconstructed cascade effective mass, under Omega- hyp.
	185	TH1F *fHistMassOmegaPlus; //! reconstructed cascade effective mass, under Omega+ hyp.
	186
	187	TH1F *fHistMassWithCombPIDXiMinus; //! reconstructed Xi- effective mass, with bach. comb PID
	188	TH1F *fHistMassWithCombPIDXiPlus; //! reconstructed Xi+ effective mass, with bach. comb PID
	189	TH1F *fHistMassWithCombPIDOmegaMinus; //! reconstructed Omega- effective mass, with bach. comb PID
	190	TH1F *fHistMassWithCombPIDOmegaPlus; //! reconstructed Omega+ effective mass, with bach. comb PID
	191
	192	// - Complements for QA
	193	TH1F *fHistXiTransvMom; //! Xi transverse momentum, around the mass peak of Xi-/+
	194	TH1F *fHistXiTotMom; //! Xi momentum norm, around the mass peak of Xi-/+
	195
	196	TH1F *fHistBachTransvMomXi; //! bachelor transverse momentum, for cand. around the mass peak of Xi-/+
	197	TH1F *fHistBachTotMomXi; //! bachelor momentum norm, for cand. around the mass peak of Xi-/+
f27a407b	198	TH1F *fHistPosTransvMomXi; //! positive daughter transverse momentum, for cand. around the mass peak of Xi-/+
f27a407b	199	TH1F *fHistNegTransvMomXi; //! negative daughter transverse momentum, for cand. around the mass peak of Xi-/+
f9a6cab5	200
	201	TH1F *fHistChargeXi; //! Charge sign of the cascade candidate
	202	TH1F *fHistV0toXiCosineOfPointingAngle; //! Cos. of Pointing angle between the V0 mom and the Xi-V0 vtx line
	203
	204	TH1F *fHistRapXi; //! rapidity of Xi candidates, around the mass peak of Xi-/+
	205	TH1F *fHistRapOmega; //! rapidity of Omega candidates, around the mass peak of Omega-/+
	206	TH1F *fHistEtaXi; //! eta distrib. of all the cascade candidates, around the mass peak of Xi-/+
	207	TH1F *fHistThetaXi; //! theta distrib. of all the cascade candidates, around the mass peak of Xi-/+
	208	TH1F *fHistPhiXi; //! phi distrib. of all the cascade candidates, around the mass peak of Xi-/+
	209
	210	TH2F *f2dHistArmenteros; //! alpha(casc. cand.) Vs PtArm(casc. cand.)
	211
	212	TH2F *f2dHistEffMassLambdaVsEffMassXiMinus; //! Xi- Eff mass Vs V0 Eff mass, under Xi- hyp.
	213	TH2F *f2dHistEffMassXiVsEffMassOmegaMinus; //! Xi- Eff mass Vs Omega- Eff mass, for negative cascades
	214	TH2F *f2dHistEffMassLambdaVsEffMassXiPlus; //! Xi+ Eff mass Vs V0 Eff mass, under Xi+ hyp.
	215	TH2F *f2dHistEffMassXiVsEffMassOmegaPlus; //! Xi+ Eff mass Vs Omega+ Eff mass, for positive cascades
	216
	217	TH2F *f2dHistXiRadiusVsEffMassXiMinus; //! transv. casc. decay radius Vs Xi- Eff mass, under Xi- hyp.
	218	TH2F *f2dHistXiRadiusVsEffMassXiPlus; //! transv. casc. decay radius Vs Xi+ Eff mass, under Xi+ hyp.
	219	TH2F *f2dHistXiRadiusVsEffMassOmegaMinus; //! transv. casc. decay radius Vs Omega- Eff mass, under Omega- hyp.
	220	TH2F *f2dHistXiRadiusVsEffMassOmegaPlus; //! transv. casc. decay radius Vs Omega+ Eff mass, under Omega+ hyp.
	221
	222	TH2F *f2dHistTPCdEdxOfCascDghters; //! TPC Bethe-Bloch curve, populated with the cascade daughters
	223
	224
	225	// PART 2 : TH3F needed for pt spectrum and yield extraction
	226	// Without any PID
	227	/* TH3F *f3dHistXiPtVsEffMassVsYXiMinus; //! casc. transv. momemtum Vs Xi- Eff mass Vs Y
	228	TH3F *f3dHistXiPtVsEffMassVsYXiPlus; //! casc. transv. momemtum Vs Xi+ Eff mass Vs Y
	229	TH3F *f3dHistXiPtVsEffMassVsYOmegaMinus; //! casc. transv. momemtum Vs Omega- Eff mass Vs Y
	230	TH3F *f3dHistXiPtVsEffMassVsYOmegaPlus; //! casc. transv. momemtum Vs Omega+ Eff mass Vs Y
	231	*/
	232	// Compilation of all PID plots (3D = casc. transv. momemtum Vs Casc Eff mass Vs Y), stored into an AliCFContainer
	233	AliCFContainer *fCFContCascadePIDXiMinus; //! for Xi- : Container to store any 3D histos with the different PID flavours
	234	AliCFContainer *fCFContCascadePIDXiPlus; //! for Xi+ : Container to store any 3D histos with the different PID flavours
	235	AliCFContainer *fCFContCascadePIDOmegaMinus; //! for Omega-: Container to store any 3D histos with the different PID flavours
	236	AliCFContainer *fCFContCascadePIDOmegaPlus; //! for Omega+: Container to store any 3D histos with the different PID flavours
	237
	238
	239
	240	// PART 3 : Towards the optimisation of topological selections / systematics
	241	AliCFContainer *fCFContCascadeCuts; //! Container meant to store all the relevant distributions corresponding to the cut variables
	242
	243
	244	// PART 4 : Azimuthal correlation study
	245	/* THnSparseF *fHnSpAngularCorrXiMinus; //! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
	246	THnSparseF *fHnSpAngularCorrXiPlus; //! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
	247	THnSparseF *fHnSpAngularCorrOmegaMinus; //! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
	248	THnSparseF *fHnSpAngularCorrOmegaPlus; //! Delta Phi(Casc,any trck) Vs Delta Eta(Casc,any trck) Vs Casc Pt Vs Pt of the tracks Vs Eff Mass
	249	*/
	250	TH1F *fV0Ampl; //! histo to check the V0 amplitude distribution
	251
	252	TH2F *fHistDcaXiDaughtersvsInvMass; //! cut variables vs inv. mass
	253	TH2F *fHistDcaBachToPrimVertexvsInvMass; //! cut variables vs inv. mass
	254	TH2F *fHistXiCosineOfPointingAnglevsInvMass; //! cut variables vs inv. mass
	255	TH2F *fHistMassLambdaAsCascDghtervsInvMass; //! cut variables vs inv. mass
	256	TH2F *fHistDcaV0DaughtersXivsInvMass; //! cut variables vs inv. mass
	257	TH2F *fHistDcaV0ToPrimVertexXivsInvMass; //! cut variables vs inv. mass
	258
6fd25d0b	259	// Control plots for reco pseudorapidity of daughter tracks (Xi- only)
	260
	261	TH1F *fHistEtaBachXiM; //! bachelor pseudorapidity
	262	TH1F *fHistEtaPosXiM; //! positive daughter pseudorapidity
	263	TH1F *fHistEtaNegXiM; //! negative daughter pseudorapidity
	264
f9a6cab5	265
	266
	267	AliAnalysisTaskCheckCascadePbPb(const AliAnalysisTaskCheckCascadePbPb&); // not implemented
	268	AliAnalysisTaskCheckCascadePbPb& operator=(const AliAnalysisTaskCheckCascadePbPb&); // not implemented
	269
6fd25d0b	270	ClassDef(AliAnalysisTaskCheckCascadePbPb, 7);
f9a6cab5	271	};
	272
	273	#endif