[u/mrichter/AliRoot.git] / PWG4 / AliAnaParticleJetFinderCorrelation.h

#ifndef AliAnaParticleJetFinderCorrelation_H
#define AliAnaParticleJetFinderCorrelation_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice     */
/* $Id: AliAnaParticleJetFinderCorrelation.h 21839 2007-10-29 13:49:42Z gustavo $ */


//_________________________________________________________________________
// Class that contains the algorithm for the analysis of particle (direct gamma) - jet 
// (standard jet found with JETAN) correlation 
// Particle and jet for correlation found by independent algorithms.
// For Example direct isolated photon found in AliAnaGammaDirect and the jet with JETAN
//
//-- Author: Gustavo Conesa (INFN-LNF)

// --- ROOT system ---
class TH2F;

//---- Analysis system ----
class AliLog;
#include "AliAnaPartCorrBaseClass.h"
     
class AliAnaParticleJetFinderCorrelation : public AliAnaPartCorrBaseClass {
       
	public: 
       
		AliAnaParticleJetFinderCorrelation() ; // default ctor
		AliAnaParticleJetFinderCorrelation(const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor
		AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ;//cpy assignment
		virtual ~AliAnaParticleJetFinderCorrelation() {;} //virtual dtor
              
		TList * GetCreateOutputObjects();
       
		void InitParameters();
       
		void Print(const Option_t * opt) const;
       
		Float_t  GetConeSize() const {return fConeSize ; }
		Float_t  GetPtThresholdInCone() const {return fPtThresholdInCone ; }	   
		Double_t GetDeltaPhiMaxCut() const {return fDeltaPhiMaxCut ; }
		Double_t GetDeltaPhiMinCut() const {return fDeltaPhiMinCut ; }
		Double_t GetRatioMaxCut() const {return fRatioMaxCut ; }
		Double_t GetRatioMinCut() const {return fRatioMinCut ; }  	   
		Bool_t	 AreJetRefTracks() const {return fUseJetRefTracks ; }
		Bool_t   IsCorrelationMadeInHistoMaker() const {return fMakeCorrelationInHistoMaker ;} 
		
		void SetConeSize(Float_t cone) {fConeSize = cone; }
		void SetPtThresholdInCone(Float_t pt){fPtThresholdInCone = pt; };	   
		void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
			{fDeltaPhiMaxCut =phimax;  fDeltaPhiMinCut =phimin;}
		void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
			{fRatioMaxCut =ratiomax;  fRatioMinCut = ratiomin ; }
		void UseJetRefTracks(Bool_t use) {fUseJetRefTracks = use ; }	
	    void SetMakeCorrelationInHistoMaker(Bool_t make) {fMakeCorrelationInHistoMaker = make ;}	
		
	 private:
    
		Int_t SelectJet(AliAODParticleCorrelation * particle) ;
	
	    void MakeAnalysisFillAOD() ;
        void MakeAnalysisFillHistograms() ;
			
  private:
       
	//selection parameters  
	Double_t   fDeltaPhiMaxCut ; //! Minimum Delta Phi Gamma-Leading
	Double_t   fDeltaPhiMinCut ; //!  Maximum Delta Phi Gamma-Leading
	Double_t   fRatioMaxCut ; //! Jet/ particle Ratio cut maximum
	Double_t   fRatioMinCut ; //! Jet/particle Ratio cut minimum
	   
	Double_t   fConeSize  ;         //! Jet cone size 
	Double_t   fPtThresholdInCone ; //! Jet pT threshold in jet cone
	Bool_t     fUseJetRefTracks ; //! Use track references from JETAN not the AOD tracks
	Bool_t	   fMakeCorrelationInHistoMaker ; //!Make particle-jet correlation in histogram maker
	   
	TH2F * fhDeltaEta; //! Difference of jet eta and trigger particle eta as function of trigger particle pT
	TH2F * fhDeltaPhi;  //! Difference of jet phi and trigger particle phi as function of trigger particle pT
	TH2F * fhDeltaPt; //! Difference of jet pT and trigger particle pT as function of trigger particle pT
	TH2F * fhPtRatio; //! Ratio of jet pT and trigger particle pT as function of trigger particle pT
	TH2F * fhPt; //! jet pT vs trigger particle pT 
		
	TH2F * fhFFz ; //! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig
	TH2F * fhFFxi; //! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet)
	TH2F * fhFFpt; //! Jet particle pt distribution in cone
	TH2F * fhNTracksInCone   ; //! jet multiplicity in cone
	   
	ClassDef(AliAnaParticleJetFinderCorrelation,1)
 } ;


#endif //AliAnaParticleJetFinderCorrelation_H
Commit	Line	Data
76518957	1	#ifndef AliAnaParticleJetFinderCorrelation_H
76518957	2	#define AliAnaParticleJetFinderCorrelation_H
d92b41ad	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
d92b41ad	4	* See cxx source for full Copyright notice */
76518957	5	/* $Id: AliAnaParticleJetFinderCorrelation.h 21839 2007-10-29 13:49:42Z gustavo $ */
76518957	6
d92b41ad	7
d92b41ad	8	//_________________________________________________________________________
76518957	9	// Class that contains the algorithm for the analysis of particle (direct gamma) - jet
	10	// (standard jet found with JETAN) correlation
	11	// Particle and jet for correlation found by independent algorithms.
	12	// For Example direct isolated photon found in AliAnaGammaDirect and the jet with JETAN
	13	//
d92b41ad	14	//-- Author: Gustavo Conesa (INFN-LNF)
d92b41ad	15
76518957	16	// --- ROOT system ---
	17	class TH2F;
	18
	19	//---- Analysis system ----
	20	class AliLog;
	21	#include "AliAnaPartCorrBaseClass.h"
d92b41ad	22
76518957	23	class AliAnaParticleJetFinderCorrelation : public AliAnaPartCorrBaseClass {
d92b41ad	24
76518957	25	public:
d92b41ad	26
76518957	27	AliAnaParticleJetFinderCorrelation() ; // default ctor
	28	AliAnaParticleJetFinderCorrelation(const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor
	29	AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ;//cpy assignment
	30	virtual ~AliAnaParticleJetFinderCorrelation() {;} //virtual dtor
d92b41ad	31
76518957	32	TList * GetCreateOutputObjects();
d92b41ad	33
76518957	34	void InitParameters();
d92b41ad	35
76518957	36	void Print(const Option_t * opt) const;
d92b41ad	37
76518957	38	Float_t GetConeSize() const {return fConeSize ; }
	39	Float_t GetPtThresholdInCone() const {return fPtThresholdInCone ; }
	40	Double_t GetDeltaPhiMaxCut() const {return fDeltaPhiMaxCut ; }
	41	Double_t GetDeltaPhiMinCut() const {return fDeltaPhiMinCut ; }
	42	Double_t GetRatioMaxCut() const {return fRatioMaxCut ; }
	43	Double_t GetRatioMinCut() const {return fRatioMinCut ; }
	44	Bool_t AreJetRefTracks() const {return fUseJetRefTracks ; }
	45	Bool_t IsCorrelationMadeInHistoMaker() const {return fMakeCorrelationInHistoMaker ;}
	46
	47	void SetConeSize(Float_t cone) {fConeSize = cone; }
	48	void SetPtThresholdInCone(Float_t pt){fPtThresholdInCone = pt; };
	49	void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
	50	{fDeltaPhiMaxCut =phimax; fDeltaPhiMinCut =phimin;}
	51	void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
	52	{fRatioMaxCut =ratiomax; fRatioMinCut = ratiomin ; }
	53	void UseJetRefTracks(Bool_t use) {fUseJetRefTracks = use ; }
	54	void SetMakeCorrelationInHistoMaker(Bool_t make) {fMakeCorrelationInHistoMaker = make ;}
	55
	56	private:
	57
	58	Int_t SelectJet(AliAODParticleCorrelation * particle) ;
	59
	60	void MakeAnalysisFillAOD() ;
	61	void MakeAnalysisFillHistograms() ;
	62
d92b41ad	63	private:
d92b41ad	64
76518957	65	//selection parameters
	66	Double_t fDeltaPhiMaxCut ; //! Minimum Delta Phi Gamma-Leading
	67	Double_t fDeltaPhiMinCut ; //! Maximum Delta Phi Gamma-Leading
	68	Double_t fRatioMaxCut ; //! Jet/ particle Ratio cut maximum
	69	Double_t fRatioMinCut ; //! Jet/particle Ratio cut minimum
	70
	71	Double_t fConeSize ; //! Jet cone size
	72	Double_t fPtThresholdInCone ; //! Jet pT threshold in jet cone
	73	Bool_t fUseJetRefTracks ; //! Use track references from JETAN not the AOD tracks
	74	Bool_t fMakeCorrelationInHistoMaker ; //!Make particle-jet correlation in histogram maker
	75
	76	TH2F * fhDeltaEta; //! Difference of jet eta and trigger particle eta as function of trigger particle pT
	77	TH2F * fhDeltaPhi; //! Difference of jet phi and trigger particle phi as function of trigger particle pT
	78	TH2F * fhDeltaPt; //! Difference of jet pT and trigger particle pT as function of trigger particle pT
	79	TH2F * fhPtRatio; //! Ratio of jet pT and trigger particle pT as function of trigger particle pT
	80	TH2F * fhPt; //! jet pT vs trigger particle pT
	81
	82	TH2F * fhFFz ; //! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig
	83	TH2F * fhFFxi; //! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet)
	84	TH2F * fhFFpt; //! Jet particle pt distribution in cone
	85	TH2F * fhNTracksInCone ; //! jet multiplicity in cone
	86
	87	ClassDef(AliAnaParticleJetFinderCorrelation,1)
d92b41ad	88	} ;
	89
	90
76518957	91	#endif //AliAnaParticleJetFinderCorrelation_H
d92b41ad	92
	93
	94