G4VElasticCollision Class Reference

#include <G4VElasticCollision.hh>

Inheritance diagram for G4VElasticCollision:

Public Member Functions
	G4VElasticCollision ()
virtual	~G4VElasticCollision ()
G4bool	operator== (const G4VElasticCollision &right) const
G4bool	operator!= (const G4VElasticCollision &right) const
virtual G4KineticTrackVector *	FinalState (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
Public Member Functions inherited from G4VCollision
	G4VCollision ()
void	establish_G4MT_TLS_G4VCollision ()
	G4VCollision (void *s1, void *s2, void *s3, void *s4, void *s5, void *s6, void *s7)
virtual	~G4VCollision ()
G4bool	operator== (const G4VCollision &right) const
G4bool	operator!= (const G4VCollision &right) const
virtual G4double	CrossSection (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const
virtual G4bool	IsInCharge (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const =0
virtual G4String	GetName () const =0
virtual void	Print () const
virtual void	Print (const G4KineticTrack &trk1, const G4KineticTrack &trk2) const

Additional Inherited Members
Protected Member Functions inherited from G4VCollision
G4int	GetNumberOfPartons (G4ParticleDefinition *aP) const
virtual const G4CollisionVector *	GetComponents () const
virtual const G4VCrossSectionSource *	GetCrossSectionSource () const =0
virtual const G4VAngularDistribution *	GetAngularDistribution () const =0
virtual const std::vector < G4String > &	GetListOfColliders (G4int whichOne) const =0

Detailed Description

Definition at line 40 of file G4VElasticCollision.hh.

Constructor & Destructor Documentation

G4VElasticCollision::G4VElasticCollision

(

)

Definition at line 45 of file G4VElasticCollision.cc.

{ 
}

G4VElasticCollision::~G4VElasticCollision ( )

virtual

Definition at line 50 of file G4VElasticCollision.cc.

{ }

Member Function Documentation

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

virtual

Implements G4VCollision.

Definition at line 54 of file G4VElasticCollision.cc.

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

{
  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;
}

G4bool G4VElasticCollision::operator!=

(

const G4VElasticCollision &

right

)

const

G4bool G4VElasticCollision::operator==

(

const G4VElasticCollision &

right

)

const

---

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

• G4VElasticCollision.hh
• G4VElasticCollision.cc