** Changed in: kaliveda
Importance: Undecided => Critical
** Changed in: kaliveda
Milestone: None => 1.8.1
--
You received this bug notification because you are a member of KaliVeda
Development Team, which is subscribed to KaliVeda.
https://bugs.launchpad.net/bugs/719251
Title:
Erreur dans calcul de Sigma Rutherford
Status in KaliVeda data analysis framework:
New
Bug description:
KV 1.8.0 au CC
Le calcul de la section efficace de diffusion Rutherford donne un résultat
incorrect.
Après examen des détails, il me semble que le problème vient du calcul de
l'angle cm (du proj)
pour un angle lab donné, le retour est thetaCM=thetaLab. Par exemple :
KVDetector *det=(KVDetector*)gIndra->GetDetectorByType(1,2,4); // 4=Si
KVTelescope*tel=det->GetTelescope();
runl=gIndraDB->GetRun(7291);
KV2Body *kin = runl->GetSystem()->GetKinematics()
KVNucleus *proj=new
KVNucleus(runl->GetSystem()->GetZproj(),runl->GetSystem()->GetAproj())
kin->CalculateKinematics();
kin->SetOutgoing(3,proj);
kin->GetXSecRuthCM(2,3)
(const Double_t)2.38437563720703074e+04
kin->GetXSecRuthLab(2,3)
(const Double_t)2.36225764072999409e+04
runl->GetSystem()->GetEproj()
(const Float_t)3.20000000000000000e+01
kin->GetThetaCMProj(2,3)
(const Double_t)2.00000000000000000e+00
Comme on le voit l'angle CM renvoyé est égal à l'angle lab, et le valeur
de
XSecRuthCM correspond bien à ThetaCM=2°. En regardant la source, on
voit que cela est possible si la valeur de K est nulle.
J'ai essayé de taper kin->SetOutgoing(3,proj) avant de calculer la
cinématique, et j'obtiens
un résultat différent, mais toujours incorrect, ThetaCM=2.47°.
La valeur correcte est ThetaCM=4.19° (calcul non relativiste)
N'ai-je pas défini assez bien mon canal de sortie ? Ou y a t-il
vraiment un bug ?
Que ceux qui ont utilisé ces fonctions me viennent en aide;
_______________________________________________
Mailing list: https://launchpad.net/~kaliveda-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~kaliveda-dev
More help : https://help.launchpad.net/ListHelp
