G4ProjectileDiffractiveChannel Class Reference

#include <G4ProjectileDiffractiveChannel.hh>

Public Member Functions
	G4ProjectileDiffractiveChannel ()
	~G4ProjectileDiffractiveChannel ()
G4double	GetFraction (G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary)
G4KineticTrackVector *	Scatter (G4Nucleus &theNucleus, const G4DynamicParticle &thePrimary)

Detailed Description

Definition at line 47 of file G4ProjectileDiffractiveChannel.hh.

Constructor & Destructor Documentation

G4ProjectileDiffractiveChannel()

G4ProjectileDiffractiveChannel::G4ProjectileDiffractiveChannel

(

)

Definition at line 49 of file G4ProjectileDiffractiveChannel.cc.

{
    theQDiffraction=G4QDiffractionRatio::GetPointer();
}

~G4ProjectileDiffractiveChannel()

G4ProjectileDiffractiveChannel::~G4ProjectileDiffractiveChannel

(

)

Definition at line 54 of file G4ProjectileDiffractiveChannel.cc.

{}

Member Function Documentation

GetFraction()

G4double G4ProjectileDiffractiveChannel::GetFraction	(	G4Nucleus &	theNucleus,
		const G4DynamicParticle &	thePrimary
	)

Definition at line 57 of file G4ProjectileDiffractiveChannel.cc.

{
    G4double ratio;
    ratio=theQDiffraction->GetRatio(thePrimary.GetTotalMomentum(),
            thePrimary.GetDefinition()->GetPDGEncoding(),
            theNucleus.GetZ_asInt(), theNucleus.GetN_asInt() );
       #ifdef debug_getFraction         
          G4cout << "G4ProjectileDiffractiveChannel::ratio " << ratio << G4endl;
       #endif
           
    return ratio;               
}

Referenced by G4TheoFSGenerator::ApplyYourself().

Scatter()

G4KineticTrackVector * G4ProjectileDiffractiveChannel::Scatter	(	G4Nucleus &	theNucleus,
		const G4DynamicParticle &	thePrimary
	)

Definition at line 71 of file G4ProjectileDiffractiveChannel.cc.

{
    
    
    G4int A=theNucleus.GetA_asInt();
    G4int Z=theNucleus.GetZ_asInt();
 
    G4ParticleDefinition * pDef=thePrimary.GetDefinition();
       #ifdef debug_scatter
         G4cout << " ProjectileDiffractive: A, Z,  proj-pdg" <<" "<< A <<" "<<Z 
            << " "<<" " << pDef->GetParticleName()<< G4endl;
       #endif
    
    G4QHadronVector * result(0);
    G4KineticTrackVector * ktv(0);
    G4bool tryAgain;
    
    do {
       tryAgain = false;
       result=theQDiffraction->ProjFragment(pDef->GetPDGEncoding(),
                    thePrimary.Get4Momentum(),Z, A-Z);
                
       if ( result && result->size() > 0 )
       {
           ktv=new G4KineticTrackVector();
           std::vector<G4QHadron *>::iterator i;
         #ifdef debug_scatter
            G4int count(0);
         #endif
           for (i=result->begin(); i!=result->end(); i++)
           {
          G4int pdgCode=(*i)->GetPDGCode();
          G4ParticleDefinition * secDef(0);
          if ( pdgCode < 80000000) {     // check if ion; should be 90Million, but 89Million is possible
              secDef= G4ParticleTable::GetParticleTable()->FindParticle(pdgCode);
          }else {
              G4int qN = pdgCode % 1000;
              G4int qZ = (pdgCode/1000) %1000;
              if ( qZ < 500 && qN < 500){   // protect for delta being coded as/in nucleus 
             secDef = G4ParticleTable::GetParticleTable()->GetIon(qZ,qN+qZ, 0);
             if ( ! secDef ) 
             {  // exceptions to the rule!
                    if ( pdgCode == 90000001 ) secDef= G4Neutron::Neutron();
                if ( pdgCode == 91000000 ) secDef= G4Lambda::Lambda();
             }
              }   
          }
          
              if  ( secDef ){
 
              G4ThreeVector pos=(*i)->GetPosition();  //GetPosition returns const &
              G4LorentzVector mom=(*i)->Get4Momentum();
 
              G4KineticTrack * sPrim=new G4KineticTrack(secDef,
                      (*i)->GetFormationTime(), pos, mom);
              ktv->push_back(sPrim);
 
            #ifdef debug_scatter
              G4cout << "G4ProjectileDiffractive sec # "  << ++count << ", "
                    << "ChipsPDGCode=" << pdgCode << ", "
                    << secDef->GetParticleName() 
                    << ", time(ns) " << (*i)->GetFormationTime()/ns
                    << ", pos " << pos
                    << ", 4mom " << (*i)->Get4Momentum()
                    << G4endl;
            #endif
               } else {
               #ifdef debug_scatter
                 G4cout << "G4ProjectileDiffractiveChannel: Particle with PDG code "<< pdgCode <<" does not converted!!!"<<G4endl;
               #endif 
               tryAgain=true;
           }
           }   
           std::for_each(result->begin(), result->end(), DeleteQHadron());
               delete result;
           result=0;
       }
        
       if ( result ) delete result;
       if ( tryAgain ) {     // QDiffraction returned a "difficult" particle in nucleus
           std::for_each(ktv->begin(), ktv->end(), DeleteKineticTrack());
           delete ktv;
           ktv=0;
       }
    } while ( ! ktv);          
 
    return ktv;                
}

Referenced by G4TheoFSGenerator::ApplyYourself().

---

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

• G4ProjectileDiffractiveChannel.hh
• G4ProjectileDiffractiveChannel.cc