54 : myA(0), myZ(0), theNucleons(250), currentNucleon(-1), theDensity(0),
55 nucleondistance(0.8*
fermi),excitationEnergy(0.),
56 places(250), momentum(250), fermiM(250), testSums(250)
66 #if defined(NON_INTEGER_A_Z)
110 for (
G4int aNucleon=0; aNucleon <
myA; aNucleon++)
181 for (
int i=0; i<
myA; i++)
183 if (
theNucleons[i].GetPosition().mag2() > maxradius2 )
222 factor * (theBeta*
theNucleons[i].GetPosition()) * theBeta;
230 if (theBoost.e() !=0 ) {
269 G4int protons=0,nucleons=0;
271 while (nucleons <
myA )
278 else if ( (nucleons-protons) < (
myA-
myZ) )
282 else G4cout <<
"G4Fancy3DNucleus::ChooseNucleons not efficient" <<
G4endl;
302 while ( (i <
myA) && (--interationsLeft>0))
309 G4RandFlat::shootArray(jr,prand);
314 aPos.set((2*arand[jx]-1.), (2*arand[jy]-1.), (2*arand[--jr]-1.));
315 }
while (aPos.mag2() > 1. );
321 std::vector<G4ThreeVector>::iterator iplace;
322 for( iplace=
places.begin(); iplace!=
places.end() && freeplace;++iplace)
324 delta = *iplace - aPos;
325 freeplace= delta.mag2() > nd2;
350 if (interationsLeft<=0) {
352 "Problem to place nucleons");
366 for (
G4int ntry=0; ntry<1 ; ntry ++ )
368 for (i=0; i <
myA; i++ )
377 if ( eMax >
theNucleons[i].GetDefinition()->GetPDGMass() )
380 fermiM[i] = std::sqrt(pmax2);
381 while ( mom.mag2() > pmax2 )
387 G4cerr <<
"G4Fancy3DNucleus: difficulty finding proton momentum" <<
G4endl;
405 for ( i=0; i<
myA ; i++ )
407 energy =
theNucleons[i].GetParticleType()->GetPDGMass()
427 if ( sum.mag() <= PFermi )
439 for (
G4int aNucleon=0; aNucleon < myA-1; aNucleon++) {
440 delta = 2.*((
momentum[aNucleon]*testDir)*testDir);
442 testSums[aNucleon].Fill(delta, delta.mag(), aNucleon);
449 while ( (sum-
testSums[--index].Vector).mag()>PFermi && index>0)
452 if ( sum.mag() > (sum-
testSums[index].Vector).mag() ) {
454 sum-=testSums[index].Vector;
458 if ( (sum-
testSums[index].Vector).mag() <= PFermi )
462 for (
G4int aNucleon=0; aNucleon<=index; aNucleon++)
467 && std::abs(
momentum[myA-1].mag() - pTry ) < pBest )
469 pBest=std::abs(
momentum[myA-1].mag() - pTry );
475 G4String text =
"G4Fancy3DNucleus.cc: Logic error in ReduceSum()";
479 momentum[myA-1]=testSums[best].Vector-sum;
489 if (
fermiM[++swapit] > PFermi )
break;
491 if (swapit == myA-1 )
return false;
static G4Pow * GetInstance()
std::vector< G4ThreeVector > momentum
void ChooseFermiMomenta()
CLHEP::Hep3Vector G4ThreeVector
G4double CoulombBarrier()
void DoTranslation(const G4ThreeVector &theShift)
std::vector< G4double > fermiM
void DoLorentzBoost(const G4LorentzVector &theBoost)
virtual const G4ThreeVector & GetPosition() const
G4double excitationEnergy
std::vector< G4Fancy3DNucleusHelper > testSums
std::vector< G4ThreeVector > places
virtual G4double GetRadius(const G4double maxRelativeDenisty) const =0
G4GLOB_DLL std::ostream G4cout
G4double Z13(G4int Z) const
std::vector< G4Nucleon > theNucleons
bool G4Fancy3DNucleusHelperForSortInZ(const G4Nucleon &nuc1, const G4Nucleon &nuc2)
void Init(G4int theA, G4int theZ)
virtual G4double GetRelativeDensity(const G4ThreeVector &aPosition) const =0
G4double GetFermiMomentum(G4double density)
static G4Proton * Proton()
G4double GetNuclearRadius()
static G4Neutron * Neutron()
const std::vector< G4Nucleon > & GetNucleons()
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
void Init(G4int anA, G4int aZ)
static const G4double factor
G4double GetPDGMass() const
G4double GetDensity(const G4ThreeVector &aPosition) const
static G4double GetBindingEnergy(const G4int A, const G4int Z)
G4double energy(const ThreeVector &p, const G4double m)
void DoLorentzContraction(const G4LorentzVector &theBoost)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4ThreeVector GetMomentum(G4double density, G4double maxMomentum=-1.)
G4double GetOuterRadius()
const G4VNuclearDensity * GetNuclearDensity() const
const G4double nucleondistance
G4Nucleon * GetNextNucleon()
static const double fermi
G4VNuclearDensity * theDensity
G4GLOB_DLL std::ostream G4cerr
CLHEP::HepLorentzVector G4LorentzVector