PWGJE/UserTasks/AliAnalysisTaskCheckSingleTrackJetRejection.h

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