]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - JETAN/AliUA1JetFinderV2.h
Adding rho's depenence on trigger track (M. Verweij)
[u/mrichter/AliRoot.git] / JETAN / AliUA1JetFinderV2.h
index b3823458681bd4725497636eba3181f4f1ddcf33..df055fc04195e284cc2f5def4017adb0d965bbd6 100644 (file)
 // manages the search for jets
 // Author: Rafael.Diaz.Valdes@cern.ch
 // (version in c++)
-// Modified to include neutral particles (magali.estienne@ires.in2p3.fr)
 //---------------------------------------------------------------------
 
 #include "AliJetFinder.h"
+
 class AliUA1JetHeaderV1;
 class TH2F;
+class TChain;
 
 class AliUA1JetFinderV2 : public AliJetFinder
 {
@@ -25,43 +26,50 @@ class AliUA1JetFinderV2 : public AliJetFinder
   ~AliUA1JetFinderV2();
 
   // others
+  void FindJetsC();
   void FindJets();
+  void RunAlgoritmC(Float_t EtbgTotal, Double_t dEtTotal, Int_t& nJets,
+                  Float_t* const etJet,Float_t* const etaJet, Float_t* const phiJet,
+                  Float_t* const etallJet, Int_t* const ncellsJet);
 
-  void RunAlgoritm(Int_t nIn, Float_t* etCell, Float_t* etaCell, Float_t* phiCell, 
-                  Int_t* flagCell, 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 RunAlgoritm(Int_t nIn, Float_t* etCell, Float_t* const etaCell, Float_t* phiCell, 
+                  Int_t* const flagCell, const Float_t* etCell2, const Float_t* etaCell2, const Float_t* phiCell2, 
+                  const Int_t* flagCell2, Float_t etbgTotal, Double_t dEtTotal, 
+                  Int_t& nJets, Float_t* const etJet, Float_t* const etaJet, Float_t* const phiJet,
+                  Float_t* const etallJet, Int_t* const ncellsJet);
+   
+  void SubtractBackgC(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* const etJet, const Float_t* etaJet, const Float_t* phiJet,
+                      Float_t* const etsigJet,Int_t* const multJet, Int_t* const injet);
 
-  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 SubtractBackg(const Int_t& nIn, const Int_t&nJ, Float_t&EtbgTotalN, const Float_t* ptT, const Int_t* vectT, 
+                    const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* cFlag2T, 
+                    const Int_t* sFlagT, Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet, 
+                    Float_t* const etsigJet, Int_t* const multJet, Int_t* const 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 SubtractBackgCone(const Int_t& nIn, const Int_t&nJ,Float_t& EtbgTotalN,
+                      const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
+                      Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
+                      Float_t* const etsigJet, Int_t* const multJet, Int_t* const 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 SubtractBackgRatio(const Int_t& nIn, const Int_t&nJ,Float_t& EtbgTotalN,
+                      const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, const Int_t* sFlagT,
+                      Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
+                      Float_t* const etsigJet, Int_t* const multJet, Int_t* const 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 SubtractBackgStat(const Int_t& nIn, const Int_t&nJ,Float_t&EtbgTotalN,
+                        const Float_t* ptT, const Float_t* etaT, const Float_t* phiT, const Int_t* cFlagT, 
+                        const Int_t* sFlagT, Float_t* const etJet, const Float_t* etaJet, const Float_t* phiJet,
+                        Float_t* const etsigJet, Int_t* const multJet, Int_t* const injet);
   void Reset();
-  void Init();
-  void WriteJHeaderToFile();
-
+  void InitTask(TChain* tree);
+  void WriteJHeaderToFile() const;
  protected:
   AliUA1JetFinderV2(const AliUA1JetFinderV2& rJetF1);
   AliUA1JetFinderV2& operator = (const AliUA1JetFinderV2& rhsf);
-  TH2F           * fLego;           //Lego Histo
-  Int_t fDebug;
-  Int_t fOpt;
+  TH2F           * fLego;           // Lego Histo
+  Int_t fOpt;                       // Detector option (charged only or charged+neutral)
 
   ClassDef(AliUA1JetFinderV2,1)
 };