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