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

#ifndef ALIANAGAMMAISOLCUT_H
#define ALIANAGAMMAISOLCUT_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice     */
/* $Id$ */

/* History of cvs commits:
 *
 * $Log$
 * Revision 1.1  2007/01/23 17:17:29  schutz
 * New Gamma package
 *
 *
 */

//_________________________________________________________________________

// Class for the analysis of gamma  (gamma-jet, 
// gamma-hadron(Arleo, TODO))
// This class make the Isolation Cut analysis for 2 methods, pt cone sum 
// and pt particle threshold and with different  cone sizes and pt thresholds.

//*-- Author: Gustavo Conesa (INFN-LNF)

// --- ROOT system ---
#include <TROOT.h>
#include <TChain.h>
#include "TTask.h"
#include "TArrayD.h"
#include "TChain.h"
#include <TH2F.h>
#include <TTree.h> 
#include <TParticle.h> 
#include "AliAnaGammaDirect.h" 

class AliESD ; 
 
class AliAnaGammaIsolCut : public AliAnaGammaDirect {

public: 

  AliAnaGammaIsolCut(const char *name) ; // default ctor
  AliAnaGammaIsolCut(const AliAnaGammaIsolCut & g) ; // cpy ctor
  virtual ~AliAnaGammaIsolCut() ; //virtual dtor
  virtual void Exec(Option_t * opt = "") ;
  virtual void Init(Option_t * opt = "");
  virtual void Terminate(Option_t * opt = "");

  void Print(const Option_t * opt)const;

  Int_t    GetNCones()                  const {return fNCones ; }
  Int_t    GetNPtThres()                const {return fNPtThres ; }
  Float_t GetConeSizes(Int_t i)      const {return fConeSizes[i] ; }
  Float_t GetPtThresholds(Int_t i)  const {return fPtThresholds[i] ; }
  
  void SetNCone(Int_t ncs)              {fNCones = ncs ; }
  void SetNPtThres(Int_t npt)        {fNPtThres = npt; }
  void SetConeSizes(Int_t i, Float_t r)         {fConeSizes[i] = r ; }
  void SetPtThresholds(Int_t i, Float_t pt)   {fPtThresholds[i] = pt; }
  
  void MakeHistos() ;
  
  
 private:

  TObjArray  *fOutputContainer ; //! output data container

  Int_t         fNCones   ; //Number of cone sizes to test
  Int_t         fNPtThres ; //Number of ptThres to test
  Float_t     fConeSizes[10] ; // Arrat with cones to test
  Float_t     fPtThresholds[10] ; // Array with pt thresholds to test

  //Histograms  
 
  TH1F * fhPtCandidate ;
  TH1F* fhPtThresIsolated[10][10] ;
  TH2F* fhPtSumIsolated[10] ;

  ClassDef(AliAnaGammaIsolCut,0)
} ;
 

#endif //ALIANAGAMMAISOLCUT_H
Commit	Line	Data
2d11ded3	1	#ifndef ALIANAGAMMAISOLCUT_H
	2	#define ALIANAGAMMAISOLCUT_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5	/* $Id$ */
	6
	7	/* History of cvs commits:
	8	*
	9	* $Log$
	10	* Revision 1.1 2007/01/23 17:17:29 schutz
	11	* New Gamma package
	12	*
	13	*
	14	*/
	15
	16	//_________________________________________________________________________
	17
	18	// Class for the analysis of gamma (gamma-jet,
	19	// gamma-hadron(Arleo, TODO))
	20	// This class make the Isolation Cut analysis for 2 methods, pt cone sum
	21	// and pt particle threshold and with different cone sizes and pt thresholds.
	22
	23	//*-- Author: Gustavo Conesa (INFN-LNF)
	24
	25	// --- ROOT system ---
	26	#include <TROOT.h>
	27	#include <TChain.h>
	28	#include "TTask.h"
	29	#include "TArrayD.h"
	30	#include "TChain.h"
	31	#include <TH2F.h>
	32	#include <TTree.h>
	33	#include <TParticle.h>
	34	#include "AliAnaGammaDirect.h"
	35
	36	class AliESD ;
	37
	38	class AliAnaGammaIsolCut : public AliAnaGammaDirect {
	39
	40	public:
	41
	42	AliAnaGammaIsolCut(const char *name) ; // default ctor
	43	AliAnaGammaIsolCut(const AliAnaGammaIsolCut & g) ; // cpy ctor
	44	virtual ~AliAnaGammaIsolCut() ; //virtual dtor
	45	virtual void Exec(Option_t * opt = "") ;
	46	virtual void Init(Option_t * opt = "");
	47	virtual void Terminate(Option_t * opt = "");
	48
	49	void Print(const Option_t * opt)const;
	50
	51	Int_t GetNCones() const {return fNCones ; }
	52	Int_t GetNPtThres() const {return fNPtThres ; }
	53	Float_t GetConeSizes(Int_t i) const {return fConeSizes[i] ; }
	54	Float_t GetPtThresholds(Int_t i) const {return fPtThresholds[i] ; }
	55
	56	void SetNCone(Int_t ncs) {fNCones = ncs ; }
	57	void SetNPtThres(Int_t npt) {fNPtThres = npt; }
	58	void SetConeSizes(Int_t i, Float_t r) {fConeSizes[i] = r ; }
	59	void SetPtThresholds(Int_t i, Float_t pt) {fPtThresholds[i] = pt; }
	60
	61	void MakeHistos() ;
	62
	63
	64	private:
65
66	TObjArray *fOutputContainer ; //! output data container
67
68	Int_t fNCones ; //Number of cone sizes to test
69	Int_t fNPtThres ; //Number of ptThres to test
70	Float_t fConeSizes[10] ; // Arrat with cones to test
71	Float_t fPtThresholds[10] ; // Array with pt thresholds to test
72
73	//Histograms
74
75	TH1F * fhPtCandidate ;
76	TH1F* fhPtThresIsolated[10][10] ;
77	TH2F* fhPtSumIsolated[10] ;
78
79	ClassDef(AliAnaGammaIsolCut,0)
80	} ;
81
82
83	#endif //ALIANAGAMMAISOLCUT_H
84
85
86