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* 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, 2); // Charged jet analysis for pA,,, increase the last parametre by 1 after a change !!!!!!!!!!!!

};
#endif