sherpa is hosted by Hepforge, IPPP Durham
close Warning: Can't synchronize with repository "(default)" (/hepforge/svn/sherpa does not appear to be a Subversion repository.). Look in the Trac log for more information.

Ticket #247: 0002-Save-some-PDF-information-for-subevents.patch

File 0002-Save-some-PDF-information-for-subevents.patch, 1.3 KB (added by Generic User (don't modify these fields), 12 years ago)

  • SHERPA/Tools/HepMC2_Interface.C 

    From d288721cfd0bbbeb5e907ea901d91318b300f1ae Mon Sep 17 00:00:00 2001
From: Valery Yundin <yuvalery@gmail.com>
Date: Thu, 14 Jun 2012 18:58:13 +0200
Subject: [PATCH 2/2] Save some PDF information for subevents.

Flavours, x1, x2 and mu2_fac.
---
 SHERPA/Tools/HepMC2_Interface.C |   11 +++++++++--
 1 file changed, 9 insertions(+), 2 deletions(-)

diff --git a/SHERPA/Tools/HepMC2_Interface.C b/SHERPA/Tools/HepMC2_Interface.C
index ac24119..d534ccd 100644
    a b bool HepMC2_Interface::Sherpa2ShortHepMC(ATOOLS::Blob_List *const blobs, 
    281281      }
    282282      subevent->add_vertex(subvertex);
    283283      // not enough info in subevents to set PDFInfo properly,
    284       // so leave it blank
    285       HepMC::PdfInfo subpdfinfo;
     284      // so set only flavours, x1, x2, and q from the Signal_Process
     285      HepMC::PdfInfo subpdfinfo(*event.pdf_info());
     286      double q = sqrt(sub->m_mu2[stp::fac]);
     287      if (q != subpdfinfo.scalePDF()) {
     288        subpdfinfo.set_scalePDF(q);
     289        subpdfinfo.set_pdf1(0.);
     290        subpdfinfo.set_pdf2(0.);
     291      }
    286292      subevent->set_pdf_info(subpdfinfo);
     293
    287294      // add weight
    288295      std::vector<double> subweight; subweight.push_back(sub->m_result);
    289296      subevent->weights()=subweight;