#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 >	fGeo

int	fNoIncomingMuons

Detailed Description

Definition at line 32 of file G4BadIdeaAction.h.

Constructor & Destructor Documentation

larg4::G4BadIdeaAction::G4BadIdeaAction ( int trkOption )

Definition at line 23 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 38 of file G4BadIdeaAction.cxx.

39 {

40 }

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.

78 {};

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.

74 {};

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().

65 {};

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.

79 {};

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

inlinevirtualinherited

Reimplemented in larg4::VisualizationAction.

Definition at line 75 of file UserAction.h.

75 {};

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().

98 { return myName; }

g4b::UserAction::myName

std::string myName

self-knowledge

Definition: UserAction.h:101

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

inlinevirtualinherited

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

Definition at line 83 of file UserAction.h.

83 {};

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.

82 {};

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().

68 {};

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().

90 { 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().

99 { myName = name; }

g4b::UserAction::myName

std::string myName

self-knowledge

Definition: UserAction.h:101

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

inlinevirtualinherited

G4UserStackingAction interfaces.

Reimplemented in altns::ExampleAction.

Definition at line 93 of file UserAction.h.

93 { return fUrgent; }

virtual void g4b::UserAction::StackNewStage ( )

inlinevirtualinherited

Reimplemented in altns::ExampleAction.

Definition at line 94 of file UserAction.h.

94 {};

virtual void g4b::UserAction::StackPrepareNewEvent ( )

inlinevirtualinherited

Reimplemented in altns::ExampleAction.

Definition at line 95 of file UserAction.h.

95 {};

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 44 of file G4BadIdeaAction.cxx.

References 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> larg4::G4BadIdeaAction::fGeo

private

Definition at line 45 of file G4BadIdeaAction.h.

int larg4::G4BadIdeaAction::fNoIncomingMuons

private

Definition at line 46 of file G4BadIdeaAction.h.

Referenced by SteppingAction().

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

larsim/v06_53_00/source/larsim/LArG4/G4BadIdeaAction.h
larsim/v06_53_00/source/larsim/LArG4/G4BadIdeaAction.cxx

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation