Changeset 8e602e5 in git for examples/ExternalFastJetHepMC.cpp

Timestamp:

Dec 11, 2014, 3:46:57 PM (10 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

a8782e8

Parents:

27eb358

Message:

move all cards from examples to cards directory, add ExternalFastJetBasic, rename StandaloneHepMC to ExternalFastJetHepMC

File:

• examples/ExternalFastJetHepMC.cpp (moved) (moved from examples/StandaloneHepMC.cpp ) (6 diffs)

examples/ExternalFastJetHepMC.cpp

@@ -46 +46 @@
 #include "fastjet/JetDefinition.hh"
 #include "fastjet/ClusterSequence.hh"
-#include "fastjet/Selector.hh"
-#include "fastjet/ClusterSequenceArea.hh"
-#include "fastjet/tools/JetMedianBackgroundEstimator.hh"
-
-#include "fastjet/plugins/SISCone/fastjet/SISConePlugin.hh"
-#include "fastjet/plugins/CDFCones/fastjet/CDFMidPointPlugin.hh"
-#include "fastjet/plugins/CDFCones/fastjet/CDFJetCluPlugin.hh"
-
-#include "fastjet/contribs/Nsubjettiness/Nsubjettiness.hh"
-#include "fastjet/contribs/Nsubjettiness/Njettiness.hh"
-#include "fastjet/contribs/Nsubjettiness/NjettinessPlugin.hh"
-#include "fastjet/contribs/Nsubjettiness/WinnerTakeAllRecombiner.hh"
+
+// #include "fastjet/contrib/Nsubjettiness.hh"
+// #include "fastjet/contrib/Njettiness.hh"
+// #include "fastjet/contrib/NjettinessPlugin.hh"
+// #include "fastjet/contrib/WinnerTakeAllRecombiner.hh"
 
 using namespace std;
 using namespace fastjet;
-using namespace fastjet::contrib;
+// using namespace fastjet::contrib;
 
 //---------------------------------------------------------------------------
@@ -76 +69 @@
 int main(int argc, char *argv[])
 {
-  char appName[] = "StandaloneHepMC";
+  char appName[] = "ExternalFastJetHepMC";
   stringstream message;
   FILE *inputFile = 0;
@@ -93 +86 @@
 
   JetDefinition *definition = 0;
-  JetDefinition::Recombiner *recomb = 0;
+//  JetDefinition::Recombiner *recomb = 0;
   vector<PseudoJet> inputList, outputList;
   PseudoJet jet;
@@ -145 +138 @@
 
     ClusterSequence::print_banner();
-    recomb = new WinnerTakeAllRecombiner();
-    definition = new JetDefinition(antikt_algorithm, 0.5, recomb, Best);
+//    recomb = new WinnerTakeAllRecombiner();
+//    definition = new JetDefinition(antikt_algorithm, 0.5, recomb, Best);
+    definition = new JetDefinition(antikt_algorithm, 0.5);
 
     inputArray = modularDelphes->ImportArray("Calorimeter/towers");
@@ -213 +207 @@
             ClusterSequence sequence(inputList, *definition);
             outputList.clear();
-            outputList = sorted_by_pt(sequence.inclusive_jets(10.0));
+            outputList = sorted_by_pt(sequence.inclusive_jets(0.0));
+
+            // tell the user what was done
+            //  - the description of the algorithm used
+            //  - show the inclusive jets as
+            //      {index, rapidity, phi, pt}
+            //----------------------------------------------------------
+            if(eventCounter == skipEvents + 1)
+            {
+              cout << "Ran " << definition->description() << endl;
+
+              // label the columns
+              printf("%5s %15s %15s %15s\n","jet #", "rapidity", "phi", "pt");
+
+              // print out the details for each jet
+              for (unsigned int i = 0; i < outputList.size(); i++) {
+                printf("%5u %15.8f %15.8f %15.8f\n",
+                       i, outputList[i].rap(), outputList[i].phi(),
+                       outputList[i].perp());
+              }
+            }
           }
 
@@ -230 +244 @@
 
     cout << "** Exiting..." << endl;
+
+    delete definition;
+//    delete recomb;
 
     delete reader;