Changeset 01f9722 in git for classes/DelphesLHEFReader.cc

Timestamp:

Aug 13, 2018, 3:01:20 PM (6 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

1d9c62a, 9047a02

Parents:

cdbca0d

Message:

read xsec from LHE files produced with Wizard

File:

• classes/DelphesLHEFReader.cc (modified) (5 diffs)

classes/DelphesLHEFReader.cc

@@ -31 +31 @@
 #include <iostream>
 #include <sstream>
+#include <string>
 
 #include <stdio.h>
@@ -54 +55 @@
 DelphesLHEFReader::DelphesLHEFReader() :
   fInputFile(0), fBuffer(0), fPDG(0),
-  fEventReady(kFALSE), fEventCounter(-1), fParticleCounter(-1)
+  fEventReady(kFALSE), fEventCounter(-1), fParticleCounter(-1),   fCrossSection(1)
+
 {
   fBuffer = new char[kBufferSize];
@@ -105 +107 @@
   if(!fgets(fBuffer, kBufferSize, fInputFile)) return kFALSE;
 
+  // read cross section string from LHE file header (meant to work on Wizard output only)
+  if(strstr(fBuffer, "<xsecinfo ") && fEventCounter < 0)
+  {
+    DelphesStream bufferStream(fBuffer);
+
+    string buf = fBuffer;    
+    string xsecstr="totxsec";
+    
+    if (buf.find(xsecstr) != string::npos)
+    {
+        buf = buf.substr(buf.find(xsecstr) + xsecstr.length()+1);
+        unsigned first = buf.find_first_of("\"")+1;
+        unsigned last = buf.find_last_of("\"");
+        xsecstr = buf.substr(first,last-first);
+        fCrossSection = stod(xsecstr);
+    }  
+  }
+  
   if(strstr(fBuffer, "<event>"))
   {
@@ -110 +130 @@
     fEventCounter = 1;
   }
+  
   else if(fEventCounter > 0)
   {
@@ -206 +227 @@
   element->ProcessID = fProcessID;
   element->Weight = fWeight;
+  element->CrossSection = fCrossSection;
+
   element->ScalePDF = fScalePDF;
   element->AlphaQED = fAlphaQED;