diff --git a/PhysicsAnalysis/NtupleDumper/src/NtupleDumperAlg.cxx b/PhysicsAnalysis/NtupleDumper/src/NtupleDumperAlg.cxx index 995b69ca66ab38c7353b83ebfd0b14b7cc4e8d6e..8b6aba7341e4c4c01aa752d8e6acbda13eb053cd 100644 --- a/PhysicsAnalysis/NtupleDumper/src/NtupleDumperAlg.cxx +++ b/PhysicsAnalysis/NtupleDumper/src/NtupleDumperAlg.cxx @@ -603,40 +603,40 @@ StatusCode NtupleDumperAlg::execute(const EventContext &ctx) const if (ipart++ > 9) break; m_truth_P.push_back(particle->p4().P()); - m_truth_px.push_back(particle->p4().X()); - m_truth_py.push_back(particle->p4().Y()); - m_truth_pz.push_back(particle->p4().Z()); - m_truth_m.push_back(particle->m()); - m_truth_pdg.push_back(particle->pdgId()); - - if ( particle->hasProdVtx()) { - m_truth_prod_x.push_back(particle->prodVtx()->x()); - m_truth_prod_y.push_back(particle->prodVtx()->y()); - m_truth_prod_z.push_back(particle->prodVtx()->z()); - } else { - m_truth_prod_x.push_back(999999); - m_truth_prod_y.push_back(999999); - m_truth_prod_z.push_back(999999); - } - - if ( particle->hasDecayVtx()) { - m_truth_dec_x.push_back(particle->decayVtx()->x()); - m_truth_dec_y.push_back(particle->decayVtx()->y()); - m_truth_dec_z.push_back(particle->decayVtx()->z()); - } else { - m_truth_dec_x.push_back(999999); - m_truth_dec_y.push_back(999999); - m_truth_dec_z.push_back(999999); - } + m_truth_px.push_back(particle->p4().X()); + m_truth_py.push_back(particle->p4().Y()); + m_truth_pz.push_back(particle->p4().Z()); + m_truth_m.push_back(particle->m()); + m_truth_pdg.push_back(particle->pdgId()); + + if ( particle->hasProdVtx()) { + m_truth_prod_x.push_back(particle->prodVtx()->x()); + m_truth_prod_y.push_back(particle->prodVtx()->y()); + m_truth_prod_z.push_back(particle->prodVtx()->z()); + } else { + m_truth_prod_x.push_back(999999); + m_truth_prod_y.push_back(999999); + m_truth_prod_z.push_back(999999); + } + + if ( particle->hasDecayVtx()) { + m_truth_dec_x.push_back(particle->decayVtx()->x()); + m_truth_dec_y.push_back(particle->decayVtx()->y()); + m_truth_dec_z.push_back(particle->decayVtx()->z()); + } else { + m_truth_dec_x.push_back(999999); + m_truth_dec_y.push_back(999999); + m_truth_dec_z.push_back(999999); + } // Find the M d0 and d1 truth information for dark photon if ( particle->barcode() == 1 || particle->barcode() == 2 || particle->barcode() == 3 ) { auto positions = m_fiducialParticleTool->getTruthPositions(particle->barcode()); if ( particle->pdgId() == 32) { // mother particle (A') - m_truthM_P.push_back(particle->p4().P()); + m_truthM_P.push_back(particle->p4().P()); m_truthM_px.push_back(particle->p4().X()); - m_truthM_py.push_back(particle->p4().Y()); - m_truthM_pz.push_back(particle->p4().Z()); + m_truthM_py.push_back(particle->p4().Y()); + m_truthM_pz.push_back(particle->p4().Z()); if ( particle->hasDecayVtx()) { // decay vertex for A' particle m_truthM_x.push_back(particle->decayVtx()->x());