larg4::G4BadIdeaAction Class Reference

#include "G4BadIdeaAction.h"

Inheritance diagram for larg4::G4BadIdeaAction:

Public Member Functions
	G4BadIdeaAction (int)
virtual	~G4BadIdeaAction ()
virtual void	SteppingAction (const G4Step *)
	G4UserSteppingAction interface. More...
virtual void	Config (fhicl::ParameterSet const &)
	Override Config() to extract any necessary parameters. More...
virtual void	PrintConfig (std::string const &)
	Override PrintConfig() to print out current configuration. More...
virtual void	BeginOfRunAction (const G4Run *)
	G4UserRunAction interfaces. More...
virtual void	EndOfRunAction (const G4Run *)
virtual void	BeginOfEventAction (const G4Event *)
	G4UserEventAction interfaces. More...
virtual void	EndOfEventAction (const G4Event *)
virtual void	PreTrackingAction (const G4Track *)
	G4UserTrackingAction interfaces. More...
virtual void	PostTrackingAction (const G4Track *)
virtual bool	ProvidesStacking ()
virtual G4ClassificationOfNewTrack	StackClassifyNewTrack (const G4Track *)
	G4UserStackingAction interfaces. More...
virtual void	StackNewStage ()
virtual void	StackPrepareNewEvent ()
std::string const &	GetName () const
void	SetName (std::string const &name)

Private Attributes
art::ServiceHandle< geo::Geometry const >	fGeo
int	fNoIncomingMuons

Detailed Description

Definition at line 29 of file G4BadIdeaAction.h.

Constructor & Destructor Documentation

larg4::G4BadIdeaAction::G4BadIdeaAction

(

int

trkOption

)

Definition at line 21 of file G4BadIdeaAction.cxx.

                                                : fNoIncomingMuons(trkOption)
  {
    // trkOption comes from LArG4's fSmartStacking
    // Negative values effect action in this routine. Positive values
    // effect action in LArStackingAction.
    mf::LogWarning("G4BadIdeaAction") << "instantiating the G4BadIdeaAction \n"
                                      << "This UserAction is only to be used with "
                                      << "Geant4 v4.9.4.p02 to solve a stepping bug.\n"
                                      << "If you are using a different version of G4, "
                                      << "remove this UserAction from your list in LArG4.cxx";
  }

larg4::G4BadIdeaAction::~G4BadIdeaAction ( )

virtual

Definition at line 35 of file G4BadIdeaAction.cxx.

{}

Member Function Documentation

virtual void g4b::UserAction::BeginOfEventAction ( const G4Event *  )

inlinevirtualinherited

G4UserEventAction interfaces.

Reimplemented in larg4::ParticleListAction, larg4::VisualizationAction, and altns::ExampleAction.

Definition at line 78 of file UserAction.h.

{};

virtual void g4b::UserAction::BeginOfRunAction ( const G4Run *  )

inlinevirtualinherited

G4UserRunAction interfaces.

The following a list of methods that correspond to the available user action classes in Geant 4.0.1 and higher.

Reimplemented in larg4::VisualizationAction.

Definition at line 74 of file UserAction.h.

{};

virtual void g4b::UserAction::Config ( fhicl::ParameterSet const &  )

inlinevirtualinherited

Override Config() to extract any necessary parameters.

Reimplemented in altns::ExampleAction.

Definition at line 65 of file UserAction.h.

Referenced by g4b::UserAction::UserAction().

{};

virtual void g4b::UserAction::EndOfEventAction ( const G4Event *  )

inlinevirtualinherited

Reimplemented in larg4::ParticleListAction, larg4::VisualizationAction, and altns::ExampleAction.

Definition at line 79 of file UserAction.h.

{};

virtual void g4b::UserAction::EndOfRunAction ( const G4Run *  )

inlinevirtualinherited

Reimplemented in larg4::VisualizationAction.

Definition at line 75 of file UserAction.h.

{};

std::string const& g4b::UserAction::GetName ( ) const

inlineinherited

Definition at line 98 of file UserAction.h.

References g4b::UserAction::myName.

Referenced by g4b::UserActionManager::PrintActionList().

{ return myName; }

virtual void g4b::UserAction::PostTrackingAction ( const G4Track *  )

inlinevirtualinherited

Reimplemented in larg4::ParticleListAction, and altns::ExampleAction.

Definition at line 83 of file UserAction.h.

{};

virtual void g4b::UserAction::PreTrackingAction ( const G4Track *  )

inlinevirtualinherited

G4UserTrackingAction interfaces.

Reimplemented in larg4::ParticleListAction, and altns::ExampleAction.

Definition at line 82 of file UserAction.h.

{};

virtual void g4b::UserAction::PrintConfig ( std::string const &  )

inlinevirtualinherited

Override PrintConfig() to print out current configuration.

Reimplemented in altns::ExampleAction.

Definition at line 68 of file UserAction.h.

