- drach09's home page
- Posts
- 2022
- 2020
- June (1)
- 2019
- 2018
- 2017
- 2016
- 2015
- 2014
- December (13)
- November (2)
- October (5)
- September (2)
- August (8)
- July (9)
- June (7)
- May (5)
- April (4)
- March (4)
- February (1)
- January (2)
- 2013
- December (2)
- November (8)
- October (5)
- September (12)
- August (5)
- July (2)
- June (3)
- May (4)
- April (8)
- March (10)
- February (9)
- January (11)
- 2012
- 2011
- October (1)
- My blog
- Post new blog entry
- All blogs
Run-11 Transverse Jets: Embedding Needs (version 1)
Here, I begin to calculate my needs for an embedding request. I begin with information from Run-6 Monte Carlo as outlined in Pibero's TAMU webpage. I note that the linked webpage does not reflect the final request, and thus following this model may be problematic. It is nonetheless a place to start.
Table 1
partonic pT [GeV] | Nevents | cross section [mb] | weight |
---|---|---|---|
3-4 | 690000 | 1.287000e+00 | 1.548130e+08 |
4-5 | 686000 | 3.117000e-01 | 3.771297e+07 |
5-7 | 398000 | 1.360000e-01 | 2.836181e+07 |
7-9 | 420000 | 2.305000e-02 | 4.555119e+06 |
9-11 | 414307 | 5.494000e-03 | 1.100638e+06 |
11-15 | 420000 | 2.228000e-03 | 4.402952e+05 |
15-25 | 397200 | 3.895000e-04 | 8.139098e+04 |
25-35 | 400000 | 1.016000e-05 | 2.108200e+03 |
35-45 | 120000 | 5.299000e-07 | 3.665142e+02 |
45-55 | 118000 | 2.851000e-08 | 2.005364e+01 |
55-65 | 120000 | 1.449000e-09 | 1.002225e+00 |
Table 1 shows the 2006 simulation request used as the starting point for estimating the 2009 inclusive jet request. For the current request I need to verify that I do not need to go below 3 GeV/c in partonic pT to cover my needs for the 5-6 GeV/c reconstructed jet bin. If I need to go to 2 GeV/c or below, the numbers could change quite a bit.
Table 2
Reconstructed pT [GeV/c] | |||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Partonic pT [GeV/c] | 6-7.1 | 7.1-8.4 | 8.4-9.9 | 9.9-11.7 | 11.7-13.8 | 13.8-16.3 | 16.3-19.2 | 19.2-22.7 | 22.7-26.8 | 26.8-31.6 | 31.6-37.3 | Partonic Sum | ⟨Njets⟩ |
3-4 | 506 | 393 | 156 | 62 | 12 | 0 | 0 | 0 | 0 | 0 | 0 | 1129 | 0.001636232 |
4-5 | 1071 | 1003 | 555 | 217 | 45 | 10 | 6 | 0 | 0 | 0 | 0 | 2907 | 0.004237609 |
5-7 | 1703 | 1654 | 1163 | 574 | 169 | 31 | 1 | 0 | 0 | 0 | 0 | 5295 | 0.01330402 |
7-9 | 4504 | 6074 | 5468 | 3221 | 1265 | 335 | 51 | 6 | 1 | 0 | 0 | 20925 | 0.049821429 |
9-11 | 6747 | 12269 | 14532 | 11767 | 5460 | 1868 | 355 | 82 | 5 | 1 | 0 | 53086 | 0.128132037 |
11-15 | 7000 | 15432 | 24895 | 30465 | 22994 | 10828 | 2957 | 593 | 70 | 2 | 0 | 115236 | 0.274371429 |
15-25 | 4803 | 10916 | 20063 | 36499 | 51895 | 53818 | 34852 | 15119 | 3776 | 533 | 54 | 232328 | 0.584914401 |
25-35 | 3003 | 5564 | 7568 | 14407 | 26674 | 47155 | 74148 | 103715 | 95576 | 45069 | 9709 | 432588 | 1.08147 |
35-45 | 759 | 1250 | 1459 | 2361 | 3950 | 6701 | 11100 | 19031 | 30158 | 40383 | 36218 | 153370 | 1.278083333 |
45-55 | 700 | 1021 | 1092 | 1847 | 2656 | 4076 | 6394 | 10394 | 17112 | 27954 | 41368 | 114614 | 0.971305085 |
55-65 | 637 | 860 | 983 | 1477 | 2162 | 3259 | 4608 | 6937 | 10976 | 18038 | 29064 | 79001 | 0.658341667 |
Reco. Sum | 31433 | 56436 | 77934 | 102897 | 117282 | 128081 | 134472 | 155877 | 157674 | 131980 | 116413 |
Table 2 shows the raw jet yields for each bin of partonic and reconstructed pT. The basic idea is to estimate the uncertainty on the mean "true" pT for each reconstructed pT bin given the statistics from Run-6. I can, then, scale the partonic pT bins to achieve the desired statistics and use the result to estimate the number of events I will need for my study.
ROOT uses the following expression to calculate the uncertainty on the mean:
For the time being I will assume an RMS of 7 GeV/c for a conservative estimate of all reconstructed pT bins. Table 3 shows the estimated uncertainty of the mean under the current assumption.
Table 3
6-7.1 | 7.1-8.4 | 8.4-9.9 | 9.9-11.7 | 11.7-13.8 | 13.8-16.3 | 16.3-19.2 | 19.2-22.7 | 22.7-26.8 | 26.8-31.6 | 31.6-37.3 |
---|---|---|---|---|---|---|---|---|---|---|
0.039482544 | 0.029465915 | 0.025074631 | 0.021822107 | 0.020440069 | 0.019559407 | 0.019088953 | 0.01772993 | 0.017628607 | 0.019268326 | 0.020516218 |
In each bin, the estimated uncertainty appears to be smaller than my current needs. Thus, it may be possible to request even fewer events than the Run-6 request.
Table 4
Reconstructed pT [GeV/c] | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Partonic pT [GeV/c] | 6-7.1 | 7.1-8.4 | 8.4-9.9 | 9.9-11.7 | 11.7-13.8 | 13.8-16.3 | 16.3-19.2 | 19.2-22.7 | 22.7-26.8 | 26.8-31.6 | 31.6-37.3 | Partonic Sum |
3-4 | 7.58E+10 | 5.84E+10 | 1.43E+10 | 4.18E+09 | 6.89E+08 | 0 | 0 | 0 | 0 | 0 | 0 | 1.53E+11 |
4-5 | 3.88E+10 | 3.63E+10 | 1.22E+10 | 3.60E+09 | 5.87E+08 | 1.38E+08 | 4.16E+07 | 0 | 0 | 0 | 0 | 9.17E+10 |
5-7 | 4.58E+10 | 4.52E+10 | 1.88E+10 | 7.10E+09 | 1.78E+09 | 2.56E+08 | 1.78E+06 | 0 | 0 | 0 | 0 | 1.19E+11 |
7-9 | 1.97E+10 | 2.64E+10 | 1.42E+10 | 6.54E+09 | 2.30E+09 | 5.88E+08 | 1.05E+08 | 2.14E+07 | 1.04E+05 | 0 | 0 | 6.99E+10 |
9-11 | 7.10E+09 | 1.29E+10 | 8.96E+09 | 5.84E+09 | 2.45E+09 | 9.02E+08 | 2.23E+08 | 5.68E+07 | 3.07E+06 | 8.62E+05 | 0 | 3.85E+10 |
11-15 | 2.95E+09 | 6.49E+09 | 6.00E+09 | 6.15E+09 | 4.43E+09 | 2.29E+09 | 7.36E+08 | 1.76E+08 | 2.38E+07 | 7.39E+05 | 0 | 2.93E+10 |
15-25 | 3.74E+08 | 8.49E+08 | 8.57E+08 | 1.39E+09 | 2.07E+09 | 2.45E+09 | 1.86E+09 | 9.18E+08 | 2.51E+08 | 3.82E+07 | 3.92E+06 | 1.11E+10 |
25-35 | 6.02E+06 | 1.12E+07 | 7.25E+06 | 1.50E+07 | 3.24E+07 | 6.71E+07 | 1.19E+08 | 1.81E+08 | 1.75E+08 | 8.47E+07 | 1.89E+07 | 7.17E+08 |
35-45 | 2.64E+05 | 4.40E+05 | 2.12E+05 | 4.44E+05 | 9.08E+05 | 1.81E+06 | 3.31E+06 | 6.04E+06 | 9.96E+06 | 1.35E+07 | 1.25E+07 | 4.94E+07 |
45-55 | 1.33E+04 | 1.93E+04 | 8.07E+03 | 1.92E+04 | 3.49E+04 | 6.20E+04 | 1.06E+05 | 1.83E+05 | 3.11E+05 | 5.13E+05 | 7.81E+05 | 2.05E+06 |
55-65 | 6.16E+02 | 8.22E+02 | 3.61E+02 | 7.79E+02 | 1.43E+03 | 2.52E+03 | 3.93E+03 | 6.21E+03 | 1.01E+04 | 1.68E+04 | 2.77E+04 | 7.12E+04 |
Table 4 contains the weighted jet yields according to Pibero's webpage. These are from weighting the yields and combing JP1 and JP2 with relative fractions from the data. This obviously will be different from Run-11 triggers, but perhaps it will suffice to be informative for my request.
Table 5
6-7.1 | 7.1-8.4 | 8.4-9.9 | 9.9-11.7 | 11.7-13.8 | 13.8-16.3 | 16.3-19.2 | 19.2-22.7 | 22.7-26.8 | 26.8-31.6 | 31.6-37.3 |
---|---|---|---|---|---|---|---|---|---|---|
0.6235485 | 0.347296052 | 0.251494855 | 0.190481744 | 0.167118569 | 0.153028162 | 0.14575525 | 0.125740167 | 0.124307115 | 0.14850735 | 0.168366076 |
With a goal of σ(Mean) = 0.05, Table 5 shows the required scales for the reconstructed jet yields.
Table 6
Reconstructed pT [GeV/c] | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Partonic pT [GeV/c] | 6-7.1 | 7.1-8.4 | 8.4-9.9 | 9.9-11.7 | 11.7-13.8 | 13.8-16.3 | 16.3-19.2 | 19.2-22.7 | 22.7-26.8 | 26.8-31.6 | 31.6-37.3 | Partonic Sum |
3-4 | 4.73E+10 | 2.03E+10 | 3.59E+09 | 7.96E+08 | 1.15E+08 | 0 | 0 | 0 | 0 | 0 | 0 | 7.20E+10 |
4-5 | 2.42E+10 | 1.26E+10 | 3.07E+09 | 6.85E+08 | 9.82E+07 | 2.12E+07 | 6.06E+06 | 0 | 0 | 0 | 0 | 4.07E+10 |
5-7 | 2.86E+10 | 1.57E+10 | 4.72E+09 | 1.35E+09 | 2.97E+08 | 3.91E+07 | 2.59E+05 | 0 | 0 | 0 | 0 | 5.07E+10 |
7-9 | 1.23E+10 | 9.18E+09 | 3.57E+09 | 1.25E+09 | 3.84E+08 | 9.00E+07 | 1.53E+07 | 2.69E+06 | 1.29E+04 | 0 | 0 | 2.68E+10 |
9-11 | 4.43E+09 | 4.48E+09 | 2.25E+09 | 1.11E+09 | 4.10E+08 | 1.38E+08 | 3.25E+07 | 7.14E+06 | 3.82E+05 | 1.28E+05 | 0 | 1.29E+10 |
11-15 | 1.84E+09 | 2.26E+09 | 1.51E+09 | 1.17E+09 | 7.40E+08 | 3.50E+08 | 1.07E+08 | 2.22E+07 | 2.96E+06 | 1.10E+05 | 0 | 8.00E+09 |
15-25 | 2.33E+08 | 2.95E+08 | 2.15E+08 | 2.65E+08 | 3.46E+08 | 3.74E+08 | 2.70E+08 | 1.15E+08 | 3.12E+07 | 5.68E+06 | 6.60E+05 | 2.15E+09 |
25-35 | 3.76E+06 | 3.87E+06 | 1.82E+06 | 2.86E+06 | 5.41E+06 | 1.03E+07 | 1.73E+07 | 2.28E+07 | 2.17E+07 | 1.26E+07 | 3.18E+06 | 1.06E+08 |
35-45 | 1.64E+05 | 1.53E+05 | 5.33E+04 | 8.45E+04 | 1.52E+05 | 2.76E+05 | 4.82E+05 | 7.60E+05 | 1.24E+06 | 2.01E+06 | 2.10E+06 | 7.47E+06 |
45-55 | 8.32E+03 | 6.71E+03 | 2.03E+03 | 3.66E+03 | 5.83E+03 | 9.48E+03 | 1.55E+04 | 2.30E+04 | 3.86E+04 | 7.63E+04 | 1.31E+05 | 3.21E+05 |
55-65 | 3.84E+02 | 2.86E+02 | 9.08E+01 | 1.48E+02 | 2.40E+02 | 3.86E+02 | 5.72E+02 | 7.80E+02 | 1.26E+03 | 2.49E+03 | 4.66E+03 | 1.13E+04 |
Table 6 shows the weighted jet yields, now, scaled by the values in Table 5. The far-right-hand column reflects the jet yields summed across reconstructed jet bins.
Table 7
Partonic pT [GeV/c] | RATIO | Event Estimate |
---|---|---|
3-4 | 4.70E-01 | 3.24E+05 |
4-5 | 4.44E-01 | 3.04E+05 |
5-7 | 4.26E-01 | 1.70E+05 |
7-9 | 3.83E-01 | 1.61E+05 |
9-11 | 3.35E-01 | 1.39E+05 |
11-15 | 2.73E-01 | 1.15E+05 |
15-25 | 1.95E-01 | 7.73E+04 |
25-35 | 1.47E-01 | 5.89E+04 |
35-45 | 1.51E-01 | 1.82E+04 |
45-55 | 1.56E-01 | 1.85E+04 |
55-65 | 1.59E-01 | 1.90E+04 |
Estimated Total | 1.40E+06 |
Table 7 shows the ratio of scaled to weighted partonic jet yields summed across reconstructed jet bins. The estimated number of required events is then the product of the ratio and the original Run-6 request. The overall total is shown in blue. Again, this table does not include a 2-3 GeV/c bin. One such bin was included for the 2009 analysis, and it would be no surprise if it is needed, here.
Figure 1
In Fig. 1 I show the normalized distribution of reconstructed jet pT for each bin of partonic pT. The black points are the 2006 jets as described in Table 4, normalized to sum to unity. The red crosses are the existing embedding data for the 2011 analysis (generated for 2009 dijets), again, normalized to sum to unity. The shapes in the overlapping bins appear to be quite similar, beginning to diverge noticeably by the 55-65 GeV/c bin.
- drach09's blog
- Login or register to post comments