<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Microsoft Yi Baiti";
panose-1:3 0 5 0 0 0 0 0 0 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:"Calibri","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoPlainText">Dear Rivet experts,<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Sorry to pester you again, I'm becoming a problem customer :)<o:p></o:p></p>
<p class="MsoPlainText">I have changed my paths to use the new MCGenerators_hepmc2.06.08/... executables but I now get the following error (just after the fastjets header) when I try to use my rivet analysis on any hemc2g file:<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><i><span style="color:#C00000">python: ../../include/Rivet/Math/Vector4.hh:380: double Rivet::FourMomentum::mass() const: Assertion `Rivet::isZero(mass2()) || mass2() > 0' failed.<o:p></o:p></span></i></p>
<p class="MsoPlainText"><i><span style="color:#C00000">Aborted<o:p></o:p></span></i></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Do you have any suggestions as to how I can solves this issue?<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Thanks again,<o:p></o:p></p>
<p class="MsoPlainText">James<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><span lang="EN-US" style="mso-fareast-language:EN-GB">-----Original Message-----<br>
From: Andy Buckley [mailto:andy.buckley@ed.ac.uk] <br>
Sent: 18 June 2012 13:19<br>
To: Frank Siegert<br>
Cc: James Henderson; rivet@projects.hepforge.org<br>
Subject: Re: [Rivet] HepMV version issue</span></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Is three concurrent answers enough for you, James? ;-)<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Andy<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">On 18/06/12 13:13, Frank Siegert wrote:<o:p></o:p></p>
<p class="MsoPlainText">> Hi James,<o:p></o:p></p>
<p class="MsoPlainText">> <o:p></o:p></p>
<p class="MsoPlainText">> I think this might be due to the 2.06.05 event file containing
<o:p></o:p></p>
<p class="MsoPlainText">> additional lines which make HepMC 2.03.11 fail. I thought this was
<o:p></o:p></p>
<p class="MsoPlainText">> fixed in 2.03.10, but seems like not.<o:p></o:p></p>
<p class="MsoPlainText">> <o:p></o:p></p>
<p class="MsoPlainText">> Would it be an option for you, to use HepMC 2.06.05 in your setup,
<o:p></o:p></p>
<p class="MsoPlainText">> i.e. use the following AFS location?<o:p></o:p></p>
<p class="MsoPlainText">> /afs/cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.05/<o:p></o:p></p>
<p class="MsoPlainText">> <o:p></o:p></p>
<p class="MsoPlainText">> Cheers,<o:p></o:p></p>
<p class="MsoPlainText">> Frank<o:p></o:p></p>
<p class="MsoPlainText">> <o:p></o:p></p>
<p class="MsoPlainText">> On 18 June 2012 14:01, James Henderson <<a href="mailto:James.Henderson@physics.ox.ac.uk"><span style="color:windowtext;text-decoration:none">James.Henderson@physics.ox.ac.uk</span></a>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>> Hello,<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> I have been asked to use my rivet code to analyse some hepmc2g files
<o:p></o:p></p>
<p class="MsoPlainText">>> that a colleague has produced from Sherpa 1.4.02. My issue is that my
<o:p></o:p></p>
<p class="MsoPlainText">>> rivet analysis won't run over his samples but, when I copy his run
<o:p></o:p></p>
<p class="MsoPlainText">>> card and use this to produce my own hepmc2g file the rivet analysis works fine.<o:p></o:p></p>
<p class="MsoPlainText">>> I can see that the difference between the hepmc files is the version
<o:p></o:p></p>
<p class="MsoPlainText">>> number; he has used 2.06.05 and my files are 2.03.11.<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> The error message I recieve is:<o:p></o:p></p>
<p class="MsoPlainText">>> User-supplied cross-section = 2.551110e+01 pb Rivet 1.8.0 running on
<o:p></o:p></p>
<p class="MsoPlainText">>> machine pplxwn34.physics.ox.ac.uk (x86_64) CommonIO::read_heavy_ion
<o:p></o:p></p>
<p class="MsoPlainText">>> setting badbit.<o:p></o:p></p>
<p class="MsoPlainText">>> CommonIO::read_pdf_info setting badbit.<o:p></o:p></p>
<p class="MsoPlainText">>> CommonIO::read_vertex setting badbit.<o:p></o:p></p>
<p class="MsoPlainText">>> CommonIO::read_vertex setting badbit.<o:p></o:p></p>
<p class="MsoPlainText">>> ...<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> I belive this is due to the 2.06.05 hepmc file not having a heavy-ion
<o:p></o:p></p>
<p class="MsoPlainText">>> line where it is expected. Could you tell me if there is a way to
<o:p></o:p></p>
<p class="MsoPlainText">>> make rivet run over the 2.06.05 hepmc files?<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> My setup is as follows:<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> Rivet 1.8.0 setup from:<o:p></o:p></p>
<p class="MsoPlainText">>> /afs/cern.ch/sw/lcg/external/MCGenerators/rivet/1.8.0/x86_64-slc5-gcc<o:p></o:p></p>
<p class="MsoPlainText">>> 43-opt/rivetenv.sh<o:p></o:p></p>
<p class="MsoPlainText">>> Agile 1.2.0 setup from:<o:p></o:p></p>
<p class="MsoPlainText">>> /afs/cern.ch/sw/lcg/external/MCGenerators/agile/1.2.0/x86_64-slc5-gcc<o:p></o:p></p>
<p class="MsoPlainText">>> 43-opt/agileenv.sh<o:p></o:p></p>
<p class="MsoPlainText">>> Sherpa 1.4.0.2 from:<o:p></o:p></p>
<p class="MsoPlainText">>> /afs/cern.ch/sw/lcg/external/MCGenerators/Sherpa/1.4.0.2/x86_64-slc5-<o:p></o:p></p>
<p class="MsoPlainText">>> gcc43-opt/<o:p></o:p></p>
<p class="MsoPlainText">>> - used to create hepmc files<o:p></o:p></p>
<p class="MsoPlainText">>> HepMC 2.03.11 setup from both the rivetenv and agileenv scripts.<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> Thank you for your time,<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> James Henderson<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> _______________________________________________<o:p></o:p></p>
<p class="MsoPlainText">>> Rivet mailing list<o:p></o:p></p>
<p class="MsoPlainText">>> <a href="mailto:Rivet@projects.hepforge.org"><span style="color:windowtext;text-decoration:none">Rivet@projects.hepforge.org</span></a><o:p></o:p></p>
<p class="MsoPlainText">>> <a href="http://www.hepforge.org/lists/listinfo/rivet">
<span style="color:windowtext;text-decoration:none">http://www.hepforge.org/lists/listinfo/rivet</span></a><o:p></o:p></p>
<p class="MsoPlainText">> _______________________________________________<o:p></o:p></p>
<p class="MsoPlainText">> Rivet mailing list<o:p></o:p></p>
<p class="MsoPlainText">> <a href="mailto:Rivet@projects.hepforge.org"><span style="color:windowtext;text-decoration:none">Rivet@projects.hepforge.org</span></a><o:p></o:p></p>
<p class="MsoPlainText">> <a href="http://www.hepforge.org/lists/listinfo/rivet">
<span style="color:windowtext;text-decoration:none">http://www.hepforge.org/lists/listinfo/rivet</span></a><o:p></o:p></p>
<p class="MsoPlainText">> <o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">--<o:p></o:p></p>
<p class="MsoPlainText">Dr Andy Buckley, SUPA Advanced Research Fellow Particle Physics Expt Group, University of Edinburgh<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
</div>
</body>
</html>