G4ElementaryParticleCollider Class Reference

#include <G4ElementaryParticleCollider.hh>

Inheritance diagram for G4ElementaryParticleCollider:

Public Member Functions
	G4ElementaryParticleCollider ()
virtual	~G4ElementaryParticleCollider ()
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)

---

Detailed Description

Definition at line 60 of file G4ElementaryParticleCollider.hh.

---

Constructor & Destructor Documentation

G4ElementaryParticleCollider::G4ElementaryParticleCollider

(

)

Definition at line 110 of file G4ElementaryParticleCollider.cc.

  : G4CascadeColliderBase("G4ElementaryParticleCollider") {}

virtual G4ElementaryParticleCollider::~G4ElementaryParticleCollider

(

)

[inline, virtual]

Definition at line 63 of file G4ElementaryParticleCollider.hh.

{};

---

Member Function Documentation

void G4ElementaryParticleCollider::collide	(	G4InuclParticle *	bullet,
		G4InuclParticle *	target,
		G4CollisionOutput &	output
	)			[virtual]

Implements G4VCascadeCollider.

Definition at line 115 of file G4ElementaryParticleCollider.cc.

References G4CollisionOutput::addOutgoingParticles(), G4LorentzConvertor::backToTheLab(), G4cerr, G4cout, G4endl, G4InuclParticle::getDefinition(), G4LorentzConvertor::getKinEnergyInTheTRS(), G4InuclParticle::getMomModule(), G4ParticleDefinition::GetParticleName(), G4LorentzConvertor::getSCMMomentum(), G4CascadeChannelTables::GetTable(), G4LorentzConvertor::getTotalSCMEnergy(), G4InteractionCase::hadrons(), G4CascadeColliderBase::interCase, G4InuclElementaryParticle::isPhoton(), G4InuclElementaryParticle::nucleon(), G4InuclElementaryParticle::pion(), G4InuclElementaryParticle::quasi_deutron(), G4InteractionCase::set(), G4LorentzConvertor::setBullet(), G4LorentzConvertor::setTarget(), G4LorentzConvertor::setVerbose(), G4LorentzConvertor::toTheCenterOfMass(), G4CascadeColliderBase::useEPCollider(), G4CascadeColliderBase::validateOutput(), and G4VCascadeCollider::verboseLevel.

Referenced by G4InuclCollider::collide(), and G4NucleiModel::generateParticleFate().

