kaliveda-dev team mailing list archive

[John Frankland <1033907@xxxxxxxxxxxxxxxxxx>]

Pour le premier point:
le problème est que tu dépasses les limites de validité nominale du calcul:
CsI->GetEmaxValid(4,7) :  1.75000000000000000e+03
CsI->GetEmaxValid(6,12) :   3.00000000000000000e+03
etc.
en fait, pour tous les Z>3 l'énergie maximum pour laquelle les tables sont nominalement
correctes est de 250 MeV/nucléon, et d'habitude elles ne sont utilisées que pour des
énergies <= 100 MeV/nucléon (i.e. celles du Ganil)

tu peux t'en affranchir en faisant:

KVedaLoss*VL = (KVedaLoss*)CsI->GetRangeTable()
TF1* eloss = VL->GetMaterial("CsI")->GetDeltaEFunction(CsI->GetThickness()*VL->GetMaterial("CsI")->GetDensity(),z,a)//z=6,a=12
eloss->SetRange(0,12*1000); // Emax = 1GeV/nucleon
TF1* range = VL->GetMaterial("CsI")->GetRangeFunction(CsI->GetThickness()*VL->GetMaterial("CsI")->GetDensity(),z,a)//z=6,a=12
range->SetRange(0,12*1000)
eloss->GetMaximumX()
(const Double_t)1.20000000000000000e+04

peut-être que ça serait mieux d'enlever les limites et d'afficher un warning ?
ou même de proposer un mode "débridé" (no limits!) ?
que pense-t-il ??

-- 
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:
  New

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