<div dir="ltr"><div class="gmail_default" style="font-size:large">Hi Adam,</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">There isn't a convenient mechanism for this, because a) it's unobservable, and b) it depends on what you consider to be the hard-scattering process.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">To elaborate on the second of those, let's say you define s_hat from the 4-vectors of the quarks entering the hard-scattering vertex in the event record. Now use that definition on two events with the same kinematics: a ttj event in which the extra jet comes from the matrix element, and an equivalent event in which the jet comes from the initial-state parton shower. You'd see a different value even though the physics is the same (in this idealised example), just because the calculation has been broken up differently: that would clearly not be a good physical definition. I'm not sure if a definition that used the partons "emerging" from the initial-state protons would be any less problematic.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Rivet's main aim is to emphasise well-defined observables rather than intermediate quantities in calculations, so this kind of quantity is not "automatically" available -- you need to manually access the event graph to compute whatever definition you choose. You can get to that via the Event::genEvent() function. </div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Andy</div><div class="gmail_default" style="font-size:large"> </div></div><div class="gmail_extra"><br><div class="gmail_quote">On 4 December 2017 at 12:39, Adam Makda <span dir="ltr"><<a href="mailto:adam.makda@durham.ac.uk" target="_blank">adam.makda@durham.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello,<div><br></div><div>I'm a student relatively new to Rivet and C++ in general.</div><div><br></div><div>I'm currently using Rivet to analyse pp -> t tbar (stable top quark pair production at LHC13). I want to be able to project out the partons' (within the proton) centre of mass energy. This is usually denoted as $\sqrt{\hat{s}}$ in LaTeX. </div><div><br></div><div>I'm aware of a projector called 'Beams' but this only projects out $\sqrt{s}$ ie. in my simulations, 13 TeV.  </div><div><br></div><div>I was wondering if there is an example within the existing analyses I could use as a template in order to project out the initial state partons in an event and extract their energies?</div><div><br></div><div>Any help will be appreciated. If this isn't the purpose of this contact email address, is there anyone who could help?</div><div><br></div><div>Thanks</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Adam Makda</div></font></span></div>
<br>______________________________<wbr>_________________<br>
Rivet mailing list<br>
<a href="mailto:Rivet@projects.hepforge.org">Rivet@projects.hepforge.org</a><br>
<a href="https://www.hepforge.org/lists/listinfo/rivet" rel="noreferrer" target="_blank">https://www.hepforge.org/<wbr>lists/listinfo/rivet</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Dr Andy Buckley, Lecturer / Royal Society University Research Fellow<br>Particle Physics Expt Group, University of Glasgow<br></div>
</div>