Operating mode changes: Regular/Light/SuperLight Cascade output storage for Pb-Pb

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

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

#include "AliAnalysisEtReconstructed.h"
class TParticle;
class TParticlePDG;
class AliESDEvent;
class AliESDtrack;
class AliEMCALTrack;
class TVector3;
class AliEMCALGeometry;
class AliExternalTrackParam;
class AliStack;

class AliAnalysisEmEtReconstructed : public AliAnalysisEtReconstructed
{

public:
   
  AliAnalysisEmEtReconstructed();
  virtual ~AliAnalysisEmEtReconstructed();

    virtual Int_t AnalyseEvent(AliVEvent* event);

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

protected:
	
	Double_t GetTrackPID(const AliESDtrack *track) const;
	

protected:

	Double_t fResCut;//track matching residual cut
	
	Double_t fAllRectotETDep;//total ET deposited - all particles
	Double_t fElectronMatchtotETDep;// total Et deposited - matched electrons
	Double_t fNeutralRectotET;// total ET - neutral particles
	Double_t fTotEMRectotET;// total electromagnetic ET
	Double_t fMuonMatchtotETDep, fPionMatchtotETDep, fKaonMatchtotETDep, fProtonMatchtotETDep;// total ET deposited - hadrons
	Double_t fTotChargedMatchtotETDep;// total Et deposited - matched chagred particles
	Double_t fTotalRectotETDep;// total ET deposited
	
	AliESDEvent *fESD;//!// ESD object
/*   	AliEMCALGeometry *fGeoUt;//!// EMCal geometry object */

	// *******************
	// all ET
	// *******************
	THnSparseF* fHistAllRecETDep;//!// ET deposited - all particles
	THnSparseF* fHistAllRec;//!// mutliplicity - all particles
	TH1F *fHistAllRectotETDep;//!// total ET deposited - all particles
	
	// *******************
	// electron ET reconstructed in EMCal
	// *******************
	THnSparseF* fHistElectronRecETDep;//!// Et deposited - matched electrons
	THnSparseF* fHistElectronRec;//!// multiplicity - matched electrons
	TH1F *fHistElectronMatchtotETDep;//!// total Et deposited - matched electrons
	
	TH2F *fHistElectronRecdEdxP;//!// electron dEdx vs p

	// *******************
	// Neutral ET reconstructed in EMCal
	// *******************
	TH1F *fHistNeutralRectotET;//!// total ET - neutral particles 

	// *******************
	// total EM ET reconstructed in EMCal
	// *******************
	TH1F *fHistTotEMRectotET;//!// total electromagnetic ET

	// *******************
	// muon ET (+ and -)
	// *******************
	THnSparseF* fHistMuonRecETDep;//!// Et deposited
	THnSparseF* fHistMuonRec;//!// multiplicity
	TH1F *fHistMuonMatchtotETDep;//!// total Et deposited

	TH2F *fHistMuonRecdEdxP;//!// dEdx vs p
	
	// *******************
	// pion ET (+ and -)
	// *******************
	THnSparseF* fHistPionRecETDep;//!// Et deposited 
	THnSparseF* fHistPionRec;//!// multiplicity
	TH1F *fHistPionMatchtotETDep;//!// total Et deposited

	TH2F *fHistPionRecdEdxP;//!// dEdx vs p

	// *******************
	// charged kaon (+ and -) ET
	// *******************
	THnSparseF* fHistKaonRecETDep;//!// Et deposited
	THnSparseF* fHistKaonRec;//!// multiplicity
	TH1F *fHistKaonMatchtotETDep;//!// total Et deposited

	TH2F *fHistKaonRecdEdxP;//!// dEdx vs p
	
	// *******************
	// proton (anti) ET
	// *******************
	THnSparseF* fHistProtonRecETDep;//!// Et deposited
	THnSparseF* fHistProtonRec;//!// multiplicity
	TH1F *fHistProtonMatchtotETDep;//!// total Et deposited

	TH2F *fHistProtonRecdEdxP;//!// dEdx vs p
	
	// *******************
	// total charged ET
	// *******************
	TH1F *fHistTotChargedMatchtotETDep;//!// total Et deposited - all charged particles 
	
	// *******************
	// total ET
	// *******************
	TH1F *fHistTotalRectotETDep;//!// total Et deposited - all particles
	
	//few checks
	TH2F *fHistDeltaRZ;//!// track-cluster matching residual
	
 private:
  //Declare it private to avoid compilation warning
    AliAnalysisEmEtReconstructed & operator = (const AliAnalysisEmEtReconstructed & g) ;//cpy assignment
    AliAnalysisEmEtReconstructed(const AliAnalysisEmEtReconstructed & g) ; // cpy ctor
    ClassDef(AliAnalysisEmEtReconstructed, 1);
};

#endif // ALIANALYSISEMETRECONSTRUCTED