set minimum pT for dicing (M. Verweij)
[u/mrichter/AliRoot.git] / JETAN / AliSISConeJetHeader.cxx
CommitLineData
2551af9d 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// SISCone (FastJet v2.3.4) finder algorithm interface
19// Finder Header Class
20// Author: swensy.jangal@ires.in2p3.fr
21//---------------------------------------------------------------------
22
23#include <Riostream.h>
24#include <TMath.h>
25
2551af9d 26#include "fastjet/ClusterSequenceArea.hh"
8838ab7a 27#include "fastjet/AreaDefinition.hh"
2551af9d 28#include "fastjet/JetDefinition.hh"
29
30#include "AliSISConeJetHeader.h"
31
3a7af7bd 32using std::cout;
33using std::endl;
2551af9d 34ClassImp(AliSISConeJetHeader)
35
36////////////////////////////////////////////////////////////////////////
37
38AliSISConeJetHeader::AliSISConeJetHeader():
39 AliJetHeader("AliSISConeJetHeader"),
40 fActiveAreaRepeats(1),
8838ab7a 41 fAreaTypeNumber(4),
42 fBGAlgo(1),
b19ec4fe 43 fNHardJets(2),
8838ab7a 44 fBGMode(1),
2551af9d 45 fCaching(0),
8838ab7a 46 fConeRadius(0.7),
2551af9d 47 fEffectiveRFact(1),
8838ab7a 48 fGhostEtaMax(4.0),
2551af9d 49 fGhostArea(0.05),
8838ab7a 50 fGridScatter(1),
51 fKtScatter(0.1),
52 fMeanGhostKt(1e-100),
53 fMinJetPt(2),
2551af9d 54 fNPassMax(0),
8838ab7a 55 fOverlapThreshold(0.75),
2551af9d 56 fPhiMax(TMath::TwoPi()),
57 fPhiMin(0),
2551af9d 58 fPtProtoJetMin(2),
8838ab7a 59 fRapMax(0.9),
60 fRapMin(-0.9),
f92ba4be 61 fRRho(0.5),
2551af9d 62 fSplitMergeScaleNumber(0),
8838ab7a 63 fSplitMergeStoppingScale(0)
2551af9d 64{
8838ab7a 65 // Constructor
2551af9d 66}
67
68////////////////////////////////////////////////////////////////////////
69
70void AliSISConeJetHeader::PrintParameters() const
71{
72 // prints out parameters of jet algorithm
73
74 cout << "SISConeJet algorithm parameters:"<<endl;
75
76 cout<<"Cone Radius = "<<fConeRadius<<endl;
77 cout<<"Overlap parameter = "<<fOverlapThreshold<<endl;
78 cout<<"Maximum number of runs = "<<fNPassMax<<endl;
79 cout<<"Pt min of protojets = "<<fPtProtoJetMin<<endl;
80 cout<<"Do we record cones of these events ? (0 = no, 1 = yes) = "<<fCaching<<endl;
81
82 cout << "Background subtraction parameters :" <<endl;
8838ab7a 83 if (fAreaTypeNumber == 1) cout<<"Kind of area used = Active area"<<endl;
84 if (fAreaTypeNumber == 2) cout<<"Kind of area used = Active area explicit ghosts"<<endl;
85 if (fAreaTypeNumber == 3) cout<<"Kind of area used = One ghost passive area"<<endl;
86 if (fAreaTypeNumber == 4) cout<<"Kind of area used = Passive area"<<endl;
87 if (fAreaTypeNumber == 5) cout<<"Kind of area used = Voronoi"<<endl;
88 if (fBGAlgo == 0) cout<<"Algorithm for rho calculus = kT"<<endl;
89 if (fBGAlgo == 1) cout<<"Algorithm for rho calculus = Cambridge"<<endl;
2551af9d 90 cout<<"Eta max in which ghosts wil be generated = "<<fGhostEtaMax<<endl;
91 cout<<"Ghost area = "<<fGhostArea<<endl;
92 cout<<"Background will be studied in ["<<fRapMin<<","<<fRapMax<<"] in eta and ["<<fPhiMin<<","<<fPhiMax<<"] in phi"<<endl;
8838ab7a 93 cout<<"Kind of recombination for split/merge procedure = SM_pttilde"<<endl;
94 cout<<"Stopping scale for split/merge procedure = "<<fSplitMergeStoppingScale<<endl;
2551af9d 95 cout<<"Do we repeat active area calculus? (0 = no, 1 = yes) = "<<fActiveAreaRepeats<<endl;
8838ab7a 96 cout<<"Fractional random fluctuations of the position of the ghosts on the y-phi grid = "<<fGridScatter<<endl;
97 cout<<"Fractional random fluctuations of the tranverse momentum of the ghosts on the y-phi grid = "<<fKtScatter<<endl;
98 cout<<"Average transverse momentum of the ghosts = "<<fMeanGhostKt<<endl;
2551af9d 99
100 cout<<"Jets PtMin = "<<fMinJetPt<<endl;
8838ab7a 101
2551af9d 102}