]>
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 | ||
32 | ClassImp(AliSISConeJetHeader) | |
33 | ||
34 | //////////////////////////////////////////////////////////////////////// | |
35 | ||
36 | AliSISConeJetHeader::AliSISConeJetHeader(): | |
37 | AliJetHeader("AliSISConeJetHeader"), | |
38 | fActiveAreaRepeats(1), | |
8838ab7a | 39 | fAreaTypeNumber(4), |
40 | fBGAlgo(1), | |
41 | fBGMode(1), | |
2551af9d | 42 | fCaching(0), |
8838ab7a | 43 | fConeRadius(0.7), |
44 | fDebug(0), | |
2551af9d | 45 | fEffectiveRFact(1), |
8838ab7a | 46 | fGhostEtaMax(4.0), |
2551af9d | 47 | fGhostArea(0.05), |
8838ab7a | 48 | fGridScatter(1), |
49 | fKtScatter(0.1), | |
50 | fMeanGhostKt(1e-100), | |
51 | fMinJetPt(2), | |
2551af9d | 52 | fNPassMax(0), |
8838ab7a | 53 | fOverlapThreshold(0.75), |
2551af9d | 54 | fPhiMax(TMath::TwoPi()), |
55 | fPhiMin(0), | |
2551af9d | 56 | fPtProtoJetMin(2), |
8838ab7a | 57 | fRapMax(0.9), |
58 | fRapMin(-0.9), | |
2551af9d | 59 | fSplitMergeScaleNumber(0), |
8838ab7a | 60 | fSplitMergeStoppingScale(0) |
2551af9d | 61 | { |
8838ab7a | 62 | // Constructor |
2551af9d | 63 | } |
64 | ||
65 | //////////////////////////////////////////////////////////////////////// | |
66 | ||
67 | void AliSISConeJetHeader::PrintParameters() const | |
68 | { | |
69 | // prints out parameters of jet algorithm | |
70 | ||
71 | cout << "SISConeJet algorithm parameters:"<<endl; | |
72 | ||
73 | cout<<"Cone Radius = "<<fConeRadius<<endl; | |
74 | cout<<"Overlap parameter = "<<fOverlapThreshold<<endl; | |
75 | cout<<"Maximum number of runs = "<<fNPassMax<<endl; | |
76 | cout<<"Pt min of protojets = "<<fPtProtoJetMin<<endl; | |
77 | cout<<"Do we record cones of these events ? (0 = no, 1 = yes) = "<<fCaching<<endl; | |
78 | ||
79 | cout << "Background subtraction parameters :" <<endl; | |
8838ab7a | 80 | if (fAreaTypeNumber == 1) cout<<"Kind of area used = Active area"<<endl; |
81 | if (fAreaTypeNumber == 2) cout<<"Kind of area used = Active area explicit ghosts"<<endl; | |
82 | if (fAreaTypeNumber == 3) cout<<"Kind of area used = One ghost passive area"<<endl; | |
83 | if (fAreaTypeNumber == 4) cout<<"Kind of area used = Passive area"<<endl; | |
84 | if (fAreaTypeNumber == 5) cout<<"Kind of area used = Voronoi"<<endl; | |
85 | if (fBGAlgo == 0) cout<<"Algorithm for rho calculus = kT"<<endl; | |
86 | if (fBGAlgo == 1) cout<<"Algorithm for rho calculus = Cambridge"<<endl; | |
2551af9d | 87 | cout<<"Eta max in which ghosts wil be generated = "<<fGhostEtaMax<<endl; |
88 | cout<<"Ghost area = "<<fGhostArea<<endl; | |
89 | cout<<"Background will be studied in ["<<fRapMin<<","<<fRapMax<<"] in eta and ["<<fPhiMin<<","<<fPhiMax<<"] in phi"<<endl; | |
8838ab7a | 90 | cout<<"Kind of recombination for split/merge procedure = SM_pttilde"<<endl; |
91 | cout<<"Stopping scale for split/merge procedure = "<<fSplitMergeStoppingScale<<endl; | |
2551af9d | 92 | cout<<"Do we repeat active area calculus? (0 = no, 1 = yes) = "<<fActiveAreaRepeats<<endl; |
8838ab7a | 93 | cout<<"Fractional random fluctuations of the position of the ghosts on the y-phi grid = "<<fGridScatter<<endl; |
94 | cout<<"Fractional random fluctuations of the tranverse momentum of the ghosts on the y-phi grid = "<<fKtScatter<<endl; | |
95 | cout<<"Average transverse momentum of the ghosts = "<<fMeanGhostKt<<endl; | |
2551af9d | 96 | |
97 | cout<<"Jets PtMin = "<<fMinJetPt<<endl; | |
8838ab7a | 98 | |
2551af9d | 99 | } |