80 if (targetMass < 500*MeV) targetMass =
89 theNeutron.
Lorentz(theNeutron, theTarget);
100 thePhotons = theMF6FinalState.
Sample( eKinetic );
102 thePhotons = theFinalStatePhotons.
GetPhotons(eKinetic);
104 if ( thePhotons == NULL ) {
105 throw G4HadronicException(__FILE__, __LINE__,
"Final state data for photon is not properly allocated");
119 for(
auto it=products->cbegin(); it!=products->cend(); ++it)
123 if ( (*it)->GetParticleDefinition() != 0 )
131 if ( (*it)->GetMomentum().mag() > 10*MeV)
134 if ( (*it)->GetExcitationEnergy() > 1.0e-2*eV) {
135 G4double ex = (*it)->GetExcitationEnergy();
138 aPhoton->
SetMomentum( (*it)->GetMomentum().vect().unit() * ex );
140 thePhotons->push_back(aPhoton);
143 theOne->
SetMomentum( (*it)->GetMomentum().vect() * ( (*it)->GetMomentum().t() - (*it)->GetExcitationEnergy() ) / (*it)->GetMomentum().t() ) ;
144 thePhotons->push_back(theOne);
152 nPhotons=(
G4int)thePhotons->size();
164 G4double theta = std::acos(costheta);
167 G4ThreeVector direction(sinth*std::cos(phi), sinth*std::sin(phi), costheta);
169 thePhotons->push_back(theOne);
175 if ( nPhotons == 1 && thePhotons->operator[](0)->GetDefinition()->GetBaryonNumber() == 0 )
182 thePhotons->operator[](0)->SetMomentum( Q*direction );
188 for(i=0; i<nPhotons; i++)
190 thePhotons->operator[](i)->Lorentz(*(thePhotons->operator[](i)), -1*theTarget);
195 if ( nPhotons == 1 && thePhotons->operator[](0)->GetDefinition()->GetBaryonNumber() == 0 )
205 -thePhotons->operator[](0)->GetMomentum();
222 for(i=0; i<nPhotons; i++)
226 theOne->
SetDefinition(thePhotons->operator[](i)->GetDefinition());
227 theOne->
SetMomentum(thePhotons->operator[](i)->GetMomentum());
229 delete thePhotons->operator[](i);
237 for ( std::size_t j = 0 ; j !=
theResult.
Get()->GetNumberOfSecondaries() ; j++ )
242 if ( residual ==
false )
246 for ( std::size_t j = 0 ; j !=
theResult.
Get()->GetNumberOfSecondaries() ; j++ )
258 if ( nPhotons - nNonZero > 0 )
261 std::vector<G4double> vRand;
262 vRand.push_back( 0.0 );
263 for (
G4int j = 0 ; j != nPhotons - nNonZero - 1 ; j++ )
267 vRand.push_back( 1.0 );
268 std::sort( vRand.begin(), vRand.end() );
270 std::vector<G4double> vEPhoton;
271 for (
G4int j = 0 ; j < (
G4int)vRand.size() - 1 ; j++ )
273 vEPhoton.push_back( deltaE * ( vRand[j+1] - vRand[j] ) );
275 std::sort( vEPhoton.begin(), vEPhoton.end() );
277 for (
G4int j = 0 ; j < nPhotons - nNonZero - 1 ; j++ )
282 G4double theta = std::acos(costheta);
286 G4ThreeVector tempVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*costheta);
326 std::stringstream ss;
327 ss << static_cast<G4int>(
Z);
331 ss << static_cast<G4int>(
A);
346 G4String filenameMF6 = dirName+
"/FSMF6/"+sZ+
"_"+sA+sM+
"_"+element_name;
349 std::istringstream theData(std::ios::in);
354 if ( theData.good() ==
true ) {
356 theMF6FinalState.
Init(theData);
371 if(!dbool || (
Z<2.5 && ( std::abs(
theBaseZ -
Z)>0.0001 || std::abs(
theBaseA -
A)>0.0001)))
std::vector< G4Fragment * > G4FragmentVector
CLHEP::HepLorentzVector G4LorentzVector
std::vector< G4ReactionProduct * > G4ReactionProductVector
void Put(const value_type &val) const
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4ParticleDefinition * GetDefinition() const
G4LorentzVector Get4Momentum() const
void SetMomentum(const G4ThreeVector &momentum)
void SetStatusChange(G4HadFinalStateStatus aS)
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
G4HadSecondary * GetSecondary(size_t i)
const G4Material * GetMaterial() const
const G4ParticleDefinition * GetDefinition() const
G4double GetKineticEnergy() const
const G4LorentzVector & Get4Momentum() const
G4DynamicParticle * GetParticle()
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
static G4IonTable * GetIonTable()
G4double GetIonMass(G4int Z, G4int A, G4int nL=0, G4int lvl=0) const
G4double GetTemperature() const
static G4Neutron * Neutron()
static G4double GetNuclearMass(const G4double A, const G4double Z)
G4ReactionProduct GetBiasedThermalNucleus(G4double aMass, G4ThreeVector aVelocity, G4double temp=-1) const
G4double GetPDGMass() const
void Init(G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *)
G4HadFinalState * ApplyYourself(const G4HadProjectile &theTrack)
void Init(std::istream &aDataFile)
void SetProjectileRP(G4ReactionProduct &aIncidentPart)
void SetTarget(G4ReactionProduct &aTarget)
G4ReactionProductVector * Sample(G4double anEnergy)
void SetAZMs(G4double anA, G4double aZ, G4int aM, G4ParticleHPDataUsed used)
G4Cache< G4HadFinalState * > theResult
static G4ParticleHPManager * GetInstance()
G4bool GetDoNotAdjustFinalState()
void GetDataStream(G4String, std::istringstream &iss)
G4ParticleHPDataUsed GetName(G4int A, G4int Z, G4String base, G4String rest, G4bool &active)
void InitEnergies(std::istream &aDataFile)
G4ReactionProductVector * GetPhotons(G4double anEnergy)
void InitAngular(std::istream &aDataFile)
G4bool InitMean(std::istream &aDataFile)
void SetICM(G4bool) override
G4FragmentVector * BreakItUp(const G4Fragment &theNucleus)
static G4int GetModelID(const G4int modelIndex)
void SetMomentum(const G4double x, const G4double y, const G4double z)
G4double GetKineticEnergy() const
G4double GetTotalEnergy() const
G4ThreeVector GetMomentum() const
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)
void SetDefinitionAndUpdateE(const G4ParticleDefinition *aParticleDefinition)
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
void SetKineticEnergy(const G4double en)