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

#ifndef ALIANALYSISTASKUE_H
#define ALIANALYSISTASKUE_H
 
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */
 
#include "AliAnalysisTaskSE.h"

class AliESDEvent;
class AliAODEvent;
class TH1F;
class TH2F;
class TH1I;
class TVector3;
       
class  AliAnalysisTaskUE : public AliAnalysisTask
{
 public:
                     AliAnalysisTaskUE(const char* name="AliAnalysisTaskUE");
  virtual           ~AliAnalysisTaskUE() {;}
  
  // Implementation of interface methods
  virtual     void   ConnectInputData(Option_t *);
  virtual     void   CreateOutputObjects();
  virtual     void   Exec(Option_t *option);
  virtual     void   Terminate(Option_t *);
  
  //  Setters
  virtual     void   SetDebugLevel(Int_t level)      { fDebug = level; }
              void   SetPtRangeInHist(Int_t bin, Double_t min, Double_t max) { 
                                          fBinsPtInHist = bin; 
                                          fMinJetPtInHist = min; 
                                          fMaxJetPtInHist = max; 
                                       }
  
              void   SetAnaTopology(Int_t val)    { fAnaType = val;    }          
              void   SetRegionType(Int_t val)     { fRegionType = val; }
              void   SetConeRadius(Double_t val)  { fConeRadius = val; }
              void   SetPtSumOrdering(Int_t val)  { fOrdering = val;   }
      // Jet cuts    
              void   SetJet1EtaCut(Double_t val)      { fJet1EtaCut = val; }
              void   SetJet2DeltaPhiCut(Double_t val) { fJet2DeltaPhiCut = val; }
              void   SetJet2RatioPtCut(Double_t val)  { fJet2RatioPtCut = val; }
              void   SetJet3PtCut(Double_t val)       { fJet3PtCut = val; }
      // track cuts
              void   SetTrackPtCut(Double_t val)  { fTrackPtCut = val; }
              void   SetTrackEtaCut(Double_t val) { fTrackEtaCut = val; }

private:
                     AliAnalysisTaskUE(const  AliAnalysisTaskUE &det);
AliAnalysisTaskUE&   operator=(const  AliAnalysisTaskUE &det);

              void   AnalyseUE();
             Int_t   IsTrackInsideRegion(TVector3 *jetVect, TVector3 *partVect);
              void   CreateHistos();
              void   FillSumPtRegion( Double_t leadingE, Double_t ptMax, Double_t ptMin );
              void   FillAvePartPtRegion( Double_t leadingE, Double_t ptMax, Double_t ptMin );
              void   FillMultRegion( Double_t leadingE, Double_t nTrackPtmax, Double_t nTrackPtmin, Double_t ptMin );
    
       
             Int_t   fDebug;           //  Debug flag
       AliAODEvent*  fAOD;             //! AOD Event 
             TList*  fListOfHistos;    //  Output list of histograms
      
      // Config     
             Int_t   fBinsPtInHist;     //  # bins for Pt histos range
          Double_t   fMinJetPtInHist;   //  min Jet Pt value for histo range
          Double_t   fMaxJetPtInHist;   //  max Jet Pt value for histo range
      
      // Cuts 
             Int_t   fAnaType;          // Analysis type on jet topology: 
                                        //     1=inclusive  (default) 
                                        //     2=back to back inclusive
                                        //     3=back to back exclusive
                                        //     4=gama jet (back to back) ???
                                        //  Minimum bias
                                        //     31 = Semi jet (charged leading particle jets)
                                        //     32 = Random jetcone  ?
                                        //     33 = Swiss chees   ?

             Int_t   fRegionType;       // 1 = transverse regions (default)
                                        // 2 = cone regions   
                                        
             Double_t   fConeRadius;    // if selected Cone-like region type set Radius (=0.7 default)
                                        
             Int_t   fOrdering;         //  Pt and multiplicity summation ordering:
                                        //     1=CDF-like -independent sorting according quantity to be scored: Double sorting- (default)
                                        //       if Pt summation will be scored take Pt minimum between both zones and 
                                        //          fill Pt Max. and Min. histog. accordingly
                                        //       if Multiplicity summation will be scored take Mult. minimum between both zones and 
                                        //          fill Mult Max and Min histog. accordingly
                                        //     2=Marchesini-like (Only Pt sorting: Single sorting)
                                        //          sort only according Pt summation scored, find minimum between both zones and
                                        //          fill Pt and Multiplicity Max and Min summation histog. following only this criterium
                                        //     3=User Selection sorting (NOTE: USER must implement it within cxx)
                                        
