PWGJE/UserTasks/AliAnalysisTaskCheckSingleTrackJetRejection.h

   1 #ifndef AliAnalysisTaskCheckSingleTrackJetRejection_cxx
   2 #define AliAnalysisTaskCheckSingleTrackJetRejection_cxx
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 class AliJetHeader;
   8 class AliESDEvent;
   9 class AliAODEvent;
  10 class AliAODJet;
  11 class AliGenPythiaEventHeader;
  12 class AliCFManager;
  13
  14 class TList;
  15 class TChain;
  16 class TH2F;
  17 class TH1F;
  18 class TH3F;
  19 class TProfile;
  20
  21
  22 #include "AliAnalysisTaskSE.h"
  23 #include  "THnSparse.h" // cannot forward declare ThnSparseF
  24 #include <TList.h>
  25 #include <TH1F.h>
  26 #include <TH2F.h>
  27
  28 class AliAnalysisTaskCheckSingleTrackJetRejection : public AliAnalysisTaskSE
  29 {
  30         public:
  31                 AliAnalysisTaskCheckSingleTrackJetRejection();
  32                 AliAnalysisTaskCheckSingleTrackJetRejection(const char *name);
  33                 virtual ~AliAnalysisTaskCheckSingleTrackJetRejection() {;}
  34
  35                 // Implementation of interface methods
  36                 virtual void    UserCreateOutputObjects();
  37                 virtual void    Init();
  38                 virtual Bool_t  Notify();
  39                 virtual void    UserExec(Option_t *option);
  40                 virtual void    Terminate(Option_t *);
  41                 virtual void    SetDebug(Int_t debug = 0) {fDebug = debug;}
  42                 virtual void    SetAlgorithm(TString jf="ANTIKT"){JFAlg=jf;}
  43                 virtual void    SetRadius(Float_t radius=0.4){Radius=radius;}
  44                 virtual void    SetFilterMask(UInt_t filter=256){Filtermask=filter;}
  45                 virtual void    SetBackSubMode(Int_t backM=0){BackM=backM;}
  46                 virtual void    SetTrackPtCut(Float_t tPtcut=0){TrackPtcut=tPtcut;}
  47                 virtual void    SetSkipCone(Int_t skipCone=0){SkipCone=skipCone;}
  48                 virtual void    SetMC(Bool_t ismc=kFALSE){IsMC=ismc;}
  49     virtual void    FinishTaskOutput();
  50     virtual Bool_t  PythiaInfoFromFile(const char* currFile,Float_t &fXsec,Float_t &fTrials);
  51     //virtual Float_t GetTotalEvents(const char* currFile);
  52     virtual Double_t DeltaPhi(Double_t phi1,Double_t phi2);
  53     virtual Double_t DeltaR(Double_t phi1,Double_t phi2,Double_t eta1,Double_t eta2);
  54
  55
  56                 enum {kNPTBINS=10};
  57
  58                 // 0 all jets
  59                 // 1 all jet in eta window
  60                 // 2 all jets with partner
  61                 // 3 all jets in eta window with partner
  62                 // 4 all jets with partner in eta window
  63                 enum {kStep0 = 0, kStep1, kStep2, kStep3, kStep4,kMaxStep};
  64
  65
  66
  67         private:
  68                 AliAnalysisTaskCheckSingleTrackJetRejection(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented
  69                 AliAnalysisTaskCheckSingleTrackJetRejection& operator=(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented
  70
  71     Bool_t  JetSelected(AliAODJet *jet);
  72
  73                 Bool_t        fUseAODInput; // read jets from input AOD
  74                 Bool_t        fFillAOD;     // option to fill AOD branch
  75                 TString       fNonStdFile;
  76                 TString       fJetBranch;   // jet branch to read
  77
  78                 AliAODEvent     *fAODIn;      // AOD event
  79                 AliAODEvent     *fAODOut;      // AOD event
  80                 AliAODExtension *fAODExtension;
  81                 TString         JFAlg;
  82                 Float_t         Radius;
  83                 UInt_t          Filtermask;
  84                 Int_t           BackM;
  85                 Float_t         TrackPtcut;
  86                 Int_t           SkipCone;
  87                 Bool_t          IsMC;
  88
  89                 TList        *fHistList; // Output list
  90                 Float_t       fxsec;
  91                 Float_t       ftrial;
  92                 Float_t       fJetRecEtaWindow;       // eta window for rec jets
  93                 Float_t       fMinJetPt;              // limits the jet p_T in addition to what already is done in the jet finder, this is important for jet matching for JF with lo threshold
  94
  95                 TProfile     *fH1Xsec;
  96                 TH1F         *fH1Trials;
  97                 TH1F         *fH1Events;
  98
  99                 TH1F         *fH1jetMCAKT04_pt         [6];
 100                 TH2F         *fH2jetMCAKT04_Jetpt_maxpt;
 101                 TH1F         *fH1jetAKT04_pt           [6];
 102                 TH2F         *fH2jetAKT04_Jetpt_maxpt  ;
 103
 104                 TH2F         *fH2jetMCAKT04_Eratio       [6];
 105                 TH1F         *fH1jetMCAKT04_ptmatch        [6];
 106
 107
 108
 109                 ClassDef(AliAnalysisTaskCheckSingleTrackJetRejection, 13); // Analysis task for standard dijet analysis
 110 };
 111
 112 #endif