salut jerem
j'ai commité une nouvelle version de KVedaLoss sous bazaar
pour "débrider" les tables de pertes il suffit d'appeler la méthode statique
KVedaLoss::SetIgnoreEnergyLimits()
avant la création du premier matériau.
on ne peut pas laisser sans limites par contre, parce que pour certains ions
dans certains matériaux au-dessus d'une certaine énergie ça diverge sec!
donc je recalcule des limites (en principe plus élevées) pour chaque ion et
chaque matériau.
au pire, on garde la limite d'avant (<400AMeV pour Z<4, <250AMeV pour Z>3).
au maximum maximal du plus max on va jusqu'à 1 GeV/nucléon.
pour ton problème du début, ça donne:
root [0] KVedaLoss::SetIgnoreEnergyLimits()
root [1] KVMaterial *CsI = new KVMaterial("CsI",10.)
Info in <KVedaLoss::init_materials>: Initialising KVedaLoss...
*************************************************************************
* VEDALOSS STOPPING POWER & RANGE TABLES
*
*
*
* 1. Si Silicon Z=14 A= 28.1 rho= 2.330 g/cm**3
*
* 2. Myl Mylar Z= 4 A= 8.7 rho= 1.395 g/cm**3
*
* 3. NE102 Plastic Z= 3 A= 6.5 rho= 1.032 g/cm**3
*
* 4. Ni Nickel Z=28 A= 58.7 rho= 8.902 g/cm**3
*
* 5. C3F8 Octofluoropropane Z= 8 A= 17.1 rho= 0.008 g/cm**3
*
* 6. C Carbon Z= 6 A= 12.0 rho= 1.900 g/cm**3
*
* 7. Ag Silver Z=47 A=107.9 rho=10.500 g/cm**3
*
* 8. Sn Tin Z=50 A=118.7 rho= 5.750 g/cm**3
*
* 9. CsI CesiumIodide Z=54 A=129.9 rho= 4.510 g/cm**3
*
* 10. Au Gold Z=79 A=197.0 rho=19.300 g/cm**3
*
* 11. U Uranium Z=92 A=238.0 rho=18.950 g/cm**3
*
* 12. Air Air Z= 7 A= 14.4 rho= 0.001 g/cm**3
*
* 13. Nb Nobium Z=41 A= 92.9 rho= 8.570 g/cm**3
*
* 14. Ta Tantalum Z=73 A=180.9 rho=16.654 g/cm**3
*
* 15. Al Aluminium Z=13 A= 27.0 rho= 2.699 g/cm**3
*
* 16. KCl KCl Z=18 A= 37.3 rho= 1.987 g/cm**3
*
* 17. CF4 Tetrafluoromethane Z= 8 A= 17.6 rho= 0.004 g/cm**3
*
* 18. Ca Calcium Z=20 A= 40.1 rho= 1.550 g/cm**3
*
* 19. Ge Germanium Z=32 A= 72.6 rho= 5.323 g/cm**3
*
* 20. Cu Copper Z=29 A= 63.5 rho= 8.960 g/cm**3
*
* 21. Ti Titanium Z=22 A= 47.9 rho= 4.540 g/cm**3
*
* 22. Bi Bismuth Z=83 A=209.0 rho= 9.747 g/cm**3
*
* 23. V Vanadium Z=23 A= 50.9 rho= 6.110 g/cm**3
*
* 24. C4H10 Isobutane Z= 2 A= 4.2 rho= 0.002 g/cm**3
*
* 25. Pb Lead Z=82 A=207.2 rho=11.350 g/cm**3
*
* 26. PbS LeadSulphide Z=49 A=120.0 rho= 7.500 g/cm**3
*
* 27. Mg Magnesium Z=12 A= 24.3 rho= 1.738 g/cm**3
*
* 28. Li Lithium Z= 3 A= 6.9 rho= 0.534 g/cm**3
*
* 29. Zn Zinc Z=30 A= 65.4 rho= 7.133 g/cm**3
*
*
*
* TF1::Range::Npx = 20 TF1::EnergyLoss::Npx = 50
*
* TF1::ResidualEnergy::Npx = 20
*
*
*
* INITIALISATION COMPLETE
*
*************************************************************************
root [2] CsI->GetEIncOfMaxDeltaE(1,1)
(Double_t)1.93964834294177791e+02
root [3] CsI->GetEIncOfMaxDeltaE(2,4)
(Double_t)7.86740550904558177e+02
root [4] CsI->GetEIncOfMaxDeltaE(3,6)
(Double_t)1.47426372511762997e+03
root [5] CsI->GetEIncOfMaxDeltaE(4,7)
(Double_t)2.28176017249893493e+03
root [6] CsI->GetEIncOfMaxDeltaE(6,12)
(Double_t)5.00209433461537355e+03
ce qui est mieux qu'avant ;-)
bien sûr si tu veux savoir exactement quelle limite a été calculée
pour chaque ion/matériau:
root [7] CsI->GetEmaxValid(6,12)
(Double_t)6.08427001953125000e+03
root [8] CsI->GetEmaxValid(6,12)/12.
(double)5.07022501627604186e+02
--
You received this bug notification because you are a member of KaliVeda
Development Team, which is subscribed to KaliVeda.
https://bugs.launchpad.net/bugs/1033907
Title:
problème dans les calculs de pertes d'énergie
Status in KaliVeda data analysis framework:
Fix Committed
Status in KaliVeda 1.8 series:
Fix Committed
Status in KaliVeda 2.0 series:
In Progress
Bug description:
J'ai quelques soucis avec les calculs de pertes d'énergies. Premièrement, une
remarque sur un phénomène étrange. Lorsque je recherche l'énergie incidente qui
donne le delta E max (via la méthode KVMaterial::GetEIncOfMaxDeltaE(z,a) voici
les résultats que j'obtiens pour quelques praticules : (pour un material :
KVMaterial *CsI = new KVMaterial("CsI",10.))
z=1, a=1 -> 1.93964834288665799e+02
z=2,a=4 -> 7.86740550982209584e+02
z=3,a=6 -> 1.47426372546234097e+03 (jusque là tout va bien)
z=4,A=7 -> 1.75000000000000000e+03
Les valeurs obtenue à partir des z=4 semblent bidon car des C12 de 400MeV/A
s'arrêtent dans 10cm de CsI mais pourtant la méthode renvoie un delta E max de
250MeV/A... De plus quelque soit l'épaisseur tant qu'elle est supérieure à 4cm,
la même valeur est retournée.
Le deuxième point (peut être lié) est que je ne retrouve pas les mêmes
pertes d'énergie entre géant 4 et Kaliveda à partir d'environ 90
MeV/A. (testé sur des cartes delta e - e silicium-CsI, les grilles
d'identifications tracées ne se superposent pas aux lignes geant 4).
Ces précisions sont ajustables ou cela veut-il dire que Kaliveda n'est
plus fiable à haute énergie ?
Jérémie Dudouet
To manage notifications about this bug go to:
https://bugs.launchpad.net/kaliveda/+bug/1033907/+subscriptions
_______________________________________________
Mailing list: https://launchpad.net/~kaliveda-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~kaliveda-dev
More help : https://help.launchpad.net/ListHelp