Referenced by g4b::UserActionManager::PrintActionList().

{};

virtual bool g4b::UserAction::ProvidesStacking ( )

inlinevirtualinherited

Does this UserAction do stacking? Override to return "true" if the following interfaces are implemented

Reimplemented in altns::ExampleAction.

Definition at line 90 of file UserAction.h.

Referenced by g4b::UserActionManager::PrintActionList().

{ return false; } 

void g4b::UserAction::SetName ( std::string const &  name )

inlineinherited

Definition at line 99 of file UserAction.h.

References g4b::UserAction::myName.

Referenced by g4b::UserActionFactory::GetUserAction().

{ myName = name; }

virtual G4ClassificationOfNewTrack g4b::UserAction::StackClassifyNewTrack ( const G4Track *  )

inlinevirtualinherited

G4UserStackingAction interfaces.

Reimplemented in altns::ExampleAction.

Definition at line 93 of file UserAction.h.

{ return fUrgent; }

virtual void g4b::UserAction::StackNewStage ( )

inlinevirtualinherited

Reimplemented in altns::ExampleAction.

Definition at line 94 of file UserAction.h.

{};

virtual void g4b::UserAction::StackPrepareNewEvent ( )

inlinevirtualinherited

Reimplemented in altns::ExampleAction.

Definition at line 95 of file UserAction.h.

{};

void larg4::G4BadIdeaAction::SteppingAction ( const G4Step *  )

virtual

G4UserSteppingAction interface.

Todo:

Do not copy the code below. Contact Brian Rebel, Eric Church,

Todo:

Bill Seligman and Andrzej Szelc for reasons why not to!

Reimplemented from g4b::UserAction.

Definition at line 39 of file G4BadIdeaAction.cxx.

References util::abs(), fNoIncomingMuons, and art::ExceptionDetail::translate().

  {

    // If the step size is such that the particle appears to be "stuck"
    // in its trajectory, give it a kick.
    const double epsilon = 5000. * std::numeric_limits<double>::epsilon();
    const double stepSize = step->GetStepLength();
    G4Track* nonConstTrack = const_cast<G4Track*>(step->GetTrack());

    if (step->GetTrack()->GetCurrentStepNumber() > 5e4 && stepSize < epsilon &&
        (step->GetTrack()->GetCurrentStepNumber() % 1000) == 1) {

      // Cast away the const-ness of the pointer to G4Step.
      // This is dangerous. Don't do this at home. We're
      // only doing this because we're desperate.
      // The need to do this is the result of a bug in Geant4 v4.9.4.p02
      // We should no longer call this code when we move beyond that version
      const double kick = 0.001;

      G4ThreeVector aValue = nonConstTrack->GetPosition();

      mf::LogWarning("G4BadIdeaAction")
        << "##### In endless loop. Kicking particle by "
        << " (+0.001,+0.001,+0.001) --- "
        << " PDG and encoding " << step->GetTrack()->GetDynamicParticle()->GetPDGcode() << " "
        << step->GetTrack()->GetDynamicParticle()->GetParticleDefinition()->GetParticleName()
        << " current step number: " << step->GetTrack()->GetCurrentStepNumber()
        << " stepsize: " << stepSize << " x,y,z  " << aValue.x() << " " << aValue.y() << " "
        << aValue.z();

      G4ThreeVector translate(kick, kick, kick);
      aValue += translate;

      nonConstTrack->SetPosition(aValue);
    }

    if (fNoIncomingMuons < 0) {
      // This is for overlays of, say, rock muons, which we have
      G4StepPoint* thePrePoint = step->GetPreStepPoint();
      G4VPhysicalVolume* thePrePV = thePrePoint->GetPhysicalVolume();
      G4String thePrePVname = thePrePV->GetName();
      G4StepPoint* thePostPoint = step->GetPostStepPoint();
      G4VPhysicalVolume* thePostPV = thePostPoint->GetPhysicalVolume();
      G4String thePostPVname("null");
      if (thePostPV) thePostPVname = thePostPV->GetName();

      if (abs(step->GetTrack()->GetDynamicParticle()->GetPDGcode()) == 13 &&
          thePostPVname.contains("volTPCActive") && !thePrePVname.contains("volTPCActive"))
        ((G4Track*)nonConstTrack)->SetTrackStatus(fStopAndKill);
    }

    return;
  }

Member Data Documentation

art::ServiceHandle<geo::Geometry const> larg4::G4BadIdeaAction::fGeo

private

Definition at line 40 of file G4BadIdeaAction.h.

int larg4::G4BadIdeaAction::fNoIncomingMuons

private

Definition at line 41 of file G4BadIdeaAction.h.

Referenced by SteppingAction().

---

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

• larsim/v09_43_00/source/larsim/LegacyLArG4/G4BadIdeaAction.h
• larsim/v09_43_00/source/larsim/LegacyLArG4/G4BadIdeaAction.cxx