AliAnalysisTaskEA.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEA_H
#define ALIANALYSISTASKEA_H


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

#include <vector>
using std::vector;

#include "AliAnalysisTaskEmcalJet.h"
#include "AliEventCuts.h"
#include "AliFiducialCut.h"
#include "AliEMCALRecoUtils.h"


// ANALYSIS OF EVENT ACTIVITY WITH HIGH PT HADRON BIAS
// Author Filip Krizek   (25 OCT 2018)

class AliAnalysisTaskEA : public AliAnalysisTaskEmcalJet {
   public:

   enum Analysis {
      kDetLevel=0, 
      kPartLevel=1,
      fkTTbins = 20 //maximum number of TT bins
   };



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

   
  // ######### SETTERS/GETTERS
  void        SetMC(Int_t bMC){
                 fMC = bMC;
              }

  void        SetUseDefaultVertexCut (Bool_t val) {fUseDefaultVertexCut = val;}  
  void        SetUsePileUpCut (Bool_t val) {fUsePileUpCut = val;} 

  void        SetAcceptanceWindows(Double_t trackEta){ 
                 fTrackEtaWindow  = trackEta; 
              } 


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

  void        SetExternalRhoTaskName(const char* name) {fRhoTaskName = name;}
  void        SetExternalRhoTaskNameMC(const char* name) {fRhoTaskNameMC = name;}

  void        SetTrackContainerName(const char* name){ fMyTrackContainerName = name;}
  void        SetMCParticleContainerName(const char* name){ fMyParticleContainerName=name;}

  void        SetJetContainerName(const char* name){ fMyJetContainerName = name;}
  void        SetMCJetContainerName(const char* name){ fMyJetParticleContainerName = name;}
  void        SetHadronTT(Int_t tl,Int_t th){ fHadronTTLowPt[fnHadronTTBins] = tl; fHadronTTHighPt[fnHadronTTBins] = th; fnHadronTTBins++; }
  void        SetJetChTT(Int_t tl,Int_t th){  fJetChTTLowPt[fnJetChTTBins] = tl;   fJetChTTHighPt[fnJetChTTBins] = th;   fnJetChTTBins++;  }

  Bool_t      PassedGATrigger();
  Bool_t      PassedMinBiasTrigger();
  Bool_t      PreSelection(AliVCluster* cluster);
  Bool_t      FinalClusterCuts(AliVCluster* cluster);

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


 private:


  // ######### CHECK FUNCTIONS
  Bool_t      IsTrackInAcceptance(AliVParticle* track, Bool_t isGen=0);  
  Bool_t      IsEventInAcceptance(AliVEvent* event);     


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

  Double_t    GetExternalRho(Bool_t isMC); 
  
  // ########## 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
  

  // ########## SOURCE INFORMATION
  TString     fMyTrackContainerName;       // name of detector level track container 
  TString     fMyParticleContainerName;    // name of particle level MC particle container
  TString     fMyJetContainerName;         // name of detector level jet container 
  TString     fMyJetParticleContainerName; // name of particle level MC jet container 

  AliTrackContainer    *fTrkContainerDetLevel;        
  AliParticleContainer *fParticleContainerPartLevel;  
  AliJetContainer      *fJetContainerDetLevel;        
  AliJetContainer      *fJetContainerPartLevel;       

  TString     fRhoTaskName;           // name of rho CMS bg task for this analysis
  TString     fRhoTaskNameMC;         // MC name of rho CMS bg task for this analysis


  // ########## CENTRALITY
  TTree   *fCentralityTree;     
  AliMultSelection*  fMultSelection; 

  char    fTrigClass[1000];     
  Bool_t  fIsMinBiasTrig;       
  Bool_t  fIsEmcalTrig;         

  Float_t fCentralityV0A;       
  Float_t fCentralityV0C;       
  Float_t fCentralityCL1;       
  Float_t fCentralityZNA;       
  Float_t fCentralityZNC;       

  Double_t fxVertex;            
  Double_t fyVertex;            
  Double_t fzVertex;            
  Bool_t   fVertexer3d;         
  //
  Int_t    fNTracklets;         
  Int_t    fNClusters[2];       
  //
  Int_t    fIsV0ATriggered;     
  Int_t    fIsV0CTriggered;     
  Float_t  fMultV0A;            
  Float_t  fMultV0C;            
  Float_t  fRingMultV0[8];      
  //
  Float_t  fZEM1Energy;         
  Float_t  fZEM2Energy;         

  Float_t  fZNCtower[5];        
  Float_t  fZPCtower[5];        
  Float_t  fZNAtower[5];        
  Float_t  fZPAtower[5];        
  Float_t  fZNCtowerLG[5];      
  Float_t  fZPCtowerLG[5];      
  Float_t  fZNAtowerLG[5];      
  Float_t  fZPAtowerLG[5];      


  // ########## CUTS 
  Double_t    fTrackEtaWindow;        //gc +- window in eta for tracks  
  Double_t    fMinTrackPt;            //gc Min track pt to be accepted  

  // ########## GENERAL ////VARS
  Bool_t              fMC;                    //  real data  or MC  flag
  AliAnalysisUtils*   fHelperClass;           
  Bool_t              fInitializedLocal;      


   TH1I               *fHistEvtSelection;     

   TH1F     *fhVertexZ;  

   TH2F     *fhTrackPhiIncl;
   TH2F     *fhTrackEtaIncl;
   TH2F     *fhJetPhiIncl;
   TH2F     *fhJetEtaIncl;


   Double_t fZVertexCut; // vertex cut in z 

   Int_t    fnHadronTTBins; // number of TT bins
   Int_t    fnJetChTTBins; // number of TT bins
   Int_t    fHadronTTLowPt[fkTTbins]; // low pt TT range
   Int_t    fHadronTTHighPt[fkTTbins]; // low pt TT range
   Int_t    fJetChTTLowPt[fkTTbins]; // low pt TT range
   Int_t    fJetChTTHighPt[fkTTbins]; // low pt TT range
  

   Int_t    fHadronTT[fkTTbins]; 
   Int_t    fJetChTT[fkTTbins];  

   AliAnalysisTaskEA(const AliAnalysisTaskEA&);
   AliAnalysisTaskEA& operator=(const AliAnalysisTaskEA&);

   ClassDef(AliAnalysisTaskEA, 1); // Charged jet analysis for pAliAnalysisTaskHJetSpectra/home/fkrizek/z501.ALIC

};
#endif