QuirkWatcher Class Reference

#include <QuirkWatcher.h>

Inheritance diagram for QuirkWatcher:

Collaboration diagram for QuirkWatcher:

Public Member Functions
	QuirkWatcher ()
virtual	~QuirkWatcher ()
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &)
virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)
virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)

Private Attributes
G4VParticleChange	m_particleChange

Detailed Description

Definition at line 11 of file QuirkWatcher.h.

Constructor & Destructor Documentation

QuirkWatcher()

QuirkWatcher::QuirkWatcher

(

)

Definition at line 14 of file QuirkWatcher.cxx.

                           : G4VProcess(G4String("QuirkWatcher")) {
    enableAtRestDoIt = false;
    enableAlongStepDoIt = false;
}

~QuirkWatcher()

QuirkWatcher::~QuirkWatcher ( )

virtual

Definition at line 19 of file QuirkWatcher.cxx.

{}

Member Function Documentation

AlongStepDoIt()

virtual G4VParticleChange * QuirkWatcher::AlongStepDoIt ( const G4Track & , const G4Step &  )

inlinevirtual

Definition at line 37 of file QuirkWatcher.h.

      {return 0;}

AlongStepGetPhysicalInteractionLength()

virtual G4double QuirkWatcher::AlongStepGetPhysicalInteractionLength ( const G4Track & , G4double , G4double , G4double & , G4GPILSelection *  )

inlinevirtual

Definition at line 26 of file QuirkWatcher.h.

      {return -1.0;}

AtRestDoIt()

virtual G4VParticleChange * QuirkWatcher::AtRestDoIt ( const G4Track & , const G4Step &  )

inlinevirtual

Definition at line 41 of file QuirkWatcher.h.

      {return 0;}

AtRestGetPhysicalInteractionLength()

virtual G4double QuirkWatcher::AtRestGetPhysicalInteractionLength ( const G4Track & , G4ForceCondition *  )

inlinevirtual

Definition at line 33 of file QuirkWatcher.h.

      {return -1.0;}

PostStepDoIt()

G4VParticleChange * QuirkWatcher::PostStepDoIt ( const G4Track & track, const G4Step &  )

virtual

Definition at line 30 of file QuirkWatcher.cxx.

  {
    // Get infracolor string
    auto part_nc ATLAS_THREAD_SAFE = // track should really be non-const in Geant4 interface
      const_cast<G4ParticleDefinition*>(track.GetParticleDefinition());
    Quirk* quirkDef = dynamic_cast<Quirk*>(part_nc);
    if (quirkDef == 0) {
        G4Exception("QuirkWatcher::PostStepDoIt", "NonQuirk", FatalErrorInArgument, "QuirkWatcher run on non-quirk particle");
        std::abort();
    }
    InfracolorForce& string = quirkDef->GetStringIn();
 
    if (track.GetCurrentStepNumber() > 1 && !string.IsSourceInitialized()) {
        // if the partner quirk is still alive, clear its string state too so that it
        // doesn't try to push a string vector to the orphan quirk and crash
        // leave as just a warning so that the event can continue and the hits recorded so far are kept, like for looperKiller'd tracks
        string.Clear();
        string.GetReactionForce()->Clear();
        G4Exception(
            "QuirkWatcher::PostStepDoIt",
            "QuirkMissingPartner",
            JustWarning,
            "QuirkWatcher: missing partner for quirk; killing orphan track"
        );
        m_particleChange.Initialize(track);
        m_particleChange.ProposeTrackStatus(fStopAndKill);
        return &m_particleChange;
    }
 
    // Update track status
    G4TrackStatus stat = track.GetTrackStatus();
    if (stat == fStopButAlive) {
        stat = fAlive;
    }
    if (stat == fAlive || stat == fSuspend) {
        G4bool passControl = !string.HasNextStringVector();
        //G4bool passControl = !string.HasNextStringVector() || (string.GetSumStrings().t() < string.GetReactionForce()->GetSumStrings().t());
        if (passControl) {
            if (string.IsSourceAlive()) {
                stat = fSuspend;
            } else {
                stat = fStopAndKill;
            }
        }
    }
    if (stat == fStopAndKill || stat == fKillTrackAndSecondaries) {
        string.TrackKilled();
    }
 
    m_particleChange.Initialize(track);
    m_particleChange.ProposeTrackStatus(stat);
    return &m_particleChange;
}

PostStepGetPhysicalInteractionLength()

G4double QuirkWatcher::PostStepGetPhysicalInteractionLength ( const G4Track & track, G4double previousStepSize, G4ForceCondition * condition )

virtual

Definition at line 21 of file QuirkWatcher.cxx.

  {
    *condition = StronglyForced;
    return DBL_MAX;
}

Member Data Documentation

m_particleChange

G4VParticleChange QuirkWatcher::m_particleChange

private

Definition at line 47 of file QuirkWatcher.h.

---

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

• QuirkWatcher.h
• QuirkWatcher.cxx