kaliveda-dev team mailing list archive

[John Frankland <719251@xxxxxxxxxxxxxxxxxx>]

** Changed in: kaliveda
       Status: In Progress => Fix Committed

-- 
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:
  Fix Committed

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;