138 using namespace CLHEP;
140 const G4double G4RadioactiveDecay::levelTolerance = 2.0*
keV;
142 #ifdef G4MULTITHREADED
148 forceDecayDirection(0.,0.,0.), forceDecayHalfAngle(0.*deg), verboseLevel(0)
152 G4cout <<
"G4RadioactiveDecay constructor: processName = " << processName
174 theUserRadioactiveDataFiles.clear();
177 #ifdef G4MULTITHREADED
195 theRadioactivityTables.push_back(rTable);
205 isAllVolumesMode =
true;
212 delete theRadioactiveDecaymessenger;
213 delete theIsotopeTable;
220 if (((
const G4Ions*)(&aParticle))->GetExcitationEnergy() > 0.) {
return true;}
229 G4int A = ((
const G4Ions*) (&aParticle))->GetAtomicMass();
230 G4int Z = ((
const G4Ions*) (&aParticle))->GetAtomicNumber();
231 if (A > theNucleusLimits.
GetAMax() || A < theNucleusLimits.
GetAMin())
233 else if (Z > theNucleusLimits.
GetZMax() || Z < theNucleusLimits.
GetZMin())
241 DecayTableMap::iterator table_ptr = dkmap->find(key);
244 if (table_ptr == dkmap->end() ) {
246 (*dkmap)[key] = theDecayTable;
248 theDecayTable = table_ptr->second;
251 return theDecayTable;
260 for (
size_t i = 0; i < theLogicalVolumes->size(); i++) {
261 volume=(*theLogicalVolumes)[i];
262 if (volume->
GetName() == aVolume) {
263 ValidVolumes.push_back(aVolume);
264 std::sort(ValidVolumes.begin(), ValidVolumes.end());
270 }
else if(i == theLogicalVolumes->size()) {
271 G4cerr <<
"SelectAVolume: "<< aVolume
272 <<
" is not a valid logical volume name" <<
G4endl;
283 for (
size_t i = 0; i < theLogicalVolumes->size(); i++){
284 volume=(*theLogicalVolumes)[i];
285 if (volume->
GetName() == aVolume) {
286 std::vector<G4String>::iterator location;
287 location = std::find(ValidVolumes.begin(),ValidVolumes.end(),aVolume);
288 if (location != ValidVolumes.end()) {
289 ValidVolumes.erase(location);
290 std::sort(ValidVolumes.begin(), ValidVolumes.end());
291 isAllVolumesMode =
false;
293 G4cerr <<
" DeselectVolume:" << aVolume <<
" is not in the list "
298 G4cout <<
" DeselectVolume: " << aVolume <<
" is removed from list "
301 }
else if (i == theLogicalVolumes->size()) {
302 G4cerr <<
" DeselectVolume:" << aVolume
303 <<
"is not a valid logical volume name" <<
G4endl;
314 ValidVolumes.clear();
319 for (
size_t i = 0; i < theLogicalVolumes->size(); i++){
320 volume = (*theLogicalVolumes)[i];
321 ValidVolumes.push_back(volume->
GetName());
327 std::sort(ValidVolumes.begin(), ValidVolumes.end());
329 isAllVolumesMode=
true;
335 ValidVolumes.clear();
336 isAllVolumesMode=
false;
349 for (
size_t i = 0; i < theDecayRateTableVector.size(); i++) {
350 if (theDecayRateTableVector[i].GetIonName() == aParticleName)
return true;
363 for (
size_t i = 0; i < theDecayRateTableVector.size(); i++) {
364 if (theDecayRateTableVector[i].GetIonName() == aParticleName) {
365 theDecayRateVector = theDecayRateTableVector[i].GetItsRates();
370 G4cout <<
"The DecayRate Table for " << aParticleName <<
" is selected."
380 long double taotime = 0.L;
382 if ( t > SBin[NSourceBin]) {
386 while (t > SBin[nbin]) nbin++;
389 long double ltao = tao;
392 for (
G4int i = 0; i < nbin; i++) {
393 taotime += (
long double)SProfile[i] *
394 (std::exp(-(lt-(
long double)SBin[i+1])/ltao)-std::exp(-(lt-(
long double)SBin[i])/ltao));
397 taotime += (
long double)SProfile[nbin] * (1.L-std::exp(-(lt-(
long double)SBin[nbin])/ltao));
399 G4cout <<
" Tao time =: " <<taotime <<
" reset to zero!"<<
G4endl;
534 while ( DProfile[i] < rand) i++;
536 decaytime = DBin[i] + rand*(DBin[i+1]-DBin[i]);
548 while ( aDecayTime > DBin[i] ) i++;
573 G4cout <<
"G4RadioactiveDecay::GetMeanLifeTime() " <<
G4endl;
575 <<
" GeV, Mass: " << theParticle->
GetMass()/
GeV
576 <<
" GeV, Life time: " << theLife/
ns <<
" ns " <<
G4endl;
580 else if (theLife < 0.0) {meanlife =
DBL_MAX;}
581 else {meanlife = theLife;}
584 if (((
const G4Ions*)(theParticleDef))->GetExcitationEnergy() > 0. &&
585 meanlife ==
DBL_MAX) {meanlife = 0.;}
589 G4cout <<
" mean life time: " << meanlife/
s <<
" s " <<
G4endl;
611 G4cout <<
"G4RadioactiveDecay::GetMeanFreePath() " <<
G4endl;
613 <<
" GeV, Mass: " << aMass/
GeV <<
" GeV, tau: " << tau <<
" ns "
623 ed <<
"Ion has negative lifetime " << tau
624 <<
" but is not stable. Setting mean free path to DBL_MAX" <<
G4endl;
625 G4Exception(
"G4RadioactiveDecay::GetMeanFreePath()",
"HAD_RDM_011",
632 pathlength =
c_light*tau*betaGamma;
638 G4cout <<
"G4Decay::GetMeanFreePath: "
640 <<
" stops, kinetic energy = "
650 G4cout <<
"mean free path: "<< pathlength/
m <<
" m" <<
G4endl;
664 if (!isInitialised) {
665 isInitialised =
true;
679 #ifdef G4MULTITHREADED
689 G4int A = ((
const G4Ions*)(&theParentNucleus))->GetAtomicMass();
690 G4int Z = ((
const G4Ions*)(&theParentNucleus))->GetAtomicNumber();
691 G4int lvl = ((
const G4Ions*)(&theParentNucleus))->GetIsomerLevel();
694 #ifdef G4MULTITHREADED
695 G4AutoLock lk(&G4RadioactiveDecay::radioactiveDecayMutex);
698 DecayTableMap::iterator master_table_ptr = master_dkmap->find(key);
700 if (master_table_ptr != master_dkmap->end() ) {
701 return master_table_ptr->second;
712 if (!getenv(
"G4RADIOACTIVEDATA") ) {
713 G4cout <<
"Please setenv G4RADIOACTIVEDATA to point to the radioactive decay data files."
715 throw G4HadronicException(__FILE__, __LINE__,
" Please setenv G4RADIOACTIVEDATA to point to the radioactive decay data files.");
717 G4String dirName = getenv(
"G4RADIOACTIVEDATA");
719 std::ostringstream os;
720 os <<dirName <<
"/z" <<Z <<
".a" <<A <<
'\0';
724 std::ifstream DecaySchemeFile(file);
727 if (DecaySchemeFile) {
730 G4bool modeFirstRecord[7];
733 for (
G4int i = 0; i < nMode; i++) {
734 modeFirstRecord[i] =
true;
739 char inputChars[100]={
' '};
746 G4int levelCounter = 0;
753 while (!complete && !DecaySchemeFile.getline(inputChars, 100).eof()) {
754 inputLine = inputChars;
755 inputLine = inputLine.
strip(1);
756 if (inputChars[0] !=
'#' && inputLine.length() != 0) {
757 std::istringstream tmpStream(inputLine);
759 if (inputChars[0] ==
'P') {
762 tmpStream >> recordType >> a >>
b;
767 found = (levelCounter == lvl);
774 if (inputChars[0] ==
'W') {
778 G4cout <<
" Warning in G4RadioactiveDecay::LoadDecayTable " <<
G4endl;
784 tmpStream >> theDecayMode >> a >> b >> c >> betaType;
788 if (inputLine.length() < 80) betaType =
allowed;
792 switch (theDecayMode) {
798 (
const G4Ions*)& theParentNucleus, b);
799 anITChannel->
SetICM(applyICM);
800 anITChannel->
SetARM(applyARM);
802 theDecayTable->
Insert(anITChannel);
808 if (modeFirstRecord[1]) {
809 modeFirstRecord[1] =
false;
823 for (
G4int ptn = 0; ptn < npti; ptn++) {
825 e = 1. + e0*(ptn+0.5)/100.;
826 p = std::sqrt(e*e - 1.);
827 f = p*e*(e0-e+1)*(e0-e+1);
833 f *= corrections.
ShapeFactor(betaType, p, e0-e+1.);
840 b, c*
MeV, a*MeV, 0, FBeta, aRandomEnergy);
841 aBetaMinusChannel->
SetICM(applyICM);
842 aBetaMinusChannel->
SetARM(applyARM);
844 theDecayTable->
Insert(aBetaMinusChannel);
854 if (modeFirstRecord[2]) {
855 modeFirstRecord[2] =
false;
858 e0 = c*
MeV/0.510999 - 2.;
871 for (
G4int ptn = 0; ptn < npti; ptn++) {
873 e = 1. + e0*(ptn+0.5)/100.;
874 p = std::sqrt(e*e - 1.);
875 f = p*e*(e0-e+1)*(e0-e+1);
881 f *= corrections.
ShapeFactor(betaType, p, e0-e+1.);
887 &theParentNucleus, b,
888 (c-1.021998)*
MeV, a*
MeV, 0,
889 FBeta, aRandomEnergy);
890 aBetaPlusChannel->
SetICM(applyICM);
891 aBetaPlusChannel->
SetARM(applyARM);
893 theDecayTable->
Insert(aBetaPlusChannel);
903 if (modeFirstRecord[3]) {
904 modeFirstRecord[3] =
false;
911 aKECChannel->
SetICM(applyICM);
912 aKECChannel->
SetARM(applyARM);
914 theDecayTable->
Insert(aKECChannel);
921 if (modeFirstRecord[4]) {
922 modeFirstRecord[4] =
false;
928 aLECChannel->
SetICM(applyICM);
929 aLECChannel->
SetARM(applyARM);
931 theDecayTable->
Insert(aLECChannel);
938 if (modeFirstRecord[5]) {
939 modeFirstRecord[5] =
false;
946 aMECChannel->
SetICM(applyICM);
947 aMECChannel->
SetARM(applyARM);
949 theDecayTable->
Insert(aMECChannel);
955 if (modeFirstRecord[6]) {
956 modeFirstRecord[6] =
false;
963 anAlphaChannel->
SetICM(applyICM);
964 anAlphaChannel->
SetARM(applyARM);
966 theDecayTable->
Insert(anAlphaChannel);
977 G4Exception(
"G4RadioactiveDecay::LoadDecayTable()",
"HAD_RDM_000",
998 if (theDecayMode !=
IT) {
999 theBR = theChannel->
GetBR();
1000 theChannel->
SetBR(theBR*modeTotalBR[theDecayMode]/modeSumBR[theDecayMode]);
1004 DecaySchemeFile.close();
1006 if (!found && lvl > 0) {
1011 anITChannel->
SetICM(applyICM);
1012 anITChannel->
SetARM(applyARM);
1014 theDecayTable->
Insert(anITChannel);
1016 if (!theDecayTable) {
1020 G4cerr <<
"G4RadoactiveDecay::LoadDecayTable() : cannot find ion radioactive decay file "
1023 return theDecayTable;
1027 G4cout <<
"G4RadioactiveDecay::LoadDecayTable()" <<
G4endl;
1031 #ifdef G4MULTITHREADED
1032 (*master_dkmap)[key] = theDecayTable;
1034 return theDecayTable;
1040 if (Z < 1 || A < 2)
G4cout <<
"Z and A not valid!" <<
G4endl;
1042 std::ifstream DecaySchemeFile(filename);
1043 if (DecaySchemeFile) {
1044 G4int ID_ion = A*1000 +
Z;
1045 theUserRadioactiveDataFiles[ID_ion] = filename;
1048 G4cout <<
"The file " << filename <<
" does not exist!" <<
G4endl;
1055 G4int theG, std::vector<G4double> theRates,
1056 std::vector<G4double> theTaos)
1059 theDecayRate.
SetZ(theZ);
1060 theDecayRate.
SetA(theA);
1061 theDecayRate.
SetE(theE);
1064 theDecayRate.
SetTaos(theTaos);
1080 theDecayRateVector.clear();
1082 G4int nGeneration = 0;
1083 std::vector<G4double> rates;
1084 std::vector<G4double> taos;
1088 rates.push_back(-1.);
1091 G4int A = ((
const G4Ions*)(&theParentNucleus))->GetAtomicMass();
1092 G4int Z = ((
const G4Ions*)(&theParentNucleus))->GetAtomicNumber();
1093 G4double E = ((
const G4Ions*)(&theParentNucleus))->GetExcitationEnergy();
1095 if (tao < 0.) tao = 1e-100;
1096 taos.push_back(tao);
1103 theDecayRateVector.push_back(theDecayRate);
1124 std::vector<G4double> TP;
1125 std::vector<G4double> RP;
1142 for (j = nS; j < nT; j++) {
1143 ZP = theDecayRateVector[j].GetZ();
1144 AP = theDecayRateVector[j].GetA();
1145 EP = theDecayRateVector[j].GetE();
1146 RP = theDecayRateVector[j].GetDecayRateC();
1147 TP = theDecayRateVector[j].GetTaos();
1149 G4cout <<
"G4RadioactiveDecay::AddDecayRateTable : daughters of ("
1150 << ZP <<
", " << AP <<
", " << EP
1151 <<
") are being calculated, generation = " << nGeneration
1155 aParentNucleus = theIonTable->
GetIon(ZP,AP,EP);
1159 for (
G4int k = 0; k < 7; k++) brs[k] = 0.0;
1162 for (i = 0; i < aTempDecayTable->
entries(); i++) {
1168 AD = ((
const G4Ions*)(theDaughterNucleus))->GetAtomicMass();
1169 ZD = ((
const G4Ions*)(theDaughterNucleus))->GetAtomicNumber();
1176 if (std::abs(daughterExcitation - level->
Energy()) < levelTolerance) {
1179 if (level->
HalfLife()*
ns >= halflifethreshold){
1181 theDecayTable->
Insert(theChannel);
1183 brs[theDecayMode] += theChannel->
GetBR();
1186 brs[theDecayMode] += theChannel->
GetBR();
1189 brs[theDecayMode] += theChannel->
GetBR();
1192 brs[2] = brs[2]+brs[3]+brs[4]+brs[5];
1193 brs[3] = brs[4] =brs[5] = 0.0;
1194 for (i= 0; i<7; i++){
1200 (
const G4Ions*) aParentNucleus, brs[0]);
1201 theDecayTable->
Insert(theITChannel);
1207 brs[1], 0.*
MeV, 0.*
MeV, 1,
false, 0);
1208 theDecayTable->
Insert(theBetaMinusChannel);
1214 aParentNucleus, brs[2], 0.*
MeV, 0.*
MeV, 1,
false, 0);
1215 theDecayTable->
Insert(theBetaPlusChannel);
1223 theDecayTable->
Insert(theAlphaChannel);
1234 for (i = 0; i < theDecayTable->
entries(); i++){
1237 theBR = theChannel->
GetBR();
1241 if (theNuclearDecayChannel->
GetDecayMode() ==
IT && nGeneration == 1) {
1242 A = ((
const G4Ions*)(theDaughterNucleus))->GetAtomicMass();
1243 Z = ((
const G4Ions*)(theDaughterNucleus))->GetAtomicNumber();
1244 theDaughterNucleus=theIonTable->
GetIon(Z,A,0.);
1247 aParentNucleus != theDaughterNucleus) {
1251 if (aTempDecayTable->
entries() ) {
1252 A = ((
const G4Ions*)(theDaughterNucleus))->GetAtomicMass();
1253 Z = ((
const G4Ions*)(theDaughterNucleus))->GetAtomicNumber();
1254 E = ((
const G4Ions*)(theDaughterNucleus))->GetExcitationEnergy();
1257 if (TaoPlus <= 0.) TaoPlus = 1e-100;
1267 taos.push_back(TaoPlus);
1273 long double ta1,ta2;
1274 ta2 = (
long double)TaoPlus;
1275 for (k = 0; k < RP.size(); k++){
1276 ta1 = (
long double)TP[k];
1280 theRate = ta1/(ta1-ta2);
1282 theRate = theRate * theBR * RP[k];
1283 rates.push_back(theRate);
1290 long double aRate, aRate1;
1292 for (k = 0; k < RP.size(); k++){
1293 ta1 = (
long double)TP[k];
1297 aRate = ta2/(ta1-ta2);
1299 aRate = aRate * (
long double)(theBR * RP[k]);
1303 rates.push_back(theRate);
1305 theDecayRateVector.push_back(theDecayRate);
1315 if (nS == nT) stable =
true;
1327 theDecayRateTable.
SetItsRates(theDecayRateVector);
1330 theDecayRateTableVector.push_back(theDecayRateTable);
1345 ed <<
" Could not open file " << filename <<
G4endl;
1346 G4Exception(
"G4RadioactiveDecay::SetSourceTimeProfile()",
"HAD_RDM_001",
1352 while (infile >> bin >> flux ) {
1354 if (NSourceBin > 99) {
1355 G4Exception(
"G4RadioactiveDecay::SetSourceTimeProfile()",
"HAD_RDM_002",
1359 SBin[NSourceBin] = bin *
s;
1360 SProfile[NSourceBin] = flux;
1368 {
G4cout <<
" Source Timeprofile Nbin = " << NSourceBin <<
G4endl;}
1383 if (!infile)
G4Exception(
"G4RadioactiveDecay::SetDecayBias()",
"HAD_RDM_003",
1390 theRadioactivityTables.clear();
1393 while (infile >> bin >> flux ) {
1395 if (NDecayBin > 99) {
1396 G4Exception(
"G4RadioactiveDecay::SetDecayBias()",
"HAD_RDM_004",
1399 DBin[NDecayBin] = bin *
s;
1400 DProfile[NDecayBin] = flux;
1402 decayWindows[NDecayBin] = dWindows;
1405 theRadioactivityTables.push_back(rTable);
1409 for ( i = 1; i<= NDecayBin; i++) DProfile[i] += DProfile[i-1];
1410 for ( i = 0; i<= NDecayBin; i++) DProfile[i] /= DProfile[NDecayBin];
1418 {
G4cout <<
" Decay Bias Profile Nbin = " << NDecayBin <<
G4endl;}
1433 fParticleChangeForRadDecay.
Initialize(theTrack);
1439 if (!isAllVolumesMode) {
1440 if (!std::binary_search(ValidVolumes.begin(), ValidVolumes.end(),
1444 G4cout <<
"G4RadioactiveDecay::DecayIt : "
1446 <<
" is not selected for the RDM"<<
G4endl;
1447 G4cout <<
" There are " << ValidVolumes.size() <<
" volumes" <<
G4endl;
1449 for (
size_t i = 0; i< ValidVolumes.size(); i++)
1450 G4cout << ValidVolumes[i] << G4endl;
1459 return &fParticleChangeForRadDecay;
1469 G4cerr <<
"G4RadioactiveDecay::DecayIt : "
1471 <<
" is not a valid nucleus for the RDM"<<
G4endl;
1480 return &fParticleChangeForRadDecay;
1485 if (theDecayTable == 0 || theDecayTable->
entries() == 0) {
1490 G4cerr <<
"G4RadioactiveDecay::DecayIt : decay table not defined for ";
1500 return &fParticleChangeForRadDecay;
1522 if ( products->
entries() == 1) {
1527 return &fParticleChangeForRadDecay;
1551 if (temptime < 0.) temptime = 0.;
1552 finalGlobalTime += temptime;
1553 finalLocalTime += temptime;
1556 products->
Boost(ParentEnergy, ParentDirection);
1563 G4cout <<
"G4RadioactiveDecay::DecayIt : Decay vertex :";
1564 G4cout <<
" Time: " <<finalGlobalTime/
ns <<
"[ns]";
1569 G4cout <<
"G4Decay::DecayIt : decay products in Lab. Frame" <<
G4endl;
1574 for (index=0; index < numberOfSecondaries; index++) {
1576 finalGlobalTime, currentPosition);
1588 G4cout <<
"DecayIt: Variance Reduction version " <<
G4endl;
1605 std::vector<G4double> PT;
1606 std::vector<G4double> PR;
1608 long double decayRate;
1612 G4int numberOfSecondaries;
1613 G4int totalNumberOfSecondaries = 0;
1617 std::vector<G4DynamicParticle*> secondaryparticles;
1618 std::vector<G4double> pw;
1619 std::vector<G4double> ptime;
1624 for (
G4int n = 0;
n < NSplit;
n++) {
1633 weight1 = 1./DProfile[nbin-1]
1634 *(DBin[nbin]-DBin[nbin-1])/NSplit;
1635 }
else if (nbin > 1) {
1636 weight1 = 1./(DProfile[nbin]-DProfile[nbin-2])
1637 *(DBin[nbin]-DBin[nbin-1])/NSplit;
1645 for (i = 0; i < theDecayRateVector.size(); i++) {
1646 PZ = theDecayRateVector[i].GetZ();
1647 PA = theDecayRateVector[i].GetA();
1648 PE = theDecayRateVector[i].GetE();
1649 PT = theDecayRateVector[i].GetTaos();
1650 PR = theDecayRateVector[i].GetDecayRateC();
1660 for (j = 0; j < PT.size(); j++) {
1662 decayRate -= PR[j] * (
long double)taotime;
1675 theRadioactivityTables[decayWindows[nbin-1]]->AddIsotope(PZ,PA,PE,weight1*decayRate,theTrack.
GetWeight());
1684 parentNucleus = theIonTable->
GetIon(PZ,PA,PE);
1696 if (theDecayChannel == 0) {
1700 G4cerr <<
" G4RadioactiveDecay::DoIt : cannot determine decay channel ";
1701 G4cerr <<
" for this nucleus; decay as if no biasing active ";
1706 tempprods =
DoDecay(*parentNucleus);
1711 tempprods = theDecayChannel->DecayIt(tempmass);
1712 weight *= (theDecayChannel->GetBR())*(decayTable->
entries());
1715 tempprods =
DoDecay(*parentNucleus);
1719 numberOfSecondaries = tempprods->
entries();
1720 currentTime = finalGlobalTime + theDecayTime;
1721 for (index = 0; index < numberOfSecondaries; index++) {
1724 pw.push_back(weight);
1725 ptime.push_back(currentTime);
1726 secondaryparticles.push_back(asecondaryparticle);
1736 totalNumberOfSecondaries = pw.size();
1738 for (index=0; index < totalNumberOfSecondaries; index++) {
1740 ptime[index], currentPosition);
1760 return &fParticleChangeForRadDecay ;
1778 if (theDecayChannel == 0) {
1780 G4cerr <<
"G4RadioactiveDecay::DoIt : can not determine decay channel";
1786 G4cerr <<
"G4RadioactiveDecay::DoIt : selected decay channel addr:";
1791 products = theDecayChannel->
DecayIt(tempmass);
1804 if (origin == forceDecayDirection)
return;
1805 if (180.*
deg == forceDecayHalfAngle)
return;
1806 if (0 == products || 0 == products->
entries())
return;
1824 if (daughterType == electron || daughterType == positron ||
1825 daughterType == neutron || daughterType == gamma ||
1833 G4cout <<
"CollimateDecayProduct for daughter "
1846 if (origin == forceDecayDirection)
return origin;
1847 if (forceDecayHalfAngle == 180.*
deg)
return origin;
1852 if (forceDecayHalfAngle > 0.) {
1855 G4double cosMin = std::cos(forceDecayHalfAngle);
1864 G4cout <<
" ChooseCollimationDirection returns " << dir <<
G4endl;
void RegisterIsotopeTable(G4VIsotopeTable *table)
G4double GetDaughterExcitation()
void SetBR(G4double value)
static G4LossTableManager * Instance()
G4double GetLocalTime() const
void InitialiseAtomicDeexcitation()
void SelectAVolume(const G4String aVolume)
std::map< G4String, G4DecayTable * > DecayTableMap
std::ostringstream G4ExceptionDescription
G4double GetKineticEnergy() const
G4bool IsRateTableReady(const G4ParticleDefinition &)
const G4DynamicParticle * GetDynamicParticle() const
G4double HalfLife() const
void SetTaos(std::vector< G4double > value)
G4String strip(G4int strip_Type=trailing, char c=' ')
void SetDecayRateC(std::vector< G4double > value)
void AddUserDecayDataFile(G4int Z, G4int A, G4String filename)
G4bool GetPDGStable() const
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
static G4Alpha * Definition()
const G4ThreeVector & GetPosition() const
void AddSecondary(G4Track *aSecondary)
G4TrackStatus GetTrackStatus() const
static G4Electron * Definition()
void SetTouchableHandle(const G4TouchableHandle &apValue)
G4VDecayChannel * GetDecayChannel(G4int index) const
G4double GetMeanLifeTime(const G4Track &theTrack, G4ForceCondition *condition)
G4int GetVerboseLevel() const
void SetMomentumDirection(const G4ThreeVector &aDirection)
static G4NuclearLevelStore * GetInstance()
G4ParticleDefinition * GetDefinition() const
void SetSourceTimeProfile(G4String filename)
G4bool IsApplicable(const G4ParticleDefinition &)
void Boost(G4double totalEnergy, const G4ThreeVector &momentumDirection)
void SetGeneration(G4int value)
static G4Positron * Definition()
void ClearNumberOfInteractionLengthLeft()
#define G4MUTEX_INITIALIZER
const G4String & GetParticleName() const
void SetWeight(G4double aValue)
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
void SetItsRates(G4RadioactiveDecayRates arate)
G4double GetTotalMomentum() const
virtual void Initialize(const G4Track &)
G4IonTable * GetIonTable() const
G4GLOB_DLL std::ostream G4cout
void ProposeLocalTime(G4double t)
void SetDecayRate(G4int, G4int, G4double, G4int, std::vector< G4double >, std::vector< G4double >)
G4DecayTable * GetDecayTable(G4ParticleDefinition *)
const G4ThreeVector & GetMomentumDirection() const
void CollimateDecayProduct(G4DynamicParticle *product)
static G4LogicalVolumeStore * GetInstance()
G4ThreeVector ChooseCollimationDirection() const
G4int GetDecayTimeBin(const G4double aDecayTime)
void SetProcessSubType(G4int)
void DeselectAllVolumes()
const G4String & GetParticleType() const
G4DecayProducts * DoDecay(G4ParticleDefinition &theParticleDef)
const G4ParticleDefinition * GetParticleDefinition() const
G4double GetGlobalTime() const
G4double GetTaoTime(const G4double, const G4double)
const G4TouchableHandle & GetTouchableHandle() const
G4VDecayChannel * SelectADecayChannel()
G4NuclearLevelManager * GetManager(G4int Z, G4int A)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4ParticleDefinition * GetDaughterNucleus()
void Insert(G4VDecayChannel *aChannel)
G4double GetMeanFreePath(const G4Track &theTrack, G4double previousStepSize, G4ForceCondition *condition)
G4LogicalVolume * GetLogicalVolume() const
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
void GetDecayRateTable(const G4ParticleDefinition &)
void SetNumberOfSecondaries(G4int totSecondaries)
G4DynamicParticle * PopProducts()
G4VParticleChange * pParticleChange
G4RadioactiveDecayMode GetDecayMode()
void AddDecayRateTable(const G4ParticleDefinition &)
static G4Neutron * Definition()
G4int NumberOfLevels() const
G4VPhysicalVolume * GetVolume() const
G4double GetWeight() const
void SetIonName(G4String name)
G4double GetPDGLifeTime() const
G4RadioactiveDecay(const G4String &processName="RadioactiveDecay")
G4VAtomDeexcitation * AtomDeexcitation()
G4DecayTable * LoadDecayTable(G4ParticleDefinition &theParentNucleus)
virtual G4DecayProducts * DecayIt(G4double parentMass=-1.0)=0
void ProposeTrackStatus(G4TrackStatus status)
void BuildPhysicsTable(const G4ParticleDefinition &)
G4double ShapeFactor(const G4BetaDecayType &, const G4double &p_e, const G4double &e_nu)
void SetHLThreshold(G4double hl)
void SetAnalogueMonteCarlo(G4bool r)
const G4NuclearLevel * NearestLevel(G4double energy, G4double eDiffMax=9999.*CLHEP::GeV) const
G4double FermiFunction(const G4double &W)
void DeselectAVolume(const G4String aVolume)
void SetE(G4double value)
void AddUserDecayDataFile(G4int Z, G4int A, G4String filename)
void CollimateDecay(G4DecayProducts *products)
void SetDecayBias(G4String filename)
static G4Gamma * Definition()
void SetGoodForTrackingFlag(G4bool value=true)
G4VParticleChange * DecayIt(const G4Track &theTrack, const G4Step &theStep)
G4GLOB_DLL std::ostream G4cerr
G4int GetBaryonNumber() const