Update from Sandun
[u/mrichter/AliRoot.git] / JETAN / AliFastJetBkg.h
CommitLineData
d89b8229 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//---------------------------------------------------------------------
13fa2a91 18#ifdef __CINT__
19namespace fastjet {
20 class PsuedoJet;
21}
22#endif
d89b8229 23class TString;
24class AliFastJetInput;
25
26class 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
7aa72c24 42 Bool_t EmcalAcceptance(Float_t eta, Float_t phi, Float_t radius) const;
d89b8229 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