This class is meant to filter for events that produce a dark brem matching originating in the target and matching the following parameters. More...

#include <TargetDarkBremFilter.h>

Public Member Functions
	TargetDarkBremFilter (const std::string &name, framework::config::Parameters &parameters)
	Class constructor.

	~TargetDarkBremFilter ()=default
	Class destructor.

std::vector< simcore::TYPE >	getTypes () override
	Get the types of actions this class can do.

void	BeginOfEventAction (const G4Event *e) override
	Reset flag signaling finding of A' to false.

void	stepping (const G4Step *step) override
	Looking for A' while primary is stepping.

void	EndOfEventAction (const G4Event *event) override
	Check flag signaling finding of A', if false, abort the event so it isn't saved.

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

virtual	~UserAction ()=default
	Destructor.

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 G4ClassificationOfNewTrack	ClassifyNewTrack (const G4Track *, const G4ClassificationOfNewTrack &cl)
	Method called when a track is updated.

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 Member Functions
bool	isOutsideTargetRegion (const G4VPhysicalVolume *vol) const
	Check if the volume is outside the target region.

bool	isOutsideTargetRegion (const G4LogicalVolume *vol) const
	Check if the volume is outside the target region.

void	AbortEvent (const std::string &reason) const
	Helper to abort an event with a message.

Private Attributes
double	threshold_
	Minimum energy [MeV] that the A' should have to keep the event.

bool	found_aprime_
	flag to signal that we saw an A'

Additional Inherited Members
Public Types inherited from simcore::UserAction
using	Factory
	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.

mutable::framework::logging::logger	theLog_
	the logging channel user actions can use `ldmx_log` with

Detailed Description

This class is meant to filter for events that produce a dark brem matching originating in the target and matching the following parameters.

 threshold: minimum energy [MeV] A' needs to have

See also: TargetBremFilter This filter is designed similar to the target brem filter where we check the secondaries of the primary electron if it is stopping within the target or if it is leaving the target region.

Definition at line 50 of file TargetDarkBremFilter.h.

Constructor & Destructor Documentation

◆ TargetDarkBremFilter()

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

Class constructor.

Retrieve the necessary configuration parameters

Definition at line 18 of file TargetDarkBremFilter.cxx.

20 : simcore::UserAction(name, parameters) {

21 threshold_ = parameters.getParameter<double>("threshold");

22}

biasing::TargetDarkBremFilter::threshold_

double threshold_

Minimum energy [MeV] that the A' should have to keep the event.

Definition TargetDarkBremFilter.h:153

simcore::UserAction

Interface that defines a user action.

Definition UserAction.h:43

References threshold_.

Member Function Documentation

◆ AbortEvent()

void TargetDarkBremFilter::AbortEvent ( const std::string & reason ) const

private

Helper to abort an event with a message.

Tells the RunManger to abort the current event after displaying the input message.

Parameters

[in] reason reason for aborting the event

Definition at line 128 of file TargetDarkBremFilter.cxx.

                                                                   {
  if (G4RunManager::GetRunManager()->GetVerboseLevel() > 1) {
    std::cout << "[ TargetDarkBremFilter ]: "
              << "("
              << G4EventManager::GetEventManager()
                     ->GetConstCurrentEvent()
                     ->GetEventID()
              << ") " << reason << " Aborting event." << std::endl;
  }
  G4RunManager::GetRunManager()->AbortEvent();
  return;
}

Referenced by EndOfEventAction(), and stepping().

◆ BeginOfEventAction()

void TargetDarkBremFilter::BeginOfEventAction ( const G4Event * e )

overridevirtual

Reset flag signaling finding of A' to false.

Parameters

[in] e event being started, unused

Reimplemented from simcore::UserAction.

Definition at line 24 of file TargetDarkBremFilter.cxx.

24 {

25 found_aprime_ = false;

26}

biasing::TargetDarkBremFilter::found_aprime_

bool found_aprime_

flag to signal that we saw an A'

Definition TargetDarkBremFilter.h:158

References found_aprime_.

◆ EndOfEventAction()

void biasing::TargetDarkBremFilter::EndOfEventAction ( const G4Event * event )

overridevirtual

Check flag signaling finding of A', if false, abort the event so it isn't saved.

Parameters

[in] event being ended, check if it is already aborted

Reimplemented from simcore::UserAction.

Definition at line 122 of file TargetDarkBremFilter.cxx.

                                                                {
  if (not event->IsAborted() and not found_aprime_) {
    AbortEvent("Did not find an A' after entire event was simulated.");
  }
}

References AbortEvent(), and found_aprime_.

◆ getTypes()

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

inlineoverridevirtual

Get the types of actions this class can do.

Returns: list of action types this class does

Implements simcore::UserAction.

Definition at line 70 of file TargetDarkBremFilter.h.

70 {

71 return {simcore::TYPE::STEPPING, simcore::TYPE::EVENT};

72 }

◆ isOutsideTargetRegion() [1/2]

bool biasing::TargetDarkBremFilter::isOutsideTargetRegion ( const G4LogicalVolume * vol ) const

inlineprivate

Check if the volume is outside the target region.

Note: will return true if vol is nullptr.

Parameters

[in] vol G4LogicalVolume to check region

Returns: true if vol is outside target region or nullptr or doesn't have a region

Definition at line 124 of file TargetDarkBremFilter.h.

                                                                      {
    if (!vol) return true;
    auto target_region = simcore::g4user::ptrretrieval::getRegion("target");
    if (!target_region) {
      ldmx_log(warn) << "Region 'target' not found in Geant4 region store";
    }
    auto region = vol->GetRegion();
    return region ? (region != target_region) : true;
  }

