Handle expressions with negations (Diego)

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

diff --git a/PWGJE/AliAnalysisTaskJetProperties.h b/PWGJE/AliAnalysisTaskJetProperties.h
index 97ec372fc00b027a5aad4c4f622940760da2da2c..0f5f00cf6dcb15e745a53e21699e9f603f01c155 100644 (file)
--- a/PWGJE/AliAnalysisTaskJetProperties.h
+++ b/PWGJE/AliAnalysisTaskJetProperties.h
@@ -44,6 +44,7 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   {fMaxVertexZ = VtxZ; fNContributors = nContributors;}
   virtual void   SetTrackCuts(Float_t trackPt = 0.15, Float_t trackEtaMin = -0.9, Float_t trackEtaMax = 0.9)
   {fTrackPtCut = trackPt; fTrackEtaMin = trackEtaMin; fTrackEtaMax = trackEtaMax;}
+  virtual void   SetPileupRejection(Bool_t IsPileupReject){fRejectPileup = IsPileupReject;}
   virtual void   SetJetCuts(Float_t jetPt = 5., Float_t jetEtaMin = -0.5, Float_t jetEtaMax = 0.5) 
   {fJetPtCut = jetPt; fJetEtaMin = jetEtaMin; fJetEtaMax = jetEtaMax;}
   virtual void   SetJetRejectType(Int_t i){fJetRejectType = i;}
@@ -53,13 +54,15 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   
   enum {kTrackUndef=0, kTrackAOD, kTrackKine,kTrackAODMC};//for track selection
   enum {kNoReject=0, kReject1Track};//for jet rejection
-  
+  enum {kRejectPileup=1};//for pileup rejection
  protected:
   Int_t           GetListOfJetTracks(TList* l, const AliAODJet* j);
   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);
   
@@ -73,6 +76,7 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   TString fBranchJets;          // branch name for reconstructed jets
   Int_t   fTrackType;           // type of generated tracks
   Int_t   fJetRejectType;       // type of jets rejected
+  Bool_t  fRejectPileup;        // pileup rejection
   Bool_t  fUseAODInputJets;     // take jets from in/output - only relevant if AOD event both in input AND output and we want to use output
   UInt_t  fFilterMask;         // filter bit for selected tracks
   Bool_t  fUsePhysicsSelection; // switch for event selection
@@ -111,7 +115,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
@@ -151,8 +156,18 @@ class AliAnalysisTaskJetProperties : public AliAnalysisTaskSE {
   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);