TargetENProcessFilter.cxx

Go to the documentation of this file.

 
/*~~~~~~~~~~~~~*/
/*   Biasing   */
/*~~~~~~~~~~~~~*/
#include "Biasing/TargetENProcessFilter.h"
 
/*~~~~~~~~~~~~*/
/*   Geant4   */
/*~~~~~~~~~~~~*/
#include "G4RunManager.hh"
 
/*~~~~~~~~~~~~~*/
/*   SimCore   */
/*~~~~~~~~~~~~~*/
#include "SimCore/G4User/PtrRetrieval.h"
 
namespace biasing {
 
TargetENProcessFilter::TargetENProcessFilter(
    const std::string& name, framework::config::Parameters& parameters)
    : simcore::UserAction(name, parameters) {
  recoilEnergyThreshold_ = parameters.getParameter<double>("recoilThreshold");
}
TargetENProcessFilter::TargetENProcessFilter( {…}
 
TargetENProcessFilter::~TargetENProcessFilter() {}
 
void TargetENProcessFilter::stepping(const G4Step* step) {
  if (reactionOccurred_) return;
 
  // Get the track associated with this step.
  G4Track* track = step->GetTrack();
 
  // Only process tracks whose parent is the primary particle
  if (track->GetParentID() != 0) return;
 
  // get the PDGID of the track.
  G4int pdgID = track->GetParticleDefinition()->GetPDGEncoding();
 
  // Make sure that the particle being processed is an electron.
  if (pdgID != 11) return;  // Throw an exception
 
  // Get the volume the particle is in.
  G4VPhysicalVolume* track_volume = track->GetVolume();
  auto target_volume =
      simcore::g4user::ptrretrieval::getPhysicalVolume("target_PV");
  if (!target_volume) {
    ldmx_log(warn) << "Volume 'target_PV' not found in Geant4 volume store";
  }
  // If the particle isn't in the target, don't continue with the processing.
  if (track_volume != target_volume) return;
 
  /*std::cout << "*******************************" << std::endl;
  std::cout << "*   Step " << track->GetCurrentStepNumber() << std::endl;
  std::cout << "********************************" << std::endl;*/
 
  if (track->GetMomentum().mag() > recoilEnergyThreshold_) {
    track->SetTrackStatus(fKillTrackAndSecondaries);
    G4RunManager::GetRunManager()->AbortEvent();
    return;
  }
 
  // Get the particles daughters.
  const G4TrackVector* secondaries = step->GetSecondary();
 
  // If the brem photon doesn't undergo any reaction in the target, stop
  // processing the rest of the event.
  if (secondaries->size() == 0) {
    /*std::cout << "[ TargetENProcessFilter ]: "
                << "Electron did not interact in the target. --> Postponing
       tracks."
                << std::endl;*/
 
    track->SetTrackStatus(fKillTrackAndSecondaries);
    G4RunManager::GetRunManager()->AbortEvent();
    return;
  } else {
    G4String processName =
        secondaries->at(0)->GetCreatorProcess()->GetProcessName();
 
    /*std::cout << "[ TargetENProcessFilter ]: "
              << "Electron produced " << secondaries->size()
              << " particle via " << processName << " process."
              << std::endl;*/
 
    // Only record the process that is being biased
    if (!processName.contains(process_)) {
      /*std::cout << "[ TargetENProcessFilter ]: "
                << "Process was not " << BiasingMessenger::getProcess() << "-->
         Killing all tracks!"
                << std::endl;*/
 
      track->SetTrackStatus(fKillTrackAndSecondaries);
      G4RunManager::GetRunManager()->AbortEvent();
      return;
    }
 
    std::cout << "[ TargetENProcessFilter ]: "
              << "Electronuclear reaction resulted in " << secondaries->size()
              << " particles via " << processName << " process." << std::endl;
    // BiasingMessenger::setEventWeight(track->GetWeight());
    reactionOccurred_ = true;
  }
}
void TargetENProcessFilter::stepping(const G4Step* step) {…}
 
void TargetENProcessFilter::EndOfEventAction(const G4Event*) {
  reactionOccurred_ = false;
}
void TargetENProcessFilter::EndOfEventAction(const G4Event*) {…}
}  // namespace biasing
 
DECLARE_ACTION(biasing, TargetENProcessFilter)