kaliveda-dev team mailing list archive

[John Frankland <1033907@xxxxxxxxxxxxxxxxxx>]

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