AliAnalysisTaskEmcalTriggerPatchClusterMatch.h

Go to the documentation of this file.

#ifndef AliAnalysisTaskEmcalTriggerPatchClusterMatch_h
#define AliAnalysisTaskEmcalTriggerPatchClusterMatch_h
//#ifndef ALIANALYSISTASKEMCALTRIGGERPATCHCLUSTERMATCH_H
//#define ALIANALYSISTASKEMCALTRIGGERPATCHCLUSTERMATCH_H

//-------------------------------------------------------------------------
// 1) Analysis task to identify the cluster that fired the trigger patch
// 2) perform some QA on the patch / cluster
// 3) and pass the saved out collection of cluster(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 THnSparse;
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 "AliAnalysisTaskEmcal.h"

class AliAnalysisTaskEmcalTriggerPatchClusterMatch : public AliAnalysisTaskEmcal {
 public:
  AliAnalysisTaskEmcalTriggerPatchClusterMatch();
  AliAnalysisTaskEmcalTriggerPatchClusterMatch(const char *name);
  virtual ~AliAnalysisTaskEmcalTriggerPatchClusterMatch();

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

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

  // cuts and biases
  virtual void            SetPatchEcut(Double_t pE)             { fPatchECut = pE; }
  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            SetClusterTriggeredEventname(const char *clcol)           { fClusterTriggeredEventname = clcol; }
  TString                 GetClusterTriggeredEvent() const          { return fClusterTriggeredEventname; }

  Int_t    GetLeadingCellId(const AliVCluster *clus) const;
  Double_t GetEnergyLeadingCell(const AliVCluster *clus) const;
  Double_t GetECross(Int_t absID) 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
  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           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            fClusterTriggeredEventname; // name of cluster that triggered event collection

  AliEMCALTriggerPatchInfo      *fMaxPatch;
  MainPatchType                 fMainPatchType;

  TH1F     *fhNEvents;                         
  TProfile *fhTriggerbit;                      
  TH3F     *fh3PtEtaPhiTracks;                 
  TH3F     *fh3PtEtaPhiTracksOnEmcal;          
  TH3F     *fh3PtEtaPhiTracksToProp;           
  TH3F     *fh3PtEtaPhiTracksProp;             
  TH3F     *fh3PtEtaPhiTracksNoProp;           
  TH3F     *fh3EEtaPhiCluster;                 
  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;        
  TH1F     *fHistClusEnergy;                   

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

  TH2F     *fHistdPhidEtaPatchCluster[16];     

  TClonesArray               *fClusterTriggeredEvent;                    
  TClonesArray               *fRecalcTriggerPatches;                    

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

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

  ClassDef(AliAnalysisTaskEmcalTriggerPatchClusterMatch, 4)
};
#endif