1 #ifndef ALIANAGENERATORKINE_H
2 #define ALIANAGENERATORKINE_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 //___________________________________________________________________________
7 // Do photon/pi0 analysis for isolation and correlation
8 // at the generator level. Only for kine stack (ESDs)
11 //-- Author: Gustavo Conesa (LPSC-CNRS-Grenoble)
17 class TLorentzVector ;
20 #include "AliAnaCaloTrackCorrBaseClass.h"
22 class AliAnaGeneratorKine : public AliAnaCaloTrackCorrBaseClass {
26 AliAnaGeneratorKine() ; // default ctor
27 virtual ~AliAnaGeneratorKine() { ; } //virtual dtor
29 Bool_t CorrelateWithPartonOrJet(TLorentzVector trigger,
36 TList * GetCreateOutputObjects() ;
38 void GetPartonsAndJets() ;
40 void GetXE(TLorentzVector trigger,
47 void InitParameters() ;
49 void IsLeadingAndIsolated(TLorentzVector trigger,
55 void MakeAnalysisFillAOD() { ; }
57 void MakeAnalysisFillHistograms() ;
61 AliStack * fStack; //! access stack
63 TParticle * fParton2; //! Initial state Parton
64 TParticle * fParton3; //! Initial state Parton
66 TParticle * fParton6; //! Final state Parton
67 TParticle * fParton7; //! Final state Parton
69 TLorentzVector fJet6; //! Pythia jet close to parton in position 6
70 TLorentzVector fJet7; //! Pythia jet close to parton in position 7
72 Float_t fPtHard; //! Generated pT hard
74 TH1F * fhPtHard; //! pt of parton
75 TH1F * fhPtParton; //! pt of parton
76 TH1F * fhPtJet; //! pt of jet
78 TH2F * fhPtPartonPtHard; //! pt of parton divided to pt hard, trigger is photon
79 TH2F * fhPtJetPtHard; //! pt of jet divided to pt hard, trigger is photon
80 TH2F * fhPtJetPtParton; //! pt of parton divided to pt parton, trigger is photon
82 TH1F * fhPtPhoton; //! Input photon
83 TH1F * fhPtPi0; //! Input pi0
85 // Histograms arrays for 4 isolation options and 2 options on leading or not leading particle
87 TH1F * fhPtPhotonLeading[4]; //! Leading photon
88 TH1F * fhPtPi0Leading[4]; //! Leading pi0
90 TH1F * fhPtPhotonLeadingIsolated[4]; //! Leading photon, isolated
91 TH1F * fhPtPi0LeadingIsolated[4]; //! Leading pi0, isolated
93 TH2F * fhPtPartonTypeNearPhoton[2][4]; //! Leading photon, particle pt versus originating parton type
94 TH2F * fhPtPartonTypeNearPi0[2][4]; //! Leading pi0, particle pt versus originating parton type
95 TH2F * fhPtPartonTypeNearPhotonIsolated[2][4]; //! Leading photon, particle pt versus originating parton type
96 TH2F * fhPtPartonTypeNearPi0Isolated[2][4]; //! Leading pi0, particle pt versus originating parton type
98 TH2F * fhPtPartonTypeAwayPhoton[2][4]; //! Leading photon, particle pt versus away side parton type
99 TH2F * fhPtPartonTypeAwayPi0[2][4]; //! Leading pi0, particle pt versus away side parton type
100 TH2F * fhPtPartonTypeAwayPhotonIsolated[2][4]; //! Leading photon, isolated, particle pt versus away side parton type
101 TH2F * fhPtPartonTypeAwayPi0Isolated[2][4]; //! Leading pi0, isolated, particle pt versus away side parton type
103 TH2F * fhZHardPhoton[2][4]; //! Leading photon, zHard
104 TH2F * fhZHardPi0[2][4]; //! Leading pi0, zHard
105 TH2F * fhZHardPhotonIsolated[2][4]; //! Leading photon, isolated, zHard
106 TH2F * fhZHardPi0Isolated[2][4]; //! Leading pi0, isolated, zHard
108 TH2F * fhZPartonPhoton[2][4]; //! Leading photon, zHard
109 TH2F * fhZPartonPi0[2][4]; //! Leading pi0, zHard
110 TH2F * fhZPartonPhotonIsolated[2][4]; //! Leading photon, isolated, zHard
111 TH2F * fhZPartonPi0Isolated[2][4]; //! Leading pi0, isolated, zHard
113 TH2F * fhZJetPhoton[2][4]; //! Leading photon, zHard
114 TH2F * fhZJetPi0[2][4]; //! Leading pi0, zHard
115 TH2F * fhZJetPhotonIsolated[2][4]; //! Leading photon, isolated, zHard
116 TH2F * fhZJetPi0Isolated[2][4]; //! Leading pi0, isolated, zHard
118 TH2F * fhXEPhoton[2][4]; //! Leading photon, xE away side
119 TH2F * fhXEPi0[2][4]; //! Leading pi0, xE away side
120 TH2F * fhXEPhotonIsolated[2][4]; //! Leading photon, xE away side
121 TH2F * fhXEPi0Isolated[2][4]; //! Leading pi0, isolated, xE away side
123 TH2F * fhXEUEPhoton[2][4]; //! Leading photon, xE away side
124 TH2F * fhXEUEPi0[2][4]; //! Leading pi0, xE away side
125 TH2F * fhXEUEPhotonIsolated[2][4]; //! Leading photon, xE away side
126 TH2F * fhXEUEPi0Isolated[2][4]; //! Leading pi0, isolated, xE away side
128 AliAnaGeneratorKine (const AliAnaGeneratorKine & gk) ; // cpy ctor
129 AliAnaGeneratorKine & operator = (const AliAnaGeneratorKine & gk) ; // cpy assignment
131 ClassDef(AliAnaGeneratorKine,2)
136 #endif //ALIANAGENERATORKINE_H