1 #ifndef AliAnalysisTaskCheckSingleTrackJetRejection_cxx
2 #define AliAnalysisTaskCheckSingleTrackJetRejection_cxx
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
11 class AliGenPythiaEventHeader;
22 #include "AliAnalysisTaskSE.h"
23 #include "THnSparse.h" // cannot forward declare ThnSparseF
28 class AliAnalysisTaskCheckSingleTrackJetRejection : public AliAnalysisTaskSE
31 AliAnalysisTaskCheckSingleTrackJetRejection();
32 AliAnalysisTaskCheckSingleTrackJetRejection(const char *name);
33 virtual ~AliAnalysisTaskCheckSingleTrackJetRejection() {;}
35 // Implementation of interface methods
36 virtual void UserCreateOutputObjects();
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);
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};
68 AliAnalysisTaskCheckSingleTrackJetRejection(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented
69 AliAnalysisTaskCheckSingleTrackJetRejection& operator=(const AliAnalysisTaskCheckSingleTrackJetRejection &det); // not implemented
71 Bool_t JetSelected(AliAODJet *jet);
73 Bool_t fUseAODInput; // read jets from input AOD
74 Bool_t fFillAOD; // option to fill AOD branch
76 TString fJetBranch; // jet branch to read
78 AliAODEvent *fAODIn; // AOD event
79 AliAODEvent *fAODOut; // AOD event
80 AliAODExtension *fAODExtension;
89 TList *fHistList; // Output list
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
99 TH1F *fH1jetMCAKT04_pt [6];
100 TH2F *fH2jetMCAKT04_Jetpt_maxpt;
101 TH1F *fH1jetAKT04_pt [6];
102 TH2F *fH2jetAKT04_Jetpt_maxpt ;
104 TH2F *fH2jetMCAKT04_Eratio [6];
105 TH1F *fH1jetMCAKT04_ptmatch [6];
109 ClassDef(AliAnalysisTaskCheckSingleTrackJetRejection, 13); // Analysis task for standard dijet analysis