-#ifndef AliAnalysisTaskCheckSingleTrackJetRejection_cxx\r
-#define AliAnalysisTaskCheckSingleTrackJetRejection_cxx\r
-\r
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
- * See cxx source for full Copyright notice */\r
-\r
-class AliJetHeader;\r
-class AliESDEvent;\r
-class AliAODEvent;\r
-class AliAODJet;\r
-class AliGenPythiaEventHeader;\r
-class AliCFManager;\r
-\r
-class TList;\r
-class TChain;\r
-class TH2F;\r
-class TH1F;\r
-class TH3F;\r
-class TProfile;\r
-\r
-\r
-#include "AliAnalysisTaskSE.h"\r
-#include "THnSparse.h" // cannot forward declare ThnSparseF \r
-#include <TList.h>\r
-#include <TH1F.h>\r
-#include <TH2F.h>\r
-\r
-class AliAnalysisTaskCheckSingleTrackJetRejection : public AliAnalysisTaskSE \r
-{\r
- public:\r
- AliAnalysisTaskCheckSingleTrackJetRejection();\r
- AliAnalysisTaskCheckSingleTrackJetRejection(const char *name);\r
- virtual ~AliAnalysisTaskCheckSingleTrackJetRejection() {;}\r
-\r
- // Implementation of interface methods\r
- virtual void UserCreateOutputObjects();\r
- virtual void Init();\r
- virtual Bool_t Notify();\r
- virtual void UserExec(Option_t *option);\r
- virtual void Terminate(Option_t *);\r
- virtual void SetDebug(Int_t debug = 0) {fDebug = debug;}\r
- virtual void SetAlgorithm(TString jf="ANTIKT"){JFAlg=jf;}\r
- virtual void SetRadius(Float_t radius=0.4){Radius=radius;}\r
- virtual void SetFilterMask(UInt_t filter=256){Filtermask=filter;}\r
- virtual void SetBackSubMode(Int_t backM=0){BackM=backM;}\r
- virtual void SetTrackPtCut(Float_t tPtcut=0){TrackPtcut=tPtcut;}\r
- virtual void SetSkipCone(Int_t skipCone=0){SkipCone=skipCone;}\r
- virtual void SetMC(Bool_t ismc=kFALSE){IsMC=ismc;}\r
- virtual void FinishTaskOutput();\r
- virtual Bool_t PythiaInfoFromFile(const char* currFile,Float_t &fXsec,Float_t &fTrials);\r
- //virtual Float_t GetTotalEvents(const char* currFile);\r
- virtual Double_t DeltaPhi(Double_t phi1,Double_t phi2);\r
- virtual Double_t DeltaR(Double_t phi1,Double_t phi2,Double_t eta1,Double_t eta2);\r
-\r
-\r
- enum {kNPTBINS=10};\r
-\r
- // 0 all jets\r
- // 1 all jet in eta window\r
- // 2 all jets with partner\r
- // 3 all jets in eta window with partner\r
- // 4 all jets with partner in eta window\r
- enum {kStep0 = 0, kStep1, kStep2, kStep3, kStep4,kMaxStep};\r
-\r
-\r
-\r
- private:\r
- AliAnalysisTaskCheckSingleTrackJetRejection(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented\r
- AliAnalysisTaskCheckSingleTrackJetRejection& operator=(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented\r
-\r
- Bool_t JetSelected(AliAODJet *jet);\r
-\r
- Bool_t fUseAODInput; // read jets from input AOD\r
- Bool_t fFillAOD; // option to fill AOD branch\r
- TString fNonStdFile;\r
- TString fJetBranch; // jet branch to read\r
-\r
- AliAODEvent *fAODIn; // AOD event\r
- AliAODEvent *fAODOut; // AOD event\r
- AliAODExtension *fAODExtension;\r
- TString JFAlg;\r
- Float_t Radius;\r
- UInt_t Filtermask;\r
- Int_t BackM;\r
- Float_t TrackPtcut;\r
- Int_t SkipCone;\r
- Bool_t IsMC;\r
-\r
- TList *fHistList; // Output list\r
- Float_t fxsec;\r
- Float_t ftrial;\r
- Float_t fJetRecEtaWindow; // eta window for rec jets\r
- 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
-\r
- TProfile *fH1Xsec;\r
- TH1F *fH1Trials;\r
- TH1F *fH1Events; \r
-\r
- TH1F *fH1jetMCAKT04_pt [6];\r
- TH2F *fH2jetMCAKT04_Jetpt_maxpt;\r
- TH1F *fH1jetAKT04_pt [6]; \r
- TH2F *fH2jetAKT04_Jetpt_maxpt ;\r
-\r
- TH2F *fH2jetMCAKT04_Eratio [6];\r
- TH1F *fH1jetMCAKT04_match [6];\r
-\r
-\r
-\r
- ClassDef(AliAnalysisTaskCheckSingleTrackJetRejection, 13); // Analysis task for standard dijet analysis\r
-};\r
-\r
-#endif\r
+#ifndef AliAnalysisTaskCheckSingleTrackJetRejection_cxx
+#define AliAnalysisTaskCheckSingleTrackJetRejection_cxx
+
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice */
+
+class AliJetHeader;
+class AliESDEvent;
+class AliAODEvent;
+class AliAODJet;
+class AliGenPythiaEventHeader;
+class AliCFManager;
+
+class TList;
+class TChain;
+class TH2F;
+class TH1F;
+class TH3F;
+class TProfile;
+
+
+#include "AliAnalysisTaskSE.h"
+#include "THnSparse.h" // cannot forward declare ThnSparseF
+#include <TList.h>
+#include <TH1F.h>
+#include <TH2F.h>
+
+class AliAnalysisTaskCheckSingleTrackJetRejection : public AliAnalysisTaskSE
+{
+ public:
+ AliAnalysisTaskCheckSingleTrackJetRejection();
+ AliAnalysisTaskCheckSingleTrackJetRejection(const char *name);
+ virtual ~AliAnalysisTaskCheckSingleTrackJetRejection() {;}
+
+ // Implementation of interface methods
+ virtual void UserCreateOutputObjects();
+ virtual void Init();
+ virtual Bool_t Notify();
+ virtual void UserExec(Option_t *option);
+ virtual void Terminate(Option_t *);
+ virtual void SetDebug(Int_t debug = 0) {fDebug = debug;}
+ virtual void SetAlgorithm(TString jf="ANTIKT"){JFAlg=jf;}
+ virtual void SetRadius(Float_t radius=0.4){Radius=radius;}
+ virtual void SetFilterMask(UInt_t filter=256){Filtermask=filter;}
+ virtual void SetBackSubMode(Int_t backM=0){BackM=backM;}
+ virtual void SetTrackPtCut(Float_t tPtcut=0){TrackPtcut=tPtcut;}
+ virtual void SetSkipCone(Int_t skipCone=0){SkipCone=skipCone;}
+ virtual void SetMC(Bool_t ismc=kFALSE){IsMC=ismc;}
+ virtual void FinishTaskOutput();
+ virtual Bool_t PythiaInfoFromFile(const char* currFile,Float_t &fXsec,Float_t &fTrials);
+ //virtual Float_t GetTotalEvents(const char* currFile);
+ virtual Double_t DeltaPhi(Double_t phi1,Double_t phi2);
+ virtual Double_t DeltaR(Double_t phi1,Double_t phi2,Double_t eta1,Double_t eta2);
+
+
+ enum {kNPTBINS=10};
+
+ // 0 all jets
+ // 1 all jet in eta window
+ // 2 all jets with partner
+ // 3 all jets in eta window with partner
+ // 4 all jets with partner in eta window
+ enum {kStep0 = 0, kStep1, kStep2, kStep3, kStep4,kMaxStep};
+
+
+
+ private:
+ AliAnalysisTaskCheckSingleTrackJetRejection(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented
+ AliAnalysisTaskCheckSingleTrackJetRejection& operator=(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented
+
+ Bool_t JetSelected(AliAODJet *jet);
+
+ Bool_t fUseAODInput; // read jets from input AOD
+ Bool_t fFillAOD; // option to fill AOD branch
+ TString fNonStdFile;
+ TString fJetBranch; // jet branch to read
+
+ AliAODEvent *fAODIn; // AOD event
+ AliAODEvent *fAODOut; // AOD event
+ AliAODExtension *fAODExtension;
+ TString JFAlg;
+ Float_t Radius;
+ UInt_t Filtermask;
+ Int_t BackM;
+ Float_t TrackPtcut;
+ Int_t SkipCone;
+ Bool_t IsMC;
+
+ TList *fHistList; // Output list
+ Float_t fxsec;
+ Float_t ftrial;
+ Float_t fJetRecEtaWindow; // eta window for rec jets
+ 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
+
+ TProfile *fH1Xsec;
+ TH1F *fH1Trials;
+ TH1F *fH1Events;
+
+ TH1F *fH1jetMCAKT04_pt [6];
+ TH2F *fH2jetMCAKT04_Jetpt_maxpt;
+ TH1F *fH1jetAKT04_pt [6];
+ TH2F *fH2jetAKT04_Jetpt_maxpt ;
+
+ TH2F *fH2jetMCAKT04_Eratio [6];
+ TH1F *fH1jetMCAKT04_ptmatch [6];
+
+
+
+ ClassDef(AliAnalysisTaskCheckSingleTrackJetRejection, 13); // Analysis task for standard dijet analysis
+};
+
+#endif