Changeset 01f9722 in git

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

Location:

classes

Files:

• DelphesClasses.h (modified) (1 diff)
• DelphesLHEFReader.cc (modified) (5 diffs)
• DelphesLHEFReader.h (modified) (1 diff)

classes/DelphesClasses.h

@@ -75 +75 @@
 
   Float_t Weight; // weight for the event | hepup.XWGTUP
+  Float_t CrossSection; // cross-section (read from init, implemented only for Wizard evgen)
   Float_t ScalePDF; // scale in GeV used in the calculation of the PDFs in the event | hepup.SCALUP
   Float_t AlphaQED; // value of the QED coupling used in the event | hepup.AQEDUP

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;

classes/DelphesLHEFReader.h

@@ -79 +79 @@
 
   int fParticleCounter, fProcessID;
-  double fWeight, fScalePDF, fAlphaQCD, fAlphaQED;
+  double fCrossSection, fWeight, fScalePDF, fAlphaQCD, fAlphaQED;
 
   int fPID, fStatus, fM1, fM2, fC1, fC2;