AliAnalysisTaskGammaTriggerQA.h

Go to the documentation of this file.

#ifndef ALIANLYSISTASKGAMMATRIGGERQA_cxx
#define ALIANLYSISTASKGAMMATRIGGERQA_cxx

#include "AliAnalysisTaskSE.h"
#include "AliESDtrack.h"
#include "AliV0ReaderV1.h"
#include "AliKFConversionPhoton.h"
#include "AliGammaConversionAODBGHandler.h"
#include "AliConversionAODBGHandlerRP.h"
#include "AliCaloPhotonCuts.h"
#include "AliConvEventCuts.h"
#include "AliConversionPhotonCuts.h"
#include "AliConversionMesonCuts.h"
#include "AliAnalysisManager.h"
#include "TProfile2D.h"
#include "TH3.h"
#include "TH3F.h"
#include <vector>
#include <map>

class AliAnalysisTaskGammaTriggerQA : public AliAnalysisTaskSE {
  public:

    AliAnalysisTaskGammaTriggerQA();
    AliAnalysisTaskGammaTriggerQA(const char *name);
    virtual ~AliAnalysisTaskGammaTriggerQA();

    virtual void   UserCreateOutputObjects();
    virtual Bool_t Notify();
    virtual void   UserExec(Option_t *);
    virtual void   Terminate(const Option_t*);
    void InitBack();

    void SetV0ReaderName(TString name){fV0ReaderName=name; return;}
    void SetIsHeavyIon(Int_t flag){
      fIsHeavyIon = flag;
    }

    // base functions for selecting photon and meson candidates in reconstructed data
    void ProcessClusters();

    // MC functions
    void SetIsMC(Int_t isMC){fIsMC=isMC;}

    // switches for additional analysis streams or outputs
    void SetLightOutput(Bool_t flag){fDoLightOutput = flag;}


    // Setting the cut lists for the conversion photons
    void SetEventCutList(Int_t nCuts, TList *CutArray){
      fnCuts = nCuts;
      fEventCutArray = CutArray;
    }

    // Setting the cut lists for the calo photons
    void SetCaloCutList(Int_t nCuts, TList *CutArray){
      fnCuts = nCuts;
      fClusterCutArray = CutArray;
    }

    // Function to set correction task setting
    void SetCorrectionTaskSetting(TString setting) {fCorrTaskSetting = setting;}

  protected:
    AliV0ReaderV1*        fV0Reader;                                            // basic photon Selection Task
    TString               fV0ReaderName;
    TString               fCorrTaskSetting;
    AliVEvent*            fInputEvent;                                          // current event
    AliMCEvent*           fMCEvent;                                             // corresponding MC event
    TList**               fCutFolder;                                           // Array of lists for containers belonging to cut
    TList**               fESDList;                                             // Array of lists with histograms with reconstructed properties
    TList*                fOutputContainer;                                     // Output container
    TList*                fClusterCandidates;                                   
    TList*                fEventCutArray;                                       // List with Event Cuts
    AliConvEventCuts*     fEventCuts;                                           // EventCutObject
    TList*                fClusterCutArray;                                     // List with Cluster Cuts
    AliCaloPhotonCuts*    fCaloPhotonCuts;                                      // CaloPhotonCutObject

    // histograms for rec photon clusters
    TH1F**                fHistoClusGammaPt;                                    
    TH1F**                fHistoClusGammaE;                                     

    // event histograms
    TH1F**                fHistoNEvents;                                        
    TH1F**                fHistoNEventsWOWeight;                                
    TH1F**                fHistoNGoodESDTracks;                                 
    TH1F**                fHistoCent;                                           
    TH1F**                fHistoVertexZ;                                        
    TH1F**                fHistoNGammaCandidates;                               
    TH1F**                fHistoNGammaCandidatesBasic;                          
    TH2F**                fHistoNGoodESDTracksVsNGammaCandidates;               
    TH2F**                fHistoSPDClusterTrackletBackground;                   
    TH1F**                fHistoNV0Tracks;                                      
    TH1F**                fHistoNV0Trigger;                                     
    TH2F**                fHistoNV0TriggerTracks;                               
    TProfile**            fProfileEtaShift;                                     
    TProfile**            fProfileJetJetXSection;                               
    TH1F**                fHistoJetJetNTrials;                                  

    // additional variables
    Double_t              fEventPlaneAngle;                                     // EventPlaneAngle
    TRandom3              fRandom;                                              // random
    Int_t                 fnCuts;                                               // number of cuts to be analysed in parallel
    Int_t                 fiCut;                                                // current cut
    Int_t                 fIsHeavyIon;                                          // switch for pp = 0, PbPb = 1, pPb = 2
    Bool_t                fDoLightOutput;                                       // switch for running light output, kFALSE -> normal mode, kTRUE -> light mode
    Int_t                 fIsMC;                                                // flag for MC information
    Double_t              fWeightJetJetMC;                                      // weight for Jet-Jet MC
    Int_t                 fNCurrentClusterBasic;                                // current number of cluster without minE

  private:
    AliAnalysisTaskGammaTriggerQA(const AliAnalysisTaskGammaTriggerQA&);                  // Prevent copy-construction
    AliAnalysisTaskGammaTriggerQA &operator=(const AliAnalysisTaskGammaTriggerQA&);       // Prevent assignment

    ClassDef(AliAnalysisTaskGammaTriggerQA, 2);
};

#endif