- Timestamp:
- Dec 11, 2014, 3:46:57 PM (10 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- a8782e8
- Parents:
- 27eb358
- Location:
- examples
- Files:
-
- 1 added
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
examples/ExternalFastJetHepMC.cpp
r27eb358 r8e602e5 46 46 #include "fastjet/JetDefinition.hh" 47 47 #include "fastjet/ClusterSequence.hh" 48 #include "fastjet/Selector.hh" 49 #include "fastjet/ClusterSequenceArea.hh" 50 #include "fastjet/tools/JetMedianBackgroundEstimator.hh" 51 52 #include "fastjet/plugins/SISCone/fastjet/SISConePlugin.hh" 53 #include "fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh" 54 #include "fastjet/plugins/CDFCones/fastjet/CDFJetCluPlugin.hh" 55 56 #include "fastjet/contribs/Nsubjettiness/Nsubjettiness.hh" 57 #include "fastjet/contribs/Nsubjettiness/Njettiness.hh" 58 #include "fastjet/contribs/Nsubjettiness/NjettinessPlugin.hh" 59 #include "fastjet/contribs/Nsubjettiness/WinnerTakeAllRecombiner.hh" 48 49 // #include "fastjet/contrib/Nsubjettiness.hh" 50 // #include "fastjet/contrib/Njettiness.hh" 51 // #include "fastjet/contrib/NjettinessPlugin.hh" 52 // #include "fastjet/contrib/WinnerTakeAllRecombiner.hh" 60 53 61 54 using namespace std; 62 55 using namespace fastjet; 63 using namespace fastjet::contrib;56 // using namespace fastjet::contrib; 64 57 65 58 //--------------------------------------------------------------------------- … … 76 69 int main(int argc, char *argv[]) 77 70 { 78 char appName[] = " StandaloneHepMC";71 char appName[] = "ExternalFastJetHepMC"; 79 72 stringstream message; 80 73 FILE *inputFile = 0; … … 93 86 94 87 JetDefinition *definition = 0; 95 JetDefinition::Recombiner *recomb = 0;88 // JetDefinition::Recombiner *recomb = 0; 96 89 vector<PseudoJet> inputList, outputList; 97 90 PseudoJet jet; … … 145 138 146 139 ClusterSequence::print_banner(); 147 recomb = new WinnerTakeAllRecombiner(); 148 definition = new JetDefinition(antikt_algorithm, 0.5, recomb, Best); 140 // recomb = new WinnerTakeAllRecombiner(); 141 // definition = new JetDefinition(antikt_algorithm, 0.5, recomb, Best); 142 definition = new JetDefinition(antikt_algorithm, 0.5); 149 143 150 144 inputArray = modularDelphes->ImportArray("Calorimeter/towers"); … … 213 207 ClusterSequence sequence(inputList, *definition); 214 208 outputList.clear(); 215 outputList = sorted_by_pt(sequence.inclusive_jets(10.0)); 209 outputList = sorted_by_pt(sequence.inclusive_jets(0.0)); 210 211 // tell the user what was done 212 // - the description of the algorithm used 213 // - show the inclusive jets as 214 // {index, rapidity, phi, pt} 215 //---------------------------------------------------------- 216 if(eventCounter == skipEvents + 1) 217 { 218 cout << "Ran " << definition->description() << endl; 219 220 // label the columns 221 printf("%5s %15s %15s %15s\n","jet #", "rapidity", "phi", "pt"); 222 223 // print out the details for each jet 224 for (unsigned int i = 0; i < outputList.size(); i++) { 225 printf("%5u %15.8f %15.8f %15.8f\n", 226 i, outputList[i].rap(), outputList[i].phi(), 227 outputList[i].perp()); 228 } 229 } 216 230 } 217 231 … … 230 244 231 245 cout << "** Exiting..." << endl; 246 247 delete definition; 248 // delete recomb; 232 249 233 250 delete reader;
Note:
See TracChangeset
for help on using the changeset viewer.