|
[Rivet] turn off "Event beams mismatch error"Andy Buckley andy.buckley at cern.chThu Jun 22 16:23:12 BST 2017
Hi all, Yang, are you trying the latest version of Rivet, 2.5.4? We released this just a couple of days ago, and I believe that it includes a relaxation of the event-to-event beam check if --ignore-beams is given, for exactly this reason. Andy On 22/06/17 14:36, Chris Pollard wrote: > Hi Yang, > > I personally am not sure of the "right way" to do this, and I think > several of the other devs are at a workshop this week. I think it'd be > best to ping the list again next week. > > Cheers, > > Chris > > On Wed, Jun 14, 2017 at 6:52 PM, Yang Ting Chien <ytchien at mit.edu > <mailto:ytchien at mit.edu>> wrote: > > Hi Chris, > > Sounds good. Could you help put this feature in the new rivet > version? Presumably the combination is following a Poisson > distribution, which is handled by the simulations. Or can it be as > simple as having a rivet option to not checking beams at all? The > hepmc file is likely to be lists of events with either proton or > neutron colliding, and the numbers of each combination is not fixed > because of probability. > > Thanks, > Yang-Ting > > ------------------------------------------------------------------------ > *From:* Chris Pollard [cpollard at cern.ch <mailto:cpollard at cern.ch>] > *Sent:* Wednesday, June 14, 2017 10:47 AM > *To:* Yang Ting Chien > *Cc:* Chris Pollard; rivet at projects.hepforge.org > <mailto:rivet at projects.hepforge.org> > > *Subject:* Re: [Rivet] turn off "Event beams mismatch error" > > Hi Yang, > > I guess the simplest solution would be to define a "mixed" beam > type: then analyses could accept specified combinations of particles. > > Chris > > On Wed, Jun 14, 2017 at 3:18 PM, Yang Ting Chien <ytchien at mit.edu > <mailto:ytchien at mit.edu>> wrote: > > Hi Chris, > > Thanks! > > The second comment explains my situation. The colliding beams > are lead nuclei, so the binary collisions can be PP, PN and NN > in a heavy ion run. > > Assuming simple binary collisions with an optical Glauber model, > I can weigh PP, PN and NN with the geometric probability. > However, it will be extremely helpful to the heavy ion physics > community if rivet can have an option for analyzing heavy ion > data more easily. I appreciate the great help! > > Best, > Yang-Ting > > LHC Theory Initiative Fellow > MIT Center for Theoretical Physics > 617-253-7194 <tel:(617)%20253-7194> > ------------------------------------------------------------------------ > *From:* Chris Pollard [cpollard at cern.ch <mailto:cpollard at cern.ch>] > *Sent:* Wednesday, June 14, 2017 5:09 AM > *To:* Yang Ting Chien > *Cc:* Deepak Kar; rivet at projects.hepforge.org > <mailto:rivet at projects.hepforge.org> > > *Subject:* Re: [Rivet] turn off "Event beams mismatch error" > > Hi, > > I believe the intended behavior is this: > > 1) --ignore-beams won't require the beam in your hepmc file to > match the beam type expected by the analysis. > > 2) --ignore-beams still won't allow different beam types to be > analyzed in the same run, because this is usually quite a > strange thing to want to do. > > Unfortunately I'm not sure the "right" way to achieve what you > want. Is it possible to split the events into two separate runs, > one with P-P, one with P-N, and one with N-N collisions? > > Chris > > On Wed, Jun 14, 2017 at 5:24 AM, Yang Ting Chien > <ytchien at mit.edu <mailto:ytchien at mit.edu>> wrote: > > Dear Deepak, > > Thanks for the prompt reply! I tried adding the option > --ignore-beams when running rivet but it still did not work. > > I did > > rivet-buildplugin RivetJEWEL_test.so JEWEL_test.cc > > and then > > rivet --ignore-beams --analysis=JEWEL_test test.hepmc > > > and I still get something like > > > Event 100 (0:00:00 elapsed) > > Event 200 (0:00:00 elapsed) > > Event 300 (0:00:00 elapsed) > > Event 400 (0:00:01 elapsed) > > Event 500 (0:00:01 elapsed) > > Event 600 (0:00:01 elapsed) > > Event 700 (0:00:01 elapsed) > > Event 800 (0:00:01 elapsed) > > Event 900 (0:00:01 elapsed) > > Event 1000 (0:00:01 elapsed) > > Event 1100 (0:00:02 elapsed) > > Event 1200 (0:00:02 elapsed) > > Event 1300 (0:00:02 elapsed) > > Event 1400 (0:00:02 elapsed) > > Event 1500 (0:00:02 elapsed) > > Event beams mismatch: [PROTON, NEUTRON] @ 5020 GeV vs. first > beams [PROTON @ (2510; 0, 0, 2510) GeV, PROTON @ (2510; 0, > 0, -2510) GeV] @ 5020 GeV > > > It looks like rivet ignores the ignore-beams option... > > Thanks, > Yang-Ting > > LHC Theory Initiative Fellow > MIT Center for Theoretical Physics > 617-253-7194 <tel:(617)%20253-7194> > ------------------------------------------------------------------------ > *From:* Deepak Kar [deepak.kar at cern.ch > <mailto:deepak.kar at cern.ch>] > *Sent:* Tuesday, June 13, 2017 5:51 PM > *To:* Yang Ting Chien > *Cc:* rivet at projects.hepforge.org > <mailto:rivet at projects.hepforge.org> > *Subject:* Re: [Rivet] turn off "Event beams mismatch error" > > --ignore-beams > > On Tue, Jun 13, 2017 at 10:21 PM, Yang Ting Chien > <ytchien at mit.edu <mailto:ytchien at mit.edu>> wrote: > > Dear Rivet authors, > > I am starting to use Rivet to analyze some hepmc files I > generated using jewel, a heavy ion collision generator. > Since protons and neutrons are both present in a heavy > nuclei, I have beams of protons and neutrons present in > the hepmc event file. However, rivet detects the beam > mismatch and gives the following error message and exits: > > Event beams mismatch: [PROTON, NEUTRON] @ 5020 GeV vs. > first beams [PROTON @ (2510; 0, 0, 2510) GeV, PROTON @ > (2510; 0, 0, -2510) GeV] @ 5020 GeV > > I couldn't find a way to successfully ask rivet to skip > this beam check. Could you help me resolve this issue? I > attached an analysis I will extend to include jet > substructure studies. > > Thank a lot, > Yang-Ting > > LHC Theory Initiative Fellow > MIT Center for Theoretical Physics > 617-253-7194 <tel:(617)%20253-7194> > > _______________________________________________ > Rivet mailing list > Rivet at projects.hepforge.org > <mailto:Rivet at projects.hepforge.org> > https://www.hepforge.org/lists/listinfo/rivet > > -- > Deepak Kar > University of Witwatersrand > Room PM15, School of Physics > (0027) 011-7176958 (office) (0027) 0736944181 (mobile) > > While at CERN: > Building 1, R-016 > (0041) 0767321349 (mobile) > > While at USA: > (001) 330-998-1500 (mobile) > <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 > <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 > <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 > <https://www.hepforge.org/lists/listinfo/rivet> > > > > > _______________________________________________ > Rivet mailing list > Rivet at projects.hepforge.org > https://www.hepforge.org/lists/listinfo/rivet > -- Dr Andy Buckley, Lecturer / Royal Society University Research Fellow Particle Physics Expt Group, University of Glasgow
More information about the Rivet mailing list |