update

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

diff --git a/JETAN/AliDAJetHeader.cxx b/JETAN/AliDAJetHeader.cxx
index b6ae5fb4ede4c5df20f74eba40c6378239075324..bde89afbb67a7c4eef4c6f57943b8f16b44d1918 100644 (file)
--- a/JETAN/AliDAJetHeader.cxx
+++ b/JETAN/AliDAJetHeader.cxx
@@ -17,25 +17,40 @@
 // Deterministic Annealing Jet header class
 // Stores parameters of DA jet algorithm
 // Author: Davide Perrino (davide.perrino@ba.infn.it, davide.perrino@cern.ch)
+// 2011:
+// Adding FiducialEta/PhiMin/Max setters/getters and variables to accommodate with reader/finder splitting 
 //----------------------------------------------------------------------------
 
 #include "AliDAJetHeader.h"
 
 ClassImp(AliDAJetHeader)
 
+///////////////////////////////////////////////////////////////////////
 
-//---------------------------------------------------------------------
 AliDAJetHeader::AliDAJetHeader():
-       fDirectory("/home/perrino/events"),
-       fFileOut("jets.root"),
-       fPytOnly(kTRUE),
-       fPtCut(0.),
-       fEtaCut(.9),
-       fChgOnly(kTRUE),
-       fSelectJets(kTRUE),
-       fNclustMax(16),
-       fFixedCl(kFALSE),
-       fEtMin(10.)
+  AliJetHeader("AliDAJetHeader"),
+  fSelectJets(kTRUE),
+  fNclustMax(10),
+  fFixedCl(kFALSE),
+  fEtMin(10.),
+  fNeff(0),
+  fEtaEff(0.9),
+  fFidEtaMin(-0.9),
+  fFidEtaMax(0.9),
+  fFidPhiMin(0.),
+  fFidPhiMax(2*TMath::Pi())
+{
+  // Constructor
+}
+
+//---------------------------------------------------------------------
+void AliDAJetHeader::SetRadius(Float_t radius)
 {
-    // Constructor
+  // The radius requested is used to estimate the number of clusters
+  // to be found, in order to obtain jets with the expected area.
+  // It must not be intended as a sharp limit on the cluster radius
+    
+  Int_t nclust = (Int_t) (4.*fEtaEff/(radius*radius)) + 1;
+  SetNclust(nclust);
+
 }