AliAnalysisTaskEmcalTriggerPatchJetMatch.h

Go to the documentation of this file.

#ifndef AliAnalysisTaskEmcalTriggerPatchJetMatch_h
#define AliAnalysisTaskEmcalTriggerPatchJetMatch_h
//#ifndef ALIANALYSISTASKEMCALTRIGGERPATCHJETMATCH_H
//#define ALIANALYSISTASKEMCALTRIGGERPATCHJETMATCH_H

//-------------------------------------------------------------------------
// 1) Analysis task to identify the jet whose cluster fired the trigger patch
// 2) perform some QA on the patch / cluster / jet
// 3) and pass the saved out collection of jet(s) to other tasks
//
// currently set up for GA trigger
//
// Author: Joel Mazer (joel.mazer@cern.ch)
//-------------------------------------------------------------------------
/* Copyright(c) 1998-2014, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

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

#include <TRef.h>
#include <TBits.h>
#include <TMath.h>
#include <AliVEvent.h>
#include "AliEMCALTriggerPatchInfo.h"
#include "AliEmcalTriggerSetupInfo.h"
#include "AliAnalysisTaskEmcalJet.h"

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

  enum MainPatchType {
    kManual = 0,    //just select highest energy patch in array
    kEmcalJet = 1   //use functionality of AliAnalysisTaskEmcal
  };

//  void                        ExecOnce();
  void                        UserCreateOutputObjects();
  void                        Terminate(Option_t *option);
  Bool_t                      SelectEvent();              //decides if event is used for analysis
  void                        FillTriggerPatchHistos();
  void                        SetAttachToEvent(Bool_t a)                      {fAttachToEvent = a;}

  //Setters
  void SetDebug(Int_t d)                    { fDebug = d;}
  void SetTriggerClass(const char *n)       { fTriggerClass = n; }
  void SetNFastorPatch(Int_t i)             { fNFastOR = i;}
 
  // containers
  void SetJetContainer(Int_t c)            { fJetContainer      = c;}

  void SetMainPatchType(MainPatchType t)    { fMainPatchType      = t;}
  void SetMainTriggerTypeCat(TriggerCategory cat, Bool_t b) {fMainTrigCat = cat; fMainTrigSimple = b;}

  // cuts and biases
  virtual void            SetJetPtcut(Double_t jpt)             { fJetPtCut = jpt; } // jet pt cut
  virtual void            SetPatchEcut(Double_t pE)             { fPatchECut = pE; } // jet pt cut
  virtual void            SetTrkBias(Double_t b)                { fTrkBias    = b; }  //require a track with pt > b in jet
  virtual void            SetClusBias(Double_t b)               { fClusBias   = b; }  //require a cluster with pt > b in jet

  // give comments setter and various switches
  void                    SetdoComments(Bool_t comm)               { doComments = comm; } // give comment switch
  void                    SetUseALLrecalcPatches(Bool_t useall)    { fUseALLrecalcPatches = useall; } // use all firing (offline) patches

  virtual void            SetJetTriggeredEventname(const char *jcol)           { fJetTriggeredEventname = jcol; }
  TString                 GetJetTriggeredEvent() const          { return fJetTriggeredEventname; }
  virtual void            SetCaloClustersName(const char *cn)   { fCaloClustersName=cn; }

  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:
  void                        ExecOnce();
  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);}

  // Trigger bit - do i need?????
  void                        ExtractMainPatch();
  TH1*                        FillTriggerPatchQA(TH1* h, UInt_t t, AliEMCALTriggerPatchInfo* fPatch); // filled trigger patch QA
  TH1*                        FillEventTriggerQA(TH1* h, UInt_t t); // fill event trigger QA
  Bool_t                      CorrelateToTrigger(Double_t etaclust, Double_t phiclust, TList *triggerpatches) const;

 private:
  Bool_t             fDebug;                 // debug level
  Bool_t             fAttachToEvent;         // attach local rho to the event
  TString            fTriggerClass;          // trigger class to analyze EJ or GA patches
  Int_t              fJetContainer;          // number of container of jets
  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
  TriggerCategory    fTriggerCategory;       // trigger category
  Bool_t             fMainTrigSimple;        // use offline trigger instead of online from AliAnalysisTaskEmcal::GetMainTriggerPatch
  Double_t       fJetPtCut;              // jet pt to cut on for correlations
  Double_t           fPatchECut;             // patch energy cut
  Double_t           fTrkBias;               // track bias
  Double_t           fClusBias;              // cluster bias
  Bool_t             doComments;             // summary (debugging comments)
  Bool_t             fUseALLrecalcPatches;   // use all/ just max recalulated (offline) patches
  TString      fJetTriggeredEventname; // name of jet that triggered event collection
  TString            fCaloClustersName;      // name of Calo Cluster collection

  AliEMCALTriggerPatchInfo      *fMaxPatch;
  MainPatchType                 fMainPatchType;

  TH1F     *fhNEvents;                         
  TProfile *fhTriggerbit;                      
  TH2F     *fHistRhovsCent;                    
  TH3F     *fh3PtEtaPhiTracks;                 
  TH3F     *fh3PtEtaPhiTracksOnEmcal;          
  TH3F     *fh3PtEtaPhiTracksToProp;           
  TH3F     *fh3PtEtaPhiTracksProp;             
  TH3F     *fh3PtEtaPhiTracksNoProp;           
  TH2F     *fh2CentPtJet;                      
  TH3F     *fh3PtEtaPhiJet;                    
  TH2F     *fh2NJetsPt;                        
  TH3F     *fh3PtEtaAreaJet;                   
  TH2F     *fh2PtNConstituentsCharged;         
  TH2F     *fh2PtNConstituents;                
  TH2F     *fh2PtMeanPtConstituentsCharged;    
  TH2F     *fh2PtMeanPtConstituentsNeutral;    
  TH2F     *fh2PtNEF;                          
  TH3F     *fh3NEFEtaPhi;                      
  TH2F     *fh2NEFNConstituentsCharged;        
  TH2F     *fh2NEFNConstituentsNeutral;        
  TH2F     *fh2Ptz;                            
  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     *fh3PatchEnergyEtaPhiCenterG1;      
  TH3F     *fh3PatchEnergyEtaPhiCenterG2;      
  TH3F     *fh3PatchEnergyEtaPhiCenterG1G2;    
  TH3F     *fh3PatchADCEnergyEtaPhiCenterG1;   
  TH3F     *fh3PatchADCEnergyEtaPhiCenterG2;   
  TH3F     *fh3PatchADCEnergyEtaPhiCenterG1G2; 
  TH3F     *fh3PatchADCEnergyEtaPhiCenterAll;  
  TH3F     *fh3EEtaPhiCell;                    
  TH2F     *fh2ECellVsCent;                    
  TH2F     *fh2CellEnergyVsTime;               
  TH3F     *fh3EClusELeadingCellVsTime;        
  TH3F     *fh3JetReacCent;                    

  TH1F     *fHistClusEnergy;                   
  TH1F     *fHistClusofJetEnergy;              
  TH1F     *fHistEventSelectionQA;             
  TH1F     *fhQAinfoAllPatchesCounter;         
  TH1F     *fhQAinfoCounter;                   
  TH1F     *fhQAmaxinfoCounter;                
  TH1F     *fhRecalcGammaPatchEnergy;          
  TH1F     *fhGammaLowPatchEnergy;             
  TH1F     *fhGammaLowSimplePatchEnergy;       
  TH1F     *fhRecalcJetPatchEnergy;            
  TH1F     *fhJetLowPatchEnergy;               
  TH1F     *fhJetLowSimplePatchEnergy;         
  TH1F     *fhMainTriggerPatchEnergy;          

  TH2F     *fHistdPhidEtaPatchJetCluster[16];  

  TClonesArray               *fJetTriggeredEvent;                    
  TClonesArray               *fRecalcTriggerPatches;                 

  THnSparse                  *fhnPatchMaxClus;
  THnSparse                  *fhnPatchMatch;
  THnSparse                  *fhnPatchMatch2;
  THnSparse                  *fhnPatchMatchJetLeadClus;

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

  ClassDef(AliAnalysisTaskEmcalTriggerPatchJetMatch, 4)
};
#endif