User action plugin that filters events that don't see a hard brem from the target undergo a photo-nuclear reaction in the ECal. More...

#include <EcalProcessFilter.h>

Public Member Functions
	EcalProcessFilter (const std::string &name, framework::config::Parameters &parameters)

	~EcalProcessFilter ()
	Destructor.

void	stepping (const G4Step *step) override
	Method called after each simulation step.

G4ClassificationOfNewTrack	ClassifyNewTrack (const G4Track *aTrack, const G4ClassificationOfNewTrack &currentTrackClass) override
	Classify a new track which postpones track processing.

std::vector< simcore::TYPE >	getTypes () override
	Retrieve the type of actions this class defines.

Public Member Functions inherited from simcore::UserAction
	UserAction (const std::string &name, framework::config::Parameters &parameters)
	Constructor.

virtual	~UserAction ()=default
	Destructor.

virtual void	BeginOfEventAction (const G4Event *)
	Method called at the beginning of every event.

virtual void	EndOfEventAction (const G4Event *)
	Method called at the end of every event.

virtual void	BeginOfRunAction (const G4Run *)
	Method called at the beginning of a run.

virtual void	EndOfRunAction (const G4Run *)
	Method called at the end of a run.

virtual void	PreUserTrackingAction (const G4Track *)
	Method called before the UserTrackingAction.

virtual void	PostUserTrackingAction (const G4Track *)
	Method called after the UserTrackingAction.

virtual void	NewStage ()
	Method called at the beginning of a new stage.

virtual void	PrepareNewEvent ()
	Method called at the beginning of a new event.

Private Attributes
G4Track *	currentTrack_ {nullptr}
	Pointer to the current track being processed.

std::string	process_ {""}
	Process to filter.

Additional Inherited Members
Public Types inherited from simcore::UserAction
using	Factory = ::simcore::Factory< UserAction, std::shared_ptr< UserAction >, const std::string &, framework::config::Parameters & >
	factory for user actions

Protected Member Functions inherited from simcore::UserAction
UserEventInformation *	getEventInfo () const
	Get a handle to the event information.

const std::map< int, ldmx::SimParticle > &	getCurrentParticleMap () const
	Get the current particle map.

Protected Attributes inherited from simcore::UserAction
std::string	name_ {""}
	Name of the UserAction.

framework::config::Parameters	parameters_
	The set of parameters used to configure this class.

Detailed Description

User action plugin that filters events that don't see a hard brem from the target undergo a photo-nuclear reaction in the ECal.

Definition at line 30 of file EcalProcessFilter.h.

Constructor & Destructor Documentation

◆ EcalProcessFilter()

biasing::EcalProcessFilter::EcalProcessFilter	(	const std::string &	name,
		framework::config::Parameters &	parameters
	)

Definition at line 19 of file EcalProcessFilter.cxx.

    : simcore::UserAction(name, parameters) {
  process_ = parameters.getParameter<std::string>("process");
}

◆ ~EcalProcessFilter()

biasing::EcalProcessFilter::~EcalProcessFilter ( )

Destructor.

Definition at line 25 of file EcalProcessFilter.cxx.

25{}

Member Function Documentation

◆ ClassifyNewTrack()

G4ClassificationOfNewTrack biasing::EcalProcessFilter::ClassifyNewTrack	(	const G4Track *	aTrack,
		const G4ClassificationOfNewTrack &	currentTrackClass
	)

overridevirtual

Classify a new track which postpones track processing.

Track processing resumes normally if a target PN interaction occurred.

Parameters

aTrack	The Geant4 track.
currentTrackClass	The current track classification.

Reimplemented from simcore::UserAction.

Definition at line 27 of file EcalProcessFilter.cxx.

                                                                               {
  // Get the particle type.
  G4String particleName = track->GetParticleDefinition()->GetParticleName();
 
  if (track == currentTrack_) {
    /*
    std::cout << "[ EcalProcessFilter ]: "
        << "Putting track " << track->GetTrackID()
        << " onto waiting stack." << std::endl;
    */
    currentTrack_ = nullptr;
    return fWaiting;
  }
 
  // Use current classification by default so values from other plugins are not
  // overridden.
  G4ClassificationOfNewTrack classification = currentTrackClass;
 
  return classification;
}

◆ getTypes()

std::vector< simcore::TYPE > biasing::EcalProcessFilter::getTypes ( )

inlineoverridevirtual

Retrieve the type of actions this class defines.

Implements simcore::UserAction.

Definition at line 56 of file EcalProcessFilter.h.

                                             {
    return {simcore::TYPE::STACKING, simcore::TYPE::STEPPING};
  }

◆ stepping()

void biasing::EcalProcessFilter::stepping ( const G4Step * )

overridevirtual

Method called after each simulation step.

TYPE::STEPPING

Parameters

current Geant4 step

Check if the photon will be exiting the ecal

The 'hcal_PV' volume name is automatically constructed by Geant4's GDML parser and was found by inspecting the geometry using a visualization. This Physical Volume (PV) is associated with the hcal parent volume and so it will break if the hcal parent volume changes its name.

Reimplemented from simcore::UserAction.

