Changeset 14ae668 in git for examples/ExternalFastJetBasic.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/ExternalFastJetBasic.cpp (modified) (8 diffs)

examples/ExternalFastJetBasic.cpp

@@ -19 +19 @@
 /*
 
-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 are specified via the multidimentionnal array "EVENTS" (for an example reading 
+an hepmc file see ExternalFastJetHepMC.cpp). 
+
+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 +34 @@
 ROOT_DIR=<path to ROOT installation>
 
-Then run the following commands to build ExternalFastJetBasic:
+Then run the following commands to build the executable:
 
 DELPHES_LIB="-L$DELPHES_DIR -lDelphesNoFastJet"
@@ -38 +47 @@
 LDFLAGS="$FASTJET_LIB $ROOT_LIB $DELPHES_LIB"
 
-g++ $CXXFLAGS $LDFLAGS ExternalFastJetBasic.cpp -o ExternalFastJetBasic
+g++ $CXXFLAGS $LDFLAGS examples/ExternalFastJetBasic.cpp -o examples/ExternalFastJetBasic
+
+Then run from the main Delphes dir:
+
+./examples/ExternalFastJetBasic cards/delphes_card_CMS_NoFastJet.tcl
+
+
+########################################################################
+
 
 */
@@ -95 +112 @@
 //---------------------------------------------------------------------------
 
+
+// this function converts input event array into Delphes candidates (defined below)
+
 void ConvertInput(Int_t event, DelphesFactory *factory, TObjArray *allParticleOutputArray, TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray);
+
+
+//----------------------------------------------------------------------------------------------------------------------------
 
 int main(int argc, char *argv[])
@@ -149 +172 @@
     modularDelphes->InitTask();
 
+    
+    // fastjet definition
     ClusterSequence::print_banner();
     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();
 
@@ -159 +188 @@
     {
       modularDelphes->Clear();
+      
+      // convert EVENT input array into Delphes internal format
       ConvertInput(event, factory, allParticleOutputArray, stableParticleOutputArray, partonOutputArray);
+      
+      // run Delphes reconstruction
       modularDelphes->ProcessTask();
 
       inputList.clear();
       inputIterator->Reset();
+      
+      
+      // pass delphes candidates to fastjet clustering  
       while((candidate = static_cast<Candidate*>(inputIterator->Next())))
       {
@@ -170 +206 @@
         inputList.push_back(jet);
       }
+     
+      // run clustering 
       ClusterSequence sequence(inputList, *definition);
       outputList.clear();
@@ -204 +242 @@
   }
 }
+
+
+// ------------------------------------------------------------------------------------------------------------------------------------
 
 void ConvertInput(Int_t event, DelphesFactory *factory, TObjArray *allParticleOutputArray, TObjArray *stableParticleOutputArray, TObjArray *partonOutputArray)