#include <InfracolorForce.h>

Collaboration diagram for InfracolorForce:

Public Member Functions
	InfracolorForce ()
	~InfracolorForce ()
void	SetReactionForce (InfracolorForce *reactionForce)
const InfracolorForce *	GetReactionForce () const
InfracolorForce *	GetReactionForce ()
void	SetFirstStringLength (G4double firstStringLength)
void	StartTracking (const G4Track *dest)
void	TrackKilled ()
G4bool	IsSourceAlive () const
G4bool	IsSourceInitialized () const
G4bool	HasNextStringVector () const
const std::deque< StringVector > &	GetStringVectors () const
void	PopTo (std::deque< StringVector >::const_iterator stringPtr, G4double fracLeft)
void	PushStringVector (const StringVector &v)
void	SetStringForce (G4double stringForce)
G4double	GetStringForce () const
void	SetMaxBoost (G4double maxBoost)
G4double	GetMaxExpRapidity () const
void	SetMaxMergeT (G4double maxMergeT)
void	SetMaxMergeMag (G4double maxMergeMag)
G4int	GetNStrings () const
G4LorentzVector	GetSumStrings () const
G4ThreeVector	GetAngMomentum () const
G4ThreeVector	GetMomentOfE () const
void	Clear ()

Private Member Functions
void	CombineStringVector (const StringVector &v)

Private Attributes
G4double	m_stringForce
InfracolorForce *	m_reactionForce
G4LorentzVector	m_initU
G4bool	m_initialized
G4bool	m_killed
G4bool	m_firstStep
std::deque< StringVector >	m_stringVectors
StringVector	m_firstString
StringVector	m_borrowedString
G4double	m_firstStringLength
G4double	m_maxExpRapidity
G4double	m_maxMergeT
G4double	m_maxMergeMag

Detailed Description

Definition at line 14 of file InfracolorForce.h.

Constructor & Destructor Documentation

◆ InfracolorForce()

InfracolorForce::InfracolorForce ( )

Definition at line 10 of file InfracolorForce.cxx.

                                 :
    m_stringForce(0.0),
    m_reactionForce(0),
    m_initU(),
    m_initialized(false),
    m_killed(false),
    m_firstStep(false),
    m_stringVectors(),
    m_firstString(),
    m_borrowedString(),
    m_firstStringLength(1*CLHEP::nm),
    m_maxExpRapidity(1e-3),
    m_maxMergeT(0.0),
    m_maxMergeMag(0.0)
{}

◆ ~InfracolorForce()

InfracolorForce::~InfracolorForce ( )

Definition at line 26 of file InfracolorForce.cxx.

26{}

Member Function Documentation

◆ Clear()

void InfracolorForce::Clear ( )

Definition at line 57 of file InfracolorForce.cxx.

                            {
    m_initU.set(0,0,0,0);
    m_initialized = false;
    m_killed = false;
    m_firstStep = false;
    m_stringVectors.clear();
    m_firstString.set(0,0,0,0);
    m_borrowedString.set(0,0,0,0);
}

◆ CombineStringVector()

void InfracolorForce::CombineStringVector ( const StringVector & v )

private

Definition at line 112 of file InfracolorForce.cxx.

                                                               {
    if (m_stringVectors.empty()) {
        m_stringVectors.push_back(v);
    } else {
        G4LorentzVector sum = m_stringVectors.back().lv() + v.lv();
        if (sum.t() < m_maxMergeT && sum.m2() < sqr(m_maxMergeMag)) {
            m_stringVectors.back().set(sum.vect(), sum.m());
        } else {
            m_stringVectors.push_back(v);
        }
    }
}

◆ GetAngMomentum()

G4ThreeVector InfracolorForce::GetAngMomentum ( ) const

Definition at line 134 of file InfracolorForce.cxx.

                                                    {
    std::deque<StringVector>::const_iterator it = m_stringVectors.begin();
    G4LorentzVector x(0,0,0,0);
    G4ThreeVector L(0,0,0);
    for (; it != m_stringVectors.end(); ++it) {
        G4LorentzVector dx = it->lv();
        L += m_stringForce * x.vect().cross(dx.vect());
        x += dx;
    }
    return L;
}

◆ GetMaxExpRapidity()

G4double InfracolorForce::GetMaxExpRapidity ( ) const

inline

Definition at line 107 of file InfracolorForce.h.

                                                         {
    return m_maxExpRapidity;
}

◆ GetMomentOfE()

G4ThreeVector InfracolorForce::GetMomentOfE ( ) const

Definition at line 146 of file InfracolorForce.cxx.

                                                  {
    std::deque<StringVector>::const_iterator it = m_stringVectors.begin();
    G4LorentzVector x(0,0,0,0);
    G4ThreeVector Excm(0,0,0);
    for (; it != m_stringVectors.end(); ++it) {
        G4LorentzVector dx = it->lv();
        Excm += m_stringForce * (dx.t() * x.vect() - dx.vect() * x.t());
        x += dx;
    }
    return Excm;
}

◆ GetNStrings()

G4int InfracolorForce::GetNStrings ( ) const

inline

