Comparison of HARP tantalum 5% interaction length data against simulations using MARS 15.07

Stephen Brooks (   August 2007

See the HARP publication, the goal is to replicate figure 30 in a simulation.

The graph shows yield of π+ and π integrated over three different regions:

Incident proton momenta used in the graph are 3, 5, 8 and 12 GeV/c.

Simulations were made in MARS for each of these cases with at least 1 million incident protons and the yield of pions integrated over the above ranges.  Incoming beam was assumed to be a pencil beam (no spatial or angular spread) with a 1% RMS Gaussian spread in momentum around the nominal value.

Dimensions for the 5% interaction length tantalum target were taken from  The simulations used a cylinder of length 5.6 mm (nominal) and radius 15.07063 mm (average of measured values).

Table 1.  Integrals of pion yield in the three regions.
Proton momentum (GeV/c)Incident protons simulatedπ+ per protonπ per proton
33.453 × 1060.0065890.004240.0022440.0064660.0038750.001728
51 × 1060.0159150.0115910.0057590.0166320.0115570.005354
81 × 1060.0241460.0183040.00890.0255170.0189290.009239
121 × 1060.0273870.0208830.0101010.0309780.0229040.01094

Note that the three integration regions are listed left-to-right in decreasing order of size in these tables.  In figure 1 it should be clear which are which, as with the graph it is based on.

Table 2.  Integrals of pion "efficiency", i.e. yield divided by kinetic energy of the incident protons, in the three regions.
Proton momentum (GeV/c)π+ per proton.GeVπ per proton.GeV

I suspect the "arbitrary units" in the HARP publication's figure 30 are in fact units of 10-4 of the ones used here.  It would be helpful to compare the absolute numbers here if they are available.

[Figure 1]
Figure 1.  Data for both pion signs over each of three integration regions, plotted in the same way as the HARP paper.

Valid XHTML 1.1