PWGJE/EMCALJetTasks/AliJetResponseMaker.h

   1 #ifndef ALIJETRESPONSEMAKER_H
   2 #define ALIJETRESPONSEMAKER_H
   3
   4 // $Id$
   5
   6 class AliGenPythiaEventHeader;
   7 class TClonesArray;
   8 class TH1F;
   9 class TH2F;
  10
  11 #include "AliAnalysisTaskEmcalJet.h"
  12
  13 class AliJetResponseMaker : public AliAnalysisTaskEmcalJet {
  14  public:
  15   AliJetResponseMaker();
  16   AliJetResponseMaker(const char *name);
  17   virtual ~AliJetResponseMaker();
  18
  19   void                        UserCreateOutputObjects();
  20
  21   void                        SetMCJetsName(const char *n)       { fMCJetsName      = n; }
  22   void                        SetMCTracksName(const char *n)     { fMCTracksName    = n; }
  23   void                        SetMaxDistance(Double_t d)         { fMaxDistance     = d; }
  24   void                        SetDoWeighting(Bool_t d = kTRUE)   { fDoWeighting     = d; }
  25   void                        SetMCFiducial(Bool_t f = kTRUE)    { fMCFiducial      = f; }
  26   void                        SetEventWeightHist(Bool_t h)       { fEventWeightHist = h; }
  27   void                        SetPtHardBin(Int_t b)              { fSelectPtHardBin = b; }
  28
  29  protected:
  30   Bool_t                      IsEventSelected();
  31   Bool_t                      AcceptJet(AliEmcalJet* jet, Bool_t /*bias*/ = kFALSE, Bool_t /*upCut*/ = kFALSE)   const;
  32   void                        ExecOnce();
  33   void                        DoJetLoop(TClonesArray *jets1, TClonesArray *jets2, Bool_t mc);
  34   Bool_t                      FillHistograms();
  35   Bool_t                      RetrieveEventObjects();
  36   Bool_t                      Run();
  37
  38   TString                     fMCTracksName;              // name of MC particle collection
  39   TString                     fMCJetsName;                // name of MC jet collection
  40   Double_t                    fMaxDistance;               // maximum distance between matched particle and detector level jets
  41   Bool_t                      fDoWeighting;               // = true, weight using trials and given x section
  42   Bool_t                      fEventWeightHist;           // = true create event weight histogram
  43   Bool_t                      fMCFiducial;                // = true MC jets in fiducial acceptance
  44   Double_t                    fMCminEta;                  // MC jets minimum eta
  45   Double_t                    fMCmaxEta;                  // MC jets maximum eta
  46   Double_t                    fMCminPhi;                  // MC jets minimum phi
  47   Double_t                    fMCmaxPhi;                  // MC jets maximum phi
  48   Int_t                       fSelectPtHardBin;           // Select one pt hard bin for analysis
  49
  50   AliGenPythiaEventHeader    *fPythiaHeader;              //!event Pythia header
  51   Double_t                    fEventWeight;               //!event weight
  52   Int_t                       fPtHardBin;                 //!event pt hard bin
  53   Int_t                       fNTrials;                   //!event trials
  54   TClonesArray               *fMCTracks;                  //!MC particles
  55   TClonesArray               *fMCJets;                    //!MC jets
  56   // General histograms
  57   TH1F                       *fHistZVertex;               //!Z vertex position
  58   TH1F                       *fHistNTrials;               //!total number of trials per pt hard bin
  59   TH1F                       *fHistEvents;                //!total number of events per pt hard bin
  60   TH1F                       *fHistEventWeight[11];       //!event weight
  61   // Particle level jets
  62   TH2F                       *fHistMCJetPhiEta;           //!phi-eta distribution of jets
  63   TH1F                       *fHistMCJetsPt;              //!inclusive jet pt spectrum
  64   TH2F                       *fHistMCJetPhiEtaFiducial;   //!phi-eta distribution of jets in fiducial acceptance (plus lead hadron bias)
  65   TH1F                       *fHistMCJetsPtFiducial;      //!inclusive jet pt spectrum in fiducial acceptance (plus lead hadron bias)
  66   TH2F                       *fHistMCJetsNEFvsPt;         //!jet neutral energy fraction vs. jet pt
  67   TH2F                       *fHistMCJetsZvsPt;           //!constituent Pt over Jet Pt ratio vs. jet pt
  68   // Detector level jets
  69   TH2F                       *fHistJetPhiEta;             //!phi-eta distribution of jets
  70   TH1F                       *fHistJetsPt;                //!inclusive jet pt spectrum
  71   TH2F                       *fHistJetsNEFvsPt;           //!jet neutral energy fraction vs. jet pt
  72   TH2F                       *fHistJetsZvsPt;             //!constituent Pt over Jet Pt ratio vs. jet pt
  73   // Detector-particle level matching
  74   TH1F                       *fHistClosestDistance;       //!distance between closest particle to detector level jet
  75   TH1F                       *fHistClosestDeltaPhi;       //!delta phi between closest particle to detector level jet
  76   TH1F                       *fHistClosestDeltaEta;       //!delta eta between closest particle to detector level jet
  77   TH1F                       *fHistClosestDeltaPt;        //!delta pt between closest particle to detector level jet
  78   TH1F                       *fHistNonMatchedMCJetPt;     //!non-matched mc jet pt distribution
  79   TH1F                       *fHistNonMatchedJetPt;       //!non-matched jet pt distribution
  80   TH2F                       *fHistPartvsDetecPt;         //!particle vs detector level jet pt
  81   TH1F                       *fHistMissedMCJets;          //!mc jets not measured
  82
  83  private:
  84   AliJetResponseMaker(const AliJetResponseMaker&);            // not implemented
  85   AliJetResponseMaker &operator=(const AliJetResponseMaker&); // not implemented
  86
  87   ClassDef(AliJetResponseMaker, 6) // Jet response matrix producing task
  88 };
  89 #endif