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

#ifndef ALIANALYSISKINKESDMC_H
#define ALIANALYSISKINKESDMC_H

/*  See cxx source for full Copyright notice */

//-----------------------------------------------------------------
//                 AliAnalysisKinkESDMC 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 AliESDEvent;
class AliESDVertex;
class AliESDtrack;
class TF1;
class TH1F;
class TH2F;
class TH1D;
class TH2D;
class TList;
class AliESDtrackCuts;
class AliPhysicsSelection;

#include "AliAnalysisTaskSE.h"

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

  virtual void   UserCreateOutputObjects();
  virtual void   UserExec(Option_t *option);
  virtual void   Terminate(Option_t *);
  
  const AliESDVertex *GetEventVertex(const AliESDEvent* esd) const;
//  const AliESDVertex *GetEventVertex(AliESDEvent* esd) ;
  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;}

 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        *frad; //!radius of kinks,  MC , inside the eta nad Pt cuts 
   TH1F        *fradMC; //!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        *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        *fgenPtEtR; //!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
   TH1F        *fPtKink; //!Pt  spectrum   of all kinks  from track bank
   TH2F        *fcodeH; //!PDG code(mother)  vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake)
   TH2F        *fdcodeH; //!inks, code  vrs dcode of BG,if mother code is 321 and daughter code > 
   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        *fSignPtEtaMC;//!signPt vrs Eta  in TPC for kaons from kink sele sample
   TH1F        *fSignPtMC;//!signPt   in TPC for kaons 
   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;//! Radis  f kink vetex   for kaons from kink sele sample
   TH2F        *fTPCSgnlP;//! TPC de/dx signal      for kaons from kink sele sample
   TH2F        *fTPCSgnlPa;//! TPC de/dx signal      for  kink  sample
   TH1F        *fSignPtGen;//!signPt   in TPC for kaonsgenerated 
   TH1D        *fRpr;//! Radius of VTX at Y , X plane
   TH1D        *fZpr;//!Z distrio of main vertex                  
   TH1D        *fdcatoVxXY;//! dca to Vertex XY distr          
   TH1F        *fMCEtaKaon;//!MC eta for kaons                                     
   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        *fPtPrKink;//! pt of Primary PDG kaons inside the selected ones by the kink topology              
   TH1F        *fgenPtEtRP;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
   TH1F        *fgenPtEtRN;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
   TH1F        *fkinkKaonP;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
   TH1F        *fkinkKaonN;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
   TH1F        *frapidESDK;//!MC ESD K  rapidity distr.                                                                   
   TH1F        *frapidKMC;//!MC MC K rapidity dis
   TH1F        *fPtKPlMC;//!MC MC K rapidity dis
   TH1F        *fPtKMnMC;//!MC MC K rapidity dis
   TH1F        *fHistPtKaoP;//!MC MC K rapidity dis
   TH1F        *fHistPtKaoN;//!MC MC K rapidity dis
   TH1F        *fHiPtKPDGP;//!MC MC K rapidity dis
   TH1F        *fHiPtKPDGN;//!MC MC K rapidity dis
   TH1F        *fKinKBGP;//!MC MC K rapidity dis
   TH1F        *fKinKBGN;//!MC MC K rapidity dis
   TH1F        *fQtKMu;//!MC MC K Qt K to mu
   TH1F        *fQtKPi;//!MC MC K Qt K to mu
   TH1F        *fQtKEl;//!MC MC K Qt K to mu
   TH1F        *fFakepipi;//!MC Fake pipi
   TH1F        *fFakeKPi;//!MC Fake Kpi
   TH1F        *fDCAkink;//!MC dcs kink
   TH1F        *fDCAkinkBG;//!MC dcs kink
   TH2F        *fPosiKink;//!MC position  kink
   TH2F        *fPosiKinkK;//!MC position  kink
   TH2F        *fPosiKinKXZ;//!MC position  kink
   TH2F        *fPosiKinKYZ;//!MC position  kink
   TH2F        *fPosiKinKBgZY;//!MC position  kink
   TH2F        *fcode2;//!PDG code(mother)  vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake)
   TH2F        *fcode4;//!PDG code(mother)  vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake)
   TH2F        *fZkinkZDau;//!   z-position of kink z position of daughter bg                       )
   TH1F        *fQtKMuMC;//!MC MC K Qt K to mu
   TH1F        *fQtKElMC;//!MC MC K Qt K to mu
   TH1F        *fQtKPiMC;//!MC MC K Qt K to mu
   TH1F        *fQtK3PiP;//!EDS   K Qt K to 3Pi
   TH1F        *fQtK3PiM;//!EDS   K Qt K to 3Pi
   TH2F        *fmaxAngMomKmu; //!Decay angle vrs Mother Mom for pdg kaons, inside the selected sample
   TH2F        *fPosiKinKBgZX;//!MC position  kink
   TH2F        *fPosiKinKBgXY;//!MC position  kink
   TH1F        *fMinvPi;//!MC R life time     
   TH1F        *fMinvKa;//!MC R life time     
   TH1F        *fMinvPr;//!MC R life time     
   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 
   TH2F        *fcodeDau1;//!PDG code(mother)  vrs PDG dcode(daughter) of kinks   daughters          
   TH2F        *fcodeDau2;//!PDG code(mother)  vrs PDG dcode(daughter) of kinks   daughters , Bg study
   TH2F        *fMothKinkMomSgnlD;//kink  mother TPC momentum vrs nsigmas of dEdx                    
   TH1F        *fInvMassMuNuAll;//kinks,  Inv Mass all kinks  MuNu                                       
   TH2F        *fInvMassMuNuPt;//kinks,Invariant Mass MuNu     vs Pt                                         
   TH2F        *fInvMassMuNuPtAll;//kinks,Invariant Mass MuNu     vs Pt      22/8/2013                                   
   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        *fRadiusPtKaon;//kinks,  Radius      vs Pt    for clean kaons                                  
   TH2F        *fRadiusPtPion;//kinks,  Radius      vs Pt    for clean kaons                                  
   TH2F        *fRadiusPtFake;//kinks,  Radius      vs Pt    for clean kaons                                  
   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                                      
    TH1F        *flengthMCK;//!MC R life time
   TH1F        *flifetiMCK;//!MC R life time
   TH1F        *flifetim2;//!MC R life time
   TH1F        *fLHelESDK;//!MC R life time
   TH1F        *flifeInt;//!MC R life time
   TH1F        *flifeYuri;//!MC R life time
   TH1F        *flenYuri;//!MC R life time
   TH1F        *flenTrRef;//!MC R life time
   TH1F        *flifeSmall;//!MC R life time
   TH1F        *flifetime;//!MC R life time
   TH1F        *flifTiESDK;//!MC R life time
   TH1F        *flifeKink;//!MC R life time
   TH1F        *flenHelx;//!MC R life time
   TH3F        *fradPtRapMC;//!MC R life time
   TH3F        *fradPtRapDC;//!MC R life time
   TH3F        *fradPtRapESD;//!MC R life time
   TH2F        *fRadNclcln;//!MC R life time


   TF1         *f1;
   TF1         *f2;
  TList        *fListOfHistos; //! list of histos

 Int_t fLowMulcut; //