Definition at line 121 of file InfracolorForce.h.

                                                {
    return m_stringVectors.size();
}

◆ GetReactionForce() [1/2]

InfracolorForce * InfracolorForce::GetReactionForce ( )

inline

Definition at line 68 of file InfracolorForce.h.

                                                          {
    return m_reactionForce;
}

◆ GetReactionForce() [2/2]

const InfracolorForce * InfracolorForce::GetReactionForce ( ) const

inline

Definition at line 64 of file InfracolorForce.h.

                                                                      {
    return m_reactionForce;
}

◆ GetStringForce()

G4double InfracolorForce::GetStringForce ( ) const

inline

Definition at line 98 of file InfracolorForce.h.

                                                      {
    return m_stringForce;
}

◆ GetStringVectors()

const std::deque< StringVector > & InfracolorForce::GetStringVectors ( ) const

inline

Definition at line 89 of file InfracolorForce.h.

                                                                             {
    return m_stringVectors;
}

◆ GetSumStrings()

G4LorentzVector InfracolorForce::GetSumStrings ( ) const

Definition at line 125 of file InfracolorForce.cxx.

                                                     {
    G4LorentzVector x(0,0,0,0);
    std::deque<StringVector>::const_iterator it = m_stringVectors.begin();
    for (; it != m_stringVectors.end(); ++it) {
        x += it->lv();
    }
    return x - m_borrowedString.lv();
}

◆ HasNextStringVector()

G4bool InfracolorForce::HasNextStringVector ( ) const

inline

Definition at line 85 of file InfracolorForce.h.

                                                         {
    return !m_stringVectors.empty();
}

◆ IsSourceAlive()

G4bool InfracolorForce::IsSourceAlive ( ) const

inline

Definition at line 77 of file InfracolorForce.h.

                                                   {
    return !m_reactionForce->m_killed;
}

◆ IsSourceInitialized()

G4bool InfracolorForce::IsSourceInitialized ( ) const

inline

Definition at line 81 of file InfracolorForce.h.

                                                         {
    return m_reactionForce->m_initialized;
}

◆ PopTo()

void InfracolorForce::PopTo	(	std::deque< StringVector >::const_iterator	stringPtr,
		G4double	fracLeft )

Definition at line 72 of file InfracolorForce.cxx.

                                                                                             {
    if (fracLeft > 1.0 || fracLeft < 0.0) {
        Clear();
        m_reactionForce->Clear();
        G4Exception("InfracolorForce::PopTo", "QuirkStringBadFraction", EventMustBeAborted, "invalid fraction of string vector");
    }
    if (m_firstStep) {
        if (stringPtr == m_stringVectors.begin()) {
            m_borrowedString = (1 - fracLeft) * m_firstString;
        } else {
            m_borrowedString = m_firstString;
        }
        m_firstStep = false;
        m_stringVectors.clear();
    } else {
        std::deque<StringVector>::iterator stringPtr2 = m_stringVectors.begin() + (stringPtr - m_stringVectors.begin());
        m_stringVectors.erase(m_stringVectors.begin(), stringPtr2);
        if (fracLeft != 1.0 && !m_stringVectors.empty()) {
            m_stringVectors[0] *= fracLeft;
        }
    }
}

◆ PushStringVector()

void InfracolorForce::PushStringVector ( const StringVector & v )

Definition at line 95 of file InfracolorForce.cxx.

                                                            {
    if (v.t() == 0) return;
    if (m_borrowedString.t() == 0) {
        CombineStringVector(v);
    } else if (m_borrowedString.t() <= v.t()) {
        G4double r = m_borrowedString.t() / v.t();
        CombineStringVector((1-r) * v);
        m_borrowedString.set(0,0,0,0);
    } else {
        G4double r = v.t() / m_borrowedString.t();
        m_borrowedString *= 1-r;
        G4Exception("InfracolorForce::PushStringVector()",
                    "BorrowedStringSplit", JustWarning,
                    "Initial step too long.");
    }
}

◆ SetFirstStringLength()

void InfracolorForce::SetFirstStringLength ( G4double firstStringLength )

inline

Definition at line 72 of file InfracolorForce.h.

                                                                            {
    m_firstStringLength = firstStringLength;
    m_reactionForce->m_firstStringLength = firstStringLength;
}

◆ SetMaxBoost()

void InfracolorForce::SetMaxBoost ( G4double maxBoost )

inline

Definition at line 102 of file InfracolorForce.h.

                                                          {
    m_maxExpRapidity = std::sqrt((1.0 + maxBoost) / (1.0 - maxBoost));
    m_reactionForce->m_maxExpRapidity = m_maxExpRapidity;
}

◆ SetMaxMergeMag()

void InfracolorForce::SetMaxMergeMag ( G4double maxMergeMag )

inline

Definition at line 116 of file InfracolorForce.h.

                                                                {
    m_maxMergeMag = maxMergeMag;
    m_reactionForce->m_maxMergeMag = maxMergeMag;
}

◆ SetMaxMergeT()

void InfracolorForce::SetMaxMergeT ( G4double maxMergeT )

inline

