RecoilMissesEcalSkimmer.cxx

Go to the documentation of this file.

 
#include "Recon/Skims/RecoilMissesEcalSkimmer.h"
 
namespace recon {
 
RecoilMissesEcalSkimmer::RecoilMissesEcalSkimmer(const std::string &name,
                                                 framework::Process &process)
    : framework::Producer(name, process) {}
 
RecoilMissesEcalSkimmer::~RecoilMissesEcalSkimmer() {}
 
void RecoilMissesEcalSkimmer::produce(framework::Event &event) {
  // Get the collection of simulated particles from the event
  auto particleMap{event.getMap<int, ldmx::SimParticle>("SimParticles")};
 
  // Search for the recoil electron
  auto [recoilTrackID, recoilElectron] = Analysis::getRecoil(particleMap);
 
  // Get the collection of simulated Ecal hits from the event.
  const std::vector<ldmx::SimCalorimeterHit> ecalSimHits =
      event.getCollection<ldmx::SimCalorimeterHit>("EcalSimHits");
 
  // Loop through the Ecal hits and check if the recoil electron is
  // associated with any of them.  If there are any recoil electron hits
  // in the Ecal, drop the event.
  bool hasRecoilElectronHits = false;
  for (const ldmx::SimCalorimeterHit &simHit : ecalSimHits) {
    /*std::cout << "[ RecoilMissesEcalSkimmer ]: "
              << "Number of hit contributions: "
              << simHit->getNumberOfContribs() << std::endl;*/
 
    for (int iContrib = 0; iContrib < simHit.getNumberOfContribs();
         ++iContrib) {
      ldmx::SimCalorimeterHit::Contrib contrib = simHit.getContrib(iContrib);
 
      if (contrib.trackID == recoilTrackID) {
        /*std::cout << "[ RecoilMissesEcalSkimmer ]: "
                  << "Ecal hit associated with recoil electron." << std::endl;
         */
 
        hasRecoilElectronHits = true;
      }
    }
  }
 
  // Tell the skimmer to keep or drop the event based on whether there
  // were recoil electron hits found in the Ecal.
  if (hasRecoilElectronHits) {
    setStorageHint(framework::hint_shouldDrop);
  } else {
    setStorageHint(framework::hint_shouldKeep);
  }
}
}  // namespace recon
 
DECLARE_PRODUCER_NS(recon, RecoilMissesEcalSkimmer);