[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / 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 ----
#include "AliAnaPartCorrBaseClass.h"

class AliAnaParticleJetFinderCorrelation : public AliAnaPartCorrBaseClass {
       
 public:   
  AliAnaParticleJetFinderCorrelation() ; // default ctor
  virtual ~AliAnaParticleJetFinderCorrelation() {;} //virtual dtor
 private:
  AliAnaParticleJetFinderCorrelation(const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor
  AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ;//cpy assignment

 public:

  TList * GetCreateOutputObjects();
  
  Bool_t OnlyIsolated() const {return fSelectIsolated ; }
  void SelectIsolated(Bool_t select) {fSelectIsolated = select ; }
  
  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(AliAODPWG4Particle * particle, AliAODEvent * event) const ;
  
  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
  Bool_t     fSelectIsolated ;      // Select only trigger particles isolated
  
  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,2)
 } ;


#endif //ALIANAPARTICLEJETFINDERCORRELATION_H
Commit	Line	Data
1c5acb87	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	/* $Id: AliAnaParticleJetFinderCorrelation.h 21839 2007-10-29 13:49:42Z gustavo $ */
	6
	7
	8	//_________________________________________________________________________
	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	//
	14	//-- Author: Gustavo Conesa (INFN-LNF)
	15
	16	// --- ROOT system ---
	17	class TH2F;
	18
	19	//---- Analysis system ----
1c5acb87	20	#include "AliAnaPartCorrBaseClass.h"
1c5acb87	21
	22	class AliAnaParticleJetFinderCorrelation : public AliAnaPartCorrBaseClass {
	23
78219bac	24	public:
477d6cee	25	AliAnaParticleJetFinderCorrelation() ; // default ctor
78219bac	26	virtual ~AliAnaParticleJetFinderCorrelation() {;} //virtual dtor
78219bac	27	private:
477d6cee	28	AliAnaParticleJetFinderCorrelation(const AliAnaParticleJetFinderCorrelation & g) ; // cpy ctor
477d6cee	29	AliAnaParticleJetFinderCorrelation & operator = (const AliAnaParticleJetFinderCorrelation & g) ;//cpy assignment
78219bac	30
	31	public:
	32
477d6cee	33	TList * GetCreateOutputObjects();
	34
	35	Bool_t OnlyIsolated() const {return fSelectIsolated ; }
	36	void SelectIsolated(Bool_t select) {fSelectIsolated = select ; }
	37
	38	void InitParameters();
	39
	40	void Print(const Option_t * opt) const;
	41
a3aebfff	42	Float_t GetConeSize() const { return fConeSize ; }
	43	Float_t GetPtThresholdInCone() const { return fPtThresholdInCone ; }
	44	Double_t GetDeltaPhiMaxCut() const { return fDeltaPhiMaxCut ; }
	45	Double_t GetDeltaPhiMinCut() const { return fDeltaPhiMinCut ; }
	46	Double_t GetRatioMaxCut() const { return fRatioMaxCut ; }
	47	Double_t GetRatioMinCut() const { return fRatioMinCut ; }
	48	Bool_t AreJetRefTracks() const { return fUseJetRefTracks ; }
477d6cee	49	Bool_t IsCorrelationMadeInHistoMaker() const {return fMakeCorrelationInHistoMaker ;}
	50
	51	void SetConeSize(Float_t cone) {fConeSize = cone; }
	52	void SetPtThresholdInCone(Float_t pt){fPtThresholdInCone = pt; };
	53	void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
	54	{fDeltaPhiMaxCut =phimax; fDeltaPhiMinCut =phimin;}
	55	void SetRatioCutRange(Double_t ratiomin, Double_t ratiomax)
	56	{fRatioMaxCut =ratiomax; fRatioMinCut = ratiomin ; }
	57	void UseJetRefTracks(Bool_t use) {fUseJetRefTracks = use ; }
	58	void SetMakeCorrelationInHistoMaker(Bool_t make) {fMakeCorrelationInHistoMaker = make ;}
	59
	60	private:
	61
f37fa8d2	62	Int_t SelectJet(AliAODPWG4Particle * particle, AliAODEvent * event) const ;
477d6cee	63
	64	void MakeAnalysisFillAOD() ;
	65	void MakeAnalysisFillHistograms() ;
	66
	67	private:
	68
	69	//selection parameters
	70	Double_t fDeltaPhiMaxCut ; //! Minimum Delta Phi Gamma-Leading
	71	Double_t fDeltaPhiMinCut ; //! Maximum Delta Phi Gamma-Leading
	72	Double_t fRatioMaxCut ; //! Jet/ particle Ratio cut maximum
	73	Double_t fRatioMinCut ; //! Jet/particle Ratio cut minimum
	74
	75	Double_t fConeSize ; //! Jet cone size
	76	Double_t fPtThresholdInCone ; //! Jet pT threshold in jet cone
	77	Bool_t fUseJetRefTracks ; //! Use track references from JETAN not the AOD tracks
	78	Bool_t fMakeCorrelationInHistoMaker ; //!Make particle-jet correlation in histogram maker
	79	Bool_t fSelectIsolated ; // Select only trigger particles isolated
	80
	81	TH2F * fhDeltaEta; //! Difference of jet eta and trigger particle eta as function of trigger particle pT
	82	TH2F * fhDeltaPhi; //! Difference of jet phi and trigger particle phi as function of trigger particle pT
	83	TH2F * fhDeltaPt; //! Difference of jet pT and trigger particle pT as function of trigger particle pT
	84	TH2F * fhPtRatio; //! Ratio of jet pT and trigger particle pT as function of trigger particle pT
	85	TH2F * fhPt; //! jet pT vs trigger particle pT
	86
	87	TH2F * fhFFz ; //! Accepted reconstructed jet fragmentation function, z=ptjet/pttrig
	88	TH2F * fhFFxi; //! Accepted reconstructed jet fragmentation function, xsi = ln(pttrig/ptjet)
	89	TH2F * fhFFpt; //! Jet particle pt distribution in cone
	90	TH2F * fhNTracksInCone ; //! jet multiplicity in cone
	91
	92	ClassDef(AliAnaParticleJetFinderCorrelation,2)
1c5acb87	93	} ;
	94
	95
	96	#endif //ALIANAPARTICLEJETFINDERCORRELATION_H
	97
	98
	99