{
  if (verboseLevel > 1)
    G4cout << " >>> G4ElementaryParticleCollider::collide" << G4endl;

  if (!useEPCollider(bullet,target)) {          // Sanity check
    G4cerr << " ElementaryParticleCollider -> can collide only particle with particle " 
           << G4endl;
    return;
  }

#ifdef G4CASCADE_DEBUG_SAMPLER
  static G4bool doPrintTables = true;   // Once and only once per job
  if (doPrintTables) {
    printFinalStateTables();            // For diagnostic reporting
    doPrintTables = false;
  }
#endif

  interCase.set(bullet, target);        // To identify kind of collision

  if (verboseLevel > 1) G4cout << *bullet << G4endl << *target << G4endl;

  G4InuclElementaryParticle* particle1 =
    dynamic_cast<G4InuclElementaryParticle*>(bullet);
  G4InuclElementaryParticle* particle2 =        
    dynamic_cast<G4InuclElementaryParticle*>(target);

  if (!particle1 || !particle2) {       // Redundant with useEPCollider()
    G4cerr << " ElementaryParticleCollider -> can only collide hadrons"
           << G4endl;
    return;
  }

  // Check for available interaction, or pion+dibaryon special case
  if (!G4CascadeChannelTables::GetTable(interCase.hadrons()) &&
      !particle1->quasi_deutron() && !particle2->quasi_deutron()) {
    G4cerr << " ElementaryParticleCollider -> cannot collide " 
           << particle1->getDefinition()->GetParticleName() << " with "
           << particle2->getDefinition()->GetParticleName() << G4endl;
    return;
  }
  // Generate nucleon or pion collision with nucleon
  // or pion with quasi-deuteron

  if (particle1->nucleon() || particle2->nucleon()) { // ok
    G4LorentzConvertor convertToSCM;
    if(particle2->nucleon()) {
      convertToSCM.setBullet(particle1);
      convertToSCM.setTarget(particle2);
    } else {
      convertToSCM.setBullet(particle2);
      convertToSCM.setTarget(particle1);
    };

    convertToSCM.setVerbose(verboseLevel);

    convertToSCM.toTheCenterOfMass();
    G4double ekin = convertToSCM.getKinEnergyInTheTRS();
    G4double etot_scm = convertToSCM.getTotalSCMEnergy();
    G4double pscm = convertToSCM.getSCMMomentum();

    generateSCMfinalState(ekin, etot_scm, pscm, particle1, particle2,
                          &convertToSCM);

    if (particles.empty()) {    // No final state possible, pass bullet through
      if (verboseLevel) {
        G4cerr << " ElementaryParticleCollider -> failed to collide " 
               << particle1->getMomModule() << " GeV/c " 
               << particle1->getDefinition()->GetParticleName() << " with "
               << particle2->getDefinition()->GetParticleName() << G4endl;
      }
    } else {                     // convert back to Lab
      G4LorentzVector mom;              // Buffer to avoid memory churn
      particleIterator ipart;
      for(ipart = particles.begin(); ipart != particles.end(); ipart++) {       
        mom = convertToSCM.backToTheLab(ipart->getMomentum());
        ipart->setMomentum(mom); 
      };

      // Check conservation in multibody final state
      if (verboseLevel && !validateOutput(bullet, target, particles)) {
        G4cout << " incoming particles: \n" << *particle1 << G4endl
               << *particle2 << G4endl
               << " outgoing particles: " << G4endl;
        for(ipart = particles.begin(); ipart != particles.end(); ipart++)
          G4cout << *ipart << G4endl;

        G4cout << " <<< Non-conservation in G4ElementaryParticleCollider"
               << G4endl;
      }

      std::sort(particles.begin(), particles.end(), G4ParticleLargerEkin());
      output.addOutgoingParticles(particles);
    }
  } else {      // neither particle is nucleon: pion on quasideuteron
    if (particle1->quasi_deutron() || particle2->quasi_deutron()) {
      if (particle1->pion() || particle2->pion() ||
          particle1->isPhoton() || particle2->isPhoton()) {
        G4LorentzConvertor convertToSCM;
        if(particle2->quasi_deutron()) {        // Quasideuteron is target
          convertToSCM.setBullet(particle1);
          convertToSCM.setTarget(particle2);
        } else {
          convertToSCM.setBullet(particle2);
          convertToSCM.setTarget(particle1);
        }; 
        convertToSCM.toTheCenterOfMass(); 
        G4double etot_scm = convertToSCM.getTotalSCMEnergy();
        
        generateSCMpionAbsorption(etot_scm, particle1, particle2);

        if (particles.empty()) {        // Failed to generate final state
          if (verboseLevel) {
            G4cerr << " ElementaryParticleCollider -> failed to collide " 
                   << particle1->getMomModule() << " GeV/c " 
                   << particle1->getDefinition()->GetParticleName() << " with "
                   << particle2->getDefinition()->GetParticleName() << G4endl;
          }
        } else {                        // convert back to Lab
          G4LorentzVector mom;  // Buffer to avoid memory churn
          particleIterator ipart;
          for(ipart = particles.begin(); ipart != particles.end(); ipart++) {
            mom = convertToSCM.backToTheLab(ipart->getMomentum());
            ipart->setMomentum(mom); 
          };

          validateOutput(bullet, target, particles);    // Check conservation

          std::sort(particles.begin(), particles.end(), G4ParticleLargerEkin());
          output.addOutgoingParticles(particles);
        };
      } else {
        G4cerr << " ElementaryParticleCollider -> can only collide pions with dibaryons " 
               << G4endl;
      };
    } else {
      G4cerr << " ElementaryParticleCollider -> can only collide something with nucleon or dibaryon " 
             << G4endl;
    };
  };  
}

---

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

• G4ElementaryParticleCollider.hh
• G4ElementaryParticleCollider.cc

---