Elong-14-03-27

From HallCWiki
Jump to navigationJump to search

ptrates.f Addition Check

D & N check All Cross-Sections
2014-03-27-cs-check.png 2014-03-27-cs-check-all.png

First thing, after a number of trouble-shooting steps (mainly making sure that all of the variables pass correctly to inclusive.f since I'm putting them in from ptrates.f instead of the normal input file), I was able to get Misak's code integrated into my rates code so that when you run either the cross sections come out the same. These are at E=6.6GeV, E'=6.07 GeV, theta=9.51deg, since those were the conditions I used previously to compare Bosted's code to Misak's. (Note: For the plots to the right, it's assuming ND3 target)

Once I made sure that the deuterium and nitrogen were working correctly, I added Misak's calculations for the other cross sections.

This immediately brings to mind how previously the Bosted code was giving a smaller dilution factor than Misak's. Comparing the full Misak dilution (including N, He, & D) to Bosted's code looks like this:

2014-03-27-fdil-check.png


Rates with Misak's Cross Sections