]>
Commit | Line | Data |
---|---|---|
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 | 32 | using std::cout; |
33 | using std::endl; | |
2551af9d | 34 | ClassImp(AliSISConeJetHeader) |
35 | ||
36 | //////////////////////////////////////////////////////////////////////// | |
37 | ||
38 | AliSISConeJetHeader::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 | ||
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; | |
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 | } |