Skip to content
Snippets Groups Projects
Commit 9e5ab881 authored by Susumu Oda's avatar Susumu Oda
Browse files

Disable z0*sin(theta) cut and hange p>500 MeV cut to pT>500 MeV cut in SCTLorentzMonAlg/Tool.

parent e497d24a
No related branches found
No related tags found
No related merge requests found
......@@ -160,14 +160,14 @@ StatusCode SCTLorentzMonAlg::fillHistograms(const EventContext& ctx) const {
bool passesCuts{true};
if ((dataType() == AthMonitorAlgorithm::DataType_t::cosmics) and
(trkp->momentum().mag() > 500.) and // Pt > 500MeV
(trkp->momentum().perp() > 500.) and // Pt > 500MeV
(summary->get(Trk::numberOfSCTHits) > 6) // #SCTHits >6
) {
passesCuts = true;
} else if ((track->perigeeParameters()->parameters()[Trk::qOverP] < 0.) and // use negative track only
(std::abs(perigee->parameters()[Trk::d0]) < 1.) and // d0 < 1mm
(std::abs(perigee->parameters()[Trk::z0] * sin(perigee->parameters()[Trk::theta])) < 1.) and // d0 < 1mm
(trkp->momentum().mag() > 500.) and // Pt > 500MeV
// (std::abs(perigee->parameters()[Trk::z0] * sin(perigee->parameters()[Trk::theta])) < 1.) and // z0*sin(theta) < 1mm
(trkp->momentum().perp() > 500.) and // Pt > 500MeV
(summary->get(Trk::numberOfSCTHits) > 6)// and // #SCTHits >6
) {
passesCuts = true;
......
......@@ -217,14 +217,14 @@ SCTLorentzMonTool::fillHistograms() {
bool passesCuts{true};
if ((AthenaMonManager::dataType() == AthenaMonManager::cosmics) and
(trkp->momentum().mag() > 500.) and // Pt > 500MeV
(trkp->momentum().perp() > 500.) and // Pt > 500MeV
(summary->get(Trk::numberOfSCTHits) > 6)// and // #SCTHits >6
) {
passesCuts = true;
} else if ((track->perigeeParameters()->parameters()[Trk::qOverP] < 0.) and // use negative track only
(std::abs(perigee->parameters()[Trk::d0]) < 1.) and // d0 < 1mm
(std::abs(perigee->parameters()[Trk::z0] * sin(perigee->parameters()[Trk::theta])) < 1.) and // d0 < 1mm
(trkp->momentum().mag() > 500.) and // Pt > 500MeV
// (std::abs(perigee->parameters()[Trk::z0] * sin(perigee->parameters()[Trk::theta])) < 1.) and // z0*sin(theta) < 1mm
(trkp->momentum().perp() > 500.) and // Pt > 500MeV
(summary->get(Trk::numberOfSCTHits) > 6)// and // #SCTHits >6
) {
passesCuts = true;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment