Changeset 457 in svn

Timestamp:

Jul 10, 2009, 12:46:53 AM (15 years ago)

Author:

Xavier Rouby

Message:

new class for inputing directly rootfiles from Delphes (with GEN tree)

Location:

trunk

Files:

• CHANGELOG (modified) (1 diff)
• Delphes.cpp (modified) (2 diffs)
• interface/DelphesRootConverter.h (added)
• src/DelphesRootConverter.cc (added)

trunk/CHANGELOG

@@ -5 +5 @@
    - additional data in the logfile : all names of used input files 
      (input event files, detector card, trigger card, beam files, particle table)
+   - it is now possible to read "delphes root files" as an input for Delphes (GEN tree)
 
  /-----------------------------------------\

trunk/Delphes.cpp

@@ -51 +51 @@
 #include "STDHEPConverter.h"
 #include "LHCOConverter.h"
+#include "DelphesRootConverter.h"
 
 #include "SmearUtil.h"
@@ -232 +233 @@
     }
   else if(line.rfind(".root") < line.length())
+  // can be either a root file from h2root (i.e. with "h101" tree)
+  // or a root file from Delphes  (i.e. with "GEN" tree)
     {
-      cout <<"**                   h2root file format detected                   **"<<endl;
-      cout <<"**                  This can take several minutes                  **"<< endl;
-      HEPEVTConverter converter(inputFileList,outputfilename,DET->PDGtable,DET->NEvents);
-    }
+      TFile f(line.c_str());
+      if (f.FindKey("GEN")) {
+        cout <<"**                Delphes ROOT file format detected                **"<<endl;
+        cout <<"**                  This can take several minutes                  **"<< endl;
+        DelphesRootConverter converter(inputFileList,outputfilename,DET->NEvents);
+      }
+      else
+      if (f.FindKey("h101")) {
+        cout <<"**                   h2root file format detected                   **"<<endl;
+        cout <<"**                  This can take several minutes                  **"<< endl;
+        HEPEVTConverter converter(inputFileList,outputfilename,DET->PDGtable,DET->NEvents);
+      }
+      else {
+            cerr << left  << setw(4) <<"**  "<<""
+                 << left  << setw(63) << line.c_str() <<""
+                 << right << setw(2) <<"**"<<endl;
+            cerr <<"**         ERROR: File format not identified --  Exiting...        **"<< endl;
+            cout <<"**                                                                 **"<< endl;
+            cout <<"*********************************************************************"<< endl;
+            return -1;
+        } // not found any interesting input tree
+        f.Close();
+      } // .root file
   else {
     cerr << left  << setw(4) <<"**  "<<""