| Summary: | Incorrect range of ions (not protons) in water with Opt3 and Opt4 | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | miancortes |
| Component: | processes/electromagnetic/standard | Assignee: | Vladimir.Ivantchenko |
| Status: | RESOLVED FIXED | ||
| Severity: | major | CC: | Susanna.Guatelli |
| Priority: | P4 | ||
| Version: | 11.0 | ||
| Hardware: | All | ||
| OS: | All | ||
| Attachments: |
Li7 range in water
C12 range in water (master-20220418) O16 range in water (master-20220418) C12 range in water (master-20220512) O16 range in water (master-20220512) |
||
Miguel, can you please, repeat the test with Hello Miguel, would it be possible to perform the test again using current master version? There are several updates, which may slightly improve ranges and dedx. Vladimir Hi Vladimir. Sure, I'll queue these calculations asap. Miguel Created attachment 768 [details]
C12 range in water (master-20220418)
Hi Vladimir, I just produced the plots using master as of yesterday. Generally, option3 has improved whereas option4 remains with the same behaviour. I uploaded all results at https://uses0-my.sharepoint.com/:u:/g/personal/miancortes_us_es/Ecu_uYSrxvdLqHvj1fgrHEUBcN35hlxGxbUKRZvc0nGBMQ New ones are under results/master-20220418 directory; I kept previous ones in separate folders just in case. Deviations have now the opposite sign (i.e. G4EmStandard_option3 produces slightly larger ranges for Li-7, C-12 and O-16). It seems that the relative deviation increases a bit with Z, as I observe that O-16 reaches 1 mm deviation for about 20 cm range, C-12 at about 22 cm, and Li-7 remains below 1 mm for all energies (up to about 28 cm in water). I attach C-12 and O-16 figures, too. Cheers, Miguel Created attachment 769 [details]
O16 range in water (master-20220418)
Hi Vladimir, I verified results with master version as of 2022-05-12, and opt4 has improved. Now opt4 and opt3 give equally good results, as shown in plots attached (C12 and O16 cases). Cheers, Miguel Created attachment 776 [details]
C12 range in water (master-20220512)
Created attachment 777 [details]
O16 range in water (master-20220512)
The problem is fixed and the fix is available in 11.1beta. |
Created attachment 756 [details] Li7 range in water I ran a testing code of G4Med (LightIonBraggPeak) aimed to test range of low-Z ions beams in water (H, Li-7, C-12 and O-16). I found out that QGSP_BIC_HP_EMY and _EMZ (i.e. Opt3 and Opt4) gave significant deviations with respect to experimental data (Schardt et al.) when using 11-00-patch-01. Older tags (such as 10-05-ref-06) did pretty well. Please see attached Li-7 case. This happens to all ions tested except protons - the latter case keeps performing well. In this code, scoring cells are physically created within the water tank (i.e., no scoring mesh is being used). Also, G4_WATER material is used with a modified density to correct due to room temperature using this method: man->BuildMaterialWithNewDensity("Water_RoomTemp","G4_WATER",theWaterDensity); man->FindOrBuildMaterial("Water_RoomTemp"); I verified that both water's were properly created: Index : 3 used in the geometry : Yes Material : Water_RoomTemp Range cuts : gamma 25 um e- 25 um e+ 25 um proton 20 um Energy thresholds : gamma 1 keV e- 38.2486 keV e+ 37.8703 keV proton 2 keV Region(s) which use this couple : ROI Index : 4 used in the geometry : Yes Material : G4_WATER Range cuts : gamma 25 um e- 25 um e+ 25 um proton 20 um Energy thresholds : gamma 1 keV e- 38.3139 keV e+ 37.9365 keV proton 2 keV Region(s) which use this couple : ROI I must remark that he relative deviations are always of about 1.25-1.28 percent for all ions - maybe it is an issue related with water density?