[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / Kinks / AliAnalysisKinkESDat.h

#ifndef ALIANALYSISKINKESDat_H
#define ALIANALYSISKINKESDat_H

/*  See cxx source for full Copyright notice */

//-----------------------------------------------------------------
//                 AliAnalysisKinkESDat class
//         This task is an example of an analysis task
//                  for kink topology Study
//          Authors: Martha Spyropoulou-Stassinaki
//           and members of the Greek group at the
//          Physics Department of Athens University
//                    mspyrop@phys.uoa.gr
//-----------------------------------------------------------------
class AliPIDResponse;
class AliESDVertex;
class AliESDtrack;
class TF1;
class TH1F;
class TH2F;
class TH3F;
class TH1D;
class TH2D;
class TList;
class AliESDtrackCuts;

#include "AliAnalysisTaskSE.h"

class AliAnalysisKinkESDat : public AliAnalysisTaskSE {
 public:
  AliAnalysisKinkESDat(const char *name = "AliAnalysisKinkESDat");
  virtual ~AliAnalysisKinkESDat() {}

  virtual void   UserCreateOutputObjects();
  virtual void   UserExec(Option_t *option);
  virtual void   Terminate(Option_t *);
  
  const AliESDVertex *GetEventVertex(const AliESDEvent* esd) const;
//           Marek multiplicity bins
  void SetMulCut(Int_t low, Int_t up){fLowMulcut=low;fUpMulcut=up;}	  
  void SetKinkRadius(Float_t lRadiusKLow, Float_t lRadiusKUp)  { fKinkRadLow=lRadiusKLow; fKinkRadUp=lRadiusKUp;}	  

  void SetNClusterCut(Int_t lowCluster){fLowCluster=lowCluster;}	  
  void SetQtCut(Float_t   lowQt){fLowQt=lowQt;}	  
         void SetYKRange(Float_t  RapidityK){fRapiK=RapidityK;}
 private:
   TH1F        *fHistPtESD; //Pt spectrum of all ESD inside eta, Pt cuts
   TH1F        *fHistPt; //Pt spectrum of all ESD tracks
   TH1F        *fHistQtAll; //Qt spectrum of all kinks
   TH1F        *fHistQt1; //Qt spectrum of Kaon selected sample
   TH1F        *fHistQt2; //Qt spectrum in Qt region of kaons
   TH1F        *fHistPtKaon; //Pt Kaon spectrum of clean sample
   TH1F        *fHistPtKPDG; //Pt Kaon spectrum , confirmed by  PDG,inside kaon Qt region
   TH1F        *fHistEta; //Eta spectrum of all kinks
   TH1F        *fHistEtaK; //Eta spectrum of kaons selected by kink topology
   TH1F        *fptKMC; //Pt Kaon spectrum MC, inside eta and pt cuts 
   TH1F        *fMultiplMC; //charge multipl MC 
   TH1F        *fESDMult; //ESD charged mult
   TH1F        *fgenpt; //Pt Kaon-Kink->mu  spectrum , MC, inside eta, Pt, radius cuts
   TH1F        *frad; //radius of kinks,  MC , inside the eta nad Pt cuts 
   TH1F        *fKinkKaon; //Pt of PDG Kaons inside the selcted ones by the KInk topology 
   TH1F        *fKinKRbn; //Pt of PDG Kaons inside the selcted ones by the KInk topology 
   TH1F        *fKinkKaonBg; //Pt of the BG inside the kink-Kaon identified spectrum
   TH1F        *fM1kaon; //inv mass of kink-tracks taken as kaons decaying to  mu + neutrino
   TH1F        *fPtKink; //Pt  spectrum   of all kinks  from track bank
   TH1F        *fptKink; //Pt  spectrum of all kinks from kink bank
   TH2F        *fAngMomK; // Decay angle vrs Mother Mom for pdg kaons
   TH2F        *fAngMomPi; // Decay angle vrs Mother Mom for pdg pions
   TH2F        *fAngMomKC; //Decay angle vrs Mother Mom for pdg kaons, inside the selected sample
   TH1F        *fMultESDK; //ESD charged mult
   TH1F        *fMultMCK; //MC K charged mult
   TH2F        *fSignPtNcl;//signPt vrs number of clusters in TPC for kaons from kink sele sample
   TH2F        *fSignPtEta;//signPt vrs Eta  in TPC for kaons from kink sele sample
   TH2F        *fEtaNcl;//Eta    vrs Nclu in TPC for kaons from kink sele sample
   TH1F        *fSignPt;//signPt  in TPC for kaons from kink sele sample
   TH2F        *fChi2NclTPC;//chi2 vrs TPC Nclusters for kaons from kink sele sample
   TH1F        *fRatChi2Ncl;// Ratio chi2/ Ncl  TPC  for kaons from kink sele sample
   TH2F        *fRadiusNcl;//kink  Radius      Ncl  TPC  for kaons from kink sele sample
   TH2F        *fTPCSgnlP;//kink  Radius      Ncl  TPC  for kaons from kink sele sample
   TH2F        *fTPCSgnlPa;//kink  Radius      Ncl  TPC  for  kink sele sample
   TH1D        *fRpr; // Radius of VTX at Y , X plane              
   TH1D        *fZpr; //Z distrio of main vertex                  
   TH1D        *fdcatoVxXY; //dca to Vertex XY  distrio                   
   TH1D        *fnSigmToVx; //nSigma to Vertex  distrio of main vertex                  
   TH2F        *fKinkMothDau; //Mother vrs Daughter                                       
   TH2F        *fZvXv; //two dime of Z vrs X of vtx main           
   TH2F        *fZvYv; // two dime of Z vrs Y of vtx main           
   TH2F        *fXvYv; // two dime of X vrs Y of main tracks vtx main           
   TH1F        *fHistPtKaoP; //Pt Kaon spectrum of clean sample pos
   TH1F        *fHistPtKaoN; //Pt Kaon spectrum of clean sample neg
   TH1F        *frapiKESD;// rapidi K      
   TH1F        *flifetime;//radius of kinks,  MC , inside the eta nad Pt cuts 
   TH1F        *fradLK;//Length  of kinks,  MC , inside the eta nad Pt cuts 
   TH3F        *fradPtRpDt;//radius of kinks,  MC , inside the eta nad Pt cuts 
   TH1F        *fInvMuNuAll;//radius of kinks,  MC , inside the eta nad Pt cuts 
   TH2F        *fQtInvM;// 
   TH1F        *fDCAkink;//!MC dcs kink
   TH2F        *fPosiKink;//!MC position  kink
   TH2F        *fPosiKinkK;//!MC position  kink
   TH2F        *fPosiKinKXZ;//!MC position  kink
   TH2F        *fPosiKinKYZ;//!MC position  kink
   TH2F        *fPosiKinKBg;//!MC position  kink
   TH2F        *fQtMothP;//!qt vrs p mother  
   TH2F        *fTPCSgnlPtpc;//Kink mother moment vrs TPC signal                 
   TH2F        *fTPCMomNSgnl;//kink  mother TPC momentum vrs nsigmas of dEdx                    
   TH2F        *fMothKinkMomSgnl;//kink  mother TPC momentum vrs nsigmas of dEdx                    
   TH1F        *fNSigmTPC;//kink  mother TPC momentum vrs nsigmas of dEdx                    
   TH2F        *fTPCSgnlKinkDau;//Kink mother moment vrs TPC signal                 
   TH1F        *fPtKinkPos; //Pos K Pt  spectrum   of all kinks  from track bank, K0 bins
   TH1F        *fPtKinkNeg; //Neg K Pt  spectrum   of all kinks  from track bank, K0 bins
   TH2F        *fRadNclCln;//kink  Radius      Ncl  TPC  for kaons from kink clean sample
   TH1F        *fRatioCrossedRows; //ratio  crossed rows                                           
   TH1F        *fRatioCrossedRowsKink; //ratio  crossed rows  for kinks                                         
   TH2F        *fRadiusPt;//kinks,  Radius      vs Pt                                        
   TH2F        *fRadiusPtcln;//kinks,  Radius      vs Pt    for clean kaons                                     
   TH2F        *fInvMassMuNuPt;//kinks,Invariant Mass MuNu     vs Pt                                         
   TH2F        *fInvMassMuNuPtAll;//kinks,Invariant Mass MuNu     vs Pt                                         
   TH1F        *fPtCut1; //K Pt  spectrum   of all kinks  from track bank, K0 bins
   TH1F        *fPtCut2; //K Pt  spectrum   of all kinks  from track bank, K0 bins
   TH1F        *fPtCut3; //K Pt  spectrum   of all kinks  from track bank, K0 bins
   TH2F        *fAngMomKKinks;//kinks,  Angle vs Momentum  for K-kinks                                      


   TF1         *f1;
   TF1         *f2;
  TList        *fListOfHistos; // list of histos
//          Marek Multiplicity
Int_t fLowMulcut;  // 
Int_t fUpMulcut;
Float_t fKinkRadUp;
Float_t fKinkRadLow;
Int_t fLowCluster;
Float_t  fLowQt;       
    Float_t  fRapiK;
AliESDtrackCuts* fCutsMul;
AliESDtrackCuts* fMaxDCAtoVtxCut;  
AliPIDResponse *fPIDResponse;     //! PID response object

  AliAnalysisKinkESDat(const AliAnalysisKinkESDat&); // not implemented
  AliAnalysisKinkESDat& operator=(const AliAnalysisKinkESDat&); // not implemented

  ClassDef(AliAnalysisKinkESDat, 1); // example of analysis
};

#endif
Commit	Line	Data
828a6e0c	1	#ifndef ALIANALYSISKINKESDat_H
	2	#define ALIANALYSISKINKESDat_H
	3
	4	/* See cxx source for full Copyright notice */
	5
	6	//-----------------------------------------------------------------
	7	// AliAnalysisKinkESDat class
	8	// This task is an example of an analysis task
	9	// for kink topology Study
	10	// Authors: Martha Spyropoulou-Stassinaki
	11	// and members of the Greek group at the
	12	// Physics Department of Athens University
	13	// mspyrop@phys.uoa.gr
	14	//-----------------------------------------------------------------
fcd98e58	15	class AliPIDResponse;
828a6e0c	16	class AliESDVertex;
	17	class AliESDtrack;
	18	class TF1;
	19	class TH1F;
	20	class TH2F;
	21	class TH3F;
	22	class TH1D;
	23	class TH2D;
	24	class TList;
	25	class AliESDtrackCuts;
	26
	27	#include "AliAnalysisTaskSE.h"
	28
	29	class AliAnalysisKinkESDat : public AliAnalysisTaskSE {
	30	public:
828a6e0c	31	AliAnalysisKinkESDat(const char *name = "AliAnalysisKinkESDat");
	32	virtual ~AliAnalysisKinkESDat() {}
	33
	34	virtual void UserCreateOutputObjects();
	35	virtual void UserExec(Option_t *option);
	36	virtual void Terminate(Option_t *);
	37
	38	const AliESDVertex GetEventVertex(const AliESDEvent esd) const;
	39	// Marek multiplicity bins
	40	void SetMulCut(Int_t low, Int_t up){fLowMulcut=low;fUpMulcut=up;}
b2fb60b8	41	void SetKinkRadius(Float_t lRadiusKLow, Float_t lRadiusKUp) { fKinkRadLow=lRadiusKLow; fKinkRadUp=lRadiusKUp;}
828a6e0c	42
5220a87d	43	void SetNClusterCut(Int_t lowCluster){fLowCluster=lowCluster;}
5220a87d	44	void SetQtCut(Float_t lowQt){fLowQt=lowQt;}
144aaa1c	45	void SetYKRange(Float_t RapidityK){fRapiK=RapidityK;}
828a6e0c	46	private:
	47	TH1F *fHistPtESD; //Pt spectrum of all ESD inside eta, Pt cuts
	48	TH1F *fHistPt; //Pt spectrum of all ESD tracks
	49	TH1F *fHistQtAll; //Qt spectrum of all kinks
	50	TH1F *fHistQt1; //Qt spectrum of Kaon selected sample
	51	TH1F *fHistQt2; //Qt spectrum in Qt region of kaons
	52	TH1F *fHistPtKaon; //Pt Kaon spectrum of clean sample
	53	TH1F *fHistPtKPDG; //Pt Kaon spectrum , confirmed by PDG,inside kaon Qt region
	54	TH1F *fHistEta; //Eta spectrum of all kinks
	55	TH1F *fHistEtaK; //Eta spectrum of kaons selected by kink topology
	56	TH1F *fptKMC; //Pt Kaon spectrum MC, inside eta and pt cuts
	57	TH1F *fMultiplMC; //charge multipl MC
	58	TH1F *fESDMult; //ESD charged mult
	59	TH1F *fgenpt; //Pt Kaon-Kink->mu spectrum , MC, inside eta, Pt, radius cuts
	60	TH1F *frad; //radius of kinks, MC , inside the eta nad Pt cuts
	61	TH1F *fKinkKaon; //Pt of PDG Kaons inside the selcted ones by the KInk topology
	62	TH1F *fKinKRbn; //Pt of PDG Kaons inside the selcted ones by the KInk topology
	63	TH1F *fKinkKaonBg; //Pt of the BG inside the kink-Kaon identified spectrum
	64	TH1F *fM1kaon; //inv mass of kink-tracks taken as kaons decaying to mu + neutrino
828a6e0c	65	TH1F *fPtKink; //Pt spectrum of all kinks from track bank
828a6e0c	66	TH1F *fptKink; //Pt spectrum of all kinks from kink bank
828a6e0c	67	TH2F *fAngMomK; // Decay angle vrs Mother Mom for pdg kaons
	68	TH2F *fAngMomPi; // Decay angle vrs Mother Mom for pdg pions
	69	TH2F *fAngMomKC; //Decay angle vrs Mother Mom for pdg kaons, inside the selected sample
	70	TH1F *fMultESDK; //ESD charged mult
	71	TH1F *fMultMCK; //MC K charged mult
	72	TH2F *fSignPtNcl;//signPt vrs number of clusters in TPC for kaons from kink sele sample
	73	TH2F *fSignPtEta;//signPt vrs Eta in TPC for kaons from kink sele sample
	74	TH2F *fEtaNcl;//Eta vrs Nclu in TPC for kaons from kink sele sample
	75	TH1F *fSignPt;//signPt in TPC for kaons from kink sele sample
	76	TH2F *fChi2NclTPC;//chi2 vrs TPC Nclusters for kaons from kink sele sample
	77	TH1F *fRatChi2Ncl;// Ratio chi2/ Ncl TPC for kaons from kink sele sample
	78	TH2F *fRadiusNcl;//kink Radius Ncl TPC for kaons from kink sele sample
	79	TH2F *fTPCSgnlP;//kink Radius Ncl TPC for kaons from kink sele sample
	80	TH2F *fTPCSgnlPa;//kink Radius Ncl TPC for kink sele sample
	81	TH1D *fRpr; // Radius of VTX at Y , X plane
	82	TH1D *fZpr; //Z distrio of main vertex
	83	TH1D *fdcatoVxXY; //dca to Vertex XY distrio
	84	TH1D *fnSigmToVx; //nSigma to Vertex distrio of main vertex
	85	TH2F *fKinkMothDau; //Mother vrs Daughter
	86	TH2F *fZvXv; //two dime of Z vrs X of vtx main
	87	TH2F *fZvYv; // two dime of Z vrs Y of vtx main
	88	TH2F *fXvYv; // two dime of X vrs Y of main tracks vtx main
828a6e0c	89	TH1F *fHistPtKaoP; //Pt Kaon spectrum of clean sample pos
	90	TH1F *fHistPtKaoN; //Pt Kaon spectrum of clean sample neg
	91	TH1F *frapiKESD;// rapidi K
	92	TH1F *flifetime;//radius of kinks, MC , inside the eta nad Pt cuts
	93	TH1F *fradLK;//Length of kinks, MC , inside the eta nad Pt cuts
	94	TH3F *fradPtRpDt;//radius of kinks, MC , inside the eta nad Pt cuts
	95	TH1F *fInvMuNuAll;//radius of kinks, MC , inside the eta nad Pt cuts
	96	TH2F *fQtInvM;//
	97	TH1F *fDCAkink;//!MC dcs kink
	98	TH2F *fPosiKink;//!MC position kink
	99	TH2F *fPosiKinkK;//!MC position kink
	100	TH2F *fPosiKinKXZ;//!MC position kink
	101	TH2F *fPosiKinKYZ;//!MC position kink
	102	TH2F *fPosiKinKBg;//!MC position kink
	103	TH2F *fQtMothP;//!qt vrs p mother
ad1036b2	104	TH2F *fTPCSgnlPtpc;//Kink mother moment vrs TPC signal
	105	TH2F *fTPCMomNSgnl;//kink mother TPC momentum vrs nsigmas of dEdx
	106	TH2F *fMothKinkMomSgnl;//kink mother TPC momentum vrs nsigmas of dEdx
	107	TH1F *fNSigmTPC;//kink mother TPC momentum vrs nsigmas of dEdx
	108	TH2F *fTPCSgnlKinkDau;//Kink mother moment vrs TPC signal
	109	TH1F *fPtKinkPos; //Pos K Pt spectrum of all kinks from track bank, K0 bins
	110	TH1F *fPtKinkNeg; //Neg K Pt spectrum of all kinks from track bank, K0 bins
	111	TH2F *fRadNclCln;//kink Radius Ncl TPC for kaons from kink clean sample
8e57b720	112	TH1F *fRatioCrossedRows; //ratio crossed rows
	113	TH1F *fRatioCrossedRowsKink; //ratio crossed rows for kinks
	114	TH2F *fRadiusPt;//kinks, Radius vs Pt
	115	TH2F *fRadiusPtcln;//kinks, Radius vs Pt for clean kaons
	116	TH2F *fInvMassMuNuPt;//kinks,Invariant Mass MuNu vs Pt
f9016b0c	117	TH2F *fInvMassMuNuPtAll;//kinks,Invariant Mass MuNu vs Pt
8e57b720	118	TH1F *fPtCut1; //K Pt spectrum of all kinks from track bank, K0 bins
	119	TH1F *fPtCut2; //K Pt spectrum of all kinks from track bank, K0 bins
	120	TH1F *fPtCut3; //K Pt spectrum of all kinks from track bank, K0 bins
	121	TH2F *fAngMomKKinks;//kinks, Angle vs Momentum for K-kinks
	122
	123
828a6e0c	124
	125	TF1 *f1;
	126	TF1 *f2;
	127	TList *fListOfHistos; // list of histos
	128	// Marek Multiplicity
	129	Int_t fLowMulcut; //
	130	Int_t fUpMulcut;
c27ffa18	131	Float_t fKinkRadUp;
c27ffa18	132	Float_t fKinkRadLow;
5220a87d	133	Int_t fLowCluster;
5220a87d	134	Float_t fLowQt;
144aaa1c	135	Float_t fRapiK;
828a6e0c	136	AliESDtrackCuts* fCutsMul;
8e57b720	137	AliESDtrackCuts* fMaxDCAtoVtxCut;
fcd98e58	138	AliPIDResponse *fPIDResponse; //! PID response object
828a6e0c	139
	140	AliAnalysisKinkESDat(const AliAnalysisKinkESDat&); // not implemented
	141	AliAnalysisKinkESDat& operator=(const AliAnalysisKinkESDat&); // not implemented
	142
	143	ClassDef(AliAnalysisKinkESDat, 1); // example of analysis
	144	};
	145
ad1036b2	146	#endif