- Timestamp:
- Apr 17, 2014, 12:28:09 PM (11 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 43c646a
- Parents:
- 210c033
- Location:
- modules
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/FastJetFinder.cc
r210c033 rc6321ad 60 60 61 61 FastJetFinder::FastJetFinder() : 62 fPlugin(0), f Definition(0), fAreaDefinition(0), fItInputArray(0), fRecomb(0), fNjettinessPlugin(0)62 fPlugin(0), fRecomb(0), fNjettinessPlugin(0), fDefinition(0), fAreaDefinition(0), fItInputArray(0) 63 63 { 64 64 … … 76 76 void FastJetFinder::Init() 77 77 { 78 79 JetDefinition::Plugin *plugin = NULL; 80 JetDefinition::Recombiner *recomb = NULL; 81 NjettinessPlugin *njet_plugin = NULL; 78 JetDefinition::Plugin *plugin = 0; 79 JetDefinition::Recombiner *recomb = 0; 80 NjettinessPlugin *njetPlugin = 0; 82 81 83 82 // read eta ranges … … 135 134 { 136 135 case 1: 137 fAreaDefinition = new fastjet::AreaDefinition(active_area_explicit_ghosts, GhostedAreaSpec(fGhostEtaMax, fRepeat, fGhostArea, fGridScatter, fPtScatter, fMeanGhostPt));136 fAreaDefinition = new AreaDefinition(active_area_explicit_ghosts, GhostedAreaSpec(fGhostEtaMax, fRepeat, fGhostArea, fGridScatter, fPtScatter, fMeanGhostPt)); 138 137 break; 139 138 case 2: 140 fAreaDefinition = new fastjet::AreaDefinition(one_ghost_passive_area, GhostedAreaSpec(fGhostEtaMax, fRepeat, fGhostArea, fGridScatter, fPtScatter, fMeanGhostPt));139 fAreaDefinition = new AreaDefinition(one_ghost_passive_area, GhostedAreaSpec(fGhostEtaMax, fRepeat, fGhostArea, fGridScatter, fPtScatter, fMeanGhostPt)); 141 140 break; 142 141 case 3: 143 fAreaDefinition = new fastjet::AreaDefinition(passive_area, GhostedAreaSpec(fGhostEtaMax, fRepeat, fGhostArea, fGridScatter, fPtScatter, fMeanGhostPt));142 fAreaDefinition = new AreaDefinition(passive_area, GhostedAreaSpec(fGhostEtaMax, fRepeat, fGhostArea, fGridScatter, fPtScatter, fMeanGhostPt)); 144 143 break; 145 144 case 4: 146 fAreaDefinition = new fastjet::AreaDefinition(VoronoiAreaSpec(fEffectiveRfact));145 fAreaDefinition = new AreaDefinition(VoronoiAreaSpec(fEffectiveRfact)); 147 146 break; 148 147 case 5: 149 fAreaDefinition = new fastjet::AreaDefinition(active_area, GhostedAreaSpec(fGhostEtaMax, fRepeat, fGhostArea, fGridScatter, fPtScatter, fMeanGhostPt));148 fAreaDefinition = new AreaDefinition(active_area, GhostedAreaSpec(fGhostEtaMax, fRepeat, fGhostArea, fGridScatter, fPtScatter, fMeanGhostPt)); 150 149 break; 151 150 default: … … 158 157 { 159 158 case 1: 160 plugin = new fastjet::CDFJetCluPlugin(fSeedThreshold, fConeRadius, fAdjacencyCut, fMaxIterations, fIratch, fOverlapThreshold);161 fDefinition = new fastjet::JetDefinition(plugin);159 plugin = new CDFJetCluPlugin(fSeedThreshold, fConeRadius, fAdjacencyCut, fMaxIterations, fIratch, fOverlapThreshold); 160 fDefinition = new JetDefinition(plugin); 162 161 break; 163 162 case 2: 164 plugin = new fastjet::CDFMidPointPlugin(fSeedThreshold, fConeRadius, fConeAreaFraction, fMaxPairSize, fMaxIterations, fOverlapThreshold);165 fDefinition = new fastjet::JetDefinition(plugin);163 plugin = new CDFMidPointPlugin(fSeedThreshold, fConeRadius, fConeAreaFraction, fMaxPairSize, fMaxIterations, fOverlapThreshold); 164 fDefinition = new JetDefinition(plugin); 166 165 break; 167 166 case 3: 168 plugin = new fastjet::SISConePlugin(fConeRadius, fOverlapThreshold, fMaxIterations, fJetPTMin);169 fDefinition = new fastjet::JetDefinition(plugin);167 plugin = new SISConePlugin(fConeRadius, fOverlapThreshold, fMaxIterations, fJetPTMin); 168 fDefinition = new JetDefinition(plugin); 170 169 break; 171 170 case 4: 172 fDefinition = new fastjet::JetDefinition(fastjet::kt_algorithm, fParameterR);171 fDefinition = new JetDefinition(kt_algorithm, fParameterR); 173 172 break; 174 173 case 5: 175 fDefinition = new fastjet::JetDefinition(fastjet::cambridge_algorithm, fParameterR);174 fDefinition = new JetDefinition(cambridge_algorithm, fParameterR); 176 175 break; 177 176 default: 178 177 case 6: 179 fDefinition = new fastjet::JetDefinition(fastjet::antikt_algorithm, fParameterR);178 fDefinition = new JetDefinition(antikt_algorithm, fParameterR); 180 179 break; 181 180 case 7: 182 recomb = new fastjet::contrib::WinnerTakeAllRecombiner();183 fDefinition = new fastjet::JetDefinition(fastjet::antikt_algorithm, fParameterR, recomb, Best);181 recomb = new WinnerTakeAllRecombiner(); 182 fDefinition = new JetDefinition(antikt_algorithm, fParameterR, recomb, Best); 184 183 break; 185 184 case 8: 186 njet _plugin = new fastjet::contrib::NjettinessPlugin(fN, Njettiness::wta_kt_axes, Njettiness::unnormalized_cutoff_measure, fBeta, fRcutOff);187 fDefinition = new fastjet::JetDefinition(njet_plugin);185 njetPlugin = new NjettinessPlugin(fN, Njettiness::wta_kt_axes, Njettiness::unnormalized_cutoff_measure, fBeta, fRcutOff); 186 fDefinition = new JetDefinition(njetPlugin); 188 187 break; 189 188 } … … 191 190 fPlugin = plugin; 192 191 fRecomb = recomb; 193 fNjettinessPlugin = njet _plugin;192 fNjettinessPlugin = njetPlugin; 194 193 195 194 ClusterSequence::print_banner(); -
modules/FastJetFinder.h
r210c033 rc6321ad 25 25 class AreaDefinition; 26 26 class Selector; 27 namespace contrib { 28 class NjettinessPlugin; 29 } 27 30 } 28 31 … … 42 45 void *fPlugin; //! 43 46 void *fRecomb; //! 44 void*fNjettinessPlugin; //!47 fastjet::contrib::NjettinessPlugin *fNjettinessPlugin; //! 45 48 46 49 fastjet::JetDefinition *fDefinition; //!
Note:
See TracChangeset
for help on using the changeset viewer.