Definition at line 49 of file EcalProcessFilter.cxx.

                                                   {
  // Get the track associated with this step.
  auto track{step->GetTrack()};
 
  if (G4EventManager::GetEventManager()->GetConstCurrentEvent()->IsAborted())
    return;
 
  // Get the track info and check if this track is a brem candidate
  auto trackInfo{simcore::UserTrackInformation::get(track)};
  if ((trackInfo != nullptr) && !trackInfo->isBremCandidate()) return;
 
  // Get the particles daughters.
  auto secondaries{step->GetSecondary()};
 
  // Get the region the particle is currently in.  Continue processing
  // the particle only if it's in the calorimeter region.
  if (auto region{
          track->GetVolume()->GetLogicalVolume()->GetRegion()->GetName()};
      region.compareTo("CalorimeterRegion") != 0) {
    // If secondaries were produced outside of the volume of interest,
    // and there aren't additional brems to process, abort the
    // event.  Otherwise, suspend the track and move on to the next
    // brem.
    if (secondaries->size() != 0) {
      /*
      std::cout << "[ EcalProcessFilter ]: "
            <<
      G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID()
            << " secondaries outside ecal...";
      */
      if (getEventInfo()->bremCandidateCount() == 1) {
        // std::cout << "aborting the event." << std::endl;
        track->SetTrackStatus(fKillTrackAndSecondaries);
        G4RunManager::GetRunManager()->AbortEvent();
        currentTrack_ = nullptr;
      } else {
        /*
        std::cout << "suspending the track " << track->GetTrackID()
            << " , " << getEventInfo()->bremCandidateCount() << " brems left."
            << std::endl;
        */
        currentTrack_ = track;
        track->SetTrackStatus(fSuspend);
        getEventInfo()->decBremCandidateCount();
        trackInfo->tagBremCandidate(false);
      }
    }
    return;
  }
 
  // If the particle doesn't interact, then move on to the next step.
  if (secondaries->size() == 0) {
    if (auto volume{track->GetNextVolume()->GetName()};
        volume.compareTo("hcal_PV") == 0) {
      /*
      std::cout << "[ EcalProcessFilter ]: "
            <<
      G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID()
            << " no secondaries when leaving ecal...";
      */
      if (getEventInfo()->bremCandidateCount() == 1) {
        // std::cout << "aborting the event." << std::endl;
        track->SetTrackStatus(fKillTrackAndSecondaries);
        G4RunManager::GetRunManager()->AbortEvent();
        currentTrack_ = nullptr;
      } else {
        /*
        std::cout << "suspending the track " << track->GetTrackID()
            << " , " << getEventInfo()->bremCandidateCount() << " brems left."
            << std::endl;
        */
        currentTrack_ = track;
        track->SetTrackStatus(fSuspend);
        getEventInfo()->decBremCandidateCount();
        trackInfo->tagBremCandidate(false);
      }
    }
 
    return;
  } else {
    // If the brem gamma interacts and produces secondaries, get the
    // process used to create them.
    auto processName{secondaries->at(0)->GetCreatorProcess()->GetProcessName()};
 
    // Only record the process that is being biased
    if (!processName.contains(process_)) {
      /*
      std::cout << "[ EcalProcessFilter ]: "
            <<
      G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID()
            << " not PN products...";
      */
      if (getEventInfo()->bremCandidateCount() == 1) {
        // std::cout << "aborting the event." << std::endl;
        track->SetTrackStatus(fKillTrackAndSecondaries);
        G4RunManager::GetRunManager()->AbortEvent();
        currentTrack_ = nullptr;
      } else {
        /*
        std::cout << "suspending the track " << track->GetTrackID()
            << " , " << getEventInfo()->bremCandidateCount() << " brems left."
            << std::endl;
        */
        currentTrack_ = track;
        track->SetTrackStatus(fSuspend);
        getEventInfo()->decBremCandidateCount();
        trackInfo->tagBremCandidate(false);
      }
      return;
    }
 
    ldmx_log(debug) << "[ EcalProcessFilter ]: "
                    << G4EventManager::GetEventManager()
                           ->GetConstCurrentEvent()
                           ->GetEventID()
                    << " Brem photon produced " << secondaries->size()
                    << " particle via " << processName << " process.";
    trackInfo->tagBremCandidate(false);
    trackInfo->setSaveFlag(true);
    trackInfo->tagPNGamma();
    getEventInfo()->decBremCandidateCount();
  }
}

References simcore::UserTrackInformation::get(), and simcore::UserTrackInformation::tagBremCandidate().

Member Data Documentation

◆ currentTrack_

G4Track* biasing::EcalProcessFilter::currentTrack_ {nullptr}

private

Pointer to the current track being processed.

Definition at line 62 of file EcalProcessFilter.h.

62{nullptr};

◆ process_

std::string biasing::EcalProcessFilter::process_ {""}

private

Process to filter.

Definition at line 65 of file EcalProcessFilter.h.

65{""};

The documentation for this class was generated from the following files:

Biasing/include/Biasing/EcalProcessFilter.h
Biasing/src/Biasing/EcalProcessFilter.cxx

Public Member Functions

Private Attributes

Additional Inherited Members