B2 analysis code

[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisEmEtMonteCarlo.h

#ifndef ALIANALYSISEMETMONTECARLO_H
#define ALIANALYSISEMETMONTECARLO_H
//_________________________________________________________________________
//  Utility Class for transverse energy studies
//  Base class for MC analysis
//  - MC output
//
//*-- Author: Marcelo G. Munhoz (USP)
//_________________________________________________________________________

#include "AliAnalysisEtMonteCarlo.h"
class TParticle;
class TParticlePDG;
class AliMCParticle;
class AliESDtrack;
class AliEMCALTrack;
class TVector3;
class AliEMCALGeometry;
class AliExternalTrackParam;
class AliStack;

class AliAnalysisEmEtMonteCarlo : public AliAnalysisEtMonteCarlo
{

public:
   
  AliAnalysisEmEtMonteCarlo();
  virtual ~AliAnalysisEmEtMonteCarlo();

    virtual Int_t AnalyseEvent(AliVEvent* event);
	virtual Int_t AnalyseEvent(AliVEvent* event, AliVEvent* event2);

    virtual void Init();
    virtual void ResetEventValues();
    virtual void CreateHistograms();
    virtual void FillOutputList(TList* list);

protected:

	virtual Bool_t IsPrimary(AliStack *stack, Int_t part, TParticlePDG *pdg, Int_t partMom, TParticlePDG *pdgMom);
	virtual Bool_t IsMotherPrimaryGamma(AliStack *stack, Int_t iPartMom, TParticlePDG *pdgMom);
	virtual Bool_t IsMotherPrimaryElectron(AliStack *stack, Int_t iPartMom, TParticlePDG *pdgMom);
	virtual Bool_t IsGammaConversion(AliStack *stack, TParticle *part, TParticlePDG *pdg);
	virtual Bool_t IsInAcceptance(TParticle *part=0, TParticlePDG *pdg=0, AliExternalTrackParam* extParam=0);
	virtual Bool_t IsInAcceptance(AliMCParticle *part=0);
	
	virtual Bool_t TrackHitsCalo(AliExternalTrackParam *extParam);
	
	virtual Bool_t GetTrackProjection(AliExternalTrackParam *trackParam, TVector3 &trackPos); // project to a radius
	virtual Bool_t GetTrackProjection(AliEMCALTrack* emcTrack, TVector3 &trackPos, TVector3 clusPos); // project to a point

	AliExternalTrackParam* CreateExternalTrackParam(TParticle *part);
	
	virtual Double_t CalcET(TParticle *part, TParticlePDG *pdg);
	virtual Double_t CalcETDep(Double_t caloE, TParticle *part, TParticlePDG *pdg);
	
protected:

	Double_t fResCut; // b(fm), for Hijing; 0 otherwise
	Double_t fPrimtotET, fPrimAcctotET, fPrimRectotET, fPrimRectotETDep;//Primary particles ET

	Double_t fElectrontotET, fElectronAcctotET, fElectronRectotET;//Electron ET
	Double_t fConvElectrontotET, fConvElectronAcctotET, fConvElectronRectotET, fScatElectrontotET, fScatElectronAcctotET, fScatElectronRectotET;//Secondary electrons ET
	Double_t fTotElectrontotET, fTotElectronAcctotET, fTotElectronRectotET;//Total Electrons ET

	Double_t fGammatotET, fGammaAcctotET, fGammaRectotET;//Gamma event ET
	Double_t fAnnihGammatotET, fAnnihGammaAcctotET, fAnnihGammaRectotET, fScatGammatotET, fScatGammaAcctotET, fScatGammaRectotET;//Secondary gamma ET
	Double_t fTotGammatotET, fTotGammaAcctotET, fTotGammaRectotET;//Total gamma ET
	Double_t fConvGammatotET, fNonConvGammatotET, fConvGammaAcctotET, fNonConvGammaAcctotET, fNPPPi0GammatotET, fNPPPi0GammaRectotET;//Conversion and pi0 gamma ET

	Double_t fTotEMtotET, fTotEMAcctotET, fTotEMRectotET;//Total EM ET

	Double_t fNPPElectrontotET, fNPPElectronRectotET, fNPPGammatotET, fNPPGammaRectotET;//Non-physical primary electron ET
	Double_t fTotNPPEMtotET, fTotNPPEMRectotET;//Total Non-physical primary electron ET

	Double_t fMuontotET, fPiontotET, fKaontotET, fProtontotET;//Charged particles ET
	Double_t fMuonAcctotET, fPionAcctotET, fKaonAcctotET, fProtonAcctotET;//Charged particles acceptance ET
	Double_t fMuonRectotET, fMuonRectotETDep, fPionRectotET, fPionRectotETDep, fKaonRectotET, fKaonRectotETDep, fProtonRectotET, fProtonRectotETDep;//Charged particles reconstructed ET
	Double_t fMuonMatchtotET, fMuonMatchtotETDep, fPionMatchtotET, fPionMatchtotETDep, fKaonMatchtotET, fKaonMatchtotETDep, fProtonMatchtotET, fProtonMatchtotETDep;//Charged particles track matched ET
	Double_t fTotChargedtotET, fTotChargedAcctotET, fTotChargedRectotET, fTotChargedRectotETDep, fTotChargedMatchtotET, fTotChargedMatchtotETDep;//Total charged particles ET

	Double_t fNeutrontotET, fNeutronAcctotET, fNeutronRectotET, fNeutronRectotETDep;//Neutrons ET
	Double_t fK0totET, fK0RectotET, fK0RectotETDep, fLambdatotET, fLambdaRectotET, fLambdaRectotETDep;//K0 and Lambda ET
	Double_t fTotNeutraltotET, fTotNeutralRectotET, fTotNeutralRectotETDep;//Neutral particles ET

	Double_t fTotaltotET, fTotalAcctotET, fTotalRectotET, fTotalRectotETDep;//Total ET
	
	AliEMCALGeometry *fGeoUt;//EMCal geometry object

	// *******************
	// primaries ET
	// *******************
	TH2F *fHistPrimEtaEET;//total ET - Eta vs E 
	TH2F *fHistPrimEtaPtET;//total ET - Eta vs pt 
	TH2F *fHistPrimEtaET;//total ET - Eta
	TH1F *fHistPrimtotET;//total ET distribution
	
	TH2F *fHistPrimAccEtaEET;//acceptance ET - Eta vs E 
	TH2F *fHistPrimAccEtaPtET;//acceptance ET - Eta vs pt 
	TH2F *fHistPrimAccEtaET;//acceptance ET - Eta
	TH1F *fHistPrimAcctotET;//acceptance ET distribution
	
	TH2F *fHistPrimRecEtaEET;//reconstructed ET - Eta vs E 
	TH2F *fHistPrimRecEtaPtET;//reconstructed ET - Eta vs pt 
	TH2F *fHistPrimRecEtaET;//reconstructed ET - Eta
	TH1F *fHistPrimRectotET;//reconstructed ET distribution

	TH2F *fHistPrimRecEtaEDepETDep;//deposited ET - Eta vs E deposited
	TH2F *fHistPrimRecEtaPtETDep;//deposited ET - Eta vs pt 
	TH2F *fHistPrimRecEtaETDep;//deposited ET - Eta 
	TH1F *fHistPrimRectotETDep;//deposited ET distribution
	
	// *******************
	// electron ET
	// *******************
	TH2F *fHistElectronEtaEET;// ET - Eta vs E 
	TH2F *fHistElectronEtaPtET;//ET - Eta vs pt
	TH2F *fHistElectronEtaET;// ET - Eta 
	TH2F *fHistElectronEtaE;// multiplicity - Eta vs E
	TH2F *fHistElectronEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistElectrontotET;// total ET distribution 

	TH2F *fHistConvElectronEtaEET;//ET - Eta vs E   
	TH2F *fHistConvElectronEtaPtET;//ET - Eta vs pt 
	TH2F *fHistConvElectronEtaET;//  ET - Eta
	TH2F *fHistConvElectronEtaE;//  multiplicity - Eta vs E
	TH2F *fHistConvElectronEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistConvElectrontotET;// total ET distribution 

	TH2F *fHistScatElectronEtaEET;//ET - Eta vs E  
	TH2F *fHistScatElectronEtaPtET;//ET - Eta vs pt 
	TH2F *fHistScatElectronEtaET;//  ET - Eta
	TH2F *fHistScatElectronEtaE;//  multiplicity - Eta vs E
	TH2F *fHistScatElectronEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistScatElectrontotET;//  total ET distribution
	
	// *******************
	// total electron ET
	// *******************
	TH1F *fHistTotElectrontotET;//total ET distribution
	
	// *******************
	// gamma ET
	// *******************
	TH2F *fHistGammaEtaEET;//ET - Eta vs E  
	TH2F *fHistGammaEtaPtET;//ET - Eta vs pt  
	TH2F *fHistGammaEtaET;//  ET - Eta
	TH2F *fHistGammaEtaE;//  multiplicity - Eta vs E
	TH2F *fHistGammaEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistGammatotET;//  total ET distribution
	
	TH2F *fHistAnnihGammaEtaEET;//ET - Eta vs E  
	TH2F *fHistAnnihGammaEtaPtET;//ET - Eta vs pt
	TH2F *fHistAnnihGammaEtaET;//  ET - Eta
	TH2F *fHistAnnihGammaEtaE;//  multiplicity - Eta vs E
	TH2F *fHistAnnihGammaEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistAnnihGammatotET;//  total ET distribution

	TH2F *fHistScatGammaEtaEET;//ET - Eta vs E   
	TH2F *fHistScatGammaEtaPtET;//ET - Eta vs pt  
	TH2F *fHistScatGammaEtaET;//  ET - Eta
	TH2F *fHistScatGammaEtaE;//  multiplicity - Eta vs E
	TH2F *fHistScatGammaEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistScatGammatotET;//  total ET distribution

	TH2F *fHistConvGammaEtaEET;//ET - Eta vs E  
	TH2F *fHistConvGammaEtaPtET;//ET - Eta vs pt
	TH2F *fHistConvGammaEtaET;//  ET - Eta
	TH2F *fHistConvGammaEtaE;//  multiplicity - Eta vs E
	TH2F *fHistConvGammaEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistConvGammatotET;//  total ET distribution
	
	TH2F *fHistNonConvGammaEtaEET;//ET - Eta vs E  
	TH2F *fHistNonConvGammaEtaPtET;//ET - Eta vs pt 
	TH2F *fHistNonConvGammaEtaET;//  ET - Eta
	TH2F *fHistNonConvGammaEtaE;//  multiplicity - Eta vs E
	TH2F *fHistNonConvGammaEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistNonConvGammatotET;//  total ET distribution
	
	
	
	// *******************
	// total gamma ET
	// *******************
	TH1F *fHistTotGammatotET;//total ET distribution

	// *******************
	// total electromagnetic ET
	// *******************
	TH1F *fHistTotEMtotET;//total ET distribution

	// non-primary electromagnetic ET
	TH2F *fHistNPPElectronEtaEET;//ET - Eta vs E 
	TH2F *fHistNPPElectronEtaPtET;//ET - Eta vs pt
	TH2F *fHistNPPElectronEtaET;// ET - Eta
	TH2F *fHistNPPElectronEtaE;// multiplicity - Eta vs E
	TH2F *fHistNPPElectronEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistNPPElectrontotET;// total ET distribution

	TH2F *fHistNPPGammaEtaEET;//ET - Eta vs E  
	TH2F *fHistNPPGammaEtaPtET;//ET - Eta vs pt 
	TH2F *fHistNPPGammaEtaET;// ET - Eta
	TH2F *fHistNPPGammaEtaE;// multiplicity - Eta vs E
	TH2F *fHistNPPGammaEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistNPPGammatotET;// total ET distribution

	TH1F *fHistTotNPPEMtotET;//total ET distribution

	TH2F *fHistNPPPi0GammaEtaEET;//ET - Eta vs E  
	TH2F *fHistNPPPi0GammaEtaPtET;//ET - Eta vs pt
	TH2F *fHistNPPPi0GammaEtaET;// ET - Eta
	TH2F *fHistNPPPi0GammaEtaE;// multiplicity - Eta vs E
	TH2F *fHistNPPPi0GammaEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistNPPPi0GammatotET;// total ET distribution
		
	// *******************
	// electron ET inside EMCal acceptance
	// *******************
	TH2F *fHistElectronAccEtaEET;//ET - Eta vs E  
	TH2F *fHistElectronAccEtaPtET;//ET - Eta vs pt 
	TH2F *fHistElectronAccEtaET;// ET - Eta
	TH2F *fHistElectronAccEtaE;// multiplicity - Eta vs E
	TH2F *fHistElectronAccEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistElectronAcctotET;// total ET distribution
	
	TH2F *fHistConvElectronAccEtaEET;//ET - Eta vs E   
	TH2F *fHistConvElectronAccEtaPtET;//ET - Eta vs pt  
	TH2F *fHistConvElectronAccEtaET;//  ET - Eta
	TH2F *fHistConvElectronAccEtaE;//  multiplicity - Eta vs E
	TH2F *fHistConvElectronAccEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistConvElectronAcctotET;//  total ET distribution
	
	TH2F *fHistScatElectronAccEtaEET;//ET - Eta vs E   
	TH2F *fHistScatElectronAccEtaPtET;//ET - Eta vs pt  
	TH2F *fHistScatElectronAccEtaET;//  ET - Eta
	TH2F *fHistScatElectronAccEtaE;//  multiplicity - Eta vs E
	TH2F *fHistScatElectronAccEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistScatElectronAcctotET;//  total ET distribution
	
	// *******************
	// total electron ET inside EMCal acceptance
	// *******************
	TH1F *fHistTotElectronAcctotET;//total ET distribution

	// *******************
	// gamma ET inside EMCal acceptance
	// *******************
	TH2F *fHistGammaAccEtaEET;//ET - Eta vs E   
	TH2F *fHistGammaAccEtaPtET;//ET - Eta vs pt  
	TH2F *fHistGammaAccEtaET;//  ET - Eta
	TH2F *fHistGammaAccEtaE;//  multiplicity - Eta vs E
	TH2F *fHistGammaAccEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistGammaAcctotET;//  total ET distribution
	
	TH2F *fHistAnnihGammaAccEtaEET;//ET - Eta vs E 
	TH2F *fHistAnnihGammaAccEtaPtET;//ET - Eta vs pt  
	TH2F *fHistAnnihGammaAccEtaET;//  ET - Eta
	TH2F *fHistAnnihGammaAccEtaE;//  multiplicity - Eta vs E
	TH2F *fHistAnnihGammaAccEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistAnnihGammaAcctotET;//  total ET distribution
	
	TH2F *fHistScatGammaAccEtaEET;//ET - Eta vs E  
	TH2F *fHistScatGammaAccEtaPtET;//ET - Eta vs pt  
	TH2F *fHistScatGammaAccEtaET;//  ET - Eta
	TH2F *fHistScatGammaAccEtaE;//  multiplicity - Eta vs E
	TH2F *fHistScatGammaAccEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistScatGammaAcctotET;//  total ET distribution
	
	TH2F *fHistConvGammaAccEtaEET;//ET - Eta vs E   
	TH2F *fHistConvGammaAccEtaPtET;//ET - Eta vs pt 
	TH2F *fHistConvGammaAccEtaET;//  ET - Eta
	TH2F *fHistConvGammaAccEtaE;//  multiplicity - Eta vs E
	TH2F *fHistConvGammaAccEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistConvGammaAcctotET;//  total ET distribution
	
	TH2F *fHistNonConvGammaAccEtaEET;//ET - Eta vs E   
	TH2F *fHistNonConvGammaAccEtaPtET;//ET - Eta vs pt 
	TH2F *fHistNonConvGammaAccEtaET;//  ET - Eta
	TH2F *fHistNonConvGammaAccEtaE;//  multiplicity - Eta vs E
	TH2F *fHistNonConvGammaAccEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistNonConvGammaAcctotET;//  total ET distribution
	
	
	// *******************
	// total gamma ET inside EMCal acceptance
	// *******************
	TH1F *fHistTotGammaAcctotET;//total ET distribution

	// *******************
	// total electromagnetic ET inside EMCal acceptance
	// *******************
	TH1F *fHistTotEMAcctotET;//total ET distribution

	// non-primary electromagnetic ET
	TH2F *fHistNPPElectronAccEtaEET;//ET - Eta vs E 
	TH2F *fHistNPPElectronAccEtaPtET;//ET - Eta vs pt 
	TH2F *fHistNPPElectronAccEtaE;// multiplicity - Eta vs E
	TH2F *fHistNPPElectronAccEtaPt;// multiplicity - Eta vs pt
	
	TH2F *fHistNPPGammaAccEtaEET;//ET - Eta vs E  
	TH2F *fHistNPPGammaAccEtaPtET;//ET - Eta vs pt 
	TH2F *fHistNPPGammaAccEtaE;// multiplicity - Eta vs E
	TH2F *fHistNPPGammaAccEtaPt;// 	multiplicity - Eta vs pt
	
	// *******************
	// electron ET reconstructed in EMCal
	// *******************
	TH2F *fHistElectronRecEtaEET;//ET - Eta vs E  
	TH2F *fHistElectronRecEtaPtET;//ET - Eta vs pt
	TH2F *fHistElectronRecEtaET;// ET - Eta
	TH2F *fHistElectronRecEtaE;// multiplicity - Eta vs E
	TH2F *fHistElectronRecEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistElectronRectotET;// total ET distribution
	
	TH2F *fHistConvElectronRecEtaEET;//ET - Eta vs E   
	TH2F *fHistConvElectronRecEtaPtET;//ET - Eta vs pt 
	TH2F *fHistConvElectronRecEtaET;//  ET - Eta
	TH2F *fHistConvElectronRecEtaE;//  multiplicity - Eta vs E
	TH2F *fHistConvElectronRecEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistConvElectronRectotET;//  total ET distribution
	
	TH2F *fHistScatElectronRecEtaEET;//ET - Eta vs E   
	TH2F *fHistScatElectronRecEtaPtET;//ET - Eta vs pt  
	TH2F *fHistScatElectronRecEtaET;//  ET - Eta
	TH2F *fHistScatElectronRecEtaE;//  multiplicity - Eta vs E
	TH2F *fHistScatElectronRecEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistScatElectronRectotET;//  total ET distribution
	
	// *******************
	// total Electron ET reconstructed in EMCal
	// *******************
	TH1F *fHistTotElectronRectotET;//total ET distribution

	// *******************
	// gamma ET reconstructed in EMCal
	// *******************
	TH2F *fHistGammaRecEtaEET;//ET - Eta vs E   
	TH2F *fHistGammaRecEtaPtET;//ET - Eta vs pt  
	TH2F *fHistGammaRecEtaET;//  ET - Eta
	TH2F *fHistGammaRecEtaE;//  multiplicity - Eta vs E
	TH2F *fHistGammaRecEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistGammaRectotET;//  total ET distribution
	
	
	
	TH2F *fHistGammaRecResEET;// ET - track matching residual vs E_Added CAIO
	TH2F *fHistGammaRecResPtET;// ET - track matching residual vs pt_Added CAIO
	TH2F *fHistGammaRecResE;// multiplicity - track matching residual vs E_Added CAIO
	TH2F *fHistGammaRecResPt;// multiplicity - track matching residual vs pt_Added CAIO 
	TH2F *fHistGammaRecResEDepETDep;// ET deposited - track matching residual vs E deposited_Added CAIO
	TH2F *fHistGammaRecResPtETDep;// ET deposited - track matching residual vs pt_Added Caio

	
	TH2F *fHistAnnihGammaRecEtaEET;//ET - Eta vs E   
	TH2F *fHistAnnihGammaRecEtaPtET;//ET - Eta vs pt  
	TH2F *fHistAnnihGammaRecEtaET;//  ET - Eta
	TH2F *fHistAnnihGammaRecEtaE;//  multiplicity - Eta vs E
	TH2F *fHistAnnihGammaRecEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistAnnihGammaRectotET;//  total ET distribution
	
	TH2F *fHistScatGammaRecEtaEET;//ET - Eta vs E   
	TH2F *fHistScatGammaRecEtaPtET;//ET - Eta vs pt  
	TH2F *fHistScatGammaRecEtaET;//  ET - Eta
	TH2F *fHistScatGammaRecEtaE;//  multiplicity - Eta vs E
	TH2F *fHistScatGammaRecEtaPt;//  multiplicity - Eta vs pt
	TH1F *fHistScatGammaRectotET;//  total ET distribution
	
	
	// *******************
	// total gamma ET reconstructed in EMCal
	// *******************
	TH1F *fHistTotGammaRectotET;//total ET distribution

	// *******************
	// total EM ET reconstructed in EMCal
	// *******************
	TH1F *fHistTotEMRectotET;//total ET distribution

	// non-primary electromagnetic ET
	TH2F *fHistNPPElectronRecEtaEET;//ET - Eta vs E  
	TH2F *fHistNPPElectronRecEtaPtET;//ET - Eta vs pt 
	TH2F *fHistNPPElectronRecEtaET;// ET - Eta
	TH2F *fHistNPPElectronRecEtaE;// multiplicity - Eta vs E
	TH2F *fHistNPPElectronRecEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistNPPElectronRectotET;// total ET distribution
	
	TH2F *fHistNPPGammaRecEtaEET;//ET - Eta vs E  
	TH2F *fHistNPPGammaRecEtaPtET;//ET - Eta vs pt 
	TH2F *fHistNPPGammaRecEtaET;// ET - Eta
	TH2F *fHistNPPGammaRecEtaE;// multiplicity - Eta vs E
	TH2F *fHistNPPGammaRecEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistNPPGammaRectotET;// total ET distribution
	
	TH1F *fHistTotNPPEMRectotET;//total ET distribution

	TH2F *fHistNPPPi0GammaRecEtaEET;//ET - Eta vs E  
	TH2F *fHistNPPPi0GammaRecEtaPtET;//ET - Eta vs pt 
	TH2F *fHistNPPPi0GammaRecEtaET;// ET - Eta
	TH2F *fHistNPPPi0GammaRecEtaE;// multiplicity - Eta vs E
	TH2F *fHistNPPPi0GammaRecEtaPt;// multiplicity - Eta vs pt
	TH1F *fHistNPPPi0GammaRectotET;// total ET distribution
	
	// *******************
	// muon ET (+ and -)
	// *******************
	TH2F *fHistMuonEtaEET;//ET - Eta vs E  
	TH2F *fHistMuonAccEtaEET;//ET - Eta vs E 
	TH2F *fHistMuonRecEtaEET;//ET - Eta vs E 
	TH2F *fHistMuonMatchEtaEET;//ET - Eta vs E 

	TH2F *fHistMuonEtaPtET;// ET - Eta vs pt
	TH2F *fHistMuonAccEtaPtET;// ET - Eta vs pt
	TH2F *fHistMuonRecEtaPtET;// ET - Eta vs pt
	TH2F *fHistMuonMatchEtaPtET;// ET - Eta vs pt

	TH2F *fHistMuonEtaET;// ET - Eta
	TH2F *fHistMuonAccEtaET;// ET - Eta
	TH2F *fHistMuonRecEtaET;// ET - Eta
	TH2F *fHistMuonMatchEtaET;// ET - Eta
	
	TH2F *fHistMuonEtaE;// multiplicity - Eta vs E
	TH2F *fHistMuonAccEtaE;// multiplicity - Eta vs E
	TH2F *fHistMuonRecEtaE;// multiplicity - Eta vs E
	TH2F *fHistMuonMatchEtaE;// multiplicity - Eta vs E
	
	TH2F *fHistMuonEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistMuonAccEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistMuonRecEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistMuonMatchEtaPt;// multiplicity - Eta vs pt
	
	TH1F *fHistMuontotET;// total ET distribution
	TH1F *fHistMuonAcctotET;// total ET distribution
	TH1F *fHistMuonRectotET;// total ET distribution
	TH1F *fHistMuonMatchtotET;// total ET distribution
	
	TH1F *fHistMuonRectotETDep;//total deposited ET distribution
	TH1F *fHistMuonMatchtotETDep;// total deposited ET distribution
	
	TH2F *fHistMuonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
	TH2F *fHistMuonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited

	TH2F *fHistMuonRecEtaPtETDep;// ET deposited - Eta vs pt
	TH2F *fHistMuonMatchEtaPtETDep;// ET deposited - Eta vs pt
	
	TH2F *fHistMuonRecEtaETDep;// ET deposited - Eta
	TH2F *fHistMuonMatchEtaETDep;// ET deposited - Eta

	TH2F *fHistMuonRecResEET;// ET - track matching residual vs E
	TH2F *fHistMuonRecResPtET;// ET - track matching residual vs pt
	TH2F *fHistMuonRecResE;// multiplicity - track matching residual vs E
	TH2F *fHistMuonRecResPt;// multiplicity - track matching residual vs pt 
	TH2F *fHistMuonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
	TH2F *fHistMuonRecResPtETDep;// ET deposited - track matching residual vs pt
	
	// *******************
	// pion ET (+ and -)
	// *******************
	
	TH2F *fHistPionEtaEET;//ET - Eta vs E  
	TH2F *fHistPionAccEtaEET;//ET - Eta vs E 
	TH2F *fHistPionRecEtaEET;//ET - Eta vs E 
	TH2F *fHistPionMatchEtaEET;//ET - Eta vs E 
	
	TH2F *fHistPionEtaPtET;// ET - Eta vs pt
	TH2F *fHistPionAccEtaPtET;// ET - Eta vs pt
	TH2F *fHistPionRecEtaPtET;// ET - Eta vs pt
	TH2F *fHistPionMatchEtaPtET;// ET - Eta vs pt
	
	TH2F *fHistPionEtaET;// ET - Eta
	TH2F *fHistPionAccEtaET;// ET - Eta
	TH2F *fHistPionRecEtaET;// ET - Eta
	TH2F *fHistPionMatchEtaET;// ET - Eta
	
	TH2F *fHistPionEtaE;// multiplicity - Eta vs E
	TH2F *fHistPionAccEtaE;// multiplicity - Eta vs E
	TH2F *fHistPionRecEtaE;// multiplicity - Eta vs E
	TH2F *fHistPionMatchEtaE;// multiplicity - Eta vs E
	
	TH2F *fHistPionEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistPionAccEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistPionRecEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistPionMatchEtaPt;// multiplicity - Eta vs pt
	
	TH1F *fHistPiontotET;// total ET distribution
	TH1F *fHistPionAcctotET;// total ET distribution
	TH1F *fHistPionRectotET;// total ET distribution
	TH1F *fHistPionMatchtotET;// total ET distribution
	
	TH1F *fHistPionRectotETDep;// total deposited ET distribution
	TH1F *fHistPionMatchtotETDep;// total deposited ET distribution
	
	TH2F *fHistPionRecEtaEDepETDep;// ET deposited - Eta vs E deposited
	TH2F *fHistPionMatchEtaEDepETDep;// ET deposited - Eta vs E deposited

	TH2F *fHistPionRecEtaPtETDep;// ET deposited - Eta vs pt
	TH2F *fHistPionMatchEtaPtETDep;// ET deposited - Eta vs pt
	
	TH2F *fHistPionRecEtaETDep;// ET deposited - Eta
	TH2F *fHistPionMatchEtaETDep;// ET deposited - Eta
	
	TH2F *fHistPionRecResEET;// ET - track matching residual vs E
	TH2F *fHistPionRecResPtET;// ET - track matching residual vs pt
	TH2F *fHistPionRecResE;// multiplicity - track matching residual vs E
	TH2F *fHistPionRecResPt;// multiplicity - track matching residual vs pt
	TH2F *fHistPionRecResEDepETDep;// ET deposited - track matching residual vs E deposited
	TH2F *fHistPionRecResPtETDep;// ET deposited - track matching residual vs pt
	
	// *******************
	// charged kaon (+ and -) ET
	// *******************
	TH2F *fHistKaonEtaEET;//ET - Eta vs E 
	TH2F *fHistKaonAccEtaEET;//ET - Eta vs E 
	TH2F *fHistKaonRecEtaEET;//ET - Eta vs E 
	TH2F *fHistKaonMatchEtaEET;//ET - Eta vs E 
	
	TH2F *fHistKaonEtaPtET;// ET - Eta vs pt
	TH2F *fHistKaonAccEtaPtET;// ET - Eta vs pt
	TH2F *fHistKaonRecEtaPtET;// ET - Eta vs pt
	TH2F *fHistKaonMatchEtaPtET;// ET - Eta vs pt
	
	TH2F *fHistKaonEtaET;// ET - Eta
	TH2F *fHistKaonAccEtaET;// ET - Eta
	TH2F *fHistKaonRecEtaET;// ET - Eta
	TH2F *fHistKaonMatchEtaET;// ET - Eta
	
	TH2F *fHistKaonEtaE;// multiplicity - Eta vs E
	TH2F *fHistKaonAccEtaE;// multiplicity - Eta vs E
	TH2F *fHistKaonRecEtaE;// multiplicity - Eta vs E
	TH2F *fHistKaonMatchEtaE;// multiplicity - Eta vs E
	
	TH2F *fHistKaonEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistKaonAccEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistKaonRecEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistKaonMatchEtaPt;// multiplicity - Eta vs pt

	TH1F *fHistKaontotET;// total ET distribution
	TH1F *fHistKaonAcctotET;// total ET distribution
	TH1F *fHistKaonRectotET;// total ET distribution
	TH1F *fHistKaonMatchtotET;// total ET distribution
	
	TH1F *fHistKaonRectotETDep;// total deposited ET distribution
	TH1F *fHistKaonMatchtotETDep;// total deposited ET distribution
	
	TH2F *fHistKaonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
	TH2F *fHistKaonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited

	TH2F *fHistKaonRecEtaPtETDep;// ET deposited - Eta vs pt
	TH2F *fHistKaonMatchEtaPtETDep;// ET deposited - Eta vs pt
	
	TH2F *fHistKaonRecEtaETDep;// ET deposited - Eta
	TH2F *fHistKaonMatchEtaETDep;// ET deposited - Eta
	
	TH2F *fHistKaonRecResEET;// ET - track matching residual vs E
	TH2F *fHistKaonRecResPtET;// ET - track matching residual vs pt
	TH2F *fHistKaonRecResE;// multiplicity - track matching residual vs E
	TH2F *fHistKaonRecResPt;// multiplicity - track matching residual vs pt
	TH2F *fHistKaonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
	TH2F *fHistKaonRecResPtETDep;// ET deposited - track matching residual vs pt
	
	// *******************
	// proton (anti) ET
	// *******************
	TH2F *fHistProtonEtaEET;//ET - Eta vs E 
	TH2F *fHistProtonAccEtaEET;//ET - Eta vs E 
	TH2F *fHistProtonRecEtaEET;//ET - Eta vs E 
	TH2F *fHistProtonMatchEtaEET;//ET - Eta vs E 
	
	TH2F *fHistProtonEtaPtET;// ET - Eta vs pt
	TH2F *fHistProtonAccEtaPtET;// ET - Eta vs pt
	TH2F *fHistProtonRecEtaPtET;// ET - Eta vs pt
	TH2F *fHistProtonMatchEtaPtET;// ET - Eta vs pt
	
	TH2F *fHistProtonEtaET;// ET - Eta
	TH2F *fHistProtonAccEtaET;// ET - Eta
	TH2F *fHistProtonRecEtaET;// ET - Eta
	TH2F *fHistProtonMatchEtaET;// ET - Eta
	
	TH2F *fHistProtonEtaE;// multiplicity - Eta vs E
	TH2F *fHistProtonAccEtaE;// multiplicity - Eta vs E
	TH2F *fHistProtonRecEtaE;// multiplicity - Eta vs E
	TH2F *fHistProtonMatchEtaE;// multiplicity - Eta vs E
	
	TH2F *fHistProtonEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistProtonAccEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistProtonRecEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistProtonMatchEtaPt;// multiplicity - Eta vs pt

	TH1F *fHistProtontotET;// total ET distribution
	TH1F *fHistProtonAcctotET;// total ET distribution
	TH1F *fHistProtonRectotET;// total ET distribution
	TH1F *fHistProtonMatchtotET;// total ET distribution
	
	TH1F *fHistProtonRectotETDep;// total deposited ET distribution
	TH1F *fHistProtonMatchtotETDep;// total deposited ET distribution
	
	TH2F *fHistProtonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
	TH2F *fHistProtonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
	
	TH2F *fHistProtonRecEtaPtETDep;// ET deposited - Eta vs pt
	TH2F *fHistProtonMatchEtaPtETDep;// ET deposited - Eta vs pt
	
	TH2F *fHistProtonRecEtaETDep;// ET deposited - Eta
	TH2F *fHistProtonMatchEtaETDep;// ET deposited - Eta

	TH2F *fHistProtonRecResEET;// ET - track matching residual vs E
	TH2F *fHistProtonRecResPtET;// ET - track matching residual vs pt
	TH2F *fHistProtonRecResE;// multiplicity - track matching residual vs E
	TH2F *fHistProtonRecResPt;// multiplicity - track matching residual vs pt
	TH2F *fHistProtonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
	TH2F *fHistProtonRecResPtETDep;// ET deposited - track matching residual vs pt
	
	// *******************
	// total charged ET
	// *******************
	TH1F *fHistTotChargedtotET;//total ET distribution
	TH1F *fHistTotChargedAcctotET;//total ET distribution
	TH1F *fHistTotChargedRectotET;//total ET distribution
	TH1F *fHistTotChargedRectotETDep;//total deposited ET distribution
	TH1F *fHistTotChargedMatchtotET;//total ET distribution
	TH1F *fHistTotChargedMatchtotETDep;//total deposited ET distribution
	
	// *******************
	// neutron (anti) ET
	// *******************
	TH2F *fHistNeutronEtaEET;//ET - Eta vs E 
	TH2F *fHistNeutronAccEtaEET;//ET - Eta vs E 
	TH2F *fHistNeutronRecEtaEET;//ET - Eta vs E 
	
	TH2F *fHistNeutronEtaPtET;// ET - Eta vs pt
	TH2F *fHistNeutronAccEtaPtET;// ET - Eta vs pt
	TH2F *fHistNeutronRecEtaPtET;// ET - Eta vs pt
	
	TH2F *fHistNeutronEtaET;// ET - Eta
	TH2F *fHistNeutronAccEtaET;// ET - Eta
	TH2F *fHistNeutronRecEtaET;// ET - Eta
	
	TH2F *fHistNeutronEtaE;// multiplicity - Eta vs E
	TH2F *fHistNeutronAccEtaE;// multiplicity - Eta vs E
	TH2F *fHistNeutronRecEtaE;// multiplicity - Eta vs E
	
	TH2F *fHistNeutronEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistNeutronAccEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistNeutronRecEtaPt;// multiplicity - Eta vs pt
	
	TH1F *fHistNeutrontotET;// total ET distribution
	TH1F *fHistNeutronAcctotET;// total ET distribution
	TH1F *fHistNeutronRectotET;// total ET distribution
	TH1F *fHistNeutronRectotETDep;// total deposited ET distribution
	
	TH2F *fHistNeutronRecEtaEDepETDep;// ET deposited - Eta vs E deposited
	TH2F *fHistNeutronRecEtaETDep;// ET deposited - Eta
	
	TH2F *fHistNeutronRecEtaPtETDep;// ET deposited - Eta vs pt
		
	// *******************
	// neutral kaon ET
	// *******************
	TH2F *fHistK0EtaEET;//ET - Eta vs E 
	TH2F *fHistK0RecEtaEET;//ET - Eta vs E 
	
	TH2F *fHistK0EtaPtET;// ET - Eta vs pt
	TH2F *fHistK0RecEtaPtET;// ET - Eta vs pt
	
	TH2F *fHistK0EtaET;// ET - Eta
	TH2F *fHistK0RecEtaET;// ET - Eta
	
	TH2F *fHistK0EtaE;// multiplicity - Eta vs E
	TH2F *fHistK0RecEtaE;// multiplicity - Eta vs E
	
	TH2F *fHistK0EtaPt;// multiplicity - Eta vs pt
	TH2F *fHistK0RecEtaPt;// multiplicity - Eta vs pt

	TH1F *fHistK0totET;// total ET distribution
	TH1F *fHistK0RectotET;// total ET distribution
	
	TH1F *fHistK0RectotETDep;// total deposited ET distribution
	
	TH2F *fHistK0RecEtaEDepETDep;// ET deposited - Eta vs E deposited
	TH2F *fHistK0RecEtaETDep;// ET deposited - Eta
	
	TH2F *fHistK0RecEtaPtETDep;// ET deposited - Eta vs pt
		
	// *******************
	// Lambda(anti) ET
	// *******************
	TH2F *fHistLambdaEtaEET;//ET - Eta vs E 
	TH2F *fHistLambdaRecEtaEET;//ET - Eta vs E 
	
	TH2F *fHistLambdaEtaPtET;// ET - Eta vs pt
	TH2F *fHistLambdaRecEtaPtET;// ET - Eta vs pt
	
	TH2F *fHistLambdaEtaET;// ET - Eta
	TH2F *fHistLambdaRecEtaET;// ET - Eta
	
	TH2F *fHistLambdaEtaE;// multiplicity - Eta vs E
	TH2F *fHistLambdaRecEtaE;// multiplicity - Eta vs E
	
	TH2F *fHistLambdaEtaPt;// multiplicity - Eta vs pt
	TH2F *fHistLambdaRecEtaPt;// multiplicity - Eta vs pt
	
	TH1F *fHistLambdatotET;// total ET distribution
	TH1F *fHistLambdaRectotET;// total ET distribution
	
	TH1F *fHistLambdaRectotETDep;// total deposited ET distribution
	
	TH2F *fHistLambdaRecEtaEDepETDep;// ET deposited - Eta vs E deposited
	TH2F *fHistLambdaRecEtaETDep;// ET deposited - Eta
	
	TH2F *fHistLambdaRecEtaPtETDep;// ET deposited - Eta vs pt

	// *******************
	// total neutral ET
	// *******************
	TH1F *fHistTotNeutraltotET;//total ET distribution
	TH1F *fHistTotNeutralRectotET;//total ET distribution
	TH1F *fHistTotNeutralRectotETDep;//total deposited ET distribution
	
	// *******************
	// total ET
	// *******************
	TH1F *fHistTotaltotET;//total ET distribution
	TH1F *fHistTotalAcctotET;//total ET distribution
	TH1F *fHistTotalRectotET;//total ET distribution
	TH1F *fHistTotalRectotETDep;//total deposited ET distribution
	
	// *******************
	// some checks
	// *******************

	// check produced electrons
	TH1F *fHistElectronFirstMother;// first mother ID
	TH2F *fHistElectronFirstMotherXY;// first mother XY position 
	TH1F *fHistElectronNDaughters;// number of daughters 
	TH1F *fHistElectronDaughters;// daughters ID
	TH2F *fHistElectronDaughtersXY;// daughters XY position

	TH1F *fHistElectronFirstMotherAcc;// first mother ID
	TH2F *fHistElectronFirstMotherXYAcc;// first mother XY position 
	TH1F *fHistElectronNDaughtersAcc;// number of daughters
	TH1F *fHistElectronDaughtersAcc;// daughters ID
	TH2F *fHistElectronDaughtersXYAcc;// daughters XY position

	TH1F *fHistElectronFirstMotherRec;// first mother ID 
	TH2F *fHistElectronFirstMotherXYRec;//  first mother XY position
	TH1F *fHistElectronNDaughtersRec;// number of daughters
	TH1F *fHistElectronDaughtersRec;// daughters ID
	TH2F *fHistElectronDaughtersXYRec;// daughters XY position

	TH1F *fHistNPPElectronFirstMother;//  first mother ID 
	TH2F *fHistNPPElectronFirstMotherXY;// first mother XY position
	TH1F *fHistNPPElectronNDaughters;// number of daughters
	TH1F *fHistNPPElectronDaughters;// daughters ID
	TH2F *fHistNPPElectronDaughtersXY;// daughters XY position
	
	TH1F *fHistNPPElectronFirstMotherAcc;//  first mother ID 
	TH2F *fHistNPPElectronFirstMotherXYAcc;//  first mother XY position
	TH1F *fHistNPPElectronNDaughtersAcc;// number of daughters
	TH1F *fHistNPPElectronDaughtersAcc;// daughters ID
	TH2F *fHistNPPElectronDaughtersXYAcc;// daughters XY position
	
	TH1F *fHistNPPElectronFirstMotherRec;// first mother ID  
	TH2F *fHistNPPElectronFirstMotherXYRec;//  first mother XY position
	TH1F *fHistNPPElectronNDaughtersRec;// number of daughters
	TH1F *fHistNPPElectronDaughtersRec;// daughters ID
	TH2F *fHistNPPElectronDaughtersXYRec;// daughters XY position
	
	// check produced gammas
	TH1F *fHistGammaFirstMother;// first mother ID 
	TH2F *fHistGammaFirstMotherXY;// first mother XY position
	TH1F *fHistGammaNDaughters;// number of daughters
	TH1F *fHistGammaDaughters;// daughters ID
	TH2F *fHistGammaDaughtersXY;// daughters XY position
	TH2F *fHistConvGammaDaughtersXY;// daughters XY position
	TH2F *fHistNonConvGammaDaughtersXY;// daughters XY position
	
	TH1F *fHistGammaFirstMotherAcc;// first mother ID  
	TH2F *fHistGammaFirstMotherXYAcc;//  first mother XY position
	TH1F *fHistGammaNDaughtersAcc;// number of daughters
	TH1F *fHistGammaDaughtersAcc;// daughters ID
	TH2F *fHistGammaDaughtersXYAcc;// daughters XY position
	TH2F *fHistConvGammaDaughtersXYAcc;// daughters XY position
	TH2F *fHistNonConvGammaDaughtersXYAcc;// daughters XY position
	
	TH1F *fHistGammaFirstMotherRec;// first mother ID  
	TH2F *fHistGammaFirstMotherXYRec;//  first mother XY position
	TH1F *fHistGammaNDaughtersRec;// number of daughters
	TH1F *fHistGammaDaughtersRec;// daughters ID
	TH2F *fHistGammaDaughtersXYRec;// daughters XY position
	TH2F *fHistConvGammaDaughtersXYRec;// daughters XY position
	TH2F *fHistNonConvGammaDaughtersXYRec;// daughters XY position
	
	TH1F *fHistNPPGammaFirstMother;// first mother ID 
	TH2F *fHistNPPGammaFirstMotherXY;// first mother XY position
	TH1F *fHistNPPGammaNDaughters;// number of daughters
	TH1F *fHistNPPGammaDaughters;// daughters ID
	TH2F *fHistNPPGammaDaughtersXY;// daughters XY position
	
	TH1F *fHistNPPGammaFirstMotherAcc;// first mother ID  
	TH2F *fHistNPPGammaFirstMotherXYAcc;//  first mother XY position
	TH1F *fHistNPPGammaNDaughtersAcc;// number of daughters
	TH1F *fHistNPPGammaDaughtersAcc;// daughters ID
	TH2F *fHistNPPGammaDaughtersXYAcc;// daughters XY position
	
	TH1F *fHistNPPGammaFirstMotherRec;// first mother ID  
	TH2F *fHistNPPGammaFirstMotherXYRec;//  first mother XY position
	TH1F *fHistNPPGammaNDaughtersRec;// number of daughters
	TH1F *fHistNPPGammaDaughtersRec;// daughters ID
	TH2F *fHistNPPGammaDaughtersXYRec;// daughters XY position

	//check projections
	TH2F *fHistAllERecEMC;// E reconstructed vs E MC	
	TH2F *fHistAllPtRecPtMC;// pt reconstructed vs pt MC
	TH2F *fHistElectronERecEMC;// E reconstructed vs E MC	
	TH2F *fHistGammaERecEMC;// E reconstructed vs E MC
	
	TH2F *fHistChargedRes;// charged particle track matching residual
	TH2F *fHistChargedRes2;// charged particle track matching residual
	TH2F *fHistChargedRes3;// charged particle track matching residual
	TH2F *fHistNeutralRes;// neutral particle track matching residual
	TH2F *fHistElectronRes;// electron track matching residual
	TH2F *fHistGammaRes;// gamma track matching residual
	
	TH2F *fHistIsInAcc;// EMCal acceptance check
	
 private:

  //Declare it private to avoid compilation warning
    AliAnalysisEmEtMonteCarlo & operator = (const AliAnalysisEmEtMonteCarlo & g) ;//cpy assignment
    AliAnalysisEmEtMonteCarlo(const AliAnalysisEmEtMonteCarlo & g) ; // cpy ctor
    ClassDef(AliAnalysisEmEtMonteCarlo, 1);
};

#endif //ALIANALYSISEMETMONTECARLO_H