JETAN/DEV/AliFastJetBkg.h

   1 #ifndef ALIFASTJETBKG_H
   2 #define ALIFASTJETBKG_H
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 /* $Id$ */
   8
   9 //---------------------------------------------------------------------
  10 // Class to calculate the background per unit area
  11 // manages the search for jets
  12 // Authors: Elena Bruna elena.bruna@yale.edu
  13 //          Sevil Salur ssalur@lbl.gov
  14 //
  15 // 2011 :
  16 // renamed from AliJetBkg to AliFastJetBkg as this class uses only FASTJET based algos
  17 //---------------------------------------------------------------------
  18
  19 class TString;
  20 class AliFastJetInput;
  21
  22 class AliFastJetBkg : public TObject
  23 {
  24  public:
  25   AliFastJetBkg();
  26   AliFastJetBkg(const AliFastJetBkg &input);
  27   AliFastJetBkg& operator=(const AliFastJetBkg& source);
  28   virtual          ~AliFastJetBkg() {;}
  29   void             SetHeader(AliJetHeader *header)  {fHeader=header;}
  30   void             SetFastJetInput(AliFastJetInput *fjinput)  {fInputFJ=fjinput;}
  31   void             BkgFastJetb(Double_t& x,Double_t& y, Double_t& z);
  32   void             BkgFastJetWoHardest(Double_t& x,Double_t& y, Double_t& z);
  33   Float_t          BkgFastJet();
  34   Float_t          BkgChargedFastJet();
  35   Float_t          BkgStat();
  36   Float_t          BkgFastJetCone(TClonesArray* fAODJets);
  37
  38   Bool_t           EmcalAcceptance(const Float_t eta, const Float_t phi, const Float_t radius) const;
  39   static Double_t  BkgFunction(Double_t *x,Double_t *par);
  40
  41  private:
  42   Double_t         CalcRho(vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
  43   void             CalcRhob(Double_t& median, Double_t& sigma, Double_t& meanarea,
  44                             vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
  45   void             CalcRhoWoHardest(Double_t& median, Double_t& sigma, Double_t& meanarea,
  46                                     vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
  47
  48   AliJetHeader*    fHeader;  //! header
  49   AliFastJetInput* fInputFJ; //! input particles
  50
  51   ClassDef(AliFastJetBkg, 2)   //  Fastjet backgroud analysis
  52
  53 };
  54
  55 #endif