4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
9 //--------------------------------------------------
10 // Method implementation for background studies and background subtraction with UA1 algorithms
12 // Author: magali.estienne@subatech.in2p3.fr
13 //-------------------------------------------------
18 class AliAODJetEventBackground;
20 class AliJetBkg : public TObject
24 AliJetBkg(const AliJetBkg &input);
26 AliJetBkg& operator=(const AliJetBkg& source);
27 void SetHeader(AliJetHeader *header) {fHeader=header;}
28 void SetCalTrkEvent(AliJetCalTrkEvent *evt) {fEvent=evt;}
29 Bool_t PtCutPass(Int_t id, Int_t nTracks);
30 Bool_t SignalCutPass(Int_t id, Int_t nTracks);
31 Float_t CalcJetAreaEtaCut(Float_t radius, const Float_t etaJet);
32 void CalcJetAndBckgAreaEtaCut(Bool_t calcOutsideArea, Float_t rc, const Int_t nJ, const Float_t* etaJet, Float_t* &areaJet, Float_t &areaOut);
34 void SubtractBackg(const Int_t& nIn, const Int_t&nJ, Float_t&EtbgTotalN, Float_t&sigmaN,
35 const Float_t* ptT, const Float_t* etaT, const Float_t* phiT,
36 Float_t* etJet, const Float_t* etaJet, const Float_t* phiJet,
37 Float_t* etsigJet, Int_t* multJetT, Int_t* multJetC, Int_t* multJet,
38 Int_t* injet, Float_t* &areaJet);
40 void SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float_t& EtbgTotalN, Float_t&sigmaN,
41 const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, Float_t* etJet,
42 const Float_t* etaJet, const Float_t* phiJet, Float_t* etsigJet,
43 Int_t* multJetT, Int_t* multJetC, Int_t* multJet, Int_t* injet, Float_t* &/*areaJet*/);
45 void SubtractBackgRatio(const Int_t& nIn, const Int_t&nJ,Float_t& EtbgTotalN, Float_t&sigmaN,
46 const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, Float_t* etJet,
47 const Float_t* etaJet, const Float_t* phiJet, Float_t* etsigJet,
48 Int_t* multJetT, Int_t* multJetC, Int_t* multJet, Int_t* injet, Float_t* &/*areaJet*/);
50 void SubtractBackgStat(const Int_t& nIn, const Int_t&nJ,Float_t&EtbgTotalN, Float_t&sigmaN,
51 const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, Float_t* etJet,
52 const Float_t* etaJet, const Float_t* phiJet, Float_t* etsigJet,
53 Int_t* multJetT, Int_t* multJetC, Int_t* multJet, Int_t* injet, Float_t* &areaJet);
54 void SetDebug(Int_t debug){fDebug = debug;}
59 // Double_t CalcRho(vector<fastjet::PseudoJet> input_particles,Double_t RparamBkg,TString method);
61 AliJetCalTrkEvent* fEvent; //! reader
62 AliJetHeader* fHeader; //! header
63 Int_t fDebug; // Debug option
65 // temporary histos for background, reset for each event, no need to stream
66 TH1F* fhEtJet[kMaxJets]; //! histogram for background subtraction
67 TH1F* fhAreaJet[kMaxJets]; //! histogram for background subtraction (store global not to create it with every event
68 TH1F* fhEtBackg; //! histogram for background subtraction
69 TH1F* fhAreaBackg; //! histogram for background subtraction
71 ClassDef(AliJetBkg, 2) // background jet analysis