Lightweight cascade checks -> code implementation

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

diff --git a/JETAN/AliSISConeJetHeader.cxx b/JETAN/AliSISConeJetHeader.cxx
index 6e3ab6601d96d2b399d4b14d39749eb185e9f2fe..ea57a9c11b01d486cc9615de3471fbf50d1c2524 100644 (file)
--- a/JETAN/AliSISConeJetHeader.cxx
+++ b/JETAN/AliSISConeJetHeader.cxx
@@ -13,6 +13,7 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
+/* $Id$ */
 
 //---------------------------------------------------------------------
 // SISCone (FastJet v2.3.4) finder algorithm interface
@@ -23,41 +24,47 @@
 #include <Riostream.h>
 #include <TMath.h>
 
-#include "fastjet/AreaDefinition.hh"
-#include "fastjet/ClusterSequenceArea.hh"
-#include "fastjet/JetDefinition.hh"
-
 #include "AliSISConeJetHeader.h"
 
+using std::cout;
+using std::endl;
 ClassImp(AliSISConeJetHeader)
 
 ////////////////////////////////////////////////////////////////////////
 
 AliSISConeJetHeader::AliSISConeJetHeader():
-    AliJetHeader("AliSISConeJetHeader"),
-    fActiveAreaRepeats(1),
-    fCaching(0),
-    fConeRadius(0.4),
-    fEffectiveRFact(1),
-    fGhostArea(0.05),
-    fGhostEtaMax(2.0),
-    fMinJetPt(0),
-    fNPassMax(0),
-    fOverlapThreshold(0.5),
-    fPhiMax(TMath::TwoPi()),
-    fPhiMin(0),
-    fRapMax(-0.9),
-    fRapMin(0.9),
-    fPtProtoJetMin(2),
-    fSplitMergeScaleNumber(0),
-    fSplitMergeStoppingScale(0)
-
+  AliJetHeader("AliSISConeJetHeader"),
+  fActiveAreaRepeats(1),
+  fAreaTypeNumber(4),
+  fBGAlgo(1),
+  fCaching(0),
+  fConeRadius(0.7),
+  fEffectiveRFact(1),
+  fGhostEtaMax(4.0),
+  fGhostArea(0.05),
+  fGridScatter(1),
+  fKtScatter(0.1),
+  fMeanGhostKt(1e-100),
+  fMinJetPt(2),
+  fNPassMax(0),
+  fOverlapThreshold(0.75),
+  fPhiMax(TMath::TwoPi()),
+  fPhiMin(0),
+  fPtProtoJetMin(2),
+  fRapMax(0.9),
+  fRapMin(-0.9),
+  fRRho(0.5),
+  fSplitMergeScaleNumber(0),
+  fSplitMergeStoppingScale(0),
+  fRparamBkg(0.4),
+  fStrategy(fastjet::Best),
+  fRecombScheme(fastjet::BIpt_scheme),
+  fkUse4VectorArea(kTRUE)
 {
-  // Constructor
+  // Constructor  
 }
 
-////////////////////////////////////////////////////////////////////////
-
+//____________________________________________________________________________
 void AliSISConeJetHeader::PrintParameters() const
 {
   // prints out parameters of jet algorithm
@@ -71,14 +78,23 @@ void AliSISConeJetHeader::PrintParameters() const
   cout<<"Do we record cones of these events ? (0 = no, 1 = yes) = "<<fCaching<<endl;
 
   cout << "Background subtraction parameters :" <<endl;
-  //cout<<"Kind of area used = "<<<<endl;
+  if (fAreaTypeNumber == 1) cout<<"Kind of area used = Active area"<<endl;
+  if (fAreaTypeNumber == 2) cout<<"Kind of area used = Active area explicit ghosts"<<endl;
+  if (fAreaTypeNumber == 3) cout<<"Kind of area used = One ghost passive area"<<endl;
+  if (fAreaTypeNumber == 4) cout<<"Kind of area used = Passive area"<<endl;
+  if (fAreaTypeNumber == 5) cout<<"Kind of area used = Voronoi"<<endl;
+  if (fBGAlgo == 0) cout<<"Algorithm for rho calculus = kT"<<endl;
+  if (fBGAlgo == 1) cout<<"Algorithm for rho calculus = Cambridge"<<endl;
   cout<<"Eta max in which ghosts wil be generated = "<<fGhostEtaMax<<endl;
   cout<<"Ghost area = "<<fGhostArea<<endl;
   cout<<"Background will be studied in ["<<fRapMin<<","<<fRapMax<<"] in eta and ["<<fPhiMin<<","<<fPhiMax<<"] in phi"<<endl;
-  //cout<<"Kind of recombination for split/merge procedure = "<<<<endl;
-  //cout<<"Stopping scale for split/merge procedure = "<<<<endl;
+  cout<<"Kind of recombination for split/merge procedure = SM_pttilde"<<endl;
+  cout<<"Stopping scale for split/merge procedure = "<<fSplitMergeStoppingScale<<endl;
   cout<<"Do we repeat active area calculus? (0 = no, 1 = yes) = "<<fActiveAreaRepeats<<endl;
+  cout<<"Fractional random fluctuations of the position of the ghosts on the y-phi grid = "<<fGridScatter<<endl;       
+  cout<<"Fractional random fluctuations of the tranverse momentum of the ghosts on the y-phi grid = "<<fKtScatter<<endl;         
+  cout<<"Average transverse momentum of the ghosts = "<<fMeanGhostKt<<endl;       
 
   cout<<"Jets PtMin  = "<<fMinJetPt<<endl;
-  
+
 }