Skip to content
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

Jump to bottom

Fix obvious bug in particle boosting - backport of PR #25412 in 10.2.X for actual MC generation #25424

Merged
merged 2 commits into from Dec 9, 2018
Merged

Fix obvious bug in particle boosting - backport of PR #25412 in 10.2.X for actual MC generation #25424

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
95d073a
fix obvious bug in particle boosting (recursive statement)
covarell Dec 4, 2018
d85b5a5
small code cleanup
covarell Dec 4, 2018
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
14 changes: 4 additions & 10 deletions GeneratorInterface/GenFilters/src/MCParticlePairFilter.cc
View file
Open in desktop
Expand Up @@ -140,8 +140,7 @@ bool MCParticlePairFilter::filter(edm::StreamID, edm::Event& iEvent, const edm::
double eta1 = mom.eta();
double eta2;
double deltaR;
//HepLorentzVector momentum1 = (*p)->momentum();
//HepLorentzVector totmomentum;

double tot_x= 0.;
double tot_y= 0.;
double tot_z= 0.;
Expand All @@ -155,9 +154,7 @@ bool MCParticlePairFilter::filter(edm::StreamID, edm::Event& iEvent, const edm::
tot_z=mom.pz();
tot_e=mom.e();
charge1 = charge((*p)->pdg_id());
//totmomentum = momentum1 + typeBpassed[i]->momentum();
//invmass = totmomentum.m();
HepMC::FourVector mom_i = MCFilterZboostHelper::zboost(typeBpassed[i]->momentum(),betaBoost);
HepMC::FourVector mom_i = MCFilterZboostHelper::zboost(typeBpassed[i]->momentum(),betaBoost);
tot_x += mom_i.px();
tot_y += mom_i.py();
tot_z += mom_i.pz();
Expand Down Expand Up @@ -212,8 +209,7 @@ bool MCParticlePairFilter::filter(edm::StreamID, edm::Event& iEvent, const edm::
double eta1 = mom.eta();
double eta2;
double deltaR;
//HepLorentzVector momentum1 = (*p)->momentum();
//HepLorentzVector totmomentum;

double tot_x= 0.;
double tot_y= 0.;
double tot_z= 0.;
Expand All @@ -228,9 +224,7 @@ bool MCParticlePairFilter::filter(edm::StreamID, edm::Event& iEvent, const edm::
tot_z=mom.pz();
tot_e=mom.e();
charge1 = charge((*p)->pdg_id());
//totmomentum = momentum1 + typeApassed[i]->momentum();
//invmass = totmomentum.m();
HepMC::FourVector mom_i = MCFilterZboostHelper::zboost(mom_i,betaBoost);
HepMC::FourVector mom_i = MCFilterZboostHelper::zboost(typeApassed[i]->momentum(),betaBoost);
tot_x += mom_i.px();
tot_y += mom_i.py();
tot_z += mom_i.pz();
Expand Down