Add missing histograms to estimate the UE in left and right cones at the generator...

[u/mrichter/AliRoot.git] / JETAN / AliFastJetHeaderV1.cxx

diff --git a/JETAN/AliFastJetHeaderV1.cxx b/JETAN/AliFastJetHeaderV1.cxx
index 4076dd1e25ec7c4e591e635ed1e071c4d5c924ae..da58d5990e91349ac0f802dc2513eb6279096470 100644 (file)
--- a/JETAN/AliFastJetHeaderV1.cxx
+++ b/JETAN/AliFastJetHeaderV1.cxx
@@ -13,6 +13,7 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
+/* $Id$ */
 
 //---------------------------------------------------------------------
 // FastJet v2.3.4 finder algorithm interface
@@ -23,67 +24,68 @@
 #include <Riostream.h>
 #include <TMath.h>
 
-#include "fastjet/ClusterSequenceArea.hh"
 #include "fastjet/AreaDefinition.hh"
 #include "fastjet/JetDefinition.hh"
 
 #include "AliFastJetHeaderV1.h"
 
+using std::cout;
+using std::endl;
 ClassImp(AliFastJetHeaderV1)
 
 ////////////////////////////////////////////////////////////////////////
 
 AliFastJetHeaderV1::AliFastJetHeaderV1():
-    AliJetHeader("AliFastJetHeaderV1"),
-    fRparam(1.0), 
-    fAlgorithm(fastjet::kt_algorithm),
-    fStrategy(fastjet::Best),
-    fRecombScheme(fastjet::BIpt_scheme),
-    fGhostEtaMax(2.0),
-    fGhostArea(0.05),
-    fActiveAreaRepeats(1),
-    fAreaType(fastjet::active_area), 
-    fPtMin(5.0),
-    fPhiMax(TMath::TwoPi()),
-    fPhiMin(0),
-    fDebug(0),
-    fBGMode(0)
+  AliJetHeader("AliFastJetHeaderV1"),
+  fRparam(0.4),
+  fRparamBkg(0.4),  
+  fAlgorithm(fastjet::kt_algorithm),
+  fBGAlgorithm(fastjet::kt_algorithm),
+  fStrategy(fastjet::Best),
+  fRecombScheme(fastjet::BIpt_scheme),
+  fGhostEtaMax(2.0),
+  fGhostArea(0.01),
+  fActiveAreaRepeats(1),
+  fAreaType(fastjet::active_area), 
+  fPtMin(5.0),
+  fMinCellEt(0.),
+  fRapMax(0.),
+  fRapMin(0.),
+  fPhiMax(TMath::TwoPi()),
+  fPhiMin(0),
+  fUse4VectorArea(kFALSE),
+  fkBkgFastJetb(kFALSE),
+  fktBkgFastJetWoHardest(kFALSE)
 {
-  // Constructor
-  
+  // Default constructor
   Double_t rapmax = fGhostEtaMax - fRparam;
   Double_t rapmin = -fGhostEtaMax + fRparam;
   SetRapRange(rapmin, rapmax);
   
 }
 
-////////////////////////////////////////////////////////////////////////
-
+//____________________________________________________________________________
 void AliFastJetHeaderV1::PrintParameters() const
 {
   // prints out parameters of jet algorithm
 
-  cout << "FastJet algorithm  parameters:" << endl;
-  
-  cout << "-- Jet Definition --- " << endl;
-  cout << "R " << fRparam << endl;
-  cout << "Jet Algorithm " << fAlgorithm << endl; 
-  cout << "Strategy " << fStrategy << endl;  
-  cout << "Recombination Scheme " << fRecombScheme << endl; 
-  
-  cout << "-- Ghosted Area Spec parameters --- " << endl;
-  cout << "Ghost Eta Max " << fGhostEtaMax << endl;
-  cout << "Ghost Area " << fGhostArea << endl;
-  cout << "Active Area Repeats " << fActiveAreaRepeats << endl;
-  
-  cout << "-- Area Definition parameters --- " << endl;
-  cout << "Area Type " << fAreaType << endl; 
-  
+  cout << "FastJet algorithm  parameters:"           << endl;
+  cout << "-- Jet Definition --- "                   << endl;
+  cout << "R "                                       << fRparam << endl;
+  cout << "Jet Algorithm "                           << fAlgorithm << endl; 
+  cout << "Bkg Jet Algorithm "                       << fBGAlgorithm << endl;
+  cout << "Strategy "                                << fStrategy << endl;  
+  cout << "Recombination Scheme "                    << fRecombScheme << endl; 
+  cout << "-- Ghosted Area Spec parameters --- "     << endl;
+  cout << "Ghost Eta Max "                           << fGhostEtaMax << endl;
+  cout << "Ghost Area "                              << fGhostArea << endl;
+  cout << "Active Area Repeats "                     << fActiveAreaRepeats << endl;
+  cout << "-- Area Definition parameters --- "       << endl;
+  cout << "Area Type "                               << fAreaType << endl; 
   cout << "-- Cluster Sequence Area parameters --- " << endl;
-  cout << "pt min " << fPtMin << endl; 
-  
-  cout << "-- Range Definition parameters --- " << endl;
-  cout << " bkg rapidity range from  " << fRapMin << " to " << fRapMax << endl;
-  cout << " bkg phi range from " << fPhiMin << " to " << fPhiMax << endl;
+  cout << "pt min "                                  << fPtMin << endl; 
+  cout << "-- Range Definition parameters --- "      << endl;
+  cout << " bkg rapidity range from  "               << fRapMin << " to " << fRapMax << endl;
+  cout << " bkg phi range from "                     << fPhiMin << " to " << fPhiMax << endl;
 
 }