Message: Incoherence on Nuclear Stopping power since last build Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Warning Incoherence on Nuclear Stopping power since last build 

Forum: Hadronic Processes
Date: 23 Mar, 2017
From: <Synthex>

Dear developers.

Still on my low energy ions application, I decided to update my installation, thanks to the Centos7 virtual machine provided by IN2P3.
My previous version was G4 10-02 (4/12/15)
The new is geant4-10-03-patch-01 [MT]   (24-February-2017)

Running the Hadr01 example in order to test PhysicsList impact on stopping power, I noted the following:
On RunAction, an EmCalculator was added to retrieve F19 stopping power in NIST_WATER.
  G4LossTableManager::Instance();
  G4NistManager* NISTman = G4NistManager::Instance();
G4Material* Water=NISTman->FindOrBuildMaterial("G4_WATER");
     G4double Thickness = DBL_MAX;
    G4double SiLET     = 0;
    G4double SiNuclearDEDX = 0;
    G4double SiElectronicDEDX = 0;
    G4double CSDAr = 0;
    G4double exDEDX=0;
   G4EmCalculator EmCalculator;
EmCalculator.SetVerbose(0);
  G4IonTable* ionTable = G4IonTable::GetIonTable();
   G4ParticleDefinition* F19 = ionTable->GetIon(9,19,0.0);
 

for (G4int i=100;i>0;i-=10)  {
  SiElectronicDEDX = EmCalculator.ComputeElectronicDEDX(i*keV,F19,Water);
  	  SiNuclearDEDX = EmCalculator.ComputeNuclearDEDX(i*keV,F19,Water);

	  CSDAr = EmCalculator.GetRangeFromRestricteDEDX(i*keV,F19,Water);

	  exDEDX = EmCalculator.GetDEDX(i*keV,F19, Water);
	  G4cout <<"At this energy: "<< i<<" keV" << G4endl;
	  G4cout <<" Elec. STOPPING POWER  : "<<SiElectronicDEDX/(MeV/mm) << "\t"<< " Nucl. STOPPING POWER  : " << SiNuclearDEDX/(MeV/mm)<<"\t"<<" GetDEDX : " << exDEDX/(MeV/mm) << "\t" << "restricted range : " << G4BestUnit(CSDAr,"Length") << G4endl;
 }


On the 10-02 build, the following were obtained with QBBC:
Energy         Ele DEDX(MeV/mm)        NucDEDX(MeV/mm)       resRange(nm)
100keV         211                     80                    717
50keV          176                     113                   460
10keV          103                     180                   197

On the 10-02 build with emstandard_opt4, the following:
Energy         Ele DEDX(MeV/mm)        NucDEDX(MeV/mm)       resRange(nm)
100keV         197                     80                    692
50keV          139                     113                   444
10keV          62                      180                   190




Results of the 10-03P1 build with the same file, QBBC outputs:
Energy         Ele DEDX(MeV/mm)        NucDEDX(MeV/mm)       resRange(nm)
100keV         211                     51.54                 717
50keV          176                     65.45                 461
10keV          103                     84.16                 197

Results of the 10-03P1 build with the same file, emstandard_opt4 outputs:
Energy         Ele DEDX(MeV/mm)        NucDEDX(MeV/mm)       resRange(nm)
100keV         197                     51.54                 717
50keV          139                     65.45                 461
10keV          62                      84.16                 197



Besides, SRIM 2013 gives
Energy         Ele DEDX(MeV/mm)        NucDEDX(MeV/mm)       resRange(nm)
100keV         141                     79                    397
50keV          122                     111                   195
10keV          48                      176                   40


I have 3 concerns:

1) The previous build gave Nuclear Stopping power in accordance with Ziegler, not the current one. Is this a bug or some features needs to be activated now?

2) Why is the GetRangeFromRestrictedDEDX not affected by the Nuclear Stopping Power change?

3) Why so much difference with Ziegler on Electronic Stopping power, even with LowEPPackage, Livermore, Penelope and _option3/4?

Thank you

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Idea: Re: Incoherence on Nuclear Stopping power since last build   (Vladimir Ivanchenko - 04 May, 2017)
 Add Message Add Message
to: "Incoherence on Nuclear Stopping power since last build"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]