[Rivet] accessing HepMC/PdfInfo.h from a rivet analysis

stefania spagnolo stefania.spagnolo at le.infn.it
Tue 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