52 G4cout <<
"G4ParticleChangeForMSC::G4ParticleChangeForMSC() " <<
G4endl;
61 G4cout <<
"G4ParticleChangeForMSC::~G4ParticleChangeForMSC() " <<
G4endl;
71 G4cout <<
"G4ParticleChangeForMSC:: copy constructor is called " <<
G4endl;
73 theMomentumDirection = right.theMomentumDirection;
74 thePosition = right.thePosition;
82 G4cout <<
"G4ParticleChangeForMSC:: assignment operator is called " <<
G4endl;
94 theMomentumDirection = right.theMomentumDirection;
95 thePosition = right.thePosition;
147 G4cout <<
" Position - x (mm) : "
148 << std::setw(20) << thePosition.
x()/
mm
150 G4cout <<
" Position - y (mm) : "
151 << std::setw(20) << thePosition.
y()/
mm
153 G4cout <<
" Position - z (mm) : "
154 << std::setw(20) << thePosition.
z()/
mm
156 G4cout <<
" Momentum Direct - x : "
157 << std::setw(20) << theMomentumDirection.
x()
159 G4cout <<
" Momentum Direct - y : "
160 << std::setw(20) << theMomentumDirection.
y()
162 G4cout <<
" Momentum Direct - z : "
163 << std::setw(20) << theMomentumDirection.
z()
172 G4bool exitWithError =
false;
179 accuracy = std::fabs(theMomentumDirection.
mag2()-1.0);
184 G4cout <<
" G4ParticleChangeForMSC::CheckIt : ";
185 G4cout <<
"the Momentum Change is not unit vector !!"
186 <<
" Difference: " << accuracy <<
G4endl;
206 "momentum direction was illegal");
211 theMomentumDirection = (1./vmag)*theMomentumDirection;