|
[Rivet] accessing HepMC/PdfInfo.h from a rivet analysisFrank Siegert frank.siegert at cern.chMon Feb 17 20:00:12 GMT 2014
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
More information about the Rivet mailing list |