Add jet mass to tagged jet correlation

[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / AliJetContainer.h

diff --git a/PWGJE/EMCALJetTasks/AliJetContainer.h b/PWGJE/EMCALJetTasks/AliJetContainer.h
index 037c6ee416bd0baa3b37db5ceb85c86d7a9485de..ec7b7905f8bf31101361db0de3455bf99bb437b0 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliJetContainer.h
+++ b/PWGJE/EMCALJetTasks/AliJetContainer.h
@@ -35,24 +35,24 @@ class AliJetContainer : public AliEmcalContainer {
   void SetEMCALGeometry(AliEMCALGeometry *p) {fGeom = p;}
   void LoadRho(AliVEvent *event);
   void LoadLocalRho(AliVEvent *event);
+  void LoadRhoMass(AliVEvent *event);
 
   void                        SetJetAcceptanceType(JetAcceptanceType type)         { fJetAcceptanceType          = type ; }
+  void                        PrintCuts();
   void                        ResetCuts();
   void                        SetJetEtaPhiEMCAL() ;
   void                        SetJetEtaPhiTPC()   ;
   void                        SetRunNumber(Int_t r)                                { fRunNumber = r;                      }
   void                        SetJetEtaLimits(Float_t min, Float_t max)            { fJetMinEta = min, fJetMaxEta = max ; }
   void                        SetJetPhiLimits(Float_t min, Float_t max)            { fJetMinPhi = min, fJetMaxPhi = max ; }
+  void                        SetJetPtCut(Float_t cut)                             { fJetPtCut       = cut              ; }
+  void                        SetJetRadius(Float_t r)                              { fJetRadius      = r                ; } 
   void                        SetJetAreaCut(Float_t cut)                           { fJetAreaCut     = cut              ; }
   void                        SetPercAreaCut(Float_t p)                            { if(fJetRadius==0.) AliWarning("JetRadius not set. Area cut will be 0"); 
                                                                                      fJetAreaCut = p*TMath::Pi()*fJetRadius*fJetRadius; }
+  void                        SetAreaEmcCut(Double_t a = 0.99)                     { fAreaEmcCut     = a                ; }
   void                        SetZLeadingCut(Float_t zemc, Float_t zch)            { fZLeadingEmcCut = zemc; fZLeadingChCut = zch ; }
   void                        SetNEFCut(Float_t min = 0., Float_t max = 1.)        { fNEFMinCut = min; fNEFMaxCut = max;  }
-  void                        SetAreaEmcCut(Double_t a = 0.99)                     { fAreaEmcCut     = a                ; }
-  void                        SetJetPtCut(Float_t cut)                             { fJetPtCut       = cut              ; }
-  void                        SetJetRadius(Float_t r)                              { fJetRadius      = r                ; } 
-  virtual void                SetRhoName(const char *n)                            { fRhoName        = n                ; }
-  virtual void                SetLocalRhoName(const char *n)                       { fLocalRhoName   = n                ; }
   void                        SetFlavourCut(Int_t myflavour)                       { fFlavourSelection = myflavour;}
   void                        SetMaxClusterPt(Float_t b)                           { fMaxClusterPt   = b                ; }
   void                        SetMaxTrackPt(Float_t b)                             { fMaxTrackPt     = b                ; }
@@ -62,8 +62,14 @@ class AliJetContainer : public AliEmcalContainer {
   void                        SetLeadingHadronType(Int_t t)                        { fLeadingHadronType = t             ; }
   void                        SetJetBitMap(UInt_t m)                               { fJetBitMap      = m                ; }
   void                        SetJetTrigger(UInt_t t=AliVEvent::kEMCEJE)           { fJetTrigger     = t                ; }
+  void                        SetTagStatus(Int_t i)                                { fTagStatus      = i                ; }
+
+  virtual void                SetRhoName(const char *n)                            { fRhoName        = n                ; }
+  virtual void                SetLocalRhoName(const char *n)                       { fLocalRhoName   = n                ; }
+  virtual void                SetRhoMassName(const char *n)                        { fRhoMassName    = n                ; }
   void                        ConnectParticleContainer(AliParticleContainer *c)    { fParticleContainer = c             ; }
   void                        ConnectClusterContainer(AliClusterContainer *c)      { fClusterContainer  = c             ; }
+
   AliEmcalJet                *GetLeadingJet(const char* opt="")          ;
   AliEmcalJet                *GetJet(Int_t i)                       const;
   AliEmcalJet                *GetAcceptJet(Int_t i)                 const;
@@ -86,6 +92,9 @@ class AliJetContainer : public AliEmcalContainer {
   const TString&              GetRhoName()                          const    {return fRhoName;}
   AliLocalRhoParameter       *GetLocalRhoParameter()                const    {return fLocalRho;}
   const TString&              GetLocalRhoName()                     const    {return fLocalRhoName;}
+  AliRhoParameter            *GetRhoMassParameter()                          {return fRhoMass;}
+  Double_t                    GetRhoMassVal()                       const    {if (fRhoMass) return fRhoMass->GetVal(); else return 0;}
+  const TString&              GetRhoMassName()                      const    {return fRhoMassName;}
   Double_t                    GetJetPtCorr(Int_t i)                 const;
   Double_t                    GetJetPtCorrLocal(Int_t i)            const;
   Float_t                     GetJetRadius()                        const    {return fJetRadius;}
@@ -97,12 +106,14 @@ class AliJetContainer : public AliEmcalContainer {
   void                        SetArray(AliVEvent *event);
   AliParticleContainer       *GetParticleContainer()                         {return fParticleContainer;}
   AliClusterContainer        *GetClusterContainer()                          {return fClusterContainer;}
+  Double_t                    GetFractionSharedPt(AliEmcalJet *jet) const;
 
  protected:
   JetAcceptanceType           fJetAcceptanceType;    //  acceptance type
   Float_t                     fJetRadius;            //  jet radius
   TString                     fRhoName;              //  Name of rho object
   TString                     fLocalRhoName;         //  Name of local rho object
+  TString                     fRhoMassName;          //  Name of rho mass object
   Int_t                       fFlavourSelection;     //  selection on jet flavour
   Float_t                     fPtBiasJetTrack;       //  select jets with a minimum pt track
   Float_t                     fPtBiasJetClus;        //  select jets with a minimum pt cluster
@@ -123,11 +134,12 @@ class AliJetContainer : public AliEmcalContainer {
   Int_t                       fNLeadingJets;         //  how many jets are to be considered the leading jet(s)
   UInt_t                      fJetBitMap;            //  bit map of accepted jets
   UInt_t                      fJetTrigger;           //  jet trigger
+  Int_t                       fTagStatus;            //  jet tag status
   AliParticleContainer       *fParticleContainer;    //  particle container (jet constituents)
   AliClusterContainer        *fClusterContainer;     //  cluster container (jet constituents)
-
   AliRhoParameter            *fRho;                  //! event rho for these jets
   AliLocalRhoParameter       *fLocalRho;             //! event local rho for these jets
+  AliRhoParameter            *fRhoMass;              //! event rho mass for these jets
   AliEMCALGeometry           *fGeom;                 //! emcal geometry
   Int_t                       fRunNumber;            //! run number
 
@@ -135,7 +147,7 @@ class AliJetContainer : public AliEmcalContainer {
   AliJetContainer(const AliJetContainer& obj); // copy constructor
   AliJetContainer& operator=(const AliJetContainer& other); // assignment
 
-  ClassDef(AliJetContainer,6);
+  ClassDef(AliJetContainer,9);
 
 };