added lines to skip FastSim decays for SM particles that cross the be…

…am pipe.
cms-sw · Nov 14, 2021 · 1066de4 · 1066de4
1 parent 092bf03
commit 1066de4
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 2 deletions.
diff --git a/FastSimulation/SimplifiedGeometryPropagator/interface/ParticleManager.h b/FastSimulation/SimplifiedGeometryPropagator/interface/ParticleManager.h
@@ -146,7 +146,9 @@ inline bool isExotic(int pdgid_) {
   return ((pdgid >= 1000000 && pdgid < 4000000 && pdgid != 3000022) ||  // SUSY, R-hadron, and technicolor particles
           pdgid == 17 ||                                                // 4th generation lepton
           pdgid == 34 ||                                                // W-prime
-          pdgid == 37);                                                 // charged Higgs
+          pdgid == 37 ||                                                // charged Higgs  
+          pdgid == 39);                                                 // bulk graviton
+
 }
 
 #endif
diff --git a/FastSimulation/SimplifiedGeometryPropagator/src/ParticleManager.cc b/FastSimulation/SimplifiedGeometryPropagator/src/ParticleManager.cc
@@ -224,11 +224,17 @@ std::unique_ptr<fastsim::Particle> fastsim::ParticleManager::nextGenParticle() {
       }
     }
 
-    // particle must not decay before it reaches the beam pipe
+    // FastSim will not make hits out of particles that decay before reaching the beam pipe
     if (endVertex && endVertex->position().perp2() * lengthUnitConversionFactor2_ < beamPipeRadius2_) {
       continue;
     }
 
+    // SM particles that descend from exotics and cross the beam pipe radius should make hits but not be decayed 
+    if (productionVertex->position().perp2() * lengthUnitConversionFactor2_ < beamPipeRadius2_ &&
+	endVertex && endVertex->position().perp2() * lengthUnitConversionFactor2_ > beamPipeRadius2_) {
+      exoticRelativesChecker(productionVertex, exoticRelativeId, 0);
+    }    
+
     // make the particle
     std::unique_ptr<Particle> newParticle(
         new Particle(particle.pdg_id(),