JETAN/AliSISConeJetHeader.h

   1 #ifndef ALISISCONEJETHEADER_H
   2 #define ALISISCONEJETHEADER_H
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 //---------------------------------------------------------------------
   8 // SISCone (FastJet v2.3.4) finder algorithm interface
   9 // Finder Header Class
  10 // Author: swensy.jangal@ires.in2p3.fr
  11 //---------------------------------------------------------------------
  12
  13
  14 #include "fastjet/AreaDefinition.hh"
  15 #include "fastjet/ClusterSequenceArea.hh"
  16 #include "fastjet/JetDefinition.hh"
  17 #include "AliJetHeader.h"
  18
  19 class AliSISConeJetHeader : public AliJetHeader
  20 {
  21  public:
  22
  23   AliSISConeJetHeader();
  24   virtual ~AliSISConeJetHeader() { }
  25
  26   // Getters
  27   Bool_t                       GetBGMode()                     const {return fBGMode;}
  28
  29   Int_t                        GetActiveAreaRepeats()          const {return fActiveAreaRepeats;}
  30   Int_t                        GetAreaTypeNumber()             const {return fAreaTypeNumber;}
  31   Int_t                        GetBGAlgorithm()                const {return fBGAlgo;}
  32   Int_t                        GetNPassMax()                   const {return fNPassMax;}
  33   Int_t                        GetNumberOfJetsToErase()        const {return fNHardJets;}
  34   Int_t                        GetSplitMergeScale()            const {return fSplitMergeScaleNumber;}
  35
  36   Double_t                     GetGhostEtaMax()                const {return fGhostEtaMax;}
  37   Double_t                     GetGhostArea()                  const {return fGhostArea;}
  38   Double_t                     GetEffectiveRFact()             const {return fEffectiveRFact;}
  39
  40   Double_t                     GetRapMax()                     const {return fRapMax;}
  41   Double_t                     GetRapMin()                     const {return fRapMin;}
  42   Double_t                     GetPhiMax()                     const {return fPhiMax;}
  43   Double_t                     GetPhiMin()                     const {return fPhiMin;}
  44   Double_t                     GetConeRadius()                 const {return fConeRadius;}
  45   Double_t                     GetOverlapThreshold()           const {return fOverlapThreshold;}
  46   Double_t                     GetPtProtojetMin()              const {return fPtProtoJetMin;}
  47   Double_t                     GetRForRho()                    const {return fRRho;}
  48   Double_t                     GetCaching()                    const {return fCaching;}
  49   Double_t                     GetSplitMergeStoppingScale()    const {return fSplitMergeStoppingScale;}
  50   Double_t                     GetMinJetPt()                   const {return fMinJetPt;}
  51   Double_t                     GetGridScatter()                const {return fGridScatter;}
  52   Double_t                     GetKtScatter()                  const {return fKtScatter;}
  53   Double_t                     GetMeanGhostKt()                const {return fMeanGhostKt;}
  54
  55   // Setters
  56   void SetBGAlgorithm(Int_t value)                     {fBGAlgo = value;}
  57   void SetBGMode(Bool_t value)                         {fBGMode = value;}
  58   void SetCaching(Bool_t value)                        {fCaching = value;}
  59   void SetComment(TString com)                         {fComment=com;}
  60   void SetComment(const char* com)                     {AliJetHeader::SetComment(com);}
  61   void SetGhostEtaMax(Double_t f)                      {fGhostEtaMax = f;}
  62   void SetGhostArea(Double_t f)                        {fGhostArea = f;}
  63   void SetActiveAreaRepeats(Int_t f)                   {fActiveAreaRepeats =f;}
  64   void SetAreaTypeNumber(Int_t f)                      {fAreaTypeNumber = f;}
  65   void SetEffectiveRFact(Double_t value)               {fEffectiveRFact = value;}
  66   void SetConeRadius(Double_t value)                   {fConeRadius = value;}
  67   void SetMinJetPt(Double_t value)                     {fMinJetPt = value;}
  68   void SetNPassMax(Int_t value)                        {fNPassMax = value;}
  69   void SetNumberOfJetsToErase(Int_t value)             {fNHardJets = value;}
  70   void SetOverlapThreshold(Double_t value)             {fOverlapThreshold = value;}
  71   void SetPhiRange(Double_t fmin, Double_t fmax)       {fPhiMin = fmin; fPhiMax = fmax;}
  72   void SetPtProtojetMin(Double_t value)                {fPtProtoJetMin = value;}
  73   void SetRapRange(Double_t fmin, Double_t fmax)       {fRapMin = fmin; fRapMax = fmax;}
  74   void SetRForRho(Double_t value)                      {fRRho = value;}
  75   void SetSplitMergeScale(Int_t value)                 {fSplitMergeScaleNumber = value;}
  76   void SetSplitMergeStoppingScale(Double_t value)      {fSplitMergeStoppingScale = value;}
  77   void SetGridScatter(Double_t value)                  {fGridScatter = value;}
  78   void SetKtScatter(Double_t value)                    {fKtScatter = value;}
  79   void SetMeanGhostKt(Double_t value)                  {fMeanGhostKt = value;}
  80
  81   // others
  82   void PrintParameters() const;
  83
  84  protected:
  85
  86
  87   Int_t    fActiveAreaRepeats;        // How many times do you want to caculate active areas?
  88   Int_t    fAreaTypeNumber;           // Kind of area
  89   Int_t    fBGAlgo;                   // Algorithm for rho calculus
  90   Int_t    fNHardJets;                // Number of hard jets not to take into account for rho estimation
  91
  92   Bool_t   fBGMode;                   // Do we subtract BG or not?
  93   Bool_t   fCaching;                  // Do we record found cones for this set of data?
  94   Double_t fConeRadius;               // Cone radius
  95   Double_t fEffectiveRFact;           // Radius for Voronoi diagram
  96   Double_t fGhostEtaMax;              // Maximum eta in which a ghost can be generated
  97   Double_t fGhostArea;                // Area of one ghost
  98   Double_t fGridScatter;              // fractional random fluctuations of the position of the ghosts on the y-phi grid
  99   Double_t fKtScatter;                // fractional random fluctuations of the tranverse momentum of the ghosts on the y-phi grid
 100   Double_t fMeanGhostKt;              // average transverse momentum of the ghosts.
 101   Double_t fMinJetPt;                 // Jet minimum energy
 102   Int_t    fNPassMax;                 // maximum number of passes
 103   Double_t fOverlapThreshold;         // overlap parameter
 104   Double_t fPhiMax, fPhiMin;          // Phi range
 105   Double_t fPtProtoJetMin;            // pT min of protojets
 106   Double_t fRapMax, fRapMin;          // Eta range
 107   Double_t fRRho;                     // Radius to determine rho
 108   Int_t    fSplitMergeScaleNumber;    // Kind of recombination in split/merge procedure, there's only one
 109   Double_t fSplitMergeStoppingScale;  // Stopping scale for split/merge procedure in case of area calculus
 110
 111   ClassDef(AliSISConeJetHeader,4)
 112 };
 113
 114 #endif