05 May

May 2010 posts

 

2010.05.03 Pythia/BFC gamma-filter accaptance vs. gamma candidate pt (after gamma-maker 3x3 cluser fix)

Related inks:

Number of generated events per partnic pt bin:
gamma-jet: 25K
QCD(2-4): 50K
QCD(4-55): 25K

Pythia filter configuration

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

BFC filter configuration

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

StEEmcSlowMaker configuration

BFC:INFO - setTowerGainSpread(): gain spread: 0; gain mean value: 1.1

GammaMaker configuration

runSimuGammaTreeMaker():: GammaMaker config: ConeRadius 0.7 ClusterEtThreshold 5.5 SeedEnergyThreshold 4.2 ClusterEnergyThreshold 5.5 BsmdRange 0.05237 EsmdR ange 20

A2Emaker configuration

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

Trigger configuration

emulate L2E-gamma trigger for Run 2006 [eemc-http-mb-l2gamma:: id 137641]
Trigger conditions:
cluster Et (2x2) = 5.2GeV
seed Et = 3.7GeV

Accept/Reject relative to the total number of offline selected events

Definition: offline selected events are events which satisfy to the following conditions:

  • Online condition (L2E-gamma trigger fired)
  • Reconstructed vetrex (|v_z|<120cm)
  • Offline condition (at least one gammaMaker candidate found)

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

Yields for each of partonic pt bins separately

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

Figure 7: Same ad Fig. 1b for partonic pt=15-25

Figure 8: Same ad Fig. 1b for partonic pt=25-55

2010.05.05 Starsim/bfc timing tests

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
 

2010.05.13 BFC and Pythia QA after Pythia Eta -> - Eta

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

Pythia filter configuration

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

BFC filter configuration

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

StEEmcSlowMaker configuration

BFC:INFO - setTowerGainSpread(): gain spread: 0; gain mean value: 1.1

GammaMaker configuration

runSimuGammaTreeMaker():: GammaMaker config: ConeRadius 0.7 ClusterEtThreshold 5.5 SeedEnergyThreshold 4.2 ClusterEnergyThreshold 5.5 BsmdRange 0.05237 EsmdR ange 20

A2Emaker configuration

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

Trigger configuration

emulate L2E-gamma trigger for Run 2006 [eemc-http-mb-l2gamma:: id 137641]
Trigger conditions:
cluster Et (2x2) = 5.2GeV
seed Et = 3.7GeV

Accept/Reject relative to the total number of Pythia generated events

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

Figure 2: False rejection

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/Reject relative to the number of Pythia filter accepted events

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

2010.05.27 Run 6 and Run 9 jet Tree production

Useful inks:


New jet tree format

Description can be found here


Jet finder Run 6 configuration

Jet tree branches

  • 12-point branch
  • 5-point branch (EEMC jets only)

General configuration

  • cone radius = 0.7
  • track/tower pT > 0.2 GeV
  • trigger IDs: See the list below
  • all primary vertices with rank > 0
  • hadronic correction: 100% subtraction scheme
  • track |DCA| < 3 cm
  • pT-dependent DCAxy cut
  • track (number of hits) / (number of possible hits) > 0.51
  • track last point radius > 125 cm (ensure we get at least one point in TPC outer sector)
  • tower E > 0
  • tower status = 1
  • tower ADC - pedestal > 3 * RMS

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 finder Run 9 configuration:

Jet tree branches

  • 12-point branch
  • Tower-only branch (jets without tracking)
  • 5-point branch (EEMC jets only)

General configuration

  • cone radius = 0.7
  • track/tower pT > 0.2 GeV
  • trigger IDs: JP1, L2JetHigh, BBCMB-Cat2 (luminosity monitor)
  • all primary vertices with rank > 0
  • hadronic correction: 100% subtraction scheme
  • track |DCA| < 3 cm
  • pT-dependent DCAxy cut (a la Run 6 with slight tuning around 0.5 < pT < 1.5 GeV)
  • track chi^2 < 4
  • track (number of hits) / (number of possible hits) > 0.51
  • track last point radius > 125 cm (ensure we get at least one point in TPC outer sector)
  • tower E > 0
  • tower status = 1
  • tower ADC - pedestal > 3 * RMS

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.0,3.0,120,-TMath::Pi(),TMath::Pi());
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);

Disk space required

1Tb (see links at the top of the page for mode details)

2010.05.28 Photon-jet candidates - Perugia vs. Tune A

Comments:

  • See no difference in yilds for different tunes for prompt photons.
  • ~ 20% shape differenbce for QCD background Monte-Carlo

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

Figure 3: QCD ratio Tune-A/Perugia

EEMC simulation spreadsheet: prompt photons and QCD

Related links:

Same thresholds for Pythia/BFC

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            

Lowered BFC threshold that Pythia

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            

 

(very rough) Run 9 estimates (with 20 sec reject and 30 sec accept bfc time)

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