|
[Rivet-svn] r3778 - trunk/src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Jun 28 12:34:24 BST 2012
Author: hoeth Date: Thu Jun 28 12:34:24 2012 New Revision: 3778 Log: use straight forward squared 3-momentum in ATLAS_2012_I1091481 Modified: trunk/src/Analyses/ATLAS_2012_I1091481.cc Modified: trunk/src/Analyses/ATLAS_2012_I1091481.cc ============================================================================== --- trunk/src/Analyses/ATLAS_2012_I1091481.cc Thu Jun 28 12:25:22 2012 (r3777) +++ trunk/src/Analyses/ATLAS_2012_I1091481.cc Thu Jun 28 12:34:24 2012 (r3778) @@ -56,14 +56,13 @@ // X_j = 0.5*E_j + sum_{k=0}^{k<j}(E_k) // // pion mass; - double m_pi = 0.1396; + double m_pi = 0.1396*GeV; std::vector<double> xj; std::vector<double> Ej; foreach (const Particle& p, part) { - double pT = p.momentum().pT(); - double theta = p.momentum().theta(); - double E_j = sqrt(pow(m_pi,2) + pow(pT/sin(theta), 2)); + double p2 = p.momentum().vector3().mod2()/(GeV*GeV); + double E_j = sqrt(pow(m_pi,2) + p2); Ej.push_back(E_j); double temp = 0.5*E_j;
More information about the Rivet-svn mailing list |