<html> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> </head> <body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""> Hi Andy, <div class=""><br class=""> </div> <div class="">We did try --ignore-beams and --nskip options before but didn't help. Now trying again with the following steps:<br class=""> <div class=""> <div class=""><br class=""> </div> <div class="">I copy one of the standard analysis in Rivet, CMS_2015_<a href="http://i1356998.cc/" class="">I1356998</a>, to my working directory and modify the info file like the following, </div> <div class=""><br class=""> </div> <div class=""> <div class="">Beams: [ANY, ANY]</div> <div class="">Energies: [13000]</div> </div> <div class=""><br class=""> </div> <div class="">compile the plug-in and run over this hepmc file "<font face="Lucida Console" class=""><a href="http://cern.ch/ssen/Rivet/epos_pPb.hepmc" class="">http://cern.ch/ssen/Rivet/epos_pPb.</a>txt"</font> which has 20 proton-lead events generated with Epos-LHC. </div> <div class=""><br class=""> </div> <div class="">I use --ignore-beams option but still get the error below. </div> <div class=""><br class=""> </div> <div class=""><i class="">Event beams mismatch: [PROTON, NEUTRON] @ 13000 GeV vs. first beams [PROTON @ (6500; 0, 0, 6500) GeV, PROTON @ (6500; 0, 0, -6500) GeV] @ 13000 GeV</i></div> <div class=""><br class=""> </div> <div class="">The output of the line below shows that there are 7 p-p events and 13 p-n events in the hepmc file.</div> <div class="">grep -e "P 2 2212" epos_pPb.txt </div> <div class="">grep -e "P 2 2112" epos_pPb.txt</div> <div class=""><br class=""> </div> <div class="">I use Rivet 2.5.3, and also try the same steps with my own analysis but the error is the same. </div> <div class=""><br class=""> </div> <div class="">Thank you for your further suggestions</div> <div class="">Sercan </div> <div class=""><br class=""> </div> </div> </div> <div class=""> <div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""> <div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""> <div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""> <div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""> <div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;"> <br class=""> </div> </div> </div> </div> </div> </div> <div> <blockquote type="cite" class=""> <div class="">On Feb 28, 2017, at 2:40 PM, Andy Buckley <<a href="mailto:a.g.buckley@gmail.com" class="">a.g.buckley@gmail.com</a>> wrote:</div> <br class="Apple-interchange-newline"> <div class=""> <div class="">Hi Hans,<br class=""> <br class=""> The rivet command-line tool has an --ignore-beams option for this sort of reason. But you haven't said what analysis you are trying to run: if it's one that you're writing yourself, just use ANY rather than PROTON in the Beams specification in the analysis' .info file.<br class=""> <br class=""> We've been upgrading tools for HI and cosmic-ray analysis over the last year, so are hoping to start building up the collection. Very happy to help with anything that's currently a blocker!<br class=""> <br class=""> Andy<br class=""> <br class=""> <br class=""> On 28/02/17 11:10, Hans Jozef H Van Haevermaet wrote:<br class=""> <blockquote type="cite" class="">Dear Rivet authors,<br class=""> <br class=""> I am a member of CMS (forward and small-x community in fact), and I’m<br class=""> trying to get Rivet running with proton-lead events coming from the EPOS<br class=""> MC generator.<br class=""> It turns out however, that this is not possible (?).<br class=""> When executing Rivet, it stops because it detected that the incoming<br class=""> beams are different for each event:<br class=""> Event beams mismatch: [PROTON, PROTON] @ 1.87656 GeV vs. first beams<br class=""> [NEUTRON @ 2511.57 GeV, NEUTRON @ 2511.57 GeV] @ 1.87914 GeV<br class=""> <br class=""> This in turn, is caused by the fact that the EPOS generator writes out<br class=""> the interacting nucleon (p or n) in the HEPMC file, and not the particle<br class=""> id code of lead itself.<br class=""> <br class=""> Could you let us know if it is possible for Rivet to get compatible with<br class=""> these kind of situations?<br class=""> One direct option (for me) is to comment out this function:<br class=""> <a href="https://rivet.hepforge.org/code/dev/AnalysisHandler_8cc_source.html" class="">https://rivet.hepforge.org/code/dev/AnalysisHandler_8cc_source.html</a><br class=""> At line 95.<br class=""> <br class=""> But that can lead to unwanted behaviour in other cases, I guess. Maybe<br class=""> one can include a parameter to turn this check on/off? So that it is the<br class=""> user’s responsibility? i.e. by default the check is done, but the user<br class=""> can turn it off manually if needed.<br class=""> <br class=""> It would be very good to solve this incompatibility, since in CMS we try<br class=""> to get Rivet routines from all our analyses, and it would be nice to be<br class=""> able to extend this to proton-lead, or even lead-lead measurements. As<br class=""> far as I can see on the web, I did not find any Heavy Ion analysis<br class=""> implemented in Rivet?<br class=""> <br class=""> Thanks!<br class=""> Best Regards,<br class=""> Hans<br class=""> <br class=""> <br class=""> <br class=""> <br class=""> _______________________________________________<br class=""> Rivet mailing list<br class=""> Rivet@projects.hepforge.org<br class=""> https://www.hepforge.org/lists/listinfo/rivet<br class=""> <br class=""> </blockquote> <br class=""> <br class=""> -- <br class=""> Dr Andy Buckley, Lecturer / Royal Society University Research Fellow<br class=""> Particle Physics Expt Group, University of Glasgow<br class=""> </div> </div> </blockquote> </div> <br class=""> </body> </html>