Reproducing a rough estimate.

 An email contained the following notes from Bill Christie regarding the BHT3 cross section, and I want to reproduce this logic to 

Absolute lumi run 10088078 , low dead time, only BHT3=31 + BBC,ZDC prescaled

• BHT3/baseZDC=1.64e-4

• ZDC cross section from Angelika=2.3 mb,  crude

• Joe's calculation of luminosityHi All,

 Here is the calculation of the run9 cross section number. So Angelika gave us 4 vernier scan luminosities

 

fill#        date       time           L(10^31 s^-1cm^-2)

10207   02/23    11:06-09  0.3

10276   03/01   11:32-36   1.8

10399   03/20   13:24-30   2.6

10415   03/24   10:02-08   2.8

 

Only the last one of those time periods had a trigger mix that seem to working correctly. I looked for the closest (in time) run number before the scan (run# 10083059) and got the ZDC rate as given by the Trigger Details link in the RunLog. That number was 71513.5. All runs near the other scans gave the ZDC rate as 1000 exactly and I figure that is not correct so I will only calculate the cross section from this last number. So starting with the standard formula

dN/dt = Rate = xsec*instant_lumi

and inverting

xsec = Rate/instant_lumi = 7.15e4/2.8e31 = 2.554e-27 (cm^2) = 2.554 mb

 

Note that this number disagreed with the number Bill remembered earlier, but provided I am getting the correct ZDC rate this should be correct. Now from run# 10088078 (that we took earlier today) we know that the ratio of the BHT3 xsec to the ZDC xsec is 1.64e-4. Recall that we restricted the run to a low rate so that dead time is not an issue in any of the determinations.

xsec(BHT3)/xsec(ZDC) = 1.64e-4

xsec(BHT3) = 1.64e-4*xsec(ZDC) = 4.18e-4 mb

 

 

 

 

The first step is to try and get the ZDC cross section.  I perform the same analysis as on the previous pages, but this time for the ZDC trigger in 10097097.  I fill the bare rate plot for the good bunches (all bunches except those in the abort gaps, or otherwise only partially filled).  This results in a clean rate plot very much like the BBC plots we've seen earlier.

 

Figure 1.  Rate as a function of time for the ZDC trigger in 10097097.  Entries in this histogram come from only the good bunch crossings where both beams had full buckets.  The vernier scan sweep can be very clearly seen.

 

Running through the same alignment and rebinning, we can then fit as before:

 

 

Figure 2.  Counts as a function of vernier scan set, with fit superimposed.  The resulting widths are a good deal smaller than I expected, differing from the BHT3 widths by roughly 20%.

 

The resulting cross section is 287nb, which we then have to correct for the prescale factor (in this case 3000 [listed] * 1 [from details] ), so the corrected cross section is 0.86mb  This is pretty far off from the 2.5mb that Bill has from a different process.  If I assume the widths are somehow aritficially small, I can correct the cross section by multiplying by the ration of the 'correct' widths to these (~23/17).  This yields 1.17mb, which is close to half of the value Angelika reports and Bill calculates.

 

For the moment I can press on.  The scaler/base_rate values for the ZDC and BHT3 triggers are 1.19e-1 and 8.78e-5, respectively.  We can calculate the BHT3 cross section as being 0.86mb* the ratio of these two values, which results in 634nb, compared to the 525nb value we get from performing the vernier scan fit directly the BHT3 data.  The value Bill gets is not significantly different from this, coming out as 418nb.

 

If I use the larger cross section (after 'correcting'), I'll get an even larger BHT3 cross section, 857nb.

If I use the ZDC cross section reported by Angelika, the result is 1696nb.

 

It's possible that the vernier scan run 10097097 is too high luminosity for the ZDC to be deadtime free, in which case we'd expect the real ZDC cross section to be larger than the value seen, which would once again shift the cross section in the wrong direction.

 

Have I missed something, or is this ~20% agreement close enough given the rough values used by Bill?