Changeset 14ae668 in git for examples/ExternalFastJetHepMC.cpp

Timestamp:

Dec 12, 2014, 1:54:46 PM (10 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

3d10d1f

Parents:

a8782e8

Message:

added comments in ExternalFast

File:

• examples/ExternalFastJetHepMC.cpp (modified) (9 diffs)

examples/ExternalFastJetHepMC.cpp

@@ -18 +18 @@
 
 /*
-
-First, you need to set DELPHES_DIR, FASTJET_DIR and ROOT_DIR:
+########################################################################
+
+
+This simple example shows how to use Delphes with an external fastjet installation.
+Events in hepmc format are read via the DelphesHepMC reader. 
+
+In order to run this example you first, you need to set the paths to your Delphes, FastJet
+and ROOT installations (DELPHES_DIR, FASTJET_DIR and ROOT_DIR):
 
 DELPHES_DIR=<path to Delphes installation>
@@ -25 +31 @@
 ROOT_DIR=<path to ROOT installation>
 
-Then run the following commands to build ExternalFastJetHepMC:
+Then run the following commands to build the executable:
 
 DELPHES_LIB="-L$DELPHES_DIR -lDelphesNoFastJet"
@@ -38 +44 @@
 LDFLAGS="$FASTJET_LIB $ROOT_LIB $DELPHES_LIB"
 
-g++ $CXXFLAGS $LDFLAGS ExternalFastJetHepMC.cpp -o ExternalFastJetHepMC
+g++ $CXXFLAGS $LDFLAGS examples/ExternalFastJetHepMC.cpp -o examples/ExternalFastJetHepMC
+
+Then run from the main Delphes dir (you need an event file in hepmc format):
+
+./examples/ExternalFastJetHepMC cards/delphes_card_CMS_NoFastJet.tcl file.hepmc
+
+
+########################################################################
 
 */
@@ -163 +176 @@
 
     ClusterSequence::print_banner();
+
 //    recomb = new WinnerTakeAllRecombiner();
 //    definition = new JetDefinition(antikt_algorithm, 0.5, recomb, Best);
+  
     definition = new JetDefinition(antikt_algorithm, 0.5);
-
-    inputArray = modularDelphes->ImportArray("Calorimeter/towers");
+    
+    
+    // Define your input candidates to fastjet (by default particle-flow objects).
+    // If you want pure calorimeter towers change "EFlowMerger/eflow" into "Calorimeter/towers":
+     
+    inputArray = modularDelphes->ImportArray("EFlowMerger/eflow");
+
     inputIterator = inputArray->MakeIterator();
+
+
+    // start reading hepmc file
 
     i = 2;
@@ -204 +227 @@
       }
 
+      
       reader->SetInputFile(inputFile);
 
@@ -214 +238 @@
         stableParticleOutputArray, partonOutputArray) && !interrupted)
       {
-        if(reader->EventReady())
+        
+         // loop over events
+        if(reader->EventReady())
         {
           ++eventCounter;
@@ -220 +246 @@
           if(eventCounter > skipEvents)
           {
-            modularDelphes->ProcessTask();
+            
+            // run delphes reconstruction
+            modularDelphes->ProcessTask();
             
             inputList.clear();
             inputIterator->Reset();
-            while((candidate = static_cast<Candidate*>(inputIterator->Next())))
+        
+        
+            // pass delphes candidates to fastjet clustering  
+             while((candidate = static_cast<Candidate*>(inputIterator->Next())))
             {
               momentum = candidate->Momentum;
@@ -230 +261 @@
               inputList.push_back(jet);
             }
-            ClusterSequence sequence(inputList, *definition);
+           
+            // run fastjet clustering 
+            ClusterSequence sequence(inputList, *definition);
             outputList.clear();
             outputList = sorted_by_pt(sequence.inclusive_jets(0.0));
 
-            // tell the user what was done
+            
+            // Prints for the first event:
             //  - the description of the algorithm used
             //  - show the inclusive jets as
             //      {index, rapidity, phi, pt}
             //----------------------------------------------------------
-            if(eventCounter == skipEvents + 1)
+           
+            if(eventCounter == skipEvents + 1)
             {
               cout << "Ran " << definition->description() << endl;
@@ -258 +293 @@
           reader->Clear();
         }
-      }
+      } // end of event loop
 
       if(inputFile != stdin) fclose(inputFile);