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