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.

Table 1: Performance of the Dijet Filter Test
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.

 

Table 2: Biases and Acceptances of Filter Test
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