46 {
47 NUP=
e->particles().size();
49 XWGTUP =
e->weights().size() ?
e->weights()[0] : 1.0;
51 IDPRUP = A_signal_process_id ? A_signal_process_id->value() : 0;
53 SCALUP = A_event_scale ? A_event_scale->value() : 0;
55 AQCDUP = A_alphaQCD ? A_alphaQCD->value() : 0;
57 AQEDUP = A_alphaQED ? A_alphaQED->value() : 0;
58 for (
int i = 0;
i <
NUP; ++
i ){
59 const auto & thisParticle =
e->particles().at(i);
60 PUP[
i][0] = thisParticle->momentum().px();
61 PUP[
i][1] = thisParticle->momentum().py();
62 PUP[
i][2] = thisParticle->momentum().pz();
63 PUP[
i][3] = thisParticle->momentum().e();
64 PUP[
i][4] = thisParticle->momentum().m();
65 IDUP[
i] = thisParticle->pdg_id();
66 const auto pv = thisParticle->production_vertex();
67 const auto ev = thisParticle->end_vertex();
70 if (thisParticle->status() == 4 || !pv )
ISTUP[
i] = -1;
77 if (pv && !
pv->particles_in().empty()) {
78 l =
pv->particles_in().front()->id();
80 for (
const auto& p :
pv->particles_in()) {
81 const int id =
p->id();
84 }
85 }
88 }
89 return true;
90 }
std::vector< std::vector< double > > PUP
std::vector< std::pair< int, int > > ICOLUP
std::vector< std::pair< int, int > > MOTHUP
const std::string event_scale
const std::string signal_process_id
const std::string alphaQED
const std::string alphaQCD
l
Printing final latex table to .tex output file.