From 0ade249ea849d38d9ae2fe4e51392d89dae4b54a Mon Sep 17 00:00:00 2001 From: Colby Thompson <52970623+ctdax@users.noreply.github.com> Date: Fri, 8 Nov 2024 11:38:08 -0500 Subject: [PATCH] Fixed FullModelHadronicProcess.cc energy deposit model A 2015 commit to FullModelHadronicProcess.cc (https://github.com/cms-sw/cmssw/commit/0f7fe9f2e7b0d63c51c023862a2eeb75d9898e1e#diff-d97c0182cf0b0204534930064c29ace096341c2b2a33ede9d1f1aac75be22d91) introduced a bug to the energy deposit model. This bug incorrectly deposited energy values close to the BSM particles mass into the simulated detector. --- SimG4Core/CustomPhysics/src/FullModelHadronicProcess.cc | 5 +---- 1 file changed, 1 insertion(+), 4 deletions(-) diff --git a/SimG4Core/CustomPhysics/src/FullModelHadronicProcess.cc b/SimG4Core/CustomPhysics/src/FullModelHadronicProcess.cc index 61412807f58e5..2a58441a5c5c9 100644 --- a/SimG4Core/CustomPhysics/src/FullModelHadronicProcess.cc +++ b/SimG4Core/CustomPhysics/src/FullModelHadronicProcess.cc @@ -354,10 +354,7 @@ G4VParticleChange* FullModelHadronicProcess::PostStepDoIt(const G4Track& aTrack, G4LorentzVector p_g_cms = gluinoMomentum; //gluino in CMS BEFORE collision p_g_cms.boost(trafo_full_cms); - double e = cloud_p4_new.e() + gluinoMomentum.e(); - if (outgoingRhadron) - e += outgoingRhadron->GetPDGMass(); - G4LorentzVector p4_new(cloud_p4_new.v() + gluinoMomentum.v(), e); + G4LorentzVector p4_new(cloud_p4_new.v() + gluinoMomentum.v(), outgoingRhadron->GetPDGMass()); // G4cout<<"P4-diff: "<<(p4_new-cloud_p4_new-gluinoMomentum)/GeV<<", magnitude: " // <<(p4_new-cloud_p4_new-gluinoMomentum).m()/MeV<<" MeV" <