-- PeterKammel - 04 Jan 2007

Diffusion Monte Carlo


  • Study P_out functions more carefully and with md+d scattering
  • Generate time and radius functions with "full kinetics" and understand them
  • Study range of mp diffusion with different scenarios and compare to Adamczak
  • Check effect of high energy pm's generated during cascade


study condition Cd (ppm) statistics output analyzed with
mudh2n mud in H2 0. 1e9 run/test/h2n.root make_h2.C->h2n.root, c_pout.C, plot_t("h2n.root", pout.C++
mudHD mud in H2+D 122. 1e8 run/test/mudHD.root make_h2.C->mudHD.root, c_pout.C, plot_t("mudHD.root", pout.C++
mupHD mup in H2+D 1.5. 1e8 run/test/mupHD mupHD.C, time, radius
mupHDC mud in H2+D 122. 1e8 run/test/mupHDC mudHD.C, time, radius

P_out function and decay rates

Extracted from mudh2n and mudHD simulation and already included in latest time analysis note.

Full kinetics range and time

mp+p hfs scattering
md+d hfs scattering
md+p scattering
mp->md transfer
ppm formation is not included, but could be.

Summary: Results follow expectations and document main features of diffusion processes. In principle, these results could be used for moment method or 3 par fit. But I leave that to Brendan, because his simulation better accounts for impact parameter and detector resolution.

Production condition: mupHD



Calibration condition: mupHDC



Special pm+p diffusion studies

Some model distributions for testing Steve's reconstruction algorithm

Technical notes

Simulation directory is /home/npl/kammel/mucap/Mucat/run/test Run Mucat on windows machine with MCK_2.bat mupHD

This creates output in directory mupHD. Analyze on the local directory. Copy make_h2 root files to test directory and continue with c_out.C etc. Completed studies saved as mupHD.tar.gz, for example.

Topic attachments
I Attachment History Action Size Date Who Comment
GIFgif mupHD.radius.gif r1 manage 17.0 K 2007-01-05 - 02:58 PeterKammel  
PDFpdf mupHD.radius.pdf r1 manage 23.9 K 2007-01-05 - 02:39 PeterKammel  
GIFgif mupHD.time.gif r1 manage 13.6 K 2007-01-05 - 03:00 PeterKammel  
PDFpdf mupHD.time.pdf r1 manage 23.9 K 2007-01-05 - 02:35 PeterKammel  
GIFgif mupHDC.radius.gif r1 manage 15.9 K 2007-01-05 - 03:07 PeterKammel  
PDFpdf mupHDC.radius.pdf r3 r2 r1 manage 25.5 K 2007-01-05 - 02:39 PeterKammel  
GIFgif mupHDC.time.gif r1 manage 14.8 K 2007-01-05 - 03:10 PeterKammel  
PDFpdf mupHDC.time.pdf r1 manage 24.0 K 2007-01-05 - 02:25 PeterKammel  

This topic: Main > TWikiUsers > PeterKammel > GeneralPK > DiffusionPK
Topic revision: r3 - 2007-01-06 - PeterKammel
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback