|
[Rivet] error when running make-plots comparing to standard published analysisRaghav Kunnawalkam Elayavalli raghav.k.e at cern.chWed Apr 27 12:13:10 BST 2016
Hi Holger, I got a chance to talk with David here and he quickly figured out what was going on. when comparing with Data, I was explicitly giving it two yoda files, one from me and one from datadir for the analysis. Due to that rivet was putting an additional tag. I just realized after David explained that it automatically takes the data points so now it works perfectly. Thanks for checking and sorry for the (silly) issue :). Cheers Raghav On Apr 25, 2016, at 10:37 PM, Holger Schulz <holger.schulz at durham.ac.uk<mailto:holger.schulz at durham.ac.uk>> wrote: Hi Raghav, I am really sorry but I cannot reproduce the error you see. Can you maybe try to reproduce it independently on another machine you have access to? I am really puzzled by now as the code that generates the YODA files with the weight vector subscripts ([WeightName]) is only existing in our Rivet 3 alpha branch. Holger On 21/04/16 17:19, Raghav Kunnawalkam Elayavalli wrote: Ive never had any issue with running rivet on my hepmc files as well. ahh i see. Here are the bootstrap scripts in a zip file Cheers Raghav On Apr 21, 2016, at 6:16 PM, Holger Schulz <holger.schulz at durham.ac.uk<mailto:holger.schulz at durham.ac.uk>> wrote: That worked, I can confirm that I can run on your hepmc files and make plots without a problem. I was trying to play with you bootstrap but that also got mangled during the mail transfer. Could you please send me a zip file of those bootstrap scripts as well? Thanks, Holger On 21/04/16 16:02, Raghav Kunnawalkam Elayavalli wrote: Hi Holger, Regarding the version difference in the HepMC, we can look into using the latest version but it I dont think its an issue. Im also able to run on the HepMC file locally pb-d-128-141-166-58:JEWEL raghav$ rivet -a JEWEL_YpJet JEWEL_7TeV_YplusJet_pp.hepmc -H JEWEL_7TeV_YplusJet_pp.yoda Rivet 2.4.1 running on machine pb-d-128-141-166-58.cern.ch<http://pb-d-128-141-166-58.cern.ch/> (x86_64) at 2016-04-21 16:59:20 Rivet.Analysis.Handler: WARN Analysis 'JEWEL_YpJet' is unvalidated: be careful, it may be broken! Reading events from 'JEWEL_7TeV_YplusJet_pp.hepmc' #-------------------------------------------------------------------------- # FastJet release 3.2.0 # M. Cacciari, G.P. Salam and G. Soyez # A software package for jet finding and analysis at colliders # <http://fastjet.fr/> http://fastjet.fr<http://fastjet.fr/> # # Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package # for scientific work and optionally PLB641(2006)57 [hep-ph/0512210]. # # FastJet is provided without warranty under the terms of the GNU GPLv2. # It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code # and 3rd party plugin jet algorithms. See COPYING file for details. #-------------------------------------------------------------------------- Finished event loop at 2016-04-21 16:59:20 Cross-section = 2.951070e+06 pb Rivet.Analysis.Handler: INFO Finalising analyses Rivet.Analysis.Handler: INFO Processed 10 events The MCnet usage guidelines apply to Rivet: see <http://www.montecarlonet.org/GUIDELINES> http://www.montecarlonet.org/GUIDELINES Please acknowledge plots made with Rivet analyses, and cite arXiv:1003.0694 (http://arxiv.org/abs/1003.0694) Rivet run completed at 2016-04-21 16:59:20, time elapsed = 0:00:00 Histograms written to /Users/raghav/WORK/JEWEL/JEWEL_7TeV_YplusJet_pp.yoda pb-d-128-141-166-58:JEWEL raghav$ probably the email made it weird. here is the zipped version. hopefully that works. Cheers Raghav On Apr 21, 2016, at 4:57 PM, Holger Schulz <<mailto:holger.schulz at durham.ac.uk>holger.schulz at durham.ac.uk<mailto:holger.schulz at durham.ac.uk>> wrote: Hi Raghav, there are some modifications in your boostrap that I am not immediately able to rule out as the problem. What I saw is that you build rivet using HepMC 2.06.09 which is fine, you events however are stored in HepMC 2.06.05 How hard would it be for you to relink JEWEL against HepMC 2.06.09 and produce some events? Also I am currently not able to process you hepmc events with rivet 2.4.1: streaming input: end of stream found setting badbit maybe some illformatting happened when mailing. Could you please send me these files again but as a tarball or zip or similar? Thank you, Holger On 21/04/16 15:35, Raghav Kunnawalkam Elayavalli wrote: Hi Holger, Since this email contains attachments, im sending it just to you (didnt want to spam the list). Here are the bootstrap scripts and a couple of hepmc files for comparing with two published analysis. local and lxplus bootstraps: hepmc files: the YplusJet is compared to ATLAS_2012_I1093738 the ZplusJet is compared to ATLAS_2013_I1230812_MU Cheers Raghav On Apr 21, 2016, at 4:27 PM, Holger Schulz <<mailto:holger.schulz at durham.ac.uk>holger.schulz at durham.ac.uk<mailto:holger.schulz at durham.ac.uk>> wrote: Ok cool, so for me to reproduce this behaviour, would you mind sending me both boostrap scripts and a few hepmc events please? Thanks, Holger On 21/04/16 15:22, Raghav Kunnawalkam Elayavalli wrote: Hi Holger, I setup rivet using the latest bootstrap version on my laptop and the lcg bootstrap in the lxplus machine. its just using the information provided here <https://rivet.hepforge.org/trac/wiki/GettingStarted> https://rivet.hepforge.org/trac/wiki/GettingStarted The events i get from my MC generator called JEWEL. thats what im comparing with published ATLAS data points. Cheers Raghav On Apr 21, 2016, at 4:19 PM, Holger Schulz <<mailto:holger.schulz at durham.ac.uk>holger.schulz at durham.ac.uk<mailto:holger.schulz at durham.ac.uk>> wrote: Hi Raghav, do you think you can figure out how rivet was set up on your system? It looks a bit like a mixture of rivet 3 alpha and the latest release. Also, where to you get your events from? Thanks, Holger On 20/04/16 14:49, Raghav Kunnawalkam Elayavalli wrote: Hi Chris, im running the latest stable version i guess. pb-d-128-141-166-190:JEWEL raghav$ rivet -v Rivet 2.4.1 running on machine pb-d-128-141-162-92.cern.ch<http://pb-d-128-141-162-92.cern.ch/> (x86_64) at 2016-04-20 15:48:54 Cheers Raghav On Apr 20, 2016, at 3:46 PM, Chris Pollard <<mailto:chris.pollard at cern.ch>chris.pollard at cern.ch<mailto:chris.pollard at cern.ch>> wrote: Hi Raghav, Can you tell us which version of rivet you are using? This looks like you might be using the Rivet-3.0alpha. I don't think you should be getting the [0]s in any full release. Chris On Wed, Apr 20, 2016 at 1:16 PM, Raghav Kunnawalkam Elayavalli <<mailto:raghav.k.e at cern.ch>raghav.k.e at cern.ch<mailto:raghav.k.e at cern.ch>> wrote: Hi, Sorry for the multiple emails but either im doing something wrong while im comparing or something weird is going on. I just remade the comparison again to ATLAS_2013_I1230812_MU and in the .dat files the data histograms have this weird [0] which is causing the problems when make-plots is run # BEGIN HISTO1D /ATLAS_2013_I1230812_MU/d05-x01-y03[0] Path=/ATLAS_2013_I1230812_MU/d05-x01-y03[0] thanks for taking a look! Cheers Raghav On Apr 20, 2016, at 11:08 AM, Raghav Kunnawalkam Elayavalli <<mailto:raghav.k.e at cern.ch>raghav.k.e at cern.ch<mailto:raghav.k.e at cern.ch>> wrote: Hi Rivet experts, I found out why ATLAS_2012_I1093738 was causing a problem when doing my comparisons. in the .dat files generated when doing compare histos command, this histogram location had a spurious [0] after the name, for example it should be /ATLAS_2012_I1093738/d01-x01-y01 when it was /ATLAS_2012_I1093738/d01-x01-y01/[0] If i removed the [0] by hand then im able to make the comparison plots. Im however not able to understand the error for the comparison with ATLAS_2014_I1319490 thought. the .dat files look fine and i cant see anything weird. Cheers Raghav On Apr 19, 2016, at 3:30 PM, Raghav Kunnawalkam Elayavalli <<mailto:raghav.k.e at cern.ch>raghav.k.e at cern.ch<mailto:raghav.k.e at cern.ch>> wrote: Im also having a similar issue when i compare with ATLAS_2014_I1319490 [rkunnawa at lxplus047 Validation]$ rivet-mkhtml myMC.yoda ATLAS_2014_I1319490.yoda Making 34 plots Traceback (most recent call last): File "/afs/cern.ch/work/r/rkunnawa/JEWEL/local/bin/make-plots<http://cern.ch/work/r/rkunnawa/JEWEL/local/bin/make-plots>", line 2607, in <module> for i, _ in enumerate(pool.imap(process_datfile, datfiles)): File "/cvmfs/sft.cern.ch/lcg/releases/LCG_84/Python/2.7.10/x86_64-slc6-gcc48-opt/lib/python2.7/multiprocessing/pool.py<http://sft.cern.ch/lcg/releases/LCG_84/Python/2.7.10/x86_64-slc6-gcc48-opt/lib/python2.7/multiprocessing/pool.py>", line 668, in next raise value ValueError: invalid literal for float(): 0$ Thanks for taking a look! Cheers Raghav On Apr 19, 2016, at 3:25 PM, Raghav Kunnawalkam Elayavalli <<mailto:raghav.k.e at cern.ch>raghav.k.e at cern.ch<mailto:raghav.k.e at cern.ch>> wrote: Hi Rivet Experts, Im getting a python error when comparing to a published analysis and it seems to be from make-plots. (rivet-cmphistos works fine). this is the published analysis ATLAS_2012_I1093738 and im comparing it with output from my MC. This is the chain of commands i followed: $ rivet -a ATLAS_2012_I1093738 myMC.hepmc -H myMC.yoda $ rivet-mkhtml myMC.yoda ATLAS_2012_I1093738.yoda Making 6 plots Traceback (most recent call last): File "/afs/cern.ch/work/r/rkunnawa/JEWEL/local/bin/make-plots<http://cern.ch/work/r/rkunnawa/JEWEL/local/bin/make-plots>", line 2607, in <module> for i, _ in enumerate(pool.imap(process_datfile, datfiles)): File "/cvmfs/sft.cern.ch/lcg/releases/LCG_84/Python/2.7.10/x86_64-slc6-gcc48-opt/lib/python2.7/multiprocessing/pool.py<http://sft.cern.ch/lcg/releases/LCG_84/Python/2.7.10/x86_64-slc6-gcc48-opt/lib/python2.7/multiprocessing/pool.py>", line 668, in next raise value KeyError: '/ATLAS_2012_I1093738/d01-x01-y01' Im getting the same error when i run it on my local computer as well. Any ideas whats going on here? thanks! Cheers Raghav _______________________________________________ Rivet mailing list <mailto:Rivet at projects.hepforge.org>Rivet at projects.hepforge.org<mailto:Rivet at projects.hepforge.org> <https://www.hepforge.org/lists/listinfo/rivet>https://www.hepforge.org/lists/listinfo/rivet _______________________________________________ Rivet mailing list Rivet at projects.hepforge.org<mailto:Rivet at projects.hepforge.org> https://www.hepforge.org/lists/listinfo/rivet -------------- next part -------------- An HTML attachment was scrubbed... URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20160427/1d482f56/attachment.html>
More information about the Rivet mailing list |