      // Jet cuts    
          Double_t   fJet1EtaCut;        // |jet1 eta| < fJet1EtaCut   (fAnaType = 1,2,3)
          Double_t   fJet2DeltaPhiCut;   // |Jet1.Phi - Jet2.Phi| < fJet2DeltaPhiCut (fAnaType = 2,3)
          Double_t   fJet2RatioPtCut;    // Jet2.Pt/Jet1Pt > fJet2RatioPtCut  (fAnaType = 2,3)
          Double_t   fJet3PtCut;         // Jet3.Pt < fJet3PtCut  (fAnaType = 3)
      // track cuts
          Double_t   fTrackPtCut;        // Pt cut of tracks in the regions
          Double_t   fTrackEtaCut;       // Eta cut on tracks in the regions (fRegionType=1)

      // Histograms    ( are owner by fListOfHistos TList )
              TH1F*  fhNJets;                  //!
              TH1F*  fhEleadingPt;             //!
              
              TH1F*  fhMinRegPtDist;
              TH1F*  fhRegionMultMin;
              TH1F*  fhMinRegAvePt; 
              TH1F*  fhMinRegSumPt;            
              TH1F*  fhMinRegMaxPtPart;
              TH1F*  fhMinRegSumPtvsMult;
              
              TH1F*  fhdNdEta_PhiDist;         //!
              TH2F*  fhFullRegPartPtDistVsEt;  //!
              TH2F*  fhTransRegPartPtDistVsEt; //!
              
              TH1F*  fhRegionSumPtMaxVsEt;     //!
              TH1I*  fhRegionMultMax;          //!
              TH1F*  fhRegionMultMaxVsEt;      //!
              TH1F*  fhRegionSumPtMinVsEt;     //!
              TH1F*  fhRegionMultMinVsEt;      //!
              TH1F*  fhRegionAveSumPtVsEt;     //!
              TH1F*  fhRegionDiffSumPtVsEt;    //!
                     
              TH1F*  fhRegionAvePartPtMaxVsEt; //!
              TH1F*  fhRegionAvePartPtMinVsEt; //!
              TH1F*  fhRegionMaxPartPtMaxVsEt; //!
              
      //        TH2F*  fhValidRegion; //! test to be canceled
         
    ClassDef( AliAnalysisTaskUE, 1) // Analysis task for Underlying Event analysis
};
 
