Elong-14-03-27
ptrates.f Addition Check
D & N check | All Cross-Sections |
---|---|
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: