HyperbolaStepper Class Reference

#include <HyperbolaStepper.h>

Inheritance diagram for HyperbolaStepper:

Collaboration diagram for HyperbolaStepper:

Public Member Functions
	HyperbolaStepper (InfracolorForce &string, const G4Track &track, const G4Field *field=0)
virtual	~HyperbolaStepper ()
virtual void	Stepper (const G4double y[], const G4double[], G4double h, G4double yout[], G4double yerr[])
virtual G4double	DistChord () const
virtual void	ComputeRightHandSide (const G4double[], G4double[])
virtual G4int	IntegratorOrder () const
G4double	GetForce () const
const G4ElectroMagneticField *	GetField () const
G4double	GetMass () const
G4double	GetCharge () const
const G4LorentzVector &	GetStartMomentum () const
G4double	GetMaxExpRapidity () const
G4double	GetMaxLength () const
void	Update (G4FieldTrack &fieldTrack, G4bool forceStep)
void	SetDebug (G4bool debug)

Private Member Functions
void	SetCurrStep (G4double length)

Private Attributes
InfracolorForce &	m_string
const G4ElectroMagneticField *const	m_field
const G4double	m_mass
const G4double	m_charge
const G4LorentzVector	m_startMomentum
const G4double	m_maxExpRapidity
std::vector< HyperbolaStep >	m_steps
HyperbolaStep	m_currStep
std::vector< HyperbolaStep >::size_type	m_nPrevSteps
G4bool	m_debug

Static Private Attributes
static const G4int	m_NUM_VARS = 9

Detailed Description

Definition at line 19 of file HyperbolaStepper.h.

Constructor & Destructor Documentation

HyperbolaStepper()

HyperbolaStepper::HyperbolaStepper	(	InfracolorForce &	string,
		const G4Track &	track,
		const G4Field *	field = 0 )

Definition at line 12 of file HyperbolaStepper.cxx.

  :
    // 9 "integration" variables:
    // x, y, z, px, py, pz, KE, t_lab, steplength (in place of t_proper)
    G4MagIntegratorStepper(new DummyEquation(), m_NUM_VARS),
    m_string(string),
    m_field(dynamic_cast<const G4ElectroMagneticField*>(field)),
    m_mass(track.GetDynamicParticle()->GetMass()),
    m_charge(track.GetDynamicParticle()->GetCharge()),
    m_startMomentum(track.GetDynamicParticle()->Get4Momentum()),
    m_maxExpRapidity(string.GetMaxExpRapidity()),
    m_steps(),
    m_currStep(this, string, track),
    m_nPrevSteps(0),
    m_debug(false)
{
    if (field != 0 && m_field == 0) {
        G4Exception("HyperbolaStepper::HyperbolaStepper", "QuirkNonEM", FatalErrorInArgument, "Non-EM fields not supported");
    }
 
    // Precompute steps
    m_currStep.PrepareNextStep();
    m_steps.push_back(m_currStep);
    m_currStep.Step(m_currStep.GetMaxLength());
    while (m_currStep.GetStringPtr() != string.GetStringVectors().end()) {
        HyperbolaStep next = m_currStep;
        next.PrepareNextStep();
        if (next.IsBoostLimited()) break; // don't allow 2nd+ step if it would be cut by boost limit
        m_steps.push_back(next);
        m_currStep = next;
        m_currStep.Step(m_currStep.GetMaxLength());
    }
    m_nPrevSteps = m_steps.size() - 1;
}

~HyperbolaStepper()

HyperbolaStepper::~HyperbolaStepper ( )

virtual

Definition at line 50 of file HyperbolaStepper.cxx.

                                    {
    delete GetEquationOfMotion();
}

Member Function Documentation

ComputeRightHandSide()

virtual void HyperbolaStepper::ComputeRightHandSide ( const G4double [], G4double [] )

inlinevirtual

Definition at line 36 of file HyperbolaStepper.h.

{}

DistChord()

G4double HyperbolaStepper::DistChord ( ) const

virtual

Definition at line 85 of file HyperbolaStepper.cxx.

                                           {
    // TODO
    return 0;
}

GetCharge()

G4double HyperbolaStepper::GetCharge ( ) const

inline

Definition at line 73 of file HyperbolaStepper.h.

{return m_charge;}

GetField()

const G4ElectroMagneticField * HyperbolaStepper::GetField ( ) const

inline

Definition at line 71 of file HyperbolaStepper.h.

{return m_field;}

GetForce()

G4double HyperbolaStepper::GetForce ( ) const

inline

Definition at line 70 of file HyperbolaStepper.h.

{return m_string.GetStringForce();}

GetMass()

G4double HyperbolaStepper::GetMass ( ) const

inline

Definition at line 72 of file HyperbolaStepper.h.

{return m_mass;}

GetMaxExpRapidity()

G4double HyperbolaStepper::GetMaxExpRapidity ( ) const

inline

Definition at line 75 of file HyperbolaStepper.h.

{return m_maxExpRapidity;}

GetMaxLength()

G4double HyperbolaStepper::GetMaxLength ( ) const

inline

Definition at line 77 of file HyperbolaStepper.h.

{return m_steps.back().GetMaxLength();}

GetStartMomentum()

