G4CollisionMesonBaryonElastic Class Reference

#include <G4CollisionMesonBaryonElastic.hh>

Inheritance diagram for G4CollisionMesonBaryonElastic:

Public Member Functions
virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
void	establish_G4MT_TLS_G4VCollision ()
virtual G4KineticTrackVector *	FinalState (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
	G4CollisionMesonBaryonElastic ()
virtual G4String	GetName () const
virtual G4bool	IsInCharge (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
G4bool	operator!= (const G4CollisionMesonBaryonElastic &right) const
G4bool	operator!= (const G4VCollision &right) const
G4bool	operator!= (const G4VElasticCollision &right) const
G4bool	operator== (const G4CollisionMesonBaryonElastic &right) const
G4bool	operator== (const G4VCollision &right) const
G4bool	operator== (const G4VElasticCollision &right) const
virtual void	Print () const
virtual void	Print (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual	~G4CollisionMesonBaryonElastic ()

Protected Member Functions
virtual const G4VAngularDistribution *	GetAngularDistribution () const
virtual const G4CollisionVector *	GetComponents () const
virtual const G4VCrossSectionSource *	GetCrossSectionSource () const
virtual const std::vector< G4String > &	GetListOfColliders (G4int) const
G4int	GetNumberOfPartons (const G4ParticleDefinition *aP) const

Private Member Functions
	G4CollisionMesonBaryonElastic (const G4CollisionMesonBaryonElastic &)
G4CollisionMesonBaryonElastic &	operator= (const G4CollisionMesonBaryonElastic &)

Private Attributes
G4VAngularDistribution *	angularDistribution
G4VCrossSectionSource *	crossSectionSource
std::vector< G4String >	dummy

Detailed Description

Definition at line 41 of file G4CollisionMesonBaryonElastic.hh.

Constructor & Destructor Documentation

G4CollisionMesonBaryonElastic() [1/2]

G4CollisionMesonBaryonElastic::G4CollisionMesonBaryonElastic

(

)

Definition at line 35 of file G4CollisionMesonBaryonElastic.cc.

{
  angularDistribution = new G4AngularDistribution(false);
  crossSectionSource = new G4XMesonBaryonElastic();
}

References angularDistribution, and crossSectionSource.

~G4CollisionMesonBaryonElastic()

G4CollisionMesonBaryonElastic::~G4CollisionMesonBaryonElastic ( )

virtual

Definition at line 41 of file G4CollisionMesonBaryonElastic.cc.

{
  delete angularDistribution;
  delete crossSectionSource;
}

References angularDistribution, and crossSectionSource.

G4CollisionMesonBaryonElastic() [2/2]

G4CollisionMesonBaryonElastic::G4CollisionMesonBaryonElastic ( const G4CollisionMesonBaryonElastic &  )

private

Member Function Documentation

CrossSection()

G4double G4VCollision::CrossSection ( const G4KineticTrack &  trk1, const G4KineticTrack &  trk2  ) const

virtualinherited

Reimplemented in G4CollisionComposite, G4CollisionMesonBaryonToResonance, and G4CollisionNN.

Definition at line 54 of file G4VCollision.cc.

{
    G4double sigma = 0.;
 
    const G4VCrossSectionSource* xSource = GetCrossSectionSource();
 
    if (xSource != 0)
    {
        // There is a cross section for this Collision
        sigma = xSource->CrossSection(aTrk1,aTrk2);
    }
    return sigma;
}

References G4VCrossSectionSource::CrossSection(), and G4VCollision::GetCrossSectionSource().

Referenced by G4Scatterer::GetCrossSection(), G4Scatterer::GetTimeToInteraction(), G4VCollision::Print(), and G4Scatterer::Scatter().

establish_G4MT_TLS_G4VCollision()

void G4VCollision::establish_G4MT_TLS_G4VCollision ( )

inherited

Definition at line 143 of file G4VCollision.cc.

{ }

Referenced by G4VScatteringCollision::establish_G4MT_TLS_G4VScatteringCollision().

FinalState()

G4KineticTrackVector * G4VElasticCollision::FinalState ( const G4KineticTrack &  trk1, const G4KineticTrack &  trk2  ) const

virtualinherited

Implements G4VCollision.

Definition at line 53 of file G4VElasticCollision.cc.

{
  const G4VAngularDistribution* angDistribution;
 
  angDistribution = GetAngularDistribution();
 
 
  G4LorentzVector pCM=trk1.Get4Momentum() + trk2.Get4Momentum();
 
  G4LorentzRotation toLabFrame(pCM.boostVector());
  G4LorentzVector Ptmp=toLabFrame.inverse() * trk1.Get4Momentum();  //trk1 in CMS
  G4LorentzRotation toZ;
  toZ.rotateZ(-Ptmp.phi());
  toZ.rotateY(-Ptmp.theta());
  toLabFrame *= toZ.inverse();
 
  G4double S = pCM.mag2();
  G4double m10 = trk1.GetDefinition()->GetPDGMass();
  G4double m20 = trk2.GetDefinition()->GetPDGMass();
  if(S-(m10+m20)*(m10+m20) < 0) return new G4KineticTrackVector;
 
  G4double m_1 = trk1.GetActualMass();
  G4double m_2 = trk2.GetActualMass();
  
  // Angles of outgoing particles
  G4double cosTheta = angDistribution->CosTheta(S,m_1,m_2);
 
   if ( (trk1.GetDefinition() == G4Proton::Proton() || trk1.GetDefinition() == G4Neutron::Neutron() )
      &&(trk2.GetDefinition() == G4Proton::Proton() || trk2.GetDefinition() == G4Neutron::Neutron() ) )
   {
      if ( trk1.GetDefinition() == trk2.GetDefinition() )
      {
      if ( trk1.GetDefinition() == G4Proton::Proton() )
      {
//        G4cout << "scatterangle pp " << cosTheta
//               << " " << typeid(*angDistribution).name() << G4endl;
      } else {
//        G4cout << "scatterangle nn " << cosTheta
//               << " " << typeid(*angDistribution).name() << G4endl;
      }
      } else {
//    G4cout << "scatterangle pn " << cosTheta
//               << " " << typeid(*angDistribution).name() << G4endl;
      }
   } else {
//      G4cout << "scatterangle other " << cosTheta
//               << " " << typeid(*angDistribution).name() << G4endl;
   }
 
  G4double phi = angDistribution->Phi();
  G4double Theta = std::acos(cosTheta);
 
  // Unit vector of three-momentum
  G4ThreeVector pFinal1(std::sin(Theta)*std::cos(phi), std::sin(Theta)*std::sin(phi), cosTheta);
  // Three momentum in cm system
  G4double pInCM = std::sqrt((S-(m10+m20)*(m10+m20))*(S-(m10-m20)*(m10-m20))/(4.*S));
  pFinal1 = pFinal1 * pInCM;
  G4ThreeVector pFinal2 = -pFinal1;
 
  G4double eFinal1 = std::sqrt(pFinal1.mag2() + m10*m10);
  G4double eFinal2 = std::sqrt(pFinal2.mag2() + m20*m20);
 
  G4LorentzVector p4Final1(pFinal1, eFinal1);
  G4LorentzVector p4Final2(pFinal2, eFinal2);
 
  // Lorentz transformation
  p4Final1 *= toLabFrame;
  p4Final2 *= toLabFrame;
 
  // Final tracks are copies of incoming ones, with modified 4-momenta
  G4KineticTrack* final1 = new G4KineticTrack(trk1);
  final1->Set4Momentum(p4Final1);
  G4KineticTrack* final2 = new G4KineticTrack(trk2);
  final2->Set4Momentum(p4Final2);
 
  G4KineticTrackVector* finalTracks = new G4KineticTrackVector;
  finalTracks->push_back(final1);
  finalTracks->push_back(final2);
 
  return finalTracks;
}

References CLHEP::HepLorentzVector::boostVector(), G4VAngularDistribution::CosTheta(), G4KineticTrack::Get4Momentum(), G4KineticTrack::GetActualMass(), G4VCollision::GetAngularDistribution(), G4KineticTrack::GetDefinition(), G4ParticleDefinition::GetPDGMass(), CLHEP::HepLorentzRotation::inverse(), CLHEP::HepLorentzVector::mag2(), CLHEP::Hep3Vector::mag2(), G4Neutron::Neutron(), CLHEP::HepLorentzVector::phi(), G4VAngularDistribution::Phi(), G4Proton::Proton(), CLHEP::HepLorentzRotation::rotateY(), CLHEP::HepLorentzRotation::rotateZ(), S(), G4KineticTrack::Set4Momentum(), and CLHEP::HepLorentzVector::theta().

GetAngularDistribution()

virtual const G4VAngularDistribution * G4CollisionMesonBaryonElastic::GetAngularDistribution ( ) const

inlineprotectedvirtual

Implements G4VCollision.

Definition at line 64 of file G4CollisionMesonBaryonElastic.hh.

  { return angularDistribution; }

References angularDistribution.

GetComponents()

virtual const G4CollisionVector * G4VCollision::GetComponents ( ) const

inlineprotectedvirtualinherited

Reimplemented in G4CollisionComposite, G4CollisionNN, G4CollisionNNToDeltaDelta1600, G4CollisionNNToDeltaDelta1620, G4CollisionNNToDeltaDelta1700, G4CollisionNNToDeltaDelta1900, G4CollisionNNToDeltaDelta1905, G4CollisionNNToDeltaDelta1910, G4CollisionNNToDeltaDelta1920, G4CollisionNNToDeltaDelta1930, G4CollisionNNToDeltaDelta1950, G4CollisionNNToNDelta, G4CollisionNNToNDelta1600, G4CollisionNNToNDelta1620, G4CollisionNNToNDelta1700, G4CollisionNNToNDelta1900, G4CollisionNNToNDelta1905, G4CollisionNNToNDelta1910, G4CollisionNNToNDelta1920, G4CollisionNNToNDelta1930, and G4CollisionNNToNDelta1950.

Definition at line 79 of file G4VCollision.hh.

{ return 0;}

Referenced by G4VCollision::Print().

GetCrossSectionSource()

virtual const G4VCrossSectionSource * G4CollisionMesonBaryonElastic::GetCrossSectionSource ( ) const

inlineprotectedvirtual

Implements G4VCollision.

Definition at line 62 of file G4CollisionMesonBaryonElastic.hh.

  { return crossSectionSource; }

References crossSectionSource.

GetListOfColliders()

const std::vector< G4String > & G4CollisionMesonBaryonElastic::GetListOfColliders ( G4int  ) const

protectedvirtual

Implements G4VCollision.

Definition at line 63 of file G4CollisionMesonBaryonElastic.cc.

 {
   throw G4HadronicException(__FILE__, __LINE__, "Called G4CollisionMesonBaryonElastic::GetListOfColliders");
   return dummy;
 }

References dummy.

GetName()

G4String G4CollisionMesonBaryonElastic::GetName ( ) const

virtual

Implements G4VCollision.

Definition at line 57 of file G4CollisionMesonBaryonElastic.cc.

 {
   return "Meson Baryon Elastic Collision";
 }

GetNumberOfPartons()

G4int G4VCollision::GetNumberOfPartons ( const G4ParticleDefinition *  aP ) const

inlineprotectedinherited

Definition at line 68 of file G4VCollision.hh.

  {
    G4int result = 0;
    for(G4int i=0; i<6; i++) 
    {
      result += aP->GetQuarkContent(i+1);
      result += aP->GetAntiQuarkContent(i+1);
    }
    return result;
  }

References G4ParticleDefinition::GetAntiQuarkContent(), and G4ParticleDefinition::GetQuarkContent().

Referenced by IsInCharge().

IsInCharge()

G4bool G4CollisionMesonBaryonElastic::IsInCharge ( const G4KineticTrack &  trk1, const G4KineticTrack &  trk2  ) const

virtual

Implements G4VCollision.

Definition at line 47 of file G4CollisionMesonBaryonElastic.cc.

 {
   G4int partons1 = GetNumberOfPartons(trk1.GetDefinition());
   G4int partons2 = GetNumberOfPartons(trk2.GetDefinition());
   G4bool result = (partons1 == 2 && partons2 ==3) ||
                   (partons2 == 2 && partons1 ==3);
   return result;
 }

References G4KineticTrack::GetDefinition(), and G4VCollision::GetNumberOfPartons().

operator!=() [1/3]

G4bool G4CollisionMesonBaryonElastic::operator!=

(

const G4CollisionMesonBaryonElastic &

right

)

const

operator!=() [2/3]

G4bool G4VCollision::operator!= ( const G4VCollision &  right ) const

inherited

Definition at line 48 of file G4VCollision.cc.

{
    return (this != (G4VCollision *) &right);
}

operator!=() [3/3]

G4bool G4VElasticCollision::operator!= ( const G4VElasticCollision &  right ) const

inherited

operator=()

G4CollisionMesonBaryonElastic & G4CollisionMesonBaryonElastic::operator= ( const G4CollisionMesonBaryonElastic &  )

private

operator==() [1/3]

G4bool G4CollisionMesonBaryonElastic::operator==

(

const G4CollisionMesonBaryonElastic &

right

)

const

operator==() [2/3]

G4bool G4VCollision::operator== ( const G4VCollision &  right ) const

inherited

Definition at line 42 of file G4VCollision.cc.

{
    return (this == (G4VCollision *) &right);
}

operator==() [3/3]

G4bool G4VElasticCollision::operator== ( const G4VElasticCollision &  right ) const

inherited

Print() [1/2]

void G4VCollision::Print ( ) const

virtualinherited

Definition at line 70 of file G4VCollision.cc.

{
    G4String name = GetName();
 
    G4cout << "---- " << name << "---- Cross section" << G4endl;
 
    const G4VCrossSectionSource* xSource = GetCrossSectionSource();
    if (xSource) xSource->Print();
 
    G4int nComponents = 0;
    const G4CollisionVector* components = GetComponents();
    if (components)
    {
        nComponents = components->size();
    }
    G4cout << "---- " << name << "---- has " << nComponents << " components" <<G4endl;
    G4int i = 0;
    G4CollisionVector::const_iterator iter;
    if (components)
    {
        for (iter = components->begin(); iter != components->end(); ++iter)
        {
            G4cout << "---- " << name << " ---- Component " << i << G4endl;
            ((*iter))->Print();
            i++;
        }
    }
 
}

References G4cout, G4endl, G4VCollision::GetComponents(), G4VCollision::GetCrossSectionSource(), G4VCollision::GetName(), G4InuclParticleNames::name(), G4VCollision::Print(), and G4VCrossSectionSource::Print().

Referenced by G4VCollision::Print().

Print() [2/2]

void G4VCollision::Print ( const G4KineticTrack &  trk1, const G4KineticTrack &  trk2  ) const

virtualinherited

Definition at line 101 of file G4VCollision.cc.

{
    G4String name = GetName();
 
    if (IsInCharge(trk1,trk2))
    {
        G4cout << "---- " << name << "is in charge ---- " << G4endl;
    }
    else
    {
        G4cout << "---- " << name << "is not in charge ---- " << G4endl;
    }
 
    G4cout << "---- " << name << "---- Cross section" << G4endl;
 
    const G4VCrossSectionSource* xSource = GetCrossSectionSource();
    if (xSource) xSource->Print();
    G4cout << "Cross section = " << CrossSection(trk1,trk2) << G4endl;
 
    G4int nComponents = 0;
    const G4CollisionVector* components = GetComponents();
    if (components)
    {
        nComponents = components->size();
    }
    G4cout << "---- " << name << "has " << nComponents << " components" <<G4endl;
 
    G4int i = 0;
    G4CollisionVector::const_iterator iter;
    if (components)
    {
        for (iter = components->begin(); iter != components->end(); ++iter)
        {
            G4cout << "Component " << i << G4endl;
            ((*iter))->Print();
            i++;
        }
    }
}

References G4VCollision::CrossSection(), G4cout, G4endl, G4VCollision::GetComponents(), G4VCollision::GetCrossSectionSource(), G4VCollision::GetName(), G4VCollision::IsInCharge(), G4InuclParticleNames::name(), G4VCollision::Print(), and G4VCrossSectionSource::Print().

Field Documentation

angularDistribution

G4VAngularDistribution* G4CollisionMesonBaryonElastic::angularDistribution

private

Definition at line 71 of file G4CollisionMesonBaryonElastic.hh.

Referenced by G4CollisionMesonBaryonElastic(), GetAngularDistribution(), and ~G4CollisionMesonBaryonElastic().

crossSectionSource

G4VCrossSectionSource* G4CollisionMesonBaryonElastic::crossSectionSource

private

Definition at line 70 of file G4CollisionMesonBaryonElastic.hh.

Referenced by G4CollisionMesonBaryonElastic(), GetCrossSectionSource(), and ~G4CollisionMesonBaryonElastic().

dummy

std::vector<G4String> G4CollisionMesonBaryonElastic::dummy

private

Definition at line 73 of file G4CollisionMesonBaryonElastic.hh.

Referenced by GetListOfColliders().

---

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

• source/processes/hadronic/models/im_r_matrix/include/G4CollisionMesonBaryonElastic.hh
• source/processes/hadronic/models/im_r_matrix/src/G4CollisionMesonBaryonElastic.cc