Merge branch 'master' of https://git.cern.ch/reps/AliRoot

[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / UserTasks / AliAnalysisTaskJetMatching.h

diff --git a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskJetMatching.h b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskJetMatching.h
index 5962afe36f51f7d21254bfbfdb49760417213c33..07bc4713cf6344c8e869f78de8c0ad595a025987 100644 (file)
--- a/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskJetMatching.h
+++ b/PWGJE/EMCALJetTasks/UserTasks/AliAnalysisTaskJetMatching.h
@@ -1,7 +1,7 @@
 #ifndef AliAnalysisTaskJetMatching_H
 #define AliAnalysisTaskJetMatching_H
 
-#include <AliAnalysisTaskEmcalJetDev.h>
+#include <AliAnalysisTaskEmcalJet.h>
 #include <AliEmcalJet.h>
 #include <AliVTrack.h>
 #include <TClonesArray.h>
@@ -16,11 +16,11 @@ class AliRhoParameter;
 class AliLocalRhoParameter;
 class TClonesArray;
 
-class AliAnalysisTaskJetMatching : public AliAnalysisTaskEmcalJetDev
+class AliAnalysisTaskJetMatching : public AliAnalysisTaskEmcalJet
 {
     public:
         // enumerators
-        enum matchingSceme      {kGeoEtaPhi, kGeoR};
+        enum matchingSceme      {kGeoEtaPhi, kGeoR, kDiJet};
         enum sourceBKG          {kNoSourceBKG, kSourceRho, kSourceLocalRho};
         enum targetBKG          {kNoTargetBKG, kTargetRho, kTargetLocalRho};
         // constructors, destructor
@@ -29,6 +29,7 @@ class AliAnalysisTaskJetMatching : public AliAnalysisTaskEmcalJetDev
         virtual                 ~AliAnalysisTaskJetMatching();
         // setting up the task and technical aspects
         void                    ExecOnce();
+        virtual Bool_t          Notify();
         virtual void            UserCreateOutputObjects();
         TH1F*                   BookTH1F(const char* name, const char* x, Int_t bins, Double_t min, Double_t max, Bool_t append = kTRUE);
         TH2F*                   BookTH2F(const char* name, const char* x, const char* y, Int_t binsx, Double_t minx, Double_t maxx, Int_t binsy, Double_t miny, Double_t maxy, Bool_t append = kTRUE);
@@ -74,24 +75,29 @@ class AliAnalysisTaskJetMatching : public AliAnalysisTaskEmcalJetDev
         void                    SetMatchEta(Float_t f)                          {fMatchEta = f;}
         void                    SetMatchPhi(Float_t f)                          {fMatchPhi = f;}
         void                    SetMatchR(Float_t f)                            {fMatchR = f;}
+        void                    SetDoDetectorResponse(Bool_t d)                 {fDoDetectorResponse = d;}
         // methods
         void                    DoGeometricMatchingEtaPhi();
         void                    DoGeometricMatchingR();
+        void                    DoDiJetMatching();
         void                    DoConstituentMatching();
         void                    GetBijection();
         void                    PostMatchedJets();
         // fill histogrmas
         void                    FillAnalysisSummaryHistogram() const;
-        void                    FillMatchedJetHistograms() const;
+        void                    FillMatchedJetHistograms();
         // setters - analysis details
         /* inline */    Bool_t PassesCuts(const AliVTrack* track) const {
             return (!track || TMath::Abs(track->Eta()) > 0.9 || track->Phi() < 0 || track->Phi() > TMath::TwoPi()) ? kFALSE : kTRUE; }
-        /* inline */    Bool_t PassesCuts(AliEmcalJet* jet) const { return (jet) ? kTRUE : kFALSE; }
+        /* inline */    Bool_t PassesCuts(AliEmcalJet* jet, Int_t c)  { if(!fUseEmcalBaseJetCuts) return (jet) ? kTRUE : kFALSE; 
+            else return AcceptJet(jet, c);}
         /* inline */    void ClearMatchedJetsCache() {
             for (Int_t i(0); i < fNoMatchedJets; i++) {
                 fMatchedJetContainer[i][0] = 0x0; fMatchedJetContainer[i][1] = 0x0; }
             fNoMatchedJets = 0;
         }
+        AliEmcalJet*            GetLeadingJet(TClonesArray* source, Int_t &leadingJetIndex, Double_t etaMin = -.7, Double_t etaMax = .7);
+        AliEmcalJet*            GetSubLeadingJet(TClonesArray* source, Int_t leadingJetIndex, Int_t &subLeadingJetIndex);
         void                    PrintInfo() const;
         virtual void            Terminate(Option_t* option);
 
@@ -119,8 +125,11 @@ class AliAnalysisTaskJetMatching : public AliAnalysisTaskEmcalJetDev
         TH1F*                   fHistUnsortedCorrelation;       //! dphi correlation of unsorted jets
         TH1F*                   fHistMatchedCorrelation;        //! dphi correlation of matched jets
         TH1F*                   fHistSourceJetPt;       //! pt of source jets
+        TH1F*                   fHistMatchedSourceJetPt;//! pt of matched source jets
         TH1F*                   fHistTargetJetPt;       //! pt of target jets
         TH1F*                   fHistMatchedJetPt;      //! pt of matched jets
+        TH2F*                   fHistSourceMatchedJetPt;//! pt or source vs matched jets
+        TH2F*                   fHistDetectorResponseProb;      //! det prob
         TH1F*                   fHistNoConstSourceJet;  //! number of constituents of source jet
         TH1F*                   fHistNoConstTargetJet;  //! number of constituents of target jet
         TH1F*                   fHistNoConstMatchJet;   //! number of constituents of matched jets
@@ -128,24 +137,34 @@ class AliAnalysisTaskJetMatching : public AliAnalysisTaskEmcalJetDev
         TProfile*               fProfFracPtJets;        //! sum pt fraction for matched jets
         TProfile*               fProfFracNoMatched;     //! no of constituents fraction found / jet
         TProfile*               fProfFracNoJets;        //! no of consstituents fraction jet / jet
+        TH2F*                   fHistDiJet;             //! matched dijet eta, phi
+        TH2F*                   fHistDiJetLeadingJet;   //! leading jet (for dijet) eta, phi
+        TH1F*                   fHistDiJetDPhi;         //! dijet dphi
+        TH1F*                   fHistDiJetDPt;          //! dijet dpt
         TH1F*                   fHistAnalysisSummary;   //! flags
         TProfile*               fProfQAMatched;         //! QA spreads of matched jets
         TProfile*               fProfQA;                //! QA spreads of source and target jets
         Int_t                   fNoMatchedJets;         //! number of matched jets
-        AliEmcalJet*            fMatchedJetContainer[100][2];   //! pointers to matched jets
+        AliEmcalJet*            fMatchedJetContainer[200][2];   //! pointers to matched jets
         // geometric matching parameters
         Float_t                 fMatchEta;              // max eta distance between centers of matched jets
         Float_t                 fMatchPhi;              // max phi distance between centers of matched jets
         Float_t                 fMatchR;                // max distance between matched jets
+        Bool_t                  fDoDetectorResponse;    // get detector response matrix (not normalized)
         Bool_t                  fMatchConstituents;     // match constituents
         Float_t                 fMinFracRecoveredConstituents;  // minimium fraction of constituents that needs to be found
         Float_t                 fMinFracRecoveredConstituentPt; // minimium fraction of constituent pt that needs to be recovered
         Bool_t                  fGetBijection;          // get bijection of source and matched jets
+        // pythia ntrials xsec bookkeeping
+        TH1F*                   fh1Trials;              //! sum of trails (in Notify)
+        TH1F*                   fh1AvgTrials;           //! average sum of trials (in Run)
+        TProfile*               fh1Xsec;                //! pythia cross section and trials
+        Float_t                 fAvgTrials;             // average number of trails per event
 
         AliAnalysisTaskJetMatching(const AliAnalysisTaskJetMatching&);                  // not implemented
         AliAnalysisTaskJetMatching& operator=(const AliAnalysisTaskJetMatching&);       // not implemented
 
-        ClassDef(AliAnalysisTaskJetMatching, 4);
+        ClassDef(AliAnalysisTaskJetMatching, 6);
 };
 
 #endif