jet cut update for matching task

[u/mrichter/AliRoot.git] / PWGJE / AliAnalysisTaskJetProperties.h

diff --git a/PWGJE/AliAnalysisTaskJetProperties.h b/PWGJE/AliAnalysisTaskJetProperties.h
index 2d44d9d43f168e7d07265d244820fa1006b6df1b..40d97bd235335b1554753e4be7738b9b937ba868 100644 (file)
--- a/PWGJE/AliAnalysisTaskJetProperties.h
+++ b/PWGJE/AliAnalysisTaskJetProperties.h
@@ -49,6 +49,7 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   virtual void   SetJetRejectType(Int_t i){fJetRejectType = i;}
   virtual void   SetFilterMask(UInt_t i) {fFilterMask = i;}
   virtual void   UsePhysicsSelection(Bool_t b) {fUsePhysicsSelection = b;}
+  virtual void   SetJetRadius(Float_t r){fJetRadius = r;}
   
   enum {kTrackUndef=0, kTrackAOD, kTrackKine,kTrackAODMC};//for track selection
   enum {kNoReject=0, kReject1Track};//for jet rejection
@@ -58,6 +59,11 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   Int_t           GetListOfJets(TList* list);
   void     FillJetProperties(TList *jetlist);
   void     FillJetShape(TList *jetlist);
+  void     FillJetShapeUE(TList *jetlist);
+  void     FillFFCorr(TList *jetlist);
+
+  void     GetTracksTiltedwrpJetAxis(Float_t alpha, TList* inputlist, TList* outputlist, const AliAODJet* jet, Double_t radius,Double_t& sumPt);
+  Int_t    GetListOfTracks(TList *list, Int_t type);
   
   AliESDEvent*     fESD;          // ESD event
   AliAODEvent*     fAOD;          // AOD event
@@ -83,9 +89,12 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   Float_t fJetEtaMin;           // jet eta cut
   Float_t fJetEtaMax;           // jet eta cut
   Float_t fAvgTrials;           // average number of trials per event
+  Float_t fJetRadius;
   
   TList            *fJetList;                //! List of jets
-  TList            *fTrackList;              //! List of tracks in a jet
+  TList            *fTrackList;              //! List of selected tracks
+  TList            *fTrackListUE;            //! List of tracks in jet cone UE
+  TList            *fTrackListJet;           //! List of tracks in a jet
   TList            *fCommonHistList;         //! List of common histos
   TH1F      *fh1EvtSelection;         //! event cuts 
   TH1F     *fh1VertexNContributors;  //! NContributors to prim vertex
@@ -104,7 +113,8 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   TH2F*     fh2EtaTrack;              //!track eta distribution
   TH2F*     fh2PhiTrack;              //!track phi distribution
   TH2F*     fh2PtTrack;               //!track pt distribution
-  TH2F*     fh2FF;                    //!fragmentation function
+  TH2F*     fh2FF;                    //!fragmentation function (z)
+  TH2F*     fh2Ksi;                   //!fragmentation function (Ksi)
   TH2F*     fh2DelEta;                //!delta eta distribution
   TH2F*     fh2DelPhi;                //!delta phi distribution
   TH2F*     fh2DelR;                  //!delta R distribution
@@ -122,12 +132,40 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   TH3F*     fh3PtDelRPtSum;           //!Pt sum vs R
   TProfile* fProDiffJetShape;         //!Diff jet shape pT=20-100
   TProfile* fProIntJetShape;          //!Int jet shape pT=20-100
-  TProfile* fProDelRNchSum[5];        //!Nch sum vs R
-  TProfile* fProDelRPtSum[5];         //!Pt sum vs R
-  TProfile* fProDiffJetShapeA[5];     //!Diff jet shape pT
-  TProfile* fProIntJetShapeA[5];      //!Int jet shape pT
+  TProfile* fProDelRNchSum[13];        //!Nch sum vs R
+  TProfile* fProDelRPtSum[13];         //!Pt sum vs R
+  TProfile* fProDiffJetShapeA[13];     //!Diff jet shape pT
+  TProfile* fProIntJetShapeA[13];      //!Int jet shape pT
 
+  TH1F*     fh1PtSumInJetConeUE;        //!pt sum of UE in cone R
+  TH2F*     fh2NtracksLeadingJetUE;     //!number of tracks in jet (UE)
+  TProfile* fProNtracksLeadingJetUE;    //!number of tracks in jet (UE)
+  TH2F*     fh2DelR80pcNchUE;           //!R containing 80% of Nch vs jet pt (UE)
+  TProfile* fProDelR80pcNchUE;          //!R containing 80% of Nch vs jet pt (UE)
+  TH2F*     fh2DelR80pcPtUE;            //!R containing 80% of pT vs jet pt (UE)
+  TProfile* fProDelR80pcPtUE;           //!R containing 80% of pT vs jet pt (UE)
+  TH2F*     fh2AreaChUE;                //!charged jet area vs jet pT (UE)
+  TProfile* fProAreaChUE;               //!charged jet area vs jet pT (UE)
+  TH3F*     fh3PtDelRNchSumUE;          //!Nch sum vs R (UE)
+  TH3F*     fh3PtDelRPtSumUE;           //!Pt sum vs R (UE)
+  TProfile* fProDiffJetShapeUE;         //!Diff jet shape pT=20-100 (UE)
+  TProfile* fProIntJetShapeUE;          //!Int jet shape pT=20-100 (UE)
+  TProfile* fProDelRNchSumUE[13];        //!Nch sum vs R (UE)
+  TProfile* fProDelRPtSumUE[13];         //!Pt sum vs R (UE)
+  TProfile* fProDiffJetShapeAUE[13];     //!Diff jet shape pT (UE)
+  TProfile* fProIntJetShapeAUE[13];      //!Int jet shape pT (UE)
 
+  TH1F* fh1CorrJetPt;                    //!jet pt in FillCorr
+  TH2F* fh2CorrPtTrack1;                 //!trk pt in FillCorr
+  TH2F* fh2CorrFF1;                      //!FF in FillCorr
+  TH2F* fh2CorrKsi1;                     //!Ksi in FillCorr
+  TH2F* fh2CorrjT1;                      //!jT in FillCorr
+  TH1F* fh1JetPtvsTrkSum;                //!QA plots to check jetPt-trkPtsum
+  TH2F* fh2CorrPt1Pt2[6];                //!pt1-pt2 in FillCorr
+  TH2F* fh2CorrZ1Z2[6];                  //!z1-z2 in FillCorr
+  TH2F* fh2CorrKsi1Ksi2[6];              //!ksi1-ksi2 in FillCorr
+  TH2F* fh2CorrjT1jT2[6];                //!jt1-jt2 in FillCorr
+  
   AliAnalysisTaskJetProperties(const AliAnalysisTaskJetProperties&);// not implemented
   AliAnalysisTaskJetProperties& operator=(const AliAnalysisTaskJetProperties&);// not implemented
   ClassDef(AliAnalysisTaskJetProperties, 2);