[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / 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     */

//_________________________________________________________________________
// 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 ----
#include "AliAnaCaloTrackCorrBaseClass.h"

class AliAnaParticleJetFinderCorrelation : public AliAnaCaloTrackCorrBaseClass {
       
 public:   
  AliAnaParticleJetFinderCorrelation() ;              // default ctor
  virtual ~AliAnaParticleJetFinderCorrelation() { ; } // virtual dtor

  // General methods
  
  void     InitParameters();
  
  TList *  GetCreateOutputObjects();

  void     MakeAnalysisFillAOD() ;
  
  void     MakeAnalysisFillHistograms() ;
  
  Int_t    SelectJet(AliAODPWG4Particle * particle, const AliAODEvent * event) const ;

  void     Print(const Option_t * opt)           const;
  
  // Settings
  
  Bool_t   OnlyIsolated()                        const { return fSelectIsolated              ; }
  void     SelectIsolated(Bool_t select)               { fSelectIsolated = select            ; }
  
  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:

  //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
  Bool_t     fSelectIsolated ;    // Select only trigger particles isolated
  
  // Histograms
  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
  
  AliAnaParticleJetFinderCorrelation(              const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor
  AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ; // cpy assignment
  
  ClassDef(AliAnaParticleJetFinderCorrelation,2)
  
 } ;

#endif //ALIANAPARTICLEJETFINDERCORRELATION_H
Commit	Line	Data
	1	#ifndef ALIANAPARTICLEJETFINDERCORRELATION_H
	2	#define ALIANAPARTICLEJETFINDERCORRELATION_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	//_________________________________________________________________________
	7	// Class that contains the algorithm for the analysis of particle (direct gamma) - jet
	8	// (standard jet found with JETAN) correlation
	9	// Particle and jet for correlation found by independent algorithms.
	10	// For Example direct isolated photon found in AliAnaGammaDirect and the jet with JETAN
	11	//
	12	//-- Author: Gustavo Conesa (INFN-LNF)
	13
	14	// --- ROOT system ---
	15	class TH2F;
	16
	17	//---- Analysis system ----
	18	#include "AliAnaCaloTrackCorrBaseClass.h"
	19
	20	class AliAnaParticleJetFinderCorrelation : public AliAnaCaloTrackCorrBaseClass {
	21
	22	public:
	23	AliAnaParticleJetFinderCorrelation() ; // default ctor
	24	virtual ~AliAnaParticleJetFinderCorrelation() { ; } // virtual dtor
	25
	26	// General methods
	27
	28	void InitParameters();
	29
	30	TList * GetCreateOutputObjects();
	31
	32	void MakeAnalysisFillAOD() ;
	33
	34	void MakeAnalysisFillHistograms() ;
	35
	36	Int_t SelectJet(AliAODPWG4Particle * particle, const AliAODEvent * event) const ;
	37
	38	void Print(const Option_t * opt) const;
	39
	40	// Settings
	41
	42	Bool_t OnlyIsolated() const { return fSelectIsolated ; }
	43	void SelectIsolated(Bool_t select) { fSelectIsolated = select ; }
	44
	45	Float_t GetConeSize() const { return fConeSize ; }
	46	Float_t GetPtThresholdInCone() const { return fPtThresholdInCone ; }
	47	Double_t GetDeltaPhiMaxCut() const { return fDeltaPhiMaxCut ; }
	48	Double_t GetDeltaPhiMinCut() const { return fDeltaPhiMinCut ; }
	49	Double_t GetRatioMaxCut() const { return fRatioMaxCut ; }
	50	Double_t GetRatioMinCut() const { return fRatioMinCut ; }
	51	Bool_t AreJetRefTracks() const { return fUseJetRefTracks ; }
	52	Bool_t IsCorrelationMadeInHistoMaker() const { return fMakeCorrelationInHistoMaker ; }
	53
	54	void SetConeSize(Float_t cone) { fConeSize = cone ; }
	55	void SetPtThresholdInCone(Float_t pt) { fPtThresholdInCone = pt ; }
	56	void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
	57	{ fDeltaPhiMaxCut =phimax; fDeltaPhiMinCut =phimin ; }
	58	void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
	59	{ fRatioMaxCut =ratiomax; fRatioMinCut = ratiomin ; }
	60	void UseJetRefTracks(Bool_t use) { fUseJetRefTracks = use ; }
	61	void SetMakeCorrelationInHistoMaker(Bool_t make) { fMakeCorrelationInHistoMaker = make ; }
	62
	63	private:
	64
	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	Bool_t fSelectIsolated ; // Select only trigger particles isolated
	76
	77	// Histograms
	78	TH2F * fhDeltaEta; //! Difference of jet eta and trigger particle eta as function of trigger particle pT
	79	TH2F * fhDeltaPhi; //! Difference of jet phi and trigger particle phi as function of trigger particle pT
	80	TH2F * fhDeltaPt; //! Difference of jet pT and trigger particle pT as function of trigger particle pT
	81	TH2F * fhPtRatio; //! Ratio of jet pT and trigger particle pT as function of trigger particle pT
	82	TH2F * fhPt; //! jet pT vs trigger particle pT
	83
	84	TH2F * fhFFz ; //! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig
	85	TH2F * fhFFxi; //! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet)
	86	TH2F * fhFFpt; //! Jet particle pt distribution in cone
	87	TH2F * fhNTracksInCone; //! jet multiplicity in cone
	88
	89	AliAnaParticleJetFinderCorrelation( const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor
	90	AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ; // cpy assignment
	91
	92	ClassDef(AliAnaParticleJetFinderCorrelation,2)
	93
	94	} ;
	95
	96	#endif //ALIANAPARTICLEJETFINDERCORRELATION_H
	97
	98
	99