Fix for the problem during PbPb run of Nov 2010 (Indra)

[u/mrichter/AliRoot.git] / JETAN / AliUA1JetFinderV1.h

diff --git a/JETAN/AliUA1JetFinderV1.h b/JETAN/AliUA1JetFinderV1.h
index a71abde7262c22becda0dfc121668a771c438c8f..f40b9f4a433d587629aa0610264e616f0ed11c35 100644 (file)
--- a/JETAN/AliUA1JetFinderV1.h
+++ b/JETAN/AliUA1JetFinderV1.h
@@ -15,6 +15,7 @@
 #include "AliJetFinder.h"
 class AliUA1JetHeaderV1;
 class TH2F;
+class TH1F;
 
 class AliUA1JetFinderV1 : public AliJetFinder
 {
@@ -23,45 +24,50 @@ class AliUA1JetFinderV1 : public AliJetFinder
   AliUA1JetFinderV1();
   ~AliUA1JetFinderV1();
 
-  // getters
-
-  // setters
-  void SetJetHeader(AliUA1JetHeaderV1* h) {fHeader= h;}
   // others
   void FindJets();
   void RunAlgoritm(Float_t EtbgTotal, Double_t dEtTotal, Int_t& nJets,
                    Float_t* etJet,Float_t* etaJet, Float_t* phiJet,
                    Float_t* etallJet, Int_t* ncellsJet);
 
-  void SubtractBackg(Int_t& nIn, Int_t&nJ, Float_t&EtbgTotalN,
-                      Float_t* ptT, Float_t* etaT, Float_t* phiT,
-                      Float_t* etJet,Float_t* etaJet, Float_t* phiJet,
-                      Float_t* etsigJet,Int_t* multJet, Int_t* injet);
-
-  void SubtractBackgCone(Int_t& nIn, Int_t&nJ,Float_t& EtbgTotalN,
-                      Float_t* ptT, Float_t* etaT, Float_t* phiT,
-                      Float_t* etJet,Float_t* etaJet, Float_t* phiJet,
-                      Float_t* etsigJet, Int_t* multJet, Int_t* injet);
-
-  void SubtractBackgRatio(Int_t& nIn, Int_t&nJ,Float_t& EtbgTotalN,
-                      Float_t* ptT, Float_t* etaT, Float_t* phiT,
-                      Float_t* etJet,Float_t* etaJet, Float_t* phiJet,
-                      Float_t* etsigJet, Int_t* multJet, Int_t* injet);
-
-  void SubtractBackgStat(Int_t& nIn, Int_t&nJ,Float_t&EtbgTotalN,
-                      Float_t* ptT, Float_t* etaT, Float_t* phiT,
-                      Float_t* etJet,Float_t* etaJet, Float_t* phiJet,
-                      Float_t* etsigJet, Int_t* multJet, Int_t* injet);
+  void SubtractBackg(const Int_t& nIn, const Int_t&nJ, Float_t&EtbgTotalN,
+                    const Float_t* ptT, const Float_t* etaT, const Float_t* phiT,
+                     Float_t* etJet, const Float_t* etaJet, const Float_t* phiJet,
+                    Float_t* etsigJet,Int_t* multJet, Int_t* injet);
+  
+  void SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float_t& EtbgTotalN,
+                        const Float_t* ptT, const Float_t* etaT, const Float_t* phiT,
+                        Float_t* etJet, const Float_t* etaJet, const Float_t* phiJet,
+                        Float_t* etsigJet, Int_t* multJet, Int_t* injet);
+  
+  void SubtractBackgRatio(const Int_t& nIn, const Int_t&nJ,Float_t& EtbgTotalN,
+                         const Float_t* ptT, const Float_t* etaT, const Float_t* phiT,
+                         Float_t* etJet, const Float_t* etaJet, const Float_t* phiJet,
+                         Float_t* etsigJet, Int_t* multJet, Int_t* injet);
+  
+  void SubtractBackgStat(const Int_t& nIn, const Int_t&nJ,Float_t&EtbgTotalN,
+                        const Float_t* ptT, const Float_t* etaT, const Float_t* phiT,
+                        Float_t* etJet, const Float_t* etaJet, const Float_t* phiJet,
+                        Float_t* etsigJet, Int_t* multJet, Int_t* injet);
   void Reset();
   void Init();
-  void WriteJHeaderToFile();
+  void WriteJHeaderToFile() const;
+  
+  enum {kMaxJets = 30};
 
  protected:
-
-  AliUA1JetHeaderV1* fHeader;         // pointer to jet header
-  TH2F           * fLego;           //! Lego Histo
-
-  ClassDef(AliUA1JetFinderV1,1)
+  AliUA1JetFinderV1(const AliUA1JetFinderV1& rJetF1);
+  AliUA1JetFinderV1& operator = (const AliUA1JetFinderV1& rhsf);
+  TH2F*          fLego;           //Lego Histo
+  // temporary histos for background, reset for each event, no need to stream
+  TH1F*          fhEtJet[kMaxJets];   //! histogram for background subtraction
+  TH1F*          fhAreaJet[kMaxJets]; //! histogram for background subtraction (store global not to create it with every event
+  TH1F*          fhEtBackg;           //! histogram for background subtraction
+  TH1F*          fhAreaBackg;           //! histogram for background subtraction
+
+  //
+
+  ClassDef(AliUA1JetFinderV1,2)
 };
 
 #endif