JETAN/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 #ifdef __CINT__
  19 namespace fastjet {
  20   class PsuedoJet;
  21 }
  22 #endif
  23 class TString;
  24 class AliFastJetInput;
  25
  26 class AliFastJetBkg : public TObject
  27 {
  28  public:
  29   AliFastJetBkg();
  30   AliFastJetBkg(const AliFastJetBkg &input);
  31   AliFastJetBkg& operator=(const AliFastJetBkg& source);
  32   virtual          ~AliFastJetBkg() {;}
  33   void             SetHeader(AliJetHeader *header)  {fHeader=header;}
  34   void             SetFastJetInput(AliFastJetInput *fjinput)  {fInputFJ=fjinput;}
  35   void             BkgFastJetb(Double_t& x,Double_t& y, Double_t& z);
  36   void             BkgFastJetWoHardest(Double_t& x,Double_t& y, Double_t& z);
  37   Float_t          BkgFastJet();
  38   Float_t          BkgChargedFastJet();
  39   Float_t          BkgStat();
  40   Float_t          BkgFastJetCone(TClonesArray* fAODJets);
  41
  42   Bool_t           EmcalAcceptance(Float_t eta, Float_t phi, Float_t radius) const;
  43   static Double_t  BkgFunction(Double_t *x,Double_t *par);
  44
  45  private:
  46   Double_t         CalcRho(vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
  47   void             CalcRhob(Double_t& median, Double_t& sigma, Double_t& meanarea,
  48                             vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
  49   void             CalcRhoWoHardest(Double_t& median, Double_t& sigma, Double_t& meanarea,
  50                                     vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
  51
  52   AliJetHeader*    fHeader;  //! header
  53   AliFastJetInput* fInputFJ; //! input particles
  54
  55   ClassDef(AliFastJetBkg, 2)   //  Fastjet backgroud analysis
  56
  57 };
  58
  59 #endif