]> git.uio.no Git - u/mrichter/AliRoot.git/blob - JETAN/AliFastJetHeaderV1.cxx
Updates needed for full jet reconstruction (charged + emcal) [Magali Estienne]
[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(1.0), 
39     fAlgorithm(fastjet::kt_algorithm),
40     fStrategy(fastjet::Best),
41     fRecombScheme(fastjet::BIpt_scheme),
42     fGhostEtaMax(2.0),
43     fGhostArea(0.05),
44     fActiveAreaRepeats(1),
45     fAreaType(fastjet::active_area), 
46     fPtMin(5.0),
47     fPhiMin(0),
48     fPhiMax(TMath::TwoPi()),
49     fDebug(0),
50     fBGMode(0)
51 {
52   // Constructor
53   
54   Double_t rapmax = fGhostEtaMax - fRparam;
55   Double_t rapmin = -fGhostEtaMax + fRparam;
56   SetRapRange(rapmin, rapmax);
57   
58 }
59
60 ////////////////////////////////////////////////////////////////////////
61
62 void AliFastJetHeaderV1::PrintParameters() const
63 {
64   // prints out parameters of jet algorithm
65
66   cout << "FastJet algorithm  parameters:" << endl;
67   
68   cout << "-- Jet Definition --- " << endl;
69   cout << "R " << fRparam << endl;
70   cout << "Jet Algorithm " << fAlgorithm << endl; 
71   cout << "Strategy " << fStrategy << endl;  
72   cout << "Recombination Scheme " << fRecombScheme << endl; 
73   
74   cout << "-- Ghosted Area Spec parameters --- " << endl;
75   cout << "Ghost Eta Max " << fGhostEtaMax << endl;
76   cout << "Ghost Area " << fGhostArea << endl;
77   cout << "Active Area Repeats " << fActiveAreaRepeats << endl;
78   
79   cout << "-- Area Definition parameters --- " << endl;
80   cout << "Area Type " << fAreaType << endl; 
81   
82   cout << "-- Cluster Sequence Area parameters --- " << endl;
83   cout << "pt min " << fPtMin << endl; 
84   
85   cout << "-- Range Definition parameters --- " << endl;
86   cout << " bkg rapidity range from  " << fRapMin << " to " << fRapMax << endl;
87   cout << " bkg phi range from " << fPhiMin << " to " << fPhiMax << endl;
88
89 }