Definition at line 111 of file InfracolorForce.h.

                                                            {
    m_maxMergeT = maxMergeT;
    m_reactionForce->m_maxMergeT = maxMergeT;
}

◆ SetReactionForce()

void InfracolorForce::SetReactionForce ( InfracolorForce * reactionForce )

inline

Definition at line 60 of file InfracolorForce.h.

                                                                            {
    m_reactionForce = reactionForce;
}

◆ SetStringForce()

void InfracolorForce::SetStringForce ( G4double stringForce )

inline

Definition at line 93 of file InfracolorForce.h.

                                                                {
    m_stringForce = stringForce;
    m_reactionForce->m_stringForce = stringForce;
}

◆ StartTracking()

void InfracolorForce::StartTracking ( const G4Track * dest )

Definition at line 28 of file InfracolorForce.cxx.

                                                       {
    if (m_reactionForce == 0) {
        G4Exception("InfracolorForce::StartTracking", "NoAntiQuirk", FatalErrorInArgument, "InfracolorForce: No antiquirk defined");
    }
    if (dest->GetCurrentStepNumber() != 0) return;
 
    // Clear old stuff
    if (m_initialized && m_reactionForce->m_initialized) {
        Clear();
        m_reactionForce->Clear();
    }
 
    // Set momentum
    m_initialized = true;
    m_initU = dest->GetDynamicParticle()->Get4Momentum();
    m_initU /= m_initU.m();
 
    // Compute first string vector
    if (m_reactionForce->m_initialized) {
        m_firstStep = true;
        G4double dot = m_initU * m_reactionForce->m_initU;
        G4LorentzVector firstString = m_reactionForce->m_initU - m_initU * (dot - std::sqrt(dot*dot - 1));
        firstString *= m_firstStringLength / firstString.t();
        m_firstString.set(firstString.vect(), 0);
        m_stringVectors.push_back(m_firstString);
        // TODO: Check oscillation period
    }
}

◆ TrackKilled()

void InfracolorForce::TrackKilled ( )

Definition at line 67 of file InfracolorForce.cxx.

                                  {
    //G4cout << "killed, and other killed is " << m_reactionForce->m_killed << G4endl;
    m_killed = true;
}

Member Data Documentation

◆ m_borrowedString

StringVector InfracolorForce::m_borrowedString

private

Definition at line 53 of file InfracolorForce.h.

◆ m_firstStep

G4bool InfracolorForce::m_firstStep

private

Definition at line 50 of file InfracolorForce.h.

◆ m_firstString

StringVector InfracolorForce::m_firstString

private

Definition at line 52 of file InfracolorForce.h.

◆ m_firstStringLength

G4double InfracolorForce::m_firstStringLength

private

Definition at line 54 of file InfracolorForce.h.

◆ m_initialized

G4bool InfracolorForce::m_initialized

private

Definition at line 48 of file InfracolorForce.h.

◆ m_initU

G4LorentzVector InfracolorForce::m_initU

private

Definition at line 47 of file InfracolorForce.h.

◆ m_killed

G4bool InfracolorForce::m_killed

private

Definition at line 49 of file InfracolorForce.h.

◆ m_maxExpRapidity

G4double InfracolorForce::m_maxExpRapidity

private

Definition at line 55 of file InfracolorForce.h.

◆ m_maxMergeMag

G4double InfracolorForce::m_maxMergeMag

private

Definition at line 57 of file InfracolorForce.h.

◆ m_maxMergeT

G4double InfracolorForce::m_maxMergeT

private

Definition at line 56 of file InfracolorForce.h.

◆ m_reactionForce

InfracolorForce* InfracolorForce::m_reactionForce

private

Definition at line 46 of file InfracolorForce.h.

◆ m_stringForce

G4double InfracolorForce::m_stringForce

private

Definition at line 45 of file InfracolorForce.h.

◆ m_stringVectors

std::deque<StringVector> InfracolorForce::m_stringVectors

private

Definition at line 51 of file InfracolorForce.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ InfracolorForce()

◆ ~InfracolorForce()

Member Function Documentation

◆ Clear()

◆ CombineStringVector()

◆ GetAngMomentum()

◆ GetMaxExpRapidity()

◆ GetMomentOfE()

◆ GetNStrings()

◆ GetReactionForce() [1/2]

◆ GetReactionForce() [2/2]

◆ GetStringForce()

◆ GetStringVectors()

◆ GetSumStrings()

◆ HasNextStringVector()

◆ IsSourceAlive()

◆ IsSourceInitialized()

◆ PopTo()

◆ PushStringVector()

◆ SetFirstStringLength()

◆ SetMaxBoost()

◆ SetMaxMergeMag()

◆ SetMaxMergeT()

◆ SetReactionForce()

◆ SetStringForce()

◆ StartTracking()

◆ TrackKilled()

Member Data Documentation

◆ m_borrowedString

◆ m_firstStep

◆ m_firstString

◆ m_firstStringLength

◆ m_initialized

◆ m_initU

◆ m_killed

◆ m_maxExpRapidity

◆ m_maxMergeMag

◆ m_maxMergeT

◆ m_reactionForce

◆ m_stringForce

◆ m_stringVectors