AliAnalysisTaskEmcalJetTriggerQA.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEMCALJETTRIGGERQA_H
#define ALIANALYSISTASKEMCALJETTRIGGERQA_H

class TH1;
class TH2;
class TH3;
class TH3F;
class TProfile;
class TClonesArray;
class TArrayI;

#include <TRef.h>
#include <TBits.h>
#include <TMath.h>

#include <AliVEvent.h>

#include "AliAnalysisTaskEmcalJet.h"

class AliAnalysisTaskEmcalJetTriggerQA : public AliAnalysisTaskEmcalJet {
 public:
  AliAnalysisTaskEmcalJetTriggerQA();
  AliAnalysisTaskEmcalJetTriggerQA(const char *name);
  virtual ~AliAnalysisTaskEmcalJetTriggerQA();

  void                        UserCreateOutputObjects();
  void                        Terminate(Option_t *option);

  Bool_t                      SelectEvent();              //decides if event is used for analysis
  void                        FillTriggerPatchHistos();

  //Setters
  void SetDebug(Int_t d)                    { fDebug = d;}
  void SetTriggerClass(const char *n)       { fTriggerClass = n; }
  void SetNFastorPatch(Int_t i)             { fNFastOR = i;}
 
  void SetContainerFull(Int_t c)            { fContainerFull      = c;}
  void SetContainerCharged(Int_t c)         { fContainerCharged   = c;}

  void SetMainTriggerTypeCat(TriggerCategory cat, Bool_t b) {fMainTrigCat = cat; fMainTrigSimple = b;}

  Int_t    GetLeadingCellId(const AliVCluster *clus) const;
  Double_t GetEnergyLeadingCell(const AliVCluster *clus) const;
  Double_t GetECross(Int_t absID) const;

  Double_t GetZ(const AliVParticle *trk, const AliEmcalJet *jet)       const;
  Double_t GetZ(Double_t trkPx, Double_t trkPy, Double_t trkPz, Double_t jetPx, Double_t jetPy, Double_t jetPz) const;

 protected:
  Bool_t                      FillHistograms()   ;
  Bool_t                      Run()              ;
  Float_t                     RelativeEP(Double_t objAng, Double_t EPAng) const;
  Bool_t                      TestFilterBit(Int_t trigBit, UInt_t bitJetTrig) const {return (Bool_t) ((trigBit & bitJetTrig) != 0);}

 private:
  Bool_t             fDebug;                 // debug level
  TString            fTriggerClass;          // trigger class to analyze EJ1 or EJ2    
  Int_t              fContainerFull;         // number of container with full jets DET
  Int_t              fContainerCharged;      // number of container with charged jets DET
  Double_t           fMaxPatchEnergy;        // energy of patch with largest energy (offline)
  Double_t           fMaxPatchADCEnergy;     // energy of patch with largest energy from online ADC
  Int_t              fTriggerType;           // trigger type
  Int_t              fNFastOR;               // size of trigger patch fNFastORxfNFastOR
  TriggerCategory    fMainTrigCat;           // trigger category for main trigger
  Bool_t             fMainTrigSimple;        // use offline trigger instead of online

  TH1F     *fhNEvents;                         
  TProfile *fhTriggerbit;                      
  TH2F     *fHistRhovsCentFull;                
  TH2F     *fHistRhovsCentCharged;             
  TH3F     *fh3PtEtaPhiTracks;                 
  TH3F     *fh3PtEtaPhiTracksOnEmcal;          
  TH3F     *fh3PtEtaPhiTracksToProp;           
  TH3F     *fh3PtEtaPhiTracksProp;             
  TH3F     *fh3PtEtaPhiTracksNoProp;           
  TH2F     *fh2CentPtJetFull;                  
  TH2F     *fh2CentPtJetCharged;               
  TH3F     *fh3PtEtaPhiJetFull;                
  TH3F     *fh3PtEtaPhiJetCharged;             
  TH2F     *fh2NJetsPtFull;                    
  TH2F     *fh2NJetsPtCharged;                 
  TH3F     *fh3PtEtaAreaJetFull;               
  TH3F     *fh3PtEtaAreaJetCharged;            
  TH2F     *fh2PtNConstituentsCharged;         
  TH2F     *fh2PtNConstituents;                
  TH2F     *fh2PtMeanPtConstituentsCharged;    
  TH2F     *fh2PtMeanPtConstituentsNeutral;    
  TH2F     *fh2PtNEF;                          
  TH3F     *fh3NEFEtaPhi;                      
  TH2F     *fh2NEFNConstituentsCharged;        
  TH2F     *fh2NEFNConstituentsNeutral;        
  TH2F     *fh2Ptz;                            
  TH2F     *fh2PtzCharged;                     
  TH2F     *fh2PtLeadJet1VsLeadJet2;           
  TH3F     *fh3EEtaPhiCluster;                 
  TH3F     *fh3PtLeadJet1VsPatchEnergy;        
  TH3F     *fh3PtLeadJet2VsPatchEnergy;        
  TH3F     *fh3PtLeadJet1PatchEnergyVZEROAmp;  
  TH3F     *fh3PtLeadJet1RawPatchEnergyVZEROAmp;  
  TH3F     *fh3PatchEnergyEtaPhiCenterJ1;      
  TH3F     *fh3PatchEnergyEtaPhiCenterJ2;      
  TH3F     *fh3PatchEnergyEtaPhiCenterJ1J2;    
  TH3F     *fh3PatchADCEnergyEtaPhiCenterJ1;   
  TH3F     *fh3PatchADCEnergyEtaPhiCenterJ2;   
  TH3F     *fh3PatchADCEnergyEtaPhiCenterJ1J2; 
  TH3F     *fh3PatchADCEnergyEtaPhiCenterAll;  
  TH3F     *fh3EEtaPhiCell;                    
  TH2F     *fh2ECellVsCent;                    
  TH2F     *fh2CellEnergyVsTime;               
  TH3F     *fh3EClusELeadingCellVsTime;        
  TH3F     *fh3JetReacCent;                    

  AliAnalysisTaskEmcalJetTriggerQA(const AliAnalysisTaskEmcalJetTriggerQA&);            // not implemented
  AliAnalysisTaskEmcalJetTriggerQA &operator=(const AliAnalysisTaskEmcalJetTriggerQA&); // not implemented

  ClassDef(AliAnalysisTaskEmcalJetTriggerQA, 15)
};
#endif