|
[Rivet] accessing HepMC/PdfInfo.h from a rivet analysisstefania spagnolo stefania.spagnolo at le.infn.itTue Feb 18 08:48:05 GMT 2014
Hi Frank, sorry for the lazy question - and thanks that's very useful- Cheers, S On Feb 17, 2014, at 9:00 PM, Frank Siegert wrote: > Hi Stefania, > > I understand you have already managed to get to the HepMC event with > the Rivet::Event::genEvent() method: > http://rivet.hepforge.org/code/dev/a00236.html#a57dbce8fdc74bca2487f624010724272 > > From there it's a matter of HepMC to get to the PdfInfo. The function > you want is HepMC::GenEvent::pdf_info(), cf. line 234ff here: > http://svnweb.cern.ch/world/wsvn/hepmc/trunk/HepMC/GenEvent.h > > Hope that helps, > Frank > > > On 17 February 2014 20:09, stefania spagnolo > <stefania.spagnolo at le.infn.it> wrote: >> >> Dear all, >> I'm running a rivet analysis on hepmc files produced with aMC at NLO and >> showered with Herwig++; I need to retrieve some extra info written in the >> lhe original file which are not propagated to the hepmc format. The problem >> I have is that the hepmc file misses a few events (parton shower failure >> while running, I guess) and therefore I need to make sure I'm associating >> the correct lhe event to the hepmc event. >> >> Therefore, I'd like to retrieve some special event info from the hepmc event >> header; in particular I've the impression that the pdfScale of the event >> might be useful ... >> >> In my hepmc file this is the 4th word in the Event header row: >> >> E 3 -1 2.8168094000000000e+01 1.4036799198303040e-01 7.6616229391653504e-03 >> 0 -1 245 10003 10007 0 1 1.0000000000000000e+00 >> >> And in the lhe file I can easily find it again in the event header: >> >> <event> >> 6 66 0.70177165E-03 0.28168094E+02 0.79577472E-01 0.11965156E+00 >> >> Cross checking the value, might allow to make sure I'm reading the same >> event from the two files (hepmc and lhe). >> >> So the question: >> How can I access the PdfInfo object from the hepmc event in the analyze >> method of a rivet analysis ? >> I've tried to browse the code without much success ... >> >> Thanks for any suggestion, >> Cheers, >> Stefania >> >> >> Stefania Spagnolo >> ---------------------------------------------------------------------------------------------- >> Sezione INFN di Lecce >> Dip. di Matematica e Fisica "Ennio DeGiorgi", Univ. del Salento >> Phone: +39 0832 297439 >> http://www.dmf.unisalento.it/~spagnolo/ >> ATLAS experiment >> ---------------------------------------------------------------------------------------------- >> >> >> >> >> _______________________________________________ >> Rivet mailing list >> Rivet at projects.hepforge.org >> http://www.hepforge.org/lists/listinfo/rivet > Stefania Spagnolo ---------------------------------------------------------------------------------------------- Sezione INFN di Lecce Dip. di Matematica e Fisica "Ennio DeGiorgi", Univ. del Salento Phone: +39 0832 297439 http://www.dmf.unisalento.it/~spagnolo/ ATLAS experiment ---------------------------------------------------------------------------------------------- -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://www.hepforge.org/lists-archive/rivet/attachments/20140218/896f9cae/attachment.html> -------------- next part -------------- A non-text attachment was scrubbed... Name: smime.p7s Type: application/pkcs7-signature Size: 1763 bytes Desc: not available URL: <http://www.hepforge.org/lists-archive/rivet/attachments/20140218/896f9cae/attachment.bin>
More information about the Rivet mailing list |