]> git.uio.no Git - u/mrichter/AliRoot.git/blob - JETAN/AliFastJetHeaderV1.cxx
Clustering with free sigma added.
[u/mrichter/AliRoot.git] / JETAN / AliFastJetHeaderV1.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16
17 //---------------------------------------------------------------------
18 // FastJet v2.3.4 finder algorithm interface
19 // Finder Header Class 
20 // Author: Rafael.Diaz.Valdes@cern.ch
21 //---------------------------------------------------------------------
22
23 #include <Riostream.h>
24 #include <TMath.h>
25
26 #include "fastjet/ClusterSequenceArea.hh"
27 #include "fastjet/AreaDefinition.hh"
28 #include "fastjet/JetDefinition.hh"
29
30 #include "AliFastJetHeaderV1.h"
31
32 ClassImp(AliFastJetHeaderV1)
33
34 ////////////////////////////////////////////////////////////////////////
35
36 AliFastJetHeaderV1::AliFastJetHeaderV1():
37     AliJetHeader("AliFastJetHeaderV1"),
38     fRparam(0.4),
39     fRparamBkg(0.4),  
40     fAlgorithm(fastjet::kt_algorithm),
41     fStrategy(fastjet::Best),
42     fRecombScheme(fastjet::BIpt_scheme),
43     fGhostEtaMax(2.0),
44     fGhostArea(0.01),
45     fActiveAreaRepeats(1),
46     fAreaType(fastjet::active_area), 
47     fPtMin(5.0),
48     fMinCellEt(0.),
49     fRapMax(0.),
50     fRapMin(0.),
51     fPhiMax(TMath::TwoPi()),
52     fPhiMin(0),
53     fBGMode(0)
54 {
55   // Constructor
56   
57   Double_t rapmax = fGhostEtaMax - fRparam;
58   Double_t rapmin = -fGhostEtaMax + fRparam;
59   SetRapRange(rapmin, rapmax);
60   
61 }
62
63 ////////////////////////////////////////////////////////////////////////
64
65 void AliFastJetHeaderV1::PrintParameters() const
66 {
67   // prints out parameters of jet algorithm
68
69   cout << "FastJet algorithm  parameters:" << endl;
70   
71   cout << "-- Jet Definition --- " << endl;
72   cout << "R " << fRparam << endl;
73   cout << "Jet Algorithm " << fAlgorithm << endl; 
74   cout << "Strategy " << fStrategy << endl;  
75   cout << "Recombination Scheme " << fRecombScheme << endl; 
76   
77   cout << "-- Ghosted Area Spec parameters --- " << endl;
78   cout << "Ghost Eta Max " << fGhostEtaMax << endl;
79   cout << "Ghost Area " << fGhostArea << endl;
80   cout << "Active Area Repeats " << fActiveAreaRepeats << endl;
81   
82   cout << "-- Area Definition parameters --- " << endl;
83   cout << "Area Type " << fAreaType << endl; 
84   
85   cout << "-- Cluster Sequence Area parameters --- " << endl;
86   cout << "pt min " << fPtMin << endl; 
87   
88   cout << "-- Range Definition parameters --- " << endl;
89   cout << " bkg rapidity range from  " << fRapMin << " to " << fRapMax << endl;
90   cout << " bkg phi range from " << fPhiMin << " to " << fPhiMax << endl;
91
92 }