Fork me on GitHub

Opened 12 years ago

Closed 12 years ago

#44 closed How to (fixed)

PDF Informations in Delphes

Reported by: muanza Owned by:
Priority: major Milestone:
Component: Delphes code Version: all recent versions (> 1.8)
Keywords: Cc: muanza@…

Description

Hello,

I produced MC samples with Herwig++ and Pythia8 and ran Delphes on their HepMC outputs.

I'd like to apply PDF reweighting to those root-tuples, so I need the usual x1,x2,flav1,flav2,Q² variables. The last one seems to be in the GEN tree whereas for the x's and flav's of the colliding partons I couldn't figure an obvious robust way to identify them for all processes.
Is there a recipe for this?

Attachments (2)

BlockClasses.h (19.0 KB ) - added by muanza 12 years ago.
Delphes_V_2.0.2/Utilities/ExRootAnalysis/interface/BlockClasses.h
HepMCConverter.cc (12.1 KB ) - added by muanza 12 years ago.
Delphes_V_2.0.2/src/HepMCConverter.cc

Download all attachments as: .zip

Change History (6)

comment:1 by Pavel Demin, 12 years ago

Hi Steve,

Is it possible to extract the x1,x2,flav1,flav2,Q² variables using information from a STDHEP common block?

Or should we calculate them from the kinematic parameters of incoming protons and partons?

Here is a basic description of the STDHEP contents:

http://cmssdt.cern.ch/SDT/lxr/source/SimDataFormats/GeneratorProducts/data/stdhep.inc

Regards,

Pavel

comment:2 by muanza, 12 years ago

Dear all,

I slightly modified Delphes_V_2.0.2 in order to add the PDF variables cited above into the GEN tree, only for the case of an HepMC input file. I retrieved these variables from the HepMC::PdfInfo class.

Please find hereby attached the 2 modified files.

Best regards, Steve

by muanza, 12 years ago

Attachment: BlockClasses.h added

Delphes_V_2.0.2/Utilities/ExRootAnalysis/interface/BlockClasses.h

by muanza, 12 years ago

Attachment: HepMCConverter.cc added

Delphes_V_2.0.2/src/HepMCConverter.cc

comment:3 by muanza, 12 years ago

Dear all,
As far as I'm concerned, this ticket can be closed.

Thanks
Regards, Steve

comment:4 by favereau, 12 years ago

Resolution: fixed
Status: newclosed
Note: See TracTickets for help on using tickets.