Related inks:
Number of generated events per partnic pt bin:
gamma-jet: 25K
QCD(2-4): 50K
QCD(4-55): 25K
StEemcGammaFilter:: running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events
StEemcGammaFilter:: mConeRadius 0.22 mSeedThreshold 3.8 mClusterThreshold 5 mEtaLow 0.95 mEtaHigh 2.1 mMaxVertex 120
StEemcGammaFilter:: mCalDepth 279.5 mMinPartEnergy 1e-05 mHadronScale 1 mFilterMode 0 mPrintLevel 1
StChain:INFO - Init() : Using gamma filter on the EEMC
StChain:INFO - Init() : EEMC Sampling Fraction = 0.05
StChain:INFO - Init() : Seed energy threshold = 3.8 GeV
StChain:INFO - Init() : Cluster eT threshold = 5 GeV
StChain:INFO - Init() : Maximum vertex = +/- 120 cm
StChain:INFO - Init() : Running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events in BFC
BFC:INFO - setTowerGainSpread(): gain spread: 0; gain mean value: 1.1
runSimuGammaTreeMaker():: GammaMaker config: ConeRadius 0.7 ClusterEtThreshold 5.5 SeedEnergyThreshold 4.2 ClusterEnergyThreshold 5.5 BsmdRange 0.05237 EsmdR ange 20
StEEmcA2EMaker *EEanalysis = new StEEmcA2EMaker("mEEanalysis");
EEanalysis->threshold(3.0, 0); // tower threshold (ped+N sigma)
EEanalysis->threshold(3.0, 1); // pre1 threshold
EEanalysis->threshold(3.0, 2); // pre2 threshold
EEanalysis->threshold(3.0, 3); // post threshold
EEanalysis->threshold(3.0, 4); // smdu threshold
EEanalysis->threshold(3.0, 5); // smdv threshold
emulate L2E-gamma trigger for Run 2006 [eemc-http-mb-l2gamma:: id 137641]
Trigger conditions:
cluster Et (2x2) = 5.2GeV
seed Et = 3.7GeV
Definition: offline selected events are events which satisfy to the following conditions:
Figure 1a:
(upper plots) Gamma candidate yields vs. candidate pt (all partonic pt bins, no pt weights)
(lower plots) False rejection: histograms in the upper panel scaled by L2E-gamma-trigger yield (blue histogram)
Figure 1b: Same ad Fig. 1a with zoom into low pt region
Figure 2: Same ad Fig. 1b for partonic pt=2-3
Figure 3: Same ad Fig. 1b for partonic pt=3-4
Figure 4: Same ad Fig. 1b for partonic pt=4-6
Figure 5: Same ad Fig. 1b for partonic pt=6-9
Figure 6: Same ad Fig. 1b for partonic pt=9-15
Related inks:
Figure 1: BFC filter processing time for accepted events
Note an extra peaks around 23/33/43 seconds for the QCD sample
(they also present but less pronounced in gamma-jet sample)
I found that these are processing times needed for the first accepted by filter events
(not always the time for the first processed event and depends on filter acceptance rate).
Processing is much longer due to time needed to initiaize additional stuff in bfc makers
and it is ignored in the total cpu time estimate
Figure 1b:
1st row: starsim time per partonic pt bin
2nd row: bfc time for accepted events per partonic pt bin
3rd row: bfc time for rejected events per partonic pt bin
Average starsim/bfc timing (ignoring times in Fig. 1 with more that 18 seconds):
Gamma-jets bfc:acc starsim bfc:rej
pt=25-55 6.49765 20.2376 0.139911
pt=15-25 4.04562 11.8268 0.095903
pt=9-15 4.84475 11.4816 0.112344
pt=6-9 6.19909 13.0528 0.143109
pt=4-6 4.86546 9.8856 0.114983
pt=3-4 5.07415 10.404 0.12363
pt=2-3 4.04254 9.04627 0.0995413
QCD bfc:acc starsim bfc:rej
pt=25-55 6.18626 14.9944 0.13752
pt=15-25 6.20077 12.9668 0.126722
pt=9-15 6.848 13.7706 0.140625
pt=6-9 5.29513 10.3708 0.110363
pt=4-6 6.29547 11.5318 0.127077
pt=3-4 4.45859 10.1538 0.0986285
pt=2-3 5.26187 13.9151 0.114372
Related inks:
Pythia filter bug details:
bug in Pythia Filter StEemcGammaFilter.cxx:
double scale = (mCalDepth-v[2]) / p[2];
for(unsigned int j = 0; j < 3; ++j) p[j] = p[j] * scale + v[j];
Should be "abs(p[2])" in the "scale" factor.
Otherwise this will transfer all - Eta -> Eta,
and such tracks will pass the consequent Endcap rapidity cut:
if(detectorV.Eta() < mEtaLow || detectorV.Eta() > mEtaHigh) continue;
what will increase Pythia filter accept rate by a factor of 2.
Number of generated events per partnic pt bin:
Note: not every jib finished due to RCF scheduler upgrade
gamma-jet: 25K
QCD(2-4): 50K
QCD(4-55): 25K
StEemcGammaFilter:: running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events
StEemcGammaFilter:: mConeRadius 0.22 mSeedThreshold 3.8 mClusterThreshold 5 mEtaLow 0.95 mEtaHigh 2.1 mMaxVertex 120
StEemcGammaFilter:: mCalDepth 279.5 mMinPartEnergy 1e-05 mHadronScale 1 mFilterMode 0 mPrintLevel 1
StChain:INFO - Init() : Using gamma filter on the EEMC
StChain:INFO - Init() : EEMC Sampling Fraction = 0.05
StChain:INFO - Init() : Seed energy threshold = 3.3 GeV
StChain:INFO - Init() : Cluster eT threshold = 4.5 GeV
StChain:INFO - Init() : Maximum vertex = +/- 120 cm
StChain:INFO - Init() : Running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events in BFC
BFC:INFO - setTowerGainSpread(): gain spread: 0; gain mean value: 1.1
runSimuGammaTreeMaker():: GammaMaker config: ConeRadius 0.7 ClusterEtThreshold 5.5 SeedEnergyThreshold 4.2 ClusterEnergyThreshold 5.5 BsmdRange 0.05237 EsmdR ange 20
StEEmcA2EMaker *EEanalysis = new StEEmcA2EMaker("mEEanalysis");
EEanalysis->threshold(3.0, 0); // tower threshold (ped+N sigma)
EEanalysis->threshold(3.0, 1); // pre1 threshold
EEanalysis->threshold(3.0, 2); // pre2 threshold
EEanalysis->threshold(3.0, 3); // post threshold
EEanalysis->threshold(3.0, 4); // smdu threshold
EEanalysis->threshold(3.0, 5); // smdv threshold
emulate L2E-gamma trigger for Run 2006 [eemc-http-mb-l2gamma:: id 137641]
Trigger conditions:
cluster Et (2x2) = 5.2GeV
seed Et = 3.7GeV
Figure 1: Fraction of accepted events
Accept rate: fract. of generated events | |||
GammaJet | |||
pt bin | Bfc Filter | Pythia Filter | L2gamma Trigger |
pt=2-3 | 0.00772455 | 0.00820359 | 0.00646707 |
pt=3-4 | 0.0276027 | 0.0311644 | 0.024726 |
pt=4-6 | 0.0924862 | 0.0980663 | 0.0743646 |
pt=6-9 | 0.165708 | 0.217554 | 0.149013 |
pt=9-15 | 0.167815 | 0.258319 | 0.162353 |
pt=15-25 | 0.123 | 0.215667 | 0.117083 |
QCD | |||
pt bin | Bfc Filter | Pythia Filter | L2gamma Trigger |
pt=2-3 | 0.000107527 | 0.000430108 | 0.000107527 |
pt=3-4 | 0.000117647 | 0.00358824 | 5.88235e-05 |
pt=4-6 | 0.00115385 | 0.0146154 | 0.000576923 |
pt=6-9 | 0.00932927 | 0.06 | 0.00530488 |
pt=9-15 | 0.0325 | 0.178846 | 0.0173077 |
pt=15-25 | 0.0860331 | 0.253802 | 0.0510331 |
False reject: fract. of generated events | |||
GammaJet | |||
pt bin | Bfc Filter | Pythia Filter | L2gamma Trigger |
pt=2-3 | 0 | 0.000359281 | 0 |
pt=3-4 | 0 | 0.000616438 | 0 |
pt=4-6 | 0.000165746 | 0.00436464 | 0 |
pt=6-9 | 0.000128755 | 0.00270386 | 0 |
pt=9-15 | 0 | 0.000210084 | 0 |
pt=15-25 | 0.00025 | 0.00025 | 0 |
QCD | |||
pt bin | Bfc Filter | Pythia Filter | L2gamma Trigger |
pt=2-3 | 0 | 0 | 0 |
pt=3-4 | 0 | 0 | 0 |
pt=4-6 | 0 | 0 | 0 |
pt=6-9 | 0 | 0.000121951 | 0 |
pt=9-15 | 0 | 0 | 0 |
pt=15-25 | 4.13223e-05 | 0.000123967 | 0 |
Accept rate: fract. of Pythia-filtered events | |||
GammaJet | |||
pt bin | Bfc Filter | Pythia Filter | L2gamma Trigger |
pt=2-3 | 0.518248 | 1 | 0.49635 |
pt=3-4 | 0.505495 | 1 | 0.525275 |
pt=4-6 | 0.663662 | 1 | 0.614648 |
pt=6-9 | 0.703492 | 1 | 0.64707 |
pt=9-15 | 0.645576 | 1 | 0.605563 |
pt=15-25 | 0.568006 | 1 | 0.510433 |
QCD | |||
pt bin | Bfc Filter | Pythia Filter | L2gamma Trigger |
pt=2-3 | 0.02857 | 1 | 0 |
pt=3-4 | 0.05161 | 1 | 0 |
pt=4-6 | 0.0526316 | 1 | 0.0394737 |
pt=6-9 | 0.126016 | 1 | 0.0731707 |
pt=9-15 | 0.173118 | 1 | 0.0935484 |
pt=15-25 | 0.334419 | 1 | 0.194074 |
Description can be found here
Jet tree branches
General configuration
Triggers
pp@200GeV
addTrigger(127611); //HTTP
137611
5
127821 //HTTP-fast
137821
137822
127212 //HT2
137213
127501 //JP0
137501
127622 //JP0-etot
137622
127221 //JP1
137221
137222
137585 //bemc-jp2
127641 // eemc-http-mb-l2gamma
137641 // eemc-http-mb-l2gamma
6 // eemc-http-mb-l2gamma
Analysis cuts
StppAnaPars* anapars = new StppAnaPars();
anapars->setFlagMin(0);
anapars->setNhits(12);
anapars->setCutPtMin(0.2);
anapars->setAbsEtaMax(2.5);
anapars->setJetPtMin(3.5);
anapars->setJetEtaMax(100.0);
anapars->setJetEtaMin(0);
anapars->setJetNmin(0);
Cone finder configuration
StConePars* cpars = new StConePars();
cpars->setGridSpacing(105, -3.0, 3.0, 120,-TMath::Pi(),TMath::Pi());
cpars->setConeRadius(0.7);
cpars->setSeedEtMin(0.5);
cpars->setAssocEtMin(0.1);
cpars->setSplitFraction(0.5);
cpars->setPerformMinimization(true);
cpars->setAddMidpoints(true);
cpars->setRequireStableMidpoints(true);
cpars->setDoSplitMerge(true);
cpars->setDebug(false);
jetMaker->addAnalyzer(anapars, cpars, bet4pMaker, "ConeJets12");
anapars->setNhits(5);
jetMaker->addAnalyzer(anapars, cpars, bet4pMaker, "ConeJets5");
Disk space required
needs to be updated
Jet tree branches
General configuration
Triggers
pp@200GeV
240410 // JP1 // lum: 0.246
240411 // JP1 // lum: 4.045
240650 // L2JetHigh // lum: 3.745
240651 // L2JetHigh // lum: 1.811
240652 // L2JetHigh // lum: 19.769
240620 // L2BGamma // lum: 23.004
240630 // L2EGamma // lum: 3.969
240631 // L2EGamma // lum: 21.592
240013 // BBCMB-Cat2 (luminosity monitor)
240113 // BBCMB-Cat2 (luminosity monitor)
240123 // BBCMB-Cat2 (luminosity monitor)
240223 // BBCMB-Cat2 (luminosity monitor)
pp@500GeV
230410 // JP1 // lum: 0.198
230411 // JP2 // lum: 8.089
230630 // L2EGamma // lum: 3.347
230013 // BBCMB-Cat2 (luminosity monitor)
Analysis cuts
StppAnaPars* anapars = new StppAnaPars;
anapars->setFlagMin(0); // track->flag() > 0
anapars->setCutPtMin(0.2); // track->pt() > 0.2
anapars->setAbsEtaMax(2.5); // abs(track->eta()) < 2.5
anapars->setJetPtMin(5.0);
anapars->setJetEtaMax(100.0);
anapars->setJetEtaMin(0);
anapars->setJetNmin(0);
Cone finder configuration
StConePars* cpars = new StConePars;
cpars->setGridSpacing(105,-3.
cpars->setSeedEtMin(0.5);
cpars->setAssocEtMin(0.1);
cpars->setSplitFraction(0.5);
cpars->setPerformMinimization(
cpars->setAddMidpoints(true);
cpars->
cpars->setDoSplitMerge(true);
cpars->setDebug(false);
Disk space required
1Tb (see links at the top of the page for mode details)
Comments:
Figure 1: Reconstructed gamma-jet candidate yield vs. photon candidate pt in the endcap
(di-jet events found with the jet finder) Prompt photon filtered Monte-Carlo, partonic pt 3-25 GeV
(a)
(b) Same as (a) on a log scale
Figure 2: Same as Fig. 1 for QCD two processes filtered Monte-Carlo, partonic pt 6-9 GeV
(a)
Same on a log scale
(b) Same as (a) on a log scale
1.4M events, 6.8 CPU years, 0.29Tb disk space
parton pt, GeV | Pythia acc | bfc acc wrt. Pythia | Total filter's acc | Sigma, pb | lumi, 1/pb | Number of filtered events to generate | Total CPU time, days | disk space, Gb | starsim CPU, sec | bfc CPU, sec | Total Starsim CPU time, days | Total bfc CPU time, days | bfc acc |
g-jets | |||||||||||||
2-3 | 0.00820 | 0.2663 | 0.00218 | 1304000 | 10.0 | 28491 | 14.60576 | 4.45 | 10.72 | 4.0 | 13.270 | 1.34 | 0.00288 |
3-4 | 0.03116 | 0.3787 | 0.01180 | 293300 | 10.0 | 34611 | 14.85782 | 5.63 | 12.19 | 4.9 | 12.901 | 1.96 | 0.01504 |
4-6 | 0.09807 | 0.5521 | 0.05414 | 126300 | 10.0 | 68382 | 20.99942 | 11.34 | 12.04 | 4.7 | 17.254 | 3.75 | 0.06824 |
6-9 | 0.21755 | 0.6688 | 0.14550 | 26090 | 10.0 | 37961 | 7.79574 | 6.48 | 9.56 | 3.4 | 6.280 | 1.52 | 0.15848 |
9-15 | 0.25832 | 0.6274 | 0.16207 | 4675 | 10.0 | 7577 | 1.95231 | 1.31 | 11.56 | 3.8 | 1.616 | 0.34 | 0.15840 |
15-25 | 0.21567 | 0.5394 | 0.11633 | 326 | 10.0 | 379 | 0.13531 | 0.07 | 14.34 | 4.2 | 0.117 | 0.02 | 0.12112 |
totals: | 1754691 | 177401 | 60.3 | 29.29 | |||||||||
0.18 | 0.16533 | years | |||||||||||
QCD | |||||||||||||
2-3 | 0.00043 | 0.0185 | 0.00001 | 8226000000 | 0.1 | 6545 | 68.81525 | 1.30 | 16.62 | 9.9 | 68.064 | 0.75 | 0.00002 |
3-4 | 0.00359 | 0.0268 | 0.00010 | 1295000000 | 0.2 | 24907 | 140.98015 | 5.32 | 12.88 | 8.3 | 138.594 | 2.39 | 0.00006 |
4-6 | 0.01462 | 0.0450 | 0.00066 | 440300000 | 1.0 | 289582 | 925.73023 | 57.50 | 12.07 | 7.9 | 899.343 | 26.39 | 0.00072 |
6-9 | 0.06000 | 0.0744 | 0.00446 | 57830000 | 2.0 | 516306 | 915.77947 | 111.13 | 10.94 | 6.2 | 878.663 | 37.12 | 0.00564 |
9-15 | 0.17885 | 0.1354 | 0.02422 | 7629000 | 2.0 | 369484 | 360.48698 | 77.18 | 10.72 | 5.1 | 338.666 | 21.82 | 0.02420 |
15-25 | 0.25380 | 0.2754 | 0.06990 | 381900 | 1.0 | 26694 | 17.99424 | 5.57 | 14.43 | 5.9 | 16.185 | 1.81 | 0.07368 |
totals: | 10027140900 | 1233518 | 2429.78632 | 258.01 | |||||||||
1.23 | 6.65695 | years | |||||||||||
number of events, x 10e6 | CPU years | disk space, Gb | total time with 50 CPU, weeks | total time with 100 CPU, weeks | |||||||||
1.41 | 6.8 | 287.3 | 7.1 | 3.6 |
2M events, 6.9 CPU years, 0.41Tb disk space
parton pt, GeV | Pythia acc | bfc acc wrt. Pythia | Total filter's acc | Sigma, pb | lumi, 1/pb | Number of filtered events to generate | Total CPU time, days | disk space, Gb | starsim CPU, sec | bfc CPU, sec | Total Starsim CPU time, days | Total bfc CPU time, days | bfc acc |
g-jets | |||||||||||||
2-3 | 0.00820 | 0.51825 | 0.00425 | 1304000 | 10.0 | 55439 | 15.74249 | 8.66 | 10.72 | 3.9 | 13.270 | 2.47 | 0.00772 |
3-4 | 0.03116 | 0.50550 | 0.01575 | 293300 | 10.0 | 46205 | 15.47012 | 7.52 | 12.19 | 4.8 | 12.901 | 2.57 | 0.02760 |
4-6 | 0.09807 | 0.66366 | 0.06508 | 126300 | 10.0 | 82200 | 21.71983 | 13.64 | 12.04 | 4.7 | 17.254 | 4.47 | 0.09249 |
6-9 | 0.21755 | 0.70349 | 0.15305 | 26090 | 10.0 | 39930 | 7.87117 | 6.82 | 9.56 | 3.4 | 6.280 | 1.59 | 0.16571 |
9-15 | 0.25832 | 0.64558 | 0.16676 | 4675 | 10.0 | 7796 | 1.96162 | 1.35 | 11.56 | 3.8 | 1.616 | 0.35 | 0.16782 |
15-25 | 0.21567 | 0.56801 | 0.12250 | 326 | 10.0 | 399 | 0.13625 | 0.07 | 14.34 | 4.2 | 0.117 | 0.02 | 0.123 |
totals: | 1754691 | 231970 | 62.9 | 38.05 | |||||||||
0.23 | 0.17233 | years | |||||||||||
QCD | |||||||||||||
2-3 | 0.00043 | 0.02857 | 0.00001 | 8226000000 | 0.1 | 10108 | 68.99113 | 2.01 | 16.62 | 7.9 | 68.064 | 0.93 | 0.00011 |
3-4 | 0.00359 | 0.05161 | 0.00019 | 1295000000 | 0.2 | 47964 | 142.14305 | 10.25 | 12.88 | 6.4 | 138.594 | 3.55 | 0.00012 |
4-6 | 0.01462 | 0.05263 | 0.00077 | 440300000 | 1.0 | 338693 | 928.66220 | 67.25 | 12.07 | 7.5 | 899.343 | 29.32 | 0.00115 |
6-9 | 0.06000 | 0.12602 | 0.00756 | 57830000 | 2.0 | 874501 | 935.22936 | 188.22 | 10.94 | 5.6 | 878.663 | 56.57 | 0.00933 |
9-15 | 0.17885 | 0.17312 | 0.03096 | 7629000 | 2.0 | 472410 | 365.62683 | 98.68 | 10.72 | 4.9 | 338.666 | 26.96 | 0.0325 |
15-25 | 0.25380 | 0.33442 | 0.08488 | 381900 | 1.0 | 32414 | 18.35070 | 6.77 | 14.43 | 5.8 | 16.185 | 2.17 | 0.08603 |
totals: | 10027140900 | 1776090 | 2459.00326 | 373.18 | |||||||||
1.78 | 6.73700 | years | |||||||||||
number of events, x 10e6 | CPU years | disk space, Gb | total time with 50 CPU, weeks | total time with 100 CPU, weeks | |||||||||
2.01 | 6.9 | 411.2 | 7.2 | 3.6 |
1.4M events, 18.3 CPU years, 0.29Tb disk space
parton pt, GeV | Pythia acc | bfc acc wrt. Pythia | Total filter's acc | Sigma, pb | lumi, 1/pb | Number of filtered events to generate | Total CPU time, days | disk space, Gb | starsim CPU, sec | bfc CPU, sec | bfc reject CPU / event | Total Starsim CPU time, days | Total bfc CPU time, days | Bfc accept CPU / event |
g-jets | ||||||||||||||
2-3 | 0.00820 | 0.2663 | 0.00218 | 1304000 | 10.0 | 28491 | 41.33076 | 4.45 | 10.72 | 85.1 | 20 | 13.270 | 28.06 | 30.00 |
3-4 | 0.03116 | 0.3787 | 0.01180 | 293300 | 10.0 | 34611 | 38.06518 | 5.63 | 12.19 | 62.8 | 20 | 12.901 | 25.16 | 30.00 |
4-6 | 0.09807 | 0.5521 | 0.05414 | 126300 | 10.0 | 68382 | 53.83970 | 11.34 | 12.04 | 46.2 | 20 | 17.254 | 36.59 | 30.00 |
6-9 | 0.21755 | 0.6688 | 0.14550 | 26090 | 10.0 | 37961 | 23.81285 | 6.48 | 9.56 | 39.9 | 20 | 6.280 | 17.53 | 30.00 |
9-15 | 0.25832 | 0.6274 | 0.16207 | 4675 | 10.0 | 7577 | 5.28863 | 1.31 | 11.56 | 41.9 | 20 | 1.616 | 3.67 | 30.00 |
15-25 | 0.21567 | 0.5394 | 0.11633 | 326 | 10.0 | 379 | 0.32331 | 0.07 | 14.34 | 47.1 | 20 | 0.117 | 0.21 | 30.00 |
totals: | 1754691 | 177401 | 162.7 | 29.29 | ||||||||||
0.18 | 0.44565 | years | ||||||||||||
QCD | ||||||||||||||
2-3 | 0.00043 | 0.0185 | 0.00001 | 8226000000 | 0.1 | 6545 | 150.72171 | 1.30 | 16.62 | 1091.1 | 20 | 68.064 | 82.66 | 30.00 |
3-4 | 0.00359 | 0.0268 | 0.00010 | 1295000000 | 0.2 | 24907 | 356.60524 | 5.32 | 12.88 | 756.3 | 20 | 138.594 | 218.01 | 30.00 |
4-6 | 0.01462 | 0.0450 | 0.00066 | 440300000 | 1.0 | 289582 | 2422.48046 | 57.50 | 12.07 | 454.4 | 20 | 899.343 | 1523.14 | 30.00 |
6-9 | 0.06000 | 0.0744 | 0.00446 | 57830000 | 2.0 | 516306 | 2544.80915 | 111.13 | 10.94 | 278.8 | 20 | 878.663 | 1666.15 | 30.00 |
9-15 | 0.17885 | 0.1354 | 0.02422 | 7629000 | 2.0 | 369484 | 1013.10425 | 77.18 | 10.72 | 157.7 | 20 | 338.666 | 674.44 | 30.00 |
15-25 | 0.25380 | 0.2754 | 0.06990 | 381900 | 1.0 | 26694 | 41.71136 | 5.57 | 14.43 | 82.6 | 20 | 16.185 | 25.53 | 30.00 |
totals: | 10027140900 | 1233518 | 6529.43217 | 258.01 | ||||||||||
1.23 | 17.88886 | years | ||||||||||||
number of events, x 10e6 | CPU years | disk space, Gb | total time with 50 CPU, weeks | total time with 100 CPU, weeks | ||||||||||
1.41 | 18.3 | 287.3 | 19.1 | 9.6 |