78 theNeutron.
Lorentz(theNeutron, theTarget);
93 if ( thePhotons == NULL ) {
94 throw G4HadronicException(__FILE__, __LINE__,
"Final state data for photon is not properly allocated");
107 G4FragmentVector::iterator
it;
109 for(it=products->begin(); it!=products->end(); it++)
113 if ( (*it)->GetParticleDefinition() != 0 )
121 if ( (*it)->GetMomentum().mag() > 10*
MeV)
124 if ( (*it)->GetExcitationEnergy() > 1.0e-2*
eV) {
125 G4double ex = (*it)->GetExcitationEnergy();
128 aPhoton->
SetMomentum( (*it)->GetMomentum().vect().unit() * ex );
130 thePhotons->push_back(aPhoton);
133 theOne->
SetMomentum( (*it)->GetMomentum().vect() * ( (*it)->GetMomentum().t() - (*it)->GetExcitationEnergy() ) / (*it)->GetMomentum().t() ) ;
134 thePhotons->push_back(theOne);
142 nPhotons=thePhotons->size();
157 G4ThreeVector direction(sinth*std::cos(phi), sinth*std::sin(phi), costheta);
159 thePhotons->push_back(theOne);
165 if ( nPhotons == 1 && thePhotons->operator[](0)->GetDefinition()->GetBaryonNumber() == 0 )
172 thePhotons->operator[](0)->SetMomentum( Q*direction );
178 for(i=0; i<nPhotons; i++)
180 thePhotons->operator[](i)->Lorentz(*(thePhotons->operator[](i)), -1*theTarget);
185 if ( nPhotons == 1 && thePhotons->operator[](0)->GetDefinition()->GetBaryonNumber() == 0 )
195 -thePhotons->operator[](0)->GetMomentum();
212 for(i=0; i<nPhotons; i++)
216 theOne->
SetDefinition(thePhotons->operator[](i)->GetDefinition());
217 theOne->
SetMomentum(thePhotons->operator[](i)->GetMomentum());
219 delete thePhotons->operator[](i);
232 if ( residual ==
false )
248 if ( nPhotons - nNonZero > 0 )
251 std::vector<G4double> vRand;
252 vRand.push_back( 0.0 );
253 for (
G4int j = 0 ; j != nPhotons - nNonZero - 1 ; j++ )
257 vRand.push_back( 1.0 );
258 std::sort( vRand.begin(), vRand.end() );
260 std::vector<G4double> vEPhoton;
261 for (
G4int j = 0 ; j < (
G4int)vRand.size() - 1 ; j++ )
263 vEPhoton.push_back( deltaE * ( vRand[j+1] - vRand[j] ) );
265 std::sort( vEPhoton.begin(), vEPhoton.end() );
267 for (
G4int j = 0 ; j < nPhotons - nNonZero - 1 ; j++ )
276 G4ThreeVector tempVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*costheta);
281 theOne->SetMomentum( tempVector );
316 std::stringstream ss;
317 ss << static_cast<G4int>(
Z);
321 ss << static_cast<G4int>(
A);
336 G4String filenameMF6 = dirName+
"/FSMF6/"+sZ+
"_"+sA+sM+
"_"+element_name;
344 if ( theData.good() ==
true ) {