Commit 5fff288e authored by Sylvester Joosten's avatar Sylvester Joosten
Browse files

Add pseudorapidity cut for truth track seeding to reduce failed tracking attempts

parent ef18290b
......@@ -52,6 +52,8 @@ namespace Jug::Reco {
Gaudi::Property<double> m_maxVertexY{this, "maxVertexY", 80. * Gaudi::Units::mm};
Gaudi::Property<double> m_maxVertexZ{this, "maxVertexZ", 200. * Gaudi::Units::mm};
Gaudi::Property<double> m_minMomentum{this, "minMomentum", 100. * Gaudi::Units::MeV};
Gaudi::Property<double> m_maxEtaForward{this, "maxEtaForward", 4.0};
Gaudi::Property<double> m_maxEtaBackward{this, "maxEtaBackward", 4.1};
SmartIF<IParticleSvc> m_pidSvc;
......@@ -114,6 +116,14 @@ namespace Jug::Reco {
continue;
}
// require minimum pseudorapidity
if (part.ps().eta() > m_maxEtaForward || part.ps().eta() < -std::abs(m_maxEtaBackward)) {
if (msgLevel(MSG::DEBUG)) {
debug() << "ignoring particle with Eta = " << part.ps().eta() << endmsg;
}
continue;
}
// get the particle charge
// note that we cannot trust the mcparticles charge, as DD4hep
// sets this value to zero! let's lookup by PDGID instead
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment