136 std::vector< G4QMDNucleus* > nucleuses;
151 G4double beta_nncm = ( - boostLABtoCM.beta() + boostLABtoNN.beta() ) / ( 1 - boostLABtoCM.beta() * boostLABtoNN.beta() ) ;
159 boostToReac = boostLABtoNN;
160 boostBackToLAB = -boostLABtoNN;
216 targ->CalEnergyAndAngularMomentumInCM();
223 for (
G4int i = 0 ; i < targ->GetTotalNumberOfParticipant() ; i++ )
279 G4int i = targ->GetTotalNumberOfParticipant();
316 if ( i / 10 * 10 == i )
341 if ( numberOfSecondary == 2 )
344 G4bool elasticLike_system =
false;
345 if ( nucleuses.size() == 2 )
349 sec_a_A = nucleuses[0]->GetMassNumber();
350 sec_b_Z = nucleuses[1]->GetAtomicNumber();
351 sec_b_A = nucleuses[1]->GetMassNumber();
353 if ( ( sec_a_Z == proj_Z && sec_a_A == proj_A && sec_b_Z == targ_Z && sec_b_A == targ_A )
354 || ( sec_a_Z == targ_Z && sec_a_A == targ_A && sec_b_Z == proj_Z && sec_b_A == proj_A ) )
356 elasticLike_system =
true;
360 else if ( nucleuses.size() == 1 )
363 sec_a_Z = nucleuses[0]->GetAtomicNumber();
364 sec_a_A = nucleuses[0]->GetMassNumber();
367 if ( ( sec_a_Z == proj_Z && sec_a_A == proj_A && sec_b_pd == targ_pd )
368 || ( sec_a_Z == targ_Z && sec_a_A == targ_A && sec_b_pd == proj_pd ) )
370 elasticLike_system =
true;
380 if ( ( sec_a_pd == proj_pd && sec_b_pd == targ_pd )
381 || ( sec_a_pd == targ_pd && sec_b_pd == proj_pd ) )
383 elasticLike_system =
true;
388 if ( elasticLike_system ==
true )
391 G4bool elasticLike_energy =
true;
393 for (
G4int i = 0 ; i < int ( nucleuses.size() ) ; i++ )
401 if ( nucleuses[i]->GetExcitationEnergy()*
GeV > 1.0*
MeV ) elasticLike_energy =
false;
406 G4bool withCollision =
true;
418 if ( elasticLike_energy ==
false ) elastic =
false;
420 if ( elasticLike_energy ==
false && withCollision ==
true ) elastic =
false;
437 if ( elastic ==
true )
440 for ( std::vector< G4QMDNucleus* >::iterator
441 it = nucleuses.begin() ; it != nucleuses.end() ; it++ )
452 for ( std::vector< G4QMDNucleus* >::iterator it
453 = nucleuses.begin() ; it != nucleuses.end() ; it++ )
474 if ( (*it)->GetAtomicNumber() == 0
475 || (*it)->GetAtomicNumber() == (*it)->GetMassNumber() )
478 for (
G4int i = 0 ; i < (*it)->GetTotalNumberOfParticipant() ; i++ )
480 G4QMDParticipant* aP =
new G4QMDParticipant( ( (*it)->GetParticipant( i ) )->GetDefinition() , ( (*it)->GetParticipant( i ) )->GetMomentum() , ( (*it)->GetParticipant( i ) )->GetPosition() );
486 G4double nucleus_e = std::sqrt ( std::pow ( (*it)->GetNuclearMass()/
GeV , 2 ) + std::pow ( (*it)->Get4Momentum().vect().mag() , 2 ) );
487 G4LorentzVector nucleus_p4CM ( (*it)->Get4Momentum().vect() , nucleus_e );
491 G4int ia = (*it)->GetMassNumber();
492 G4int iz = (*it)->GetAtomicNumber();
501 for ( G4ReactionProductVector::iterator itt
502 = rv->begin() ; itt != rv->end() ; itt++ )
510 G4LorentzVector p4_CM = CLHEP::boostOf( p4 , -nucleus_p4CM.findBoostToCM() );
525 randomized_direction = randomized_direction.unit();
531 G4LorentzVector p4_a1_CM = CLHEP::boostOf ( p4_a1_Be8 , -nucleus_p4CM.findBoostToCM() );
532 G4LorentzVector p4_a1_LAB = CLHEP::boostOf ( p4_a1_CM , boostBackToLAB );
537 G4LorentzVector p4_a2_CM = CLHEP::boostOf ( p4_a2_Be8 , -nucleus_p4CM.findBoostToCM() );
538 G4LorentzVector p4_a2_LAB = CLHEP::boostOf ( p4_a2_CM , boostBackToLAB );
574 if ( notBreak ==
true )
585 for ( G4ReactionProductVector::iterator itt
586 = rv->begin() ; itt != rv->end() ; itt++ )
618 for ( std::vector< G4QMDNucleus* >::iterator it
619 = nucleuses.begin() ; it != nucleuses.end() ; it++ )
645 G4double stot = std::sqrt ( etot*etot - ptot*ptot );
647 G4double pstt = std::sqrt ( ( stot*stot - ( mass_proj + mass_targ ) * ( mass_proj + mass_targ )
648 ) * ( stot*stot - ( mass_proj - mass_targ ) * ( mass_proj - mass_targ ) ) )
652 G4double eccm = stot - ( mass_proj + mass_targ );
675 G4double rmax = std::sqrt( rmax0*rmax0 + b*b );
678 G4double pcca = 1.0 - double ( zp * zt ) * ccoul / eccm / rmax - ( b / rmax )*( b / rmax );
688 G4double aas = 2.0 * eccm * b / double ( zp * zt ) / ccoul;
689 bbs = 1.0 / std::sqrt ( 1.0 + aas*aas );
690 aas1 = ( 1.0 + aas * b / rmax ) * bbs;
697 if ( 1.0 - aas1*aas1 <= 0 || 1.0 - bbs*bbs <= 0.0 )
704 G4double aat1 = aas1 / std::sqrt ( 1.0 - aas1*aas1 );
705 G4double aat2 = bbs / std::sqrt ( 1.0 - bbs*bbs );
707 thet1 = std::atan ( aat1 );
708 thet2 = std::atan ( aat2 );
712 cost = std::cos( theta );
713 sint = std::sin( theta );
716 G4double rzpr = -rmax * cost * ( mass_targ ) / ( mass_proj + mass_targ );
717 G4double rzta = rmax * cost * ( mass_proj ) / ( mass_proj + mass_targ );
724 G4double pzpc = pzcc * ( cost * pccf + sint * b / rmax );
725 G4double pxpr = pzcc * ( -sint * pccf + cost * b / rmax );
730 G4double epc = std::sqrt ( pzpc*pzpc + pxpr*pxpr + mass_proj*mass_proj );
731 G4double etc = std::sqrt ( pztc*pztc + pxta*pxta + mass_targ*mass_targ );
739 G4double gammacm = boostToCM.gamma();
742 pzpr = pzpc + betacm * gammacm * ( gammacm / ( 1. + gammacm ) * pzpc * betacm + epc );
743 pzta = pztc + betacm * gammacm * ( gammacm / ( 1. + gammacm ) * pztc * betacm + etc );
744 epr = gammacm * ( epc + betacm * pzpc );
745 eta = gammacm * ( etc + betacm * pztc );
750 G4double gammpr = epr / ( mass_proj );
751 G4double gammta = eta / ( mass_targ );
753 pzta = pzta / double ( at );
754 pxta = pxta / double ( at );
756 pzpr = pzpr / double ( ap );
757 pxpr = pxpr / double ( ap );
G4VCrossSectionDataSet * GetCrossSectionDataSet(const G4String &name, G4bool warning=true)
G4double coulomb_collision_rx_proj
G4ThreeVector GetPosition()
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
G4double coulomb_collision_px_proj
G4double coulomb_collision_rx_targ
G4double coulomb_collision_pz_targ
G4Evaporation * evaporation
void SetParticipant(G4QMDParticipant *particle)
CLHEP::Hep3Vector G4ThreeVector
G4PiNuclearCrossSection * piNucXS
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
G4ParticleDefinition * GetDefinition() const
const G4ParticleDefinition * GetDefinition()
void SetNucleus(G4QMDNucleus *aSystem)
static const char * Default_Name()
G4ReactionProductVector * BreakItUp(const G4Fragment &theInitialState)
const G4String & GetParticleName() const
G4int GetAtomicNumber() const
void SetStatusChange(G4HadFinalStateStatus aS)
G4ThreeVector GetMomentum()
std::vector< G4ReactionProduct * > G4ReactionProductVector
void SetMeanField(G4QMDMeanField *meanfield)
const G4ParticleDefinition * GetDefinition() const
G4double GetIonMass(G4int Z, G4int A, G4int L=0, G4int lvl=0) const
G4double GetKineticEnergy() const
G4double coulomb_collision_px_targ
static G4CrossSectionDataSetRegistry * Instance()
const G4String & GetParticleType() const
void CalKinematicsOfBinaryCollisions(G4double)
G4int GetTotalNumberOfParticipant()
G4double GetTotalPotential()
const G4LorentzVector & Get4Momentum() const
G4int GetAtomicMass() const
G4LorentzVector Get4Momentum() const
std::vector< G4QMDNucleus * > DoClusterJudgment()
G4QMDParticipant * GetParticipant(G4int i)
static G4IonTable * GetIonTable()
G4QMDMeanField * meanField
G4double coulomb_collision_rz_targ
void calcOffSetOfCollision(G4double, const G4ParticleDefinition *, const G4ParticleDefinition *, G4double, G4double, G4double, G4ThreeVector)
G4IonsShenCrossSection * shenXS
G4double GetPDGMass() const
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
G4double coulomb_collision_gamma_targ
G4LorentzVector Get4Momentum()
void DoPropagation(G4double)
void SetSystem(G4QMDSystem *aSystem)
void SetTotalPotential(G4double x)
G4double coulomb_collision_rz_proj
G4ExcitationHandler * excitationHandler
void SetEvaporation(G4VEvaporation *ptr)
G4HadFinalState theParticleChange
G4QMDCollision * collision
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
G4double coulomb_collision_gamma_proj
void CalEnergyAndAngularMomentumInCM()
static const double eplus
G4double GetPDGCharge() const
static const double fermi
G4ThreeVector GetMomentum() const
G4double coulomb_collision_pz_proj
G4double GetTotalMomentum() const
G4double elastic(Particle const *const p1, Particle const *const p2)
CLHEP::HepLorentzVector G4LorentzVector