AliAnalysisTaskJetPP.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKJETPP_H
#define ALIANALYSISTASKJETPP_H


class TH1I;
class TH1F;
class TH2F;
class TH2D;
class TH1D;
class TArrayD;
class THnSparse;
class TProfile;
class TList;
class AliEmcalList;
class TClonesArray;
class TString;
class AliEmcalJet;
class AliRhoParameter;
class AliVParticle;
class AliLog;
class AliAnalysisUtils;
class TRandom3;
class AliJetContainer;
class AliParticleContainer;
class AliClusterContainer;

#include <vector>
#include "AliAnalysisTaskEmcalJet.h"

using std::vector;



class AliAnalysisTaskJetPP : public AliAnalysisTaskEmcalJet {
   public:

   enum MyContainer {
      kContainerOne = 0, 
      kContainerTwo = 1,
      kContainerThree = 2,
      kContainerFour = 3
   };



   // ######### CONTRUCTORS/DESTRUCTORS AND STD FUNCTIONS
   AliAnalysisTaskJetPP();
   AliAnalysisTaskJetPP(const char *name);
   virtual ~AliAnalysisTaskJetPP();
   void     UserCreateOutputObjects();
   void     Terminate(Option_t *);

   
  // ######### SETTERS/GETTERS

   void        SetUseDefaultVertexCut (Bool_t val) {fUseDefaultVertexCut = val;}  
   void        SetUsePileUpCut (Bool_t val) {fUsePileUpCut = val;} 
   void        SetSignalJetMinArea(Double_t minArea) {fMinJetArea = minArea;} 
   void        SetMC(Bool_t a) {fIsMC = a;}

   void        SetAcceptanceWindows(Double_t trackEta, Double_t signalJetRadius){ 
                  fTrackEtaWindow     = trackEta; 
                  fSignalJetRadius    = signalJetRadius; 
                  fSignalJetEtaWindow = fTrackEtaWindow-fSignalJetRadius; 
               } 

   void       SetCentralityType(const char* type){
                 fCentralityType = type;    
              } 

   void       SetVertexCut(Double_t vz){   fZVertexCut = vz; }   
   void       SetMinTrackPt(Double_t mpt){ fMinTrackPt = mpt;}


   Bool_t   RetrieveEventObjects();
   Bool_t   Run();
   Bool_t   FillHistograms();

   private:

   Double_t EstimateBgKT(AliJetContainer *jetCont);  // median p/A of kt jets
   Double_t EstimateLocalBg(AliJetContainer *jetCont, AliParticleContainer *trkCont);  // local bg in a perpendicular cone

   Double_t Convert(Double_t input);
   // ######### CHECK FUNCTIONS
   Bool_t      IsTrackInAcceptance(AliVParticle* track, Bool_t isprimary=0);  
   Bool_t      IsEventInAcceptance(AliVEvent* event);     
   Bool_t      IsSignalJetInAcceptance(AliEmcalJet* jet, Bool_t suppressGhost=1); 


   // ######### STANDARD FUNCTIONS
   void      ExecOnceLocal();                    

   // ########## USAGE TRIGGERS 

   Bool_t              fUseDefaultVertexCut;   // trigger if automatic vertex cut from helper class should be done 
   Bool_t              fUsePileUpCut;          // trigger if pileup cut should be done
   Bool_t              fIsMC;                  // Use MC

   // ########## JET/DIJET/RC PROPERTIES
   Double_t            fSignalJetRadius;       // Radius for the signal jets
   // ########## CUTS 
   Double_t            fSignalJetEtaWindow;    // +- window in eta for signal jets 
   Double_t            fTrackEtaWindow;        //gc +- window in eta for tracks  
   Double_t            fMinTrackPt;            //gc Min track pt to be accepted  
   Double_t            fMinJetArea;            // Min jet area to be accepted
   TString             fCentralityType;        //gc Used centrality estimate (V0A, V0C, V0M, ...) 

   Double_t            fZVertexCut;            //zvertex cut

   // ########## GENERAL ////VARS
   AliAnalysisUtils*   fHelperClass;           
   Bool_t              fInitializedLocal;      


   TH1D*  fhJetPt;     
   TH1D*  fhKTJetPt;   
   TH1D*  fhJetPtRho;  
   TH1D*  fhJetPtConeRho; 

   TH1D*  fhAtimesRho; 
   TH1D*  fhCuts;      
   TH1D*  fhJetConstituentPt; 
   TH1D*  fhTrackPt;   

   TH1D*  fhRho;       
   TH1D*  fhConeRho;   
   TH2D*  fhJetAreaPt; 
   TH2D*  fhJetEtaPt;  

   TH2D* fhTrackEtaPt; 
   TH2D* fhJetPhiPt;   
   TH2D* fhTrackPhiPt; 
   TH2D* fhRemx;       

   TH1D* fhZVertex;    
   TH1D* fhZVertexBC;  
   TH1D* fhXVertex;    
   TH1D* fhYVertex;    

   TH1D* fhPrimGenTrkPt; 
   TH1D* fhGenJetPt;     
   TH1D* fhRecTrkPt;     
   TH1D* fhFakeTrkPt;    

   TH1D* fhMult;          
   TH2D* fhTrackPhiCG;    
   TH2D* fhTrackPhiTPCG;  
   TH2D* fhInvPtQVsPhi[2];

   TH2D* fhInvPtQVsEta[2];      
   TH2D* fhInvPtQVsPhiASide[2]; 
   TH2D* fhInvPtQVsPhiCSide[2]; 
   TH2D* fhSigmaPtOverPtVsPt[2];


   AliAnalysisTaskJetPP(const AliAnalysisTaskJetPP&);
   AliAnalysisTaskJetPP& operator=(const AliAnalysisTaskJetPP&);

   ClassDef(AliAnalysisTaskJetPP, 1); // Charged jet analysis for pA

};
#endif