◆ isOutsideTargetRegion() [2/2]

bool biasing::TargetDarkBremFilter::isOutsideTargetRegion ( const G4VPhysicalVolume * vol ) const

inlineprivate

Check if the volume is outside the target region.

Note: will return true if vol is nullptr

Parameters

[in] vol G4VPhysicalVolume to check region

Returns: true if vol is outside target region or nullptr

Definition at line 111 of file TargetDarkBremFilter.h.

111 {

112 return vol ? isOutsideTargetRegion(vol->GetLogicalVolume()) : true;

113 }

biasing::TargetDarkBremFilter::isOutsideTargetRegion

bool isOutsideTargetRegion(const G4VPhysicalVolume *vol) const

Check if the volume is outside the target region.

Definition TargetDarkBremFilter.h:111

References isOutsideTargetRegion().

Referenced by isOutsideTargetRegion(), and stepping().

◆ stepping()

void TargetDarkBremFilter::stepping ( const G4Step * step )

overridevirtual

Looking for A' while primary is stepping.

We make sure that the current track is the primary electron that is within the target region. Then if the track is either stopped or leaving the target region, we look through its secondaries for a good A'.

Parameters

[in] step current G4Step

check on first step of A' to see if it originated in correct volume this needs to be here because sometimes the electron dark brems in the tagger and misses the target completely std::cout << "[ TargetDarkBremFilter ] : " << "(" << G4EventManager::GetEventManager() ->GetConstCurrentEvent() ->GetEventID() << ") " << "A' originated in " << track->GetLogicalVolumeAtVertex()->GetName() << std::endl;

Reimplemented from simcore::UserAction.

Definition at line 28 of file TargetDarkBremFilter.cxx.

                                                      {
  // don't process if event is already aborted
  if (G4EventManager::GetEventManager()->GetConstCurrentEvent()->IsAborted())
    return;
 
  // Get the track associated with this step.
  auto track{step->GetTrack()};
 
  // Leave if track is not an electron
  auto particle_def{track->GetParticleDefinition()};
  if (particle_def != G4Electron::Electron()) {
    if (particle_def == G4APrime::APrime() and
        track->GetCurrentStepNumber() == 1) {
      if (isOutsideTargetRegion(track->GetLogicalVolumeAtVertex()))
        AbortEvent("A' was not created within target.");
      // A' was found and originated in correct region
      found_aprime_ = true;
    }  // first step of A'
    return;
  }
 
  // Leave if track is not primary
  if (track->GetParentID() != 0) return;
 
  // Leave if track is not in target region
  if (isOutsideTargetRegion(track->GetVolume())) return;
 
  if (isOutsideTargetRegion(track->GetNextVolume())  // leaving target region
      or
      track->GetTrackStatus() == fStopAndKill  // stopping within target region
      or track->GetKineticEnergy() == 0.       // stopping within target region
  ) {
    // Get the electron secondries
    const std::vector<G4Track*>* secondaries{step->GetSecondary()};
    if (not secondaries or secondaries->size() == 0) {
      AbortEvent("Primary electron did not create secondaries.");
      return;
    } else {
      // check secondaries to see if we made a dark brem
      for (auto& secondary_track : *secondaries) {
        if (secondary_track->GetParticleDefinition() == G4APrime::APrime()) {
          // we found an A', woo-hoo!
 
          if (secondary_track->GetTotalEnergy() < threshold_) {
            AbortEvent(
                "A' was not created with total energy above input threshold.");
            return;
          }
 
          if (isOutsideTargetRegion(secondary_track->GetVolume())) {
            AbortEvent("A' was not created within target.");
            return;
          }
 
          /* debug message
          std::cout << "[ TargetDarkBremFilter ] : "
              << "(" <<
          G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID()
          << ") "
              << "Found an A' weighted " << secondary_track->GetWeight()
              << " originating in " <<
          secondary_track->GetVolume()->GetLogicalVolume()->GetName()
              << " with energy " << secondary_track->GetTotalEnergy() << " MeV."
              << std::endl;
           */
 
          // we found a good A', so we can leave
          return;
        }  // this secondary is an A'
      }    // loop through secondaries
    }      // are there secondaries to loop through
 
    // got here without finding A'
    AbortEvent("Primary electron did not create A'.");
    return;
  }  // should check secondaries of primary
 
  // get here if we are the primary in the target region,
  //  but shouldn't check the secondaries yet
  return;
}

References AbortEvent(), found_aprime_, isOutsideTargetRegion(), and threshold_.

Member Data Documentation

◆ found_aprime_

bool biasing::TargetDarkBremFilter::found_aprime_

private

flag to signal that we saw an A'

Definition at line 158 of file TargetDarkBremFilter.h.

Referenced by BeginOfEventAction(), EndOfEventAction(), and stepping().

◆ threshold_

double biasing::TargetDarkBremFilter::threshold_

private

Minimum energy [MeV] that the A' should have to keep the event.

Also used by PartialEnergySorter to determine which tracks should be processed first.

Parameter Name: 'threshold'

Definition at line 153 of file TargetDarkBremFilter.h.

Referenced by stepping(), and TargetDarkBremFilter().

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

Biasing/include/Biasing/TargetDarkBremFilter.h
Biasing/src/Biasing/TargetDarkBremFilter.cxx

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TargetDarkBremFilter()

Member Function Documentation

◆ AbortEvent()

◆ BeginOfEventAction()

◆ EndOfEventAction()

◆ getTypes()

◆ isOutsideTargetRegion() [1/2]

◆ isOutsideTargetRegion() [2/2]

◆ stepping()

Member Data Documentation

◆ found_aprime_

◆ threshold_