]> git.uio.no Git - u/mrichter/AliRoot.git/blame - JETAN/FASTJETAN/AliFastJetBkg.h
FASTJETAN module
[u/mrichter/AliRoot.git] / JETAN / FASTJETAN / 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//---------------------------------------------------------------------
f88804ed 18#include "fastjet/PseudoJet.hh"
13fa2a91 19#ifdef __CINT__
20namespace fastjet {
21 class PsuedoJet;
22}
23#endif
d89b8229 24class TString;
25class AliFastJetInput;
26
27class AliFastJetBkg : public TObject
28{
29 public:
30 AliFastJetBkg();
31 AliFastJetBkg(const AliFastJetBkg &input);
32 AliFastJetBkg& operator=(const AliFastJetBkg& source);
33 virtual ~AliFastJetBkg() {;}
34 void SetHeader(AliJetHeader *header) {fHeader=header;}
35 void SetFastJetInput(AliFastJetInput *fjinput) {fInputFJ=fjinput;}
36 void BkgFastJetb(Double_t& x,Double_t& y, Double_t& z);
37 void BkgFastJetWoHardest(Double_t& x,Double_t& y, Double_t& z);
38 Float_t BkgFastJet();
39 Float_t BkgChargedFastJet();
40 Float_t BkgStat();
41 Float_t BkgFastJetCone(TClonesArray* fAODJets);
42
7aa72c24 43 Bool_t EmcalAcceptance(Float_t eta, Float_t phi, Float_t radius) const;
d89b8229 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& meanarea,
49 vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
50 void CalcRhoWoHardest(Double_t& median, Double_t& sigma, Double_t& meanarea,
51 vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
52
53 AliJetHeader* fHeader; //! header
54 AliFastJetInput* fInputFJ; //! input particles
55
56 ClassDef(AliFastJetBkg, 2) // Fastjet backgroud analysis
57
58};
59
60#endif