Skip to content

Adding Xi-compatible tracks to lostTracks collection #33891

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
cmsbuild merged 9 commits into from
Jun 3, 2021
Merged

Adding Xi-compatible tracks to lostTracks collection #33891

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
4e8135d
Adding tracks compatible with a Xi to lost tracks
May 28, 2021
28116a9
Code checks and code formatting
May 28, 2021
40e810e
Lowering cut to 1.5
May 29, 2021
4f044f8
Moving to XYZTLorentz vectors
May 31, 2021
93bcbbe
Header cleanup
May 31, 2021
592c219
Code format
May 31, 2021
743f900
Adding sqrt for energy mass conversion
May 31, 2021
830b6f4
Moving the NOTUSED cut above
May 31, 2021
807d8e0
Code format
May 31, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
23 changes: 22 additions & 1 deletion PhysicsTools/PatAlgos/plugins/PATLostTracks.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -68,6 +68,8 @@ namespace pat {
const double maxDxyForNotReconstructedPrimary_;
const double maxDxySigForNotReconstructedPrimary_;
const bool useLegacySetup_;
const bool xiSelection_;
const double xiMassCut_;
};
} // namespace pat

Expand Down Expand Up @@ -100,7 +102,9 @@ pat::PATLostTracks::PATLostTracks(const edm::ParameterSet& iConfig)
.getParameter<double>("maxDxyForNotReconstructedPrimary")),
maxDxySigForNotReconstructedPrimary_(iConfig.getParameter<edm::ParameterSet>("pvAssignment")
.getParameter<double>("maxDxySigForNotReconstructedPrimary")),
useLegacySetup_(iConfig.getParameter<bool>("useLegacySetup")) {
useLegacySetup_(iConfig.getParameter<bool>("useLegacySetup")),
xiSelection_(iConfig.getParameter<bool>("xiSelection")),
xiMassCut_(iConfig.getParameter<double>("xiMassCut")) {
std::vector<std::string> trkQuals(iConfig.getParameter<std::vector<std::string>>("qualsToAutoAccept"));
std::transform(
trkQuals.begin(), trkQuals.end(), std::back_inserter(qualsToAutoAccept_), reco::TrackBase::qualityByName);
Expand Down Expand Up @@ -190,10 +194,27 @@ void pat::PATLostTracks::produce(edm::StreamID, edm::Event& iEvent, const edm::E
}
}
for (const auto& v0 : *lambdas) {
double protonCharge = 0;
for (size_t dIdx = 0; dIdx < v0.numberOfDaughters(); dIdx++) {
size_t key = (dynamic_cast<const reco::RecoChargedCandidate*>(v0.daughter(dIdx)))->track().key();
if (trkStatus[key] == TrkStatus::NOTUSED)
trkStatus[key] = TrkStatus::VTX;
protonCharge += v0.daughter(dIdx)->charge() * v0.daughter(dIdx)->momentum().mag2();
}
if (xiSelection_) {
// selecting potential Xi- -> Lambda pi candidates
math::XYZTLorentzVector p4Lambda(v0.px(), v0.py(), v0.pz(), sqrt(v0.momentum().mag2() + v0.mass() * v0.mass()));

for (unsigned int trkIndx = 0; trkIndx < tracks->size(); trkIndx++) {
reco::TrackRef trk(tracks, trkIndx);
if ((*trk).charge() * protonCharge < 0 || trkStatus[trkIndx] != TrkStatus::NOTUSED)
continue;

math::XYZTLorentzVector p4pi(
trk->px(), trk->py(), trk->pz(), sqrt(trk->momentum().mag2() + 0.01947995518)); // pion mass ^2
if ((p4Lambda + p4pi).M() < xiMassCut_)
trkStatus[trkIndx] = TrkStatus::VTX; // selecting potential Xi- candidates
}
}
}
std::vector<int> mapping(tracks->size(), -1);
Expand Down
4 changes: 3 additions & 1 deletion PhysicsTools/PatAlgos/python/slimming/lostTracks_cfi.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,7 +22,9 @@
minPtToStoreLowQualityProps = cms.double(0.0),
passThroughCut = cms.string("pt>2"),
pvAssignment = primaryVertexAssociation.assignment,
useLegacySetup = cms.bool(False) #When True: check only if track used to fit vertex[0] and do not store track detailed info for Pt between 0.5 and minPtToStoreProps GeV
useLegacySetup = cms.bool(False), #When True: check only if track used to fit vertex[0] and do not store track detailed info for Pt between 0.5 and minPtToStoreProps GeV
xiSelection = cms.bool(True),
xiMassCut = cms.double(1.5)
)
from Configuration.Eras.Modifier_phase1Pixel_cff import phase1Pixel
phase1Pixel.toModify(lostTracks, covarianceVersion =1 )
Expand Down