Int_t fUpMulcut;
Int_t fKinkRadUp;
Int_t fKinkRadLow;
    Int_t fLowCluster;
Float_t  fLowQt;

AliESDtrackCuts*  fCutsMul;

     AliESDtrackCuts* fMaxDCAtoVtxCut;
 AliPIDResponse *fPIDResponse;     //! PID response object

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

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

#endif
Commit	Line	Data
10eaad41	1	#ifndef ALIANALYSISKINKESDMC_H
	2	#define ALIANALYSISKINKESDMC_H
	3
	4	/* See cxx source for full Copyright notice */
	5
	6	//-----------------------------------------------------------------
	7	// AliAnalysisKinkESDMC 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	//-----------------------------------------------------------------
	15
a0271296	16	class AliPIDResponse;
a0271296	17	class AliESDEvent;
be1a7181	18	class AliESDVertex;
894840ad	19	class AliESDtrack;
894840ad	20	class TF1;
10eaad41	21	class TH1F;
	22	class TH2F;
	23	class TH1D;
894840ad	24	class TH2D;
92adf4f6	25	class TList;
a0271296	26	class AliESDtrackCuts;
a0271296	27	class AliPhysicsSelection;
85f91829	28
894840ad	29	#include "AliAnalysisTaskSE.h"
894840ad	30
10eaad41	31	class AliAnalysisKinkESDMC : public AliAnalysisTaskSE {
10eaad41	32	public:
92adf4f6	33	AliAnalysisKinkESDMC(const char *name = "AliAnalysisKinkESDMC");
10eaad41	34	virtual ~AliAnalysisKinkESDMC() {}
92adf4f6	35
	36	virtual void UserCreateOutputObjects();
	37	virtual void UserExec(Option_t *option);
10eaad41	38	virtual void Terminate(Option_t *);
92adf4f6	39
10eaad41	40	const AliESDVertex GetEventVertex(const AliESDEvent esd) const;
a0271296	41	// const AliESDVertex GetEventVertex(AliESDEvent esd) ;
	42	void SetMulCut(Int_t low, Int_t up){fLowMulcut=low;fUpMulcut=up;}
	43
	44	void SetKinkRadius(Float_t lRadiusKLow, Float_t lRadiusKUp) { fKinkRadLow=lRadiusKLow; fKinkRadUp=lRadiusKUp;}
	45
10eaad41	46
5220a87d	47	void SetNClusterCut(Int_t lowCluster){fLowCluster=lowCluster;}
	48	void SetQtCut(Float_t lowQt){fLowQt=lowQt;}
	49
10eaad41	50	private:
a0271296	51	TH1F *fHistPtESD; //!Pt spectrum of all ESD inside eta, Pt cuts
	52	TH1F *fHistPt; //!Pt spectrum of all ESD tracks
	53	TH1F *fHistQtAll; //!Qt spectrum of all kinks
	54	TH1F *fHistQt1; //!Qt spectrum of Kaon selected sample
	55	TH1F *fHistQt2; //!Qt spectrum in Qt region of kaons
	56	TH1F *fHistPtKaon; //!Pt Kaon spectrum of clean sample
	57	TH1F *fHistPtKPDG; //!Pt Kaon spectrum , confirmed by PDG,inside kaon Qt region
	58	TH1F *fHistEta; //!Eta spectrum of all kinks
	59	TH1F *fHistEtaK; //!Eta spectrum of kaons selected by kink topology
	60	TH1F *fptKMC; //!Pt Kaon spectrum MC, inside eta and pt cuts
	61	TH1F *fMultiplMC; //!charge multipl MC
	62	TH1F *fESDMult; //!ESD charged mult
	63	TH1F *frad; //!radius of kinks, MC , inside the eta nad Pt cuts
	64	TH1F *fradMC; //!radius of kinks, MC , inside the eta nad Pt cuts
	65	TH1F *fKinkKaon; //!Pt of PDG Kaons inside the selcted ones by the KInk topology
	66	TH1F *fKinkKaonBg; //!Pt of the BG inside the kink-Kaon identified spectrum
	67	TH1F *fM1kaon; //!inv mass of kink-tracks taken as kaons decaying to mu + neutrino
	68	TH1F *fgenPtEtR; //!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
	69	TH1F *fPtKink; //!Pt spectrum of all kinks from track bank
	70	TH2F *fcodeH; //!PDG code(mother) vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake)
	71	TH2F *fdcodeH; //!inks, code vrs dcode of BG,if mother code is 321 and daughter code >
	72	TH2F *fAngMomK; //! Decay angle vrs Mother Mom for pdg kaons
	73	TH2F *fAngMomPi; //! Decay angle vrs Mother Mom for pdg pions
	74	TH2F *fAngMomKC; //!Decay angle vrs Mother Mom for pdg kaons, inside the selected sample
	75	TH1F *fMultESDK; //!ESD charged mult
	76	TH1F *fMultMCK; //!MC K charged mult
	77	TH2F *fSignPtNcl;//!signPt vrs number of clusters in TPC for kaons from kink sele sample
	78	TH2F *fSignPtEta;//!signPt vrs Eta in TPC for kaons from kink sele sample
	79	TH2F *fSignPtEtaMC;//!signPt vrs Eta in TPC for kaons from kink sele sample
	80	TH1F *fSignPtMC;//!signPt in TPC for kaons
	81	TH2F *fEtaNcl;//!Eta vrs Nclu in TPC for kaons from kink sele sample
	82	TH1F *fSignPt;//!signPt in TPC for kaons from kink sele sample
	83	TH2F *fChi2NclTPC;//!chi2 vrs TPC Nclusters for kaons from kink sele sample
	84	TH1F *fRatChi2Ncl;//! Ratio chi2/ Ncl TPC for kaons from kink sele sample
	85	TH2F *fRadiusNcl;//! Radis f kink vetex for kaons from kink sele sample
	86	TH2F *fTPCSgnlP;//! TPC de/dx signal for kaons from kink sele sample
	87	TH2F *fTPCSgnlPa;//! TPC de/dx signal for kink sample
	88	TH1F *fSignPtGen;//!signPt in TPC for kaonsgenerated
	89	TH1D *fRpr;//! Radius of VTX at Y , X plane
	90	TH1D *fZpr;//!Z distrio of main vertex
	91	TH1D *fdcatoVxXY;//! dca to Vertex XY distr
	92	TH1F *fMCEtaKaon;//!MC eta for kaons
	93	TH2F *fZvXv;//! two dime of Z vrs X of vtx main
	94	TH2F *fZvYv;//! two dime of Z vrs Y of vtx main
	95	TH2F *fXvYv;//! two dime of X vrs Y of main tracks vtx main
	96	TH1F *fPtPrKink;//! pt of Primary PDG kaons inside the selected ones by the kink topology
	97	TH1F *fgenPtEtRP;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
	98	TH1F *fgenPtEtRN;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
	99	TH1F *fkinkKaonP;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
	100	TH1F *fkinkKaonN;//!MC Pt spectrum of kaons decaying to muon+neutrino and pi +pi, inside eta,Pt,Rad cuts
	101	TH1F *frapidESDK;//!MC ESD K rapidity distr.
	102	TH1F *frapidKMC;//!MC MC K rapidity dis
	103	TH1F *fPtKPlMC;//!MC MC K rapidity dis
	104	TH1F *fPtKMnMC;//!MC MC K rapidity dis
	105	TH1F *fHistPtKaoP;//!MC MC K rapidity dis
	106	TH1F *fHistPtKaoN;//!MC MC K rapidity dis
	107	TH1F *fHiPtKPDGP;//!MC MC K rapidity dis
	108	TH1F *fHiPtKPDGN;//!MC MC K rapidity dis
	109	TH1F *fKinKBGP;//!MC MC K rapidity dis
	110	TH1F *fKinKBGN;//!MC MC K rapidity dis
	111	TH1F *fQtKMu;//!MC MC K Qt K to mu
	112	TH1F *fQtKPi;//!MC MC K Qt K to mu
	113	TH1F *fQtKEl;//!MC MC K Qt K to mu
	114	TH1F *fFakepipi;//!MC Fake pipi
115	TH1F *fFakeKPi;//!MC Fake Kpi
116	TH1F *fDCAkink;//!MC dcs kink
117	TH1F *fDCAkinkBG;//!MC dcs kink
118	TH2F *fPosiKink;//!MC position kink
119	TH2F *fPosiKinkK;//!MC position kink
120	TH2F *fPosiKinKXZ;//!MC position kink
121	TH2F *fPosiKinKYZ;//!MC position kink
122	TH2F *fPosiKinKBgZY;//!MC position kink
123	TH2F *fcode2;//!PDG code(mother) vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake)
124	TH2F *fcode4;//!PDG code(mother) vrs PDG dcode(daughter) of kinks with Qt <0.12 (fake)
125	TH2F *fZkinkZDau;//! z-position of kink z position of daughter bg )
126	TH1F *fQtKMuMC;//!MC MC K Qt K to mu
127	TH1F *fQtKElMC;//!MC MC K Qt K to mu
128	TH1F *fQtKPiMC;//!MC MC K Qt K to mu
129	TH1F *fQtK3PiP;//!EDS K Qt K to 3Pi
130	TH1F *fQtK3PiM;//!EDS K Qt K to 3Pi
131	TH2F *fmaxAngMomKmu; //!Decay angle vrs Mother Mom for pdg kaons, inside the selected sample
132	TH2F *fPosiKinKBgZX;//!MC position kink
133	TH2F *fPosiKinKBgXY;//!MC position kink
134	TH1F *fMinvPi;//!MC R life time
135	TH1F *fMinvKa;//!MC R life time
136	TH1F *fMinvPr;//!MC R life time
137	TH2F *fTPCSgnlPtpc;//Kink mother moment vrs TPC signal
138	TH2F *fTPCMomNSgnl;//kink mother TPC momentum vrs nsigmas of dEdx
139	TH2F *fMothKinkMomSgnl;//kink mother TPC momentum vrs nsigmas of dEdx
140	TH1F *fNSigmTPC;//kink mother TPC momentum vrs nsigmas of dEdx
141	TH2F *fTPCSgnlKinkDau;//Kink mother moment vrs TPC signal
142	TH2F *fcodeDau1;//!PDG code(mother) vrs PDG dcode(daughter) of kinks daughters
143	TH2F *fcodeDau2;//!PDG code(mother) vrs PDG dcode(daughter) of kinks daughters , Bg study
144	TH2F *fMothKinkMomSgnlD;//kink mother TPC momentum vrs nsigmas of dEdx
145	TH1F *fInvMassMuNuAll;//kinks, Inv Mass all kinks MuNu
146	TH2F *fInvMassMuNuPt;//kinks,Invariant Mass MuNu vs Pt
89db1464	147	TH2F *fInvMassMuNuPtAll;//kinks,Invariant Mass MuNu vs Pt 22/8/2013
a0271296	148	TH1F *fRatioCrossedRows; //ratio crossed rows
	149	TH1F *fRatioCrossedRowsKink; //ratio crossed rows for kinks
	150	TH2F *fRadiusPt;//kinks, Radius vs Pt
	151	TH2F *fRadiusPtcln;//kinks, Radius vs Pt for clean kaons
5220a87d	152	TH2F *fRadiusPtKaon;//kinks, Radius vs Pt for clean kaons
	153	TH2F *fRadiusPtPion;//kinks, Radius vs Pt for clean kaons
	154	TH2F *fRadiusPtFake;//kinks, Radius vs Pt for clean kaons
a0271296	155	TH1F *fPtCut1; //K Pt spectrum of all kinks from track bank, K0 bins
	156	TH1F *fPtCut2; //K Pt spectrum of all kinks from track bank, K0 bins
	157	TH1F *fPtCut3; //K Pt spectrum of all kinks from track bank, K0 bins
	158	TH2F *fAngMomKKinks;//kinks, Angle vs Momentum for K-kinks
	159	TH1F *flengthMCK;//!MC R life time
	160	TH1F *flifetiMCK;//!MC R life time
	161	TH1F *flifetim2;//!MC R life time
	162	TH1F *fLHelESDK;//!MC R life time
	163	TH1F *flifeInt;//!MC R life time
	164	TH1F *flifeYuri;//!MC R life time
	165	TH1F *flenYuri;//!MC R life time
	166	TH1F *flenTrRef;//!MC R life time
	167	TH1F *flifeSmall;//!MC R life time
	168	TH1F *flifetime;//!MC R life time
	169	TH1F *flifTiESDK;//!MC R life time
	170	TH1F *flifeKink;//!MC R life time
	171	TH1F *flenHelx;//!MC R life time
	172	TH3F *fradPtRapMC;//!MC R life time
	173	TH3F *fradPtRapDC;//!MC R life time
	174	TH3F *fradPtRapESD;//!MC R life time
	175	TH2F *fRadNclcln;//!MC R life time
	176
	177
	178
	179	TF1 *f1;
	180	TF1 *f2;
	181	TList *fListOfHistos; //! list of histos
	182
	183	Int_t fLowMulcut; //
	184	Int_t fUpMulcut;
	185	Int_t fKinkRadUp;
	186	Int_t fKinkRadLow;
5220a87d	187	Int_t fLowCluster;
	188	Float_t fLowQt;
	189
a0271296	190	AliESDtrackCuts* fCutsMul;
	191
	192	AliESDtrackCuts* fMaxDCAtoVtxCut;
	193	AliPIDResponse *fPIDResponse; //! PID response object
10eaad41	194
	195	AliAnalysisKinkESDMC(const AliAnalysisKinkESDMC&); // not implemented
	196	AliAnalysisKinkESDMC& operator=(const AliAnalysisKinkESDMC&); // not implemented
	197
	198	ClassDef(AliAnalysisKinkESDMC, 1); // example of analysis
	199	};
	200
	201	#endif