G4SampleResonance Class Reference

#include <G4SampleResonance.hh>

Collaboration diagram for G4SampleResonance:

Public Types
typedef std::map< const G4ParticleDefinition *, G4double, std::less< const G4ParticleDefinition * > >::const_iterator	minMassMapIterator
typedef std::map< const G4ParticleDefinition *, G4double, std::less< const G4ParticleDefinition * > >	minMassMapType

Public Member Functions
G4double	GetMinimumMass (const G4ParticleDefinition *p) const
G4double	SampleMass (const G4double poleMass, const G4double gamma, const G4double minMass, const G4double maxMass) const
G4double	SampleMass (const G4ParticleDefinition *p, const G4double maxMass) const

Private Member Functions
G4double	BrWigInt0 (const G4double x, const G4double gamma, const G4double m0) const
G4double	BrWigInt1 (const G4double x, const G4double gamma, const G4double m0) const
G4double	BrWigInv (const G4double x, const G4double gamma, const G4double m0) const

Static Private Attributes
static G4ThreadLocal minMassMapType *	minMassCache_G4MT_TLS_ = 0

Detailed Description

Definition at line 43 of file G4SampleResonance.hh.

Member Typedef Documentation

minMassMapIterator

typedef std::map<const G4ParticleDefinition*, G4double, std::less<const G4ParticleDefinition*> >::const_iterator G4SampleResonance::minMassMapIterator

Definition at line 67 of file G4SampleResonance.hh.

minMassMapType

typedef std::map<const G4ParticleDefinition*, G4double, std::less<const G4ParticleDefinition*> > G4SampleResonance::minMassMapType

Definition at line 68 of file G4SampleResonance.hh.

Member Function Documentation

BrWigInt0()

G4double G4SampleResonance::BrWigInt0 ( const G4double  x, const G4double  gamma, const G4double  m0  ) const

inlineprivate

Definition at line 56 of file G4SampleResonance.hh.

    { return 2.0*gamma*std::atan( 2.0 * (x-m0)/ gamma  ); }

Here is the caller graph for this function:

BrWigInt1()

G4double G4SampleResonance::BrWigInt1 ( const G4double  x, const G4double  gamma, const G4double  m0  ) const

inlineprivate

Definition at line 59 of file G4SampleResonance.hh.

    { return 0.5*gamma*gamma*G4Log( (x-m0)*(x-m0)+gamma*gamma/4.0 ) + m0*BrWigInt0(x,gamma,m0); }

Here is the call graph for this function:

BrWigInv()

G4double G4SampleResonance::BrWigInv ( const G4double  x, const G4double  gamma, const G4double  m0  ) const

inlineprivate

Definition at line 62 of file G4SampleResonance.hh.

    { return 0.5*gamma*std::tan( 0.5*x/gamma )+m0; }

Here is the caller graph for this function:

GetMinimumMass()

G4double G4SampleResonance::GetMinimumMass

(

const G4ParticleDefinition *

p

)

const

Definition at line 46 of file G4SampleResonance.cc.

{  ;;;   if (!minMassCache_G4MT_TLS_) minMassCache_G4MT_TLS_ = new G4SampleResonance::minMassMapType  ; G4SampleResonance::minMassMapType &minMassCache = *minMassCache_G4MT_TLS_;  ;;;  

    G4double minResonanceMass = DBL_MAX;

    if ( p->IsShortLived() )
    {
        minMassMapIterator iter = minMassCache.find(p);
        if ( iter!=minMassCache.end() )
        {
            minResonanceMass = (*iter).second;
        }
        else
        {
            // G4cout << "--> request for " << p->GetParticleName() << G4endl;

            const G4DecayTable* theDecays = p->GetDecayTable();
            const G4int nDecays = theDecays->entries();

            for (G4int i=0; i<nDecays; i++)
            {
                const G4VDecayChannel* aDecay = theDecays->GetDecayChannel(i);
                const G4int nDaughters = aDecay->GetNumberOfDaughters();

                G4double minChannelMass = 0;

                for (G4int j=0; j<nDaughters; j++)
                {
                    const G4ParticleDefinition* aDaughter = const_cast<G4VDecayChannel*>(aDecay)->GetDaughter(j);
                    G4double minMass = GetMinimumMass(aDaughter);
                    if (!minMass) minMass = DBL_MAX; // exclude gamma channel;
                    minChannelMass+=minMass;
                }
                // G4cout << "channel mass for the above is " << minChannelMass/MeV << G4endl;
                if (minChannelMass < minResonanceMass) minResonanceMass = minChannelMass;

            }
            // replace this as soon as the compiler supports mutable!!
            G4SampleResonance* self = const_cast<G4SampleResonance*>(this);
            //Andrea Dotti (13Jan2013): Change needed for G4MT
            //(self->minMassCache)[p] = minResonanceMass;
            self->minMassCache_G4MT_TLS_->operator[](p) = minResonanceMass;

        }
    } 
    else
    {

        minResonanceMass = p->GetPDGMass();

    }
    // G4cout << "minimal mass for " << p->GetParticleName() << " is " << minResonanceMass/MeV << G4endl;

    return minResonanceMass;
} 

Here is the call graph for this function:

Here is the caller graph for this function:

SampleMass() [1/2]

G4double G4SampleResonance::SampleMass	(	const G4double	poleMass,
		const G4double	gamma,
		const G4double	minMass,
		const G4double	maxMass
	)		const

Definition at line 110 of file G4SampleResonance.cc.

{ if (!minMassCache_G4MT_TLS_) minMassCache_G4MT_TLS_ = new G4SampleResonance::minMassMapType  ;
    // Chooses a mass randomly between minMass and maxMass
    //     according to a Breit-Wigner function with constant
    //     width gamma and pole poleMass

    if ( minMass > maxMass )
    {
        throw G4HadronicException(__FILE__, __LINE__,
                "SampleResonanceMass: mass range negative (minMass>maxMass)");
    }

    G4double returnMass;

    if ( gamma < DBL_EPSILON )
    {
        returnMass = std::max(minMass, std::min(maxMass, poleMass));
    }
    else
    {
        double fmin = BrWigInt0(minMass, gamma, poleMass);
        double fmax = BrWigInt0(maxMass, gamma, poleMass);
        double f = fmin + (fmax-fmin)*G4UniformRand();
        returnMass = BrWigInv(f, gamma, poleMass);
    }

    return returnMass;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SampleMass() [2/2]

G4double G4SampleResonance::SampleMass	(	const G4ParticleDefinition *	p,
		const G4double	maxMass
	)		const

Definition at line 104 of file G4SampleResonance.cc.

{ if (!minMassCache_G4MT_TLS_) minMassCache_G4MT_TLS_ = new G4SampleResonance::minMassMapType  ;
    return SampleMass(p->GetPDGMass(), p->GetPDGWidth(), GetMinimumMass(p), maxMass);
}

Here is the call graph for this function:

Member Data Documentation

minMassCache_G4MT_TLS_

G4ThreadLocal G4SampleResonance::minMassMapType * G4SampleResonance::minMassCache_G4MT_TLS_ = 0

staticprivate

Definition at line 72 of file G4SampleResonance.hh.

---

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

• G4SampleResonance.hh
• G4SampleResonance.cc