JETAN/AliJetBkg.h

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