1 #ifndef ALIANALYSISKINKESDat_H
2 #define ALIANALYSISKINKESDat_H
4 /* See cxx source for full Copyright notice */
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 //-----------------------------------------------------------------
25 class AliESDtrackCuts;
27 #include "AliAnalysisTaskSE.h"
29 class AliAnalysisKinkESDat : public AliAnalysisTaskSE {
31 AliAnalysisKinkESDat(const char *name = "AliAnalysisKinkESDat");
32 virtual ~AliAnalysisKinkESDat() {}
34 virtual void UserCreateOutputObjects();
35 virtual void UserExec(Option_t *option);
36 virtual void Terminate(Option_t *);
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;}
43 TH1F *fHistPtESD; //Pt spectrum of all ESD inside eta, Pt cuts
44 TH1F *fHistPt; //Pt spectrum of all ESD tracks
45 TH1F *fHistQtAll; //Qt spectrum of all kinks
46 TH1F *fHistQt1; //Qt spectrum of Kaon selected sample
47 TH1F *fHistQt2; //Qt spectrum in Qt region of kaons
48 TH1F *fHistPtKaon; //Pt Kaon spectrum of clean sample
49 TH1F *fHistPtKPDG; //Pt Kaon spectrum , confirmed by PDG,inside kaon Qt region
50 TH1F *fHistEta; //Eta spectrum of all kinks
51 TH1F *fHistEtaK; //Eta spectrum of kaons selected by kink topology
52 TH1F *fptKMC; //Pt Kaon spectrum MC, inside eta and pt cuts
53 TH1F *fMultiplMC; //charge multipl MC
54 TH1F *fESDMult; //ESD charged mult
55 TH1F *fgenpt; //Pt Kaon-Kink->mu spectrum , MC, inside eta, Pt, radius cuts
56 TH1F *frad; //radius of kinks, MC , inside the eta nad Pt cuts
57 TH1F *fKinkKaon; //Pt of PDG Kaons inside the selcted ones by the KInk topology
58 TH1F *fKinKRbn; //Pt of PDG Kaons inside the selcted ones by the KInk topology
59 TH1F *fKinkKaonBg; //Pt of the BG inside the kink-Kaon identified spectrum
60 TH1F *fM1kaon; //inv mass of kink-tracks taken as kaons decaying to mu + neutrino
61 TH1F *fPtKink; //Pt spectrum of all kinks from track bank
62 TH1F *fptKink; //Pt spectrum of all kinks from kink bank
63 TH2F *fAngMomK; // Decay angle vrs Mother Mom for pdg kaons
64 TH2F *fAngMomPi; // Decay angle vrs Mother Mom for pdg pions
65 TH2F *fAngMomKC; //Decay angle vrs Mother Mom for pdg kaons, inside the selected sample
66 TH1F *fMultESDK; //ESD charged mult
67 TH1F *fMultMCK; //MC K charged mult
68 TH2F *fSignPtNcl;//signPt vrs number of clusters in TPC for kaons from kink sele sample
69 TH2F *fSignPtEta;//signPt vrs Eta in TPC for kaons from kink sele sample
70 TH2F *fEtaNcl;//Eta vrs Nclu in TPC for kaons from kink sele sample
71 TH1F *fSignPt;//signPt in TPC for kaons from kink sele sample
72 TH2F *fChi2NclTPC;//chi2 vrs TPC Nclusters for kaons from kink sele sample
73 TH1F *fRatChi2Ncl;// Ratio chi2/ Ncl TPC for kaons from kink sele sample
74 TH2F *fRadiusNcl;//kink Radius Ncl TPC for kaons from kink sele sample
75 TH2F *fTPCSgnlP;//kink Radius Ncl TPC for kaons from kink sele sample
76 TH2F *fTPCSgnlPa;//kink Radius Ncl TPC for kink sele sample
77 TH1D *fRpr; // Radius of VTX at Y , X plane
78 TH1D *fZpr; //Z distrio of main vertex
79 TH1D *fdcatoVxXY; //dca to Vertex XY distrio
80 TH1D *fnSigmToVx; //nSigma to Vertex distrio of main vertex
81 TH2F *fKinkMothDau; //Mother vrs Daughter
82 TH2F *fZvXv; //two dime of Z vrs X of vtx main
83 TH2F *fZvYv; // two dime of Z vrs Y of vtx main
84 TH2F *fXvYv; // two dime of X vrs Y of main tracks vtx main
85 TH1F *fHistPtKaoP; //Pt Kaon spectrum of clean sample pos
86 TH1F *fHistPtKaoN; //Pt Kaon spectrum of clean sample neg
87 TH1F *frapiKESD;// rapidi K
88 TH1F *flifetime;//radius of kinks, MC , inside the eta nad Pt cuts
89 TH1F *fradLK;//Length of kinks, MC , inside the eta nad Pt cuts
90 TH3F *fradPtRpDt;//radius of kinks, MC , inside the eta nad Pt cuts
91 TH1F *fInvMuNuAll;//radius of kinks, MC , inside the eta nad Pt cuts
93 TH1F *fDCAkink;//!MC dcs kink
94 TH2F *fPosiKink;//!MC position kink
95 TH2F *fPosiKinkK;//!MC position kink
96 TH2F *fPosiKinKXZ;//!MC position kink
97 TH2F *fPosiKinKYZ;//!MC position kink
98 TH2F *fPosiKinKBg;//!MC position kink
99 TH2F *fQtMothP;//!qt vrs p mother
100 TH2F *fTPCSgnlPtpc;//Kink mother moment vrs TPC signal
101 TH2F *fTPCMomNSgnl;//kink mother TPC momentum vrs nsigmas of dEdx
102 TH2F *fMothKinkMomSgnl;//kink mother TPC momentum vrs nsigmas of dEdx
103 TH1F *fNSigmTPC;//kink mother TPC momentum vrs nsigmas of dEdx
104 TH2F *fTPCSgnlKinkDau;//Kink mother moment vrs TPC signal
105 TH1F *fPtKinkPos; //Pos K Pt spectrum of all kinks from track bank, K0 bins
106 TH1F *fPtKinkNeg; //Neg K Pt spectrum of all kinks from track bank, K0 bins
107 TH2F *fRadNclCln;//kink Radius Ncl TPC for kaons from kink clean sample
111 TList *fListOfHistos; // list of histos
112 // Marek Multiplicity
115 AliESDtrackCuts* fCutsMul;
116 AliPIDResponse *fPIDResponse; //! PID response object
118 AliAnalysisKinkESDat(const AliAnalysisKinkESDat&); // not implemented
119 AliAnalysisKinkESDat& operator=(const AliAnalysisKinkESDat&); // not implemented
121 ClassDef(AliAnalysisKinkESDat, 1); // example of analysis