]>
Commit | Line | Data |
---|---|---|
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 | ||
26 | #include "fastjet/ClusterSequenceArea.hh" | |
27 | #include "fastjet/AreaDefinition.hh" | |
28 | #include "fastjet/JetDefinition.hh" | |
29 | ||
30 | #include "AliSISConeJetHeader.h" | |
31 | ||
32 | using std::cout; | |
33 | using std::endl; | |
34 | ClassImp(AliSISConeJetHeader) | |
35 | ||
36 | //////////////////////////////////////////////////////////////////////// | |
37 | ||
38 | AliSISConeJetHeader::AliSISConeJetHeader(): | |
39 | AliJetHeader("AliSISConeJetHeader"), | |
40 | fActiveAreaRepeats(1), | |
41 | fAreaTypeNumber(4), | |
42 | fBGAlgo(1), | |
43 | fNHardJets(2), | |
44 | fBGMode(1), | |
45 | fCaching(0), | |
46 | fConeRadius(0.7), | |
47 | fEffectiveRFact(1), | |
48 | fGhostEtaMax(4.0), | |
49 | fGhostArea(0.05), | |
50 | fGridScatter(1), | |
51 | fKtScatter(0.1), | |
52 | fMeanGhostKt(1e-100), | |
53 | fMinJetPt(2), | |
54 | fNPassMax(0), | |
55 | fOverlapThreshold(0.75), | |
56 | fPhiMax(TMath::TwoPi()), | |
57 | fPhiMin(0), | |
58 | fPtProtoJetMin(2), | |
59 | fRapMax(0.9), | |
60 | fRapMin(-0.9), | |
61 | fRRho(0.5), | |
62 | fSplitMergeScaleNumber(0), | |
63 | fSplitMergeStoppingScale(0) | |
64 | { | |
65 | // Constructor | |
66 | } | |
67 | ||
68 | //////////////////////////////////////////////////////////////////////// | |
69 | ||
70 | void 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; | |
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; | |
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; | |
93 | cout<<"Kind of recombination for split/merge procedure = SM_pttilde"<<endl; | |
94 | cout<<"Stopping scale for split/merge procedure = "<<fSplitMergeStoppingScale<<endl; | |
95 | cout<<"Do we repeat active area calculus? (0 = no, 1 = yes) = "<<fActiveAreaRepeats<<endl; | |
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; | |
99 | ||
100 | cout<<"Jets PtMin = "<<fMinJetPt<<endl; | |
101 | ||
102 | } |