DiJet Filter Performance at 500GeV
Using the similar filtering procedure developed by Matthew Walker (see link), I examined the performance of this filter for proton-proton collisions at a center of mass energy of 500 GeV. There are two filters implemented in this setup: a dijet filter executed at the particle level (used to create the fzd files) and a trigger filter invoked before the computationally expensive track reconstruction (makes the MuDst and geant files). Every filter performed raises concerns of biasing the data (removing events that would have reconstructed dijets) and thus possibly altering physics. However, a sucessful filter, one that doesn't give a large bias, reduces the computation time and provides a sample with more events of the desired signal.
Therefore, I adminstered a test of the filter setup to determine the characteristics needed to make a dijet filtered simulation request. It is important to note that in this test no events were removed. Instead, I record only the result if the filter were applied. For example, I print out whether or not an event at the pythia level contained the signal.
At the Pythia level a dijet filter is performed and it requires signal events to have a dijet with the following characteristics:
both jet |eta| < 1.3
delta phi > 2
max pt > 10.0
min pt > 7.0
dijet M > 10.0
Once the dijet pythia filter is completed, the TriggerFilter is then applied, which accepts events that would have fired the JP1, AJP or BHT3 triggers.
The following table displays the results of this test of the filter performance. I threw 50K events in each partonic pT bin. The first column indicates the actual events thrown and for some reason events were not reconstructed; I did not try to recover them. In column 3 and 4 we repectively see the number events that pass the dijet pythia filter and both the pythia and trigger filter. Column 5 informs us about the number of events that failed the pythia filer but had a dijet signal after the BFC was performed. However, these events did not pass the Trigger filter. Column 6 is the true indicator if there is a bias for each pT bin. It tells us the number of events that failed the pythia filter but had dijets after the BFC and passed the Trigger filter. Column 7 tells the number of events that pass the trigger filter and column 8 shows the number of events need to obatin 1pb-1 for that pT bin.
pT Bin | Events Thrown | Events Passing Pythia Filter | Events Passing Pythia Filter & Trigger | Events Failing Pythia Filter but had dijets after BFC | Events Failing Pythia Filter but had dijets after BFC & Passed the Trigger | Events Passing Trigger Filter | Events needed for 1pb-1 |
3-4 | 50000 | 348 | 0 | 0 | 0 | 0 | 5.23E+09 |
4-5 | 50000 | 845 | 7 | 0 | 0 | 17 | 1.48E+09 |
5-7 | 45058 | 2117 | 20 | 0 | 0 | 76 | 7.54E+08 |
7-9 | 49994 | 6276 | 207 | 2 | 1 | 518 | 1.50E+08 |
9-11 | 40990 | 12726 | 853 | 6 | 4 | 1640 | 4.25E+07 |
11-15 | 48955 | 18039 | 2930 | 107 | 15 | 4724 | 2.11E+07 |
15-25 | 47845 | 23708 | 6687 | 352 | 64 | 9985 | 5.59E+06 |
25-35 | 48144 | 31444 | 23324 | 225 | 110 | 31806 | 3.38E+05 |
35-45 | 49602 | 37410 | 32589 | 116 | 63 | 40614 | 4.18E+04 |
45-1000 | 47694 | 43288 | 38681 | 59 | 44 | 43370 | 9.97E+03 |
Now one can use this table to determine the biases of the filter for each pT bin. This is shown in Table 2.
pT bin | Filter Bias | Pythia Filter Acceptance | Total Filter Acceptance | Total Expected Events |
3-4 | 0 | 6.96E-03 | 2.0E-05 | 1.05E+05 |
4-5 | 0 | 0.0169 | 1.40E-04 | 2.07E+05 |
5-7 | 0 | 0.0423 | 4.00E-04 | 3.02E+05 |
7-9 | 2.00E-05 | 0.1255 | 4.13E-03 | 6.21E+05 |
9-11 | 8.00E-05 | 0.2545 | 0.0171 | 7.48E+05 |
11-15 | 3.01E-04 | 0.3608 | 0.0586 | 1.24E+06 |
15-25 | 1.28E-03 | 0.4742 | 0.134 | 7.48E+05 |
25-35 | 2.20E-03 | 0.6289 | 0.467 | 1.58E+06 |
35-45 | 1.26E-03 | 0.7482 | 0.652 | 2.72E+04 |
45-1000 | 8.80E-04 | 0.8658 | 0.774 | 7.71E+03 |
In order for a filter to be considered unbiased the values in column 2 of Table 2 should be << 10-3. One can then see that the largest bias occurs at the partonic pT bin of 25-35GeV. Therefore, using the tables above as an indicator, we propose to use apply the full filter on pT bin of 3-11GeV, the Trigger Filter of pT bin 11-25GeV, and no filter on 25GeV and above.
From the analysis meeting in March of 2011, I was informed that JP1 and JP2 trigger thresholds are higher than those at 2009 200 GeV. Therefore, the asymetric pT cuts (see above) should be increased as well. I have changed this values to be
max pt > 13.0 GeV
min pt > 10.0 GeV
pT bin | Events Thrown | Events Passing Pythia Filter |
Events Passing Pythia &Trigger Filter |
Events Failing Pythia Filter but had dijets after BFC | Events Failing Pythia Filter but had dijets after BFC & Passed the Trigger | Events Passing Trigger Filter | Events needed for 1pb-1 |
11-15 | 49966 | 17941 | 2944 | 41 | 19 | 4913 | 2.01E+07 |
15-25 | 49744 | 23694 | 9821 | 182 | 64 | 14631 | 5.57E+06 |
25-35 | 48144 | 29438 | 22218 | 36 | 18 | 32398 | 3.22E+05 |
Using this information I began calculating the time and disk space needed for the simulation request. I threw 10K events with the filters implemented as indicated above and the results can be seen on the follwing table:
pT bin | Events Thrown | Fzd File Size | average time/event | Trigger Accepted | CPU time | Real Time | MuDst File Size |
Geant File Size |
Pythia File Size |
7-9 | 79062 | 1.82 GB | 6.642 | 323 | 1.07E+06 | 1.08E+06 | 26.4MB | 73.28 MB | 51.40 MB |