#endif
Commit	Line	Data
e989c76a	1	#ifndef ALIANALYSISTASKUE_H
	2	#define ALIANALYSISTASKUE_H
	3
	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	#include "AliAnalysisTaskSE.h"
	8
	9	class AliESDEvent;
	10	class AliAODEvent;
	11	class TH1F;
	12	class TH2F;
	13	class TH1I;
	14	class TVector3;
	15
	16	class AliAnalysisTaskUE : public AliAnalysisTask
	17	{
	18	public:
	19	AliAnalysisTaskUE(const char* name="AliAnalysisTaskUE");
	20	virtual ~AliAnalysisTaskUE() {;}
	21
	22	// Implementation of interface methods
	23	virtual void ConnectInputData(Option_t *);
	24	virtual void CreateOutputObjects();
	25	virtual void Exec(Option_t *option);
	26	virtual void Terminate(Option_t *);
	27
	28	// Setters
	29	virtual void SetDebugLevel(Int_t level) { fDebug = level; }
	30	void SetPtRangeInHist(Int_t bin, Double_t min, Double_t max) {
	31	fBinsPtInHist = bin;
	32	fMinJetPtInHist = min;
	33	fMaxJetPtInHist = max;
	34	}
	35
	36	void SetAnaTopology(Int_t val) { fAnaType = val; }
	37	void SetRegionType(Int_t val) { fRegionType = val; }
	38	void SetConeRadius(Double_t val) { fConeRadius = val; }
	39	void SetPtSumOrdering(Int_t val) { fOrdering = val; }
	40	// Jet cuts
	41	void SetJet1EtaCut(Double_t val) { fJet1EtaCut = val; }
	42	void SetJet2DeltaPhiCut(Double_t val) { fJet2DeltaPhiCut = val; }
	43	void SetJet2RatioPtCut(Double_t val) { fJet2RatioPtCut = val; }
	44	void SetJet3PtCut(Double_t val) { fJet3PtCut = val; }
	45	// track cuts
	46	void SetTrackPtCut(Double_t val) { fTrackPtCut = val; }
	47	void SetTrackEtaCut(Double_t val) { fTrackEtaCut = val; }
	48
	49	private:
	50	AliAnalysisTaskUE(const AliAnalysisTaskUE &det);
	51	AliAnalysisTaskUE& operator=(const AliAnalysisTaskUE &det);
	52
	53	void AnalyseUE();
	54	Int_t IsTrackInsideRegion(TVector3 jetVect, TVector3 partVect);
	55	void CreateHistos();
	56	void FillSumPtRegion( Double_t leadingE, Double_t ptMax, Double_t ptMin );
	57	void FillAvePartPtRegion( Double_t leadingE, Double_t ptMax, Double_t ptMin );
	58	void FillMultRegion( Double_t leadingE, Double_t nTrackPtmax, Double_t nTrackPtmin, Double_t ptMin );
	59
	60
	61	Int_t fDebug; // Debug flag
	62	AliAODEvent* fAOD; //! AOD Event
	63	TList* fListOfHistos; // Output list of histograms
	64
65	// Config
66	Int_t fBinsPtInHist; // # bins for Pt histos range
67	Double_t fMinJetPtInHist; // min Jet Pt value for histo range
68	Double_t fMaxJetPtInHist; // max Jet Pt value for histo range
69
70	// Cuts
71	Int_t fAnaType; // Analysis type on jet topology:
72	// 1=inclusive (default)
73	// 2=back to back inclusive
74	// 3=back to back exclusive
75	// 4=gama jet (back to back) ???
76	// Minimum bias
77	// 31 = Semi jet (charged leading particle jets)
78	// 32 = Random jetcone ?
79	// 33 = Swiss chees ?
80
81	Int_t fRegionType; // 1 = transverse regions (default)
82	// 2 = cone regions
83
84	Double_t fConeRadius; // if selected Cone-like region type set Radius (=0.7 default)
85
86	Int_t fOrdering; // Pt and multiplicity summation ordering:
87	// 1=CDF-like -independent sorting according quantity to be scored: Double sorting- (default)
88	// if Pt summation will be scored take Pt minimum between both zones and
89	// fill Pt Max. and Min. histog. accordingly
90	// if Multiplicity summation will be scored take Mult. minimum between both zones and
91	// fill Mult Max and Min histog. accordingly
92	// 2=Marchesini-like (Only Pt sorting: Single sorting)
93	// sort only according Pt summation scored, find minimum between both zones and
94	// fill Pt and Multiplicity Max and Min summation histog. following only this criterium
95	// 3=User Selection sorting (NOTE: USER must implement it within cxx)
96
97	// Jet cuts
98	Double_t fJet1EtaCut; // \|jet1 eta\| < fJet1EtaCut (fAnaType = 1,2,3)
99	Double_t fJet2DeltaPhiCut; // \|Jet1.Phi - Jet2.Phi\| < fJet2DeltaPhiCut (fAnaType = 2,3)
100	Double_t fJet2RatioPtCut; // Jet2.Pt/Jet1Pt > fJet2RatioPtCut (fAnaType = 2,3)
101	Double_t fJet3PtCut; // Jet3.Pt < fJet3PtCut (fAnaType = 3)
102	// track cuts
103	Double_t fTrackPtCut; // Pt cut of tracks in the regions
104	Double_t fTrackEtaCut; // Eta cut on tracks in the regions (fRegionType=1)
105
106	// Histograms ( are owner by fListOfHistos TList )
107	TH1F* fhNJets; //!
108	TH1F* fhEleadingPt; //!
109
110	TH1F* fhMinRegPtDist;
111	TH1F* fhRegionMultMin;
112	TH1F* fhMinRegAvePt;
113	TH1F* fhMinRegSumPt;
114	TH1F* fhMinRegMaxPtPart;
115	TH1F* fhMinRegSumPtvsMult;
116
117	TH1F* fhdNdEta_PhiDist; //!
118	TH2F* fhFullRegPartPtDistVsEt; //!
119	TH2F* fhTransRegPartPtDistVsEt; //!
120
121	TH1F* fhRegionSumPtMaxVsEt; //!
122	TH1I* fhRegionMultMax; //!
123	TH1F* fhRegionMultMaxVsEt; //!
124	TH1F* fhRegionSumPtMinVsEt; //!
125	TH1F* fhRegionMultMinVsEt; //!
126	TH1F* fhRegionAveSumPtVsEt; //!
127	TH1F* fhRegionDiffSumPtVsEt; //!
128
129	TH1F* fhRegionAvePartPtMaxVsEt; //!
130	TH1F* fhRegionAvePartPtMinVsEt; //!
131	TH1F* fhRegionMaxPartPtMaxVsEt; //!
132
133	// TH2F* fhValidRegion; //! test to be canceled
134
cfff6259	135	ClassDef( AliAnalysisTaskUE, 1) // Analysis task for Underlying Event analysis
e989c76a	136	};
	137
	138	#endif