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;}

   void SetASPDCvsTCut(Float_t a) { fASPDCvsTCut = a; }//PILEUP
   void SetBSPDCvsTCut(Float_t b) { fBSPDCvsTCut = b; }

   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 EstimateBgKTCMS(AliJetContainer *jetCont); //CMS-style background form kT jets

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


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

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

   Bool_t              fUseDefaultVertexCut;   
   Bool_t              fUsePileUpCut;          
   Bool_t              fIsMC;          

   // ########## JET/DIJET/RC PROPERTIES
   Double_t            fSignalJetRadius;       
   // ########## CUTS 
   Double_t            fSignalJetEtaWindow;    
   Double_t            fTrackEtaWindow;        
   Double_t            fMinTrackPt;            
   Double_t            fMinJetArea;            
   TString             fCentralityType;        

   Double_t            fZVertexCut;            
   Float_t             fASPDCvsTCut;           
   Float_t             fBSPDCvsTCut;           

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


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

   TH1D* fhAtimesRhoMedian; 
   TH1D* fhAtimesRhoCone; 
   TH1D* fhAtimesRhoCMS; 
   TH1D* fhCuts; 
   TH1D* fhJetConstituentPt; 
   TH1D* fhTrackPt; 

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

   TH2D* fhAktJetEtaPhi; 
   TH2D* fhKtJetEtaPhi; 
   TH2D* fhTrackEtaPt; 
   TH2D* fhGenTrackEtaPt; 
   TH2D* fhJetPhiPt; 

   TH2D* fhTrackPhiPt; 
   TH2D* fhTrackEtaPhi; 
   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, 4); // Charged jet analysis for pA,,, increase the last parametre by 1 after a change 
};
#endif