Changes in readers/DelphesROOT.cpp [94083e7:341014c] in git

File:

• readers/DelphesROOT.cpp (modified) (12 diffs)

readers/DelphesROOT.cpp

@@ -18 +18 @@
 
 #include <algorithm>
+#include <iostream>
+#include <memory>
+#include <sstream>
 #include <stdexcept>
-#include <iostream>
-#include <sstream>
-#include <memory>
 
 #include <map>
 #include <vector>
 
-#include <stdlib.h>
 #include <signal.h>
 #include <stdio.h>
-
+#include <stdlib.h>
+
+#include "TApplication.h"
 #include "TROOT.h"
-#include "TApplication.h"
-
+
+#include "TClonesArray.h"
+#include "TDatabasePDG.h"
 #include "TFile.h"
-#include "TClonesArray.h"
+#include "TLorentzVector.h"
 #include "TObjArray.h"
+#include "TParticlePDG.h"
 #include "TStopwatch.h"
-#include "TDatabasePDG.h"
-#include "TParticlePDG.h"
-#include "TLorentzVector.h"
-
-#include "modules/Delphes.h"
-#include "classes/DelphesStream.h"
+
 #include "classes/DelphesClasses.h"
 #include "classes/DelphesFactory.h"
-
+#include "classes/DelphesStream.h"
+#include "modules/Delphes.h"
+
+#include "ExRootAnalysis/ExRootProgressBar.h"
+#include "ExRootAnalysis/ExRootTreeBranch.h"
+#include "ExRootAnalysis/ExRootTreeReader.h"
 #include "ExRootAnalysis/ExRootTreeWriter.h"
-#include "ExRootAnalysis/ExRootTreeReader.h"
-#include "ExRootAnalysis/ExRootTreeBranch.h"
-#include "ExRootAnalysis/ExRootProgressBar.h"
-
-
 
 using namespace std;
 
 //---------------------------------------------------------------------------
-
 
 //---------------------------------------------------------------------------
@@ -94 +91 @@
   if(argc < 4)
   {
-    cout << " Usage: " << appName << " config_file" << " output_file" << " input_file(s)" << endl;
+    cout << " Usage: " << appName << " config_file"
+         << " output_file"
+         << " input_file(s)" << endl;
     cout << " config_file - configuration file in Tcl format," << endl;
     cout << " output_file - output file in ROOT format," << endl;
@@ -122 +121 @@
 
     branchEvent = treeWriter->NewBranch("Event", HepMCEvent::Class());
-   
+
     confReader = new ExRootConfReader;
     confReader->ReadFile(argv[1]);
@@ -129 +128 @@
     modularDelphes->SetConfReader(confReader);
     modularDelphes->SetTreeWriter(treeWriter);
-    
+
     TChain *chain = new TChain("Delphes");
-    
+
     factory = modularDelphes->GetFactory();
     allParticleOutputArray = modularDelphes->ExportArray("allParticles");
@@ -145 +144 @@
       chain->Add(argv[i]);
       ExRootTreeReader *treeReader = new ExRootTreeReader(chain);
-      
+
       inputFile = TFile::Open(argv[i]);
 
@@ -153 +152 @@
         throw runtime_error(message.str());
       }
-      
+
       numberOfEvents = treeReader->GetEntries();
-      TClonesArray *branchParticle   = treeReader->UseBranch("Particle");
+      TClonesArray *branchParticle = treeReader->UseBranch("Particle");
       TClonesArray *branchHepMCEvent = treeReader->UseBranch("Event");
-     
+
       if(numberOfEvents <= 0) continue;
 
@@ -169 +168 @@
       for(Int_t entry = 0; entry < numberOfEvents && !interrupted; ++entry)
       {
-    
+
         treeReader->ReadEntry(entry);
 
-        // -- TBC need also to include event weights --  
-        
-        eve = (HepMCEvent*) branchHepMCEvent->At(0);
+        // -- TBC need also to include event weights --
+
+        eve = (HepMCEvent *)branchHepMCEvent->At(0);
         element = static_cast<HepMCEvent *>(branchEvent->NewEntry());
-        
+
         element->Number = eventCounter;
 
@@ -197 +196 @@
         element->ProcTime = eve->ProcTime;
 
-        for(Int_t j=0; j < branchParticle->GetEntriesFast(); j++)
-        {     
-         
-          gen = (GenParticle*) branchParticle->At(j);     
+        for(Int_t j = 0; j < branchParticle->GetEntriesFast(); j++)
+        {
+
+          gen = (GenParticle *)branchParticle->At(j);
           candidate = factory->NewCandidate();
 
           candidate->Momentum = gen->P4();
-          candidate->Position.SetXYZT(gen->X, gen->Y, gen->Z, gen->T*1.0E3*c_light);
-          
+          candidate->Position.SetXYZT(gen->X, gen->Y, gen->Z, gen->T * 1.0E3 * c_light);
+
           candidate->PID = gen->PID;
           candidate->Status = gen->Status;
-     
+
           candidate->M1 = gen->M1;
           candidate->M2 = gen->M2;
@@ -216 +215 @@
 
           candidate->Charge = gen->Charge;
-          candidate->Mass   = gen->Mass;
-    
+          candidate->Mass = gen->Mass;
+
           allParticleOutputArray->Add(candidate);
- 
+
           pdgCode = TMath::Abs(gen->PID);
 
@@ -231 +230 @@
           }
         }
-        
+
         modularDelphes->ProcessTask();
 
@@ -242 +241 @@
         ++eventCounter;
       }
- 
+
       progressBar.Update(eventCounter, eventCounter, kTRUE);
       progressBar.Finish();
 
       inputFile->Close();
-    
+
       delete treeReader;
-      
     }
 
@@ -262 +260 @@
     delete outputFile;
     delete chain;
-    
+
     return 0;
   }