one more tpcrefit added

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

#ifndef ALIUA1JETFINDERV1_H
#define ALIUA1JETFINDERV1_H

/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */


//---------------------------------------------------------------------
// UA1 Cone Algorithm Finder V1
// manages the search for jets
// Author: Rafael.Diaz.Valdes@cern.ch
// (version in c++)
//---------------------------------------------------------------------

#include "AliJetFinder.h"
class AliUA1JetHeaderV1;
class TH2F;
class TH1F;

class AliUA1JetFinderV1 : public AliJetFinder
{
 public:

  AliUA1JetFinderV1();
  ~AliUA1JetFinderV1();

  // 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(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() const;
  
  enum {kMaxJets = 30};

 protected:
  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