const G4LorentzVector & HyperbolaStepper::GetStartMomentum ( ) const

inline

Definition at line 74 of file HyperbolaStepper.h.

{return m_startMomentum;}

IntegratorOrder()

virtual G4int HyperbolaStepper::IntegratorOrder ( ) const

inlinevirtual

Definition at line 37 of file HyperbolaStepper.h.

{return 1;}

SetCurrStep()

void HyperbolaStepper::SetCurrStep ( G4double length )

private

Definition at line 54 of file HyperbolaStepper.cxx.

                                                  {
    if (m_currStep.GetLength() == length) return;
    std::vector<HyperbolaStep>::size_type i = 0;
    while (i+1 < m_steps.size() && m_steps[i+1].GetLength() < length) i++;
    m_currStep = m_steps[i];
    m_currStep.Step(length);
    m_nPrevSteps = i;
}

SetDebug()

void HyperbolaStepper::SetDebug ( G4bool debug )

inline

Definition at line 79 of file HyperbolaStepper.h.

{m_debug = debug;}

Stepper()

void HyperbolaStepper::Stepper ( const G4double y[], const G4double [], G4double h, G4double yout[], G4double yerr[] )

virtual

Definition at line 63 of file HyperbolaStepper.cxx.

  {
    G4int maxvar = GetNumberOfStateVariables();
    for (G4int i = 0; i < maxvar; i++) yout[i] = y[i];
 
    if (m_debug) G4cout << "HyperbolaStepper: asked to move " << h << G4endl;
    if (m_debug) G4cout << "HyperbolaStepper: start = " << y[0] << ", " << y[1] << ", " << y[2] << " [" << y[8] << "]" << G4endl;
 
    SetCurrStep(y[8] + h);
    m_currStep.Dump(yout);
 
    if (m_debug) G4cout << "HyperbolaStepper: end   = " << yout[0] << ", " << yout[1] << ", " << yout[2] << " [" << yout[8] << "]" << G4endl;
 
    // TODO: Errors
    for (G4int i = 0; i < m_NUM_VARS; i++) yerr[i] = 0;
}

Update()

void HyperbolaStepper::Update	(	G4FieldTrack &	fieldTrack,
		G4bool	forceStep )

Definition at line 90 of file HyperbolaStepper.cxx.

                                                                        {
    if (forceStep) {
        // Propagator refused to take step through detector geometry because length too short;
        // take step ourselves and update string state and quirk 4-momentum only (don't change position)
        m_currStep = m_steps.back();
        m_currStep.Step(m_currStep.GetMaxLength());
        m_nPrevSteps = m_steps.size() - 1;
        G4LorentzVector p = m_currStep.GetMomentum();
        fieldTrack.UpdateFourMomentum(p.t() - m_mass, p.vect().unit());
    } else {
        SetCurrStep(fieldTrack.GetProperTimeOfFlight());
    }
 
    // Update string vectors
    m_string.PopTo(m_currStep.GetStringPtr(), m_currStep.GetFracLeft());
    for (std::vector<HyperbolaStep>::size_type i = 0; i < m_nPrevSteps; i++) {
        if (m_debug) G4cout << "HyperbolaStepper: pushing vector " << m_steps[i].GetStringOut().lv() << G4endl;
        m_string.GetReactionForce()->PushStringVector(m_steps[i].GetStringOut());
    }
    if (m_debug) G4cout << "HyperbolaStepper: pushing vector " << m_currStep.GetStringOut().lv() << G4endl;
    m_string.GetReactionForce()->PushStringVector(m_currStep.GetStringOut());
}

Member Data Documentation

m_charge

const G4double HyperbolaStepper::m_charge

private

Definition at line 57 of file HyperbolaStepper.h.

m_currStep

HyperbolaStep HyperbolaStepper::m_currStep

private

Definition at line 62 of file HyperbolaStepper.h.

m_debug

G4bool HyperbolaStepper::m_debug

private

Definition at line 65 of file HyperbolaStepper.h.

m_field

const G4ElectroMagneticField* const HyperbolaStepper::m_field

private

Definition at line 55 of file HyperbolaStepper.h.

m_mass

const G4double HyperbolaStepper::m_mass

private

Definition at line 56 of file HyperbolaStepper.h.

m_maxExpRapidity

const G4double HyperbolaStepper::m_maxExpRapidity

private

Definition at line 59 of file HyperbolaStepper.h.

m_nPrevSteps

std::vector<HyperbolaStep>::size_type HyperbolaStepper::m_nPrevSteps

private

Definition at line 63 of file HyperbolaStepper.h.

m_NUM_VARS

const G4int HyperbolaStepper::m_NUM_VARS = 9

staticprivate

Definition at line 52 of file HyperbolaStepper.h.

m_startMomentum

const G4LorentzVector HyperbolaStepper::m_startMomentum

private

Definition at line 58 of file HyperbolaStepper.h.

m_steps

std::vector<HyperbolaStep> HyperbolaStepper::m_steps

private

Definition at line 61 of file HyperbolaStepper.h.

m_string

InfracolorForce& HyperbolaStepper::m_string

private

Definition at line 54 of file HyperbolaStepper.h.

---

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

• HyperbolaStepper.h
• HyperbolaStepper.cxx