01 Jan

January 2010 posts

2010.01.04 y2006 vs y2009 EEMC pre-shower migration

Monte-Carlo setup:

Throwing one photon per event
Full STAR y2006h/y2009a
Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and energy (5-35 GeV)
Using A2Emaker to get reconstructed Tower/SMD energy (no EEMC SlowSimulator in chain)
Vertex z=0
~50K/per particle type
Non-zero energy: 3 sigma above pedestal

Geometry configurations:

y6h:v6.2:10KeV y2006h after fixing TPC/EEMC overlap (CVS:2009/12/17)
y6h:10KeV y2006h Mising TPC material bug is resolved (CVS:2009/12/29)
y9a:10KeV y2009a after Mising TPC material bug is resolved (CVS:2009/12/29)

Note: results are with CVS before "15Deg rotated volume" bug being fixed

Figure 1:Pre-shower migration: y2006h (red - CVS:2009/12/17) vs. y2006h (black - CVS:2009/12/29)

Figure 2: Pre-shower migration: y2006h (red) vs. y2009a (black) all with CVS:2009/12/29

2010.01.07 EEMC response to single photons with y2006h vs y2009a geometries

Monte-Carlo setup

Throwing one photon per event
Full STAR y2006h/y2009a configurations
Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and energy (5-35 GeV)
Using A2Emaker to get reconstructed Tower/SMD energy (no EEMC SlowSimulator in chain)
Vertex z=0
~50K/per particle type
Non-zero energy: 3 sigma above pedestal

Geometry configurations:

y6h:10KeV (black) - y2006h with emc_10KeV
y9a:10KeV (red) - y2009a with emc_10KeV

STAR geometry includes the latest "15Deg rotated volume" bug bug fix

Figure 1: EEMC sampling fraction
(left) vs. thrown photon energy (with 1.2 < eta < 1.9 cut)
(right) vs. thrown photon eta

Figure 2: 2x1/3x3 clustering

Figure 3: Shower shapes (u plane)

Figure 4: Shower shape ratios (u plane)

Figure 5: Pre-shower migration (1.2 < eta < 1.9)

Figure 6: Average pre-shower1 energy vs. thown photon position in EEMC
(left) y2009a with emc_10KeV
(right) y2006h with emc_10KeV

2010.01.08 EEMC response 2006 vs. 2009: phi cuts

EEMC migration plots with cuts on TPC sector boundaries

Click here for results before phi cuts

Monte-Carlo setup

Throwing one photon per event
Full STAR y2006h/y2009a configurations
Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and energy (5-35 GeV)
Using A2Emaker to get reconstructed Tower/SMD energy (no EEMC SlowSimulator in chain)
Vertex z=0
~50K/per particle type
Non-zero energy: 3 sigma above pedestal

Geometry configurations:

y6h:10KeV (black) - y2006h with emc_10KeV option
y9a:10KeV (red) - y2009a with emc_10KeV option

Figure 1: Average pre-shower1 energy vs. thown photon position in EEMC
with cuts on TPC sector boundaries: cos(12*(phi-Pi/6.)) < -0.65 (similar plot before phi cuts)
(left) y2009a with emc_10KeV
(right) y2006h with emc_10KeV

Figure 2: Pre-shower 1 sampling fraction (E_pre1/E_thrown) vs. thrown eta

Figure 3: EEMC sampling fraction
(left) vs. thrown photon energy (with 1.2 < eta < 1.9 cut)
(right) vs. thrown photon eta

Figure 4: Pre-shower migration (1.2 < eta < 1.9)

2010.01.12 W test sample QA

All plost from second (with vertex distribution) test W-sample from Lidia/Jason.
generated files are from /star/rcf/test/Wprod_test2/

The previous sample with fixed (zero) vertex can be found was announced here:
http://www.star.bnl.gov/HyperNews-star/protected/get/starsimu/435.html

Figure 1: Electron from W decay (a) eta, (b) phi, (c) pt and (d) energy distributions
from geant record (no kinematic cuts)

Figure 2: Reconstructed vs. geant vertex (Cuts: abs(geant_eta_electron) < 1)
(left) difference, (right) ratio

Figure 3:
(left) Correlation between thrown and reconsructed energy: abs(geant_eta_electron) < 1
(right) ratio of the reconsructed to thrown energy (Bemc_Etow > 25)
Reconstructed energy is the total energy in all Barrel towers

Data base setup

The follwoing DB tables are used to read MuDst (dbMk->SetDateTime(20090325,0)):

StEmcSimulatorMaker:INFO - loaded a new bemcPed table with beginTime 2009-03-24 22:16:13 and endTime 2009-03-26 06:03:44
StEmcSimulatorMaker:INFO - loaded a new bemcStatus table with beginTime 2009-03-24 02:16:58 and endTime 2009-03-26 04:07:02
StEmcSimulatorMaker:INFO - loaded a new bemcCalib table with beginTime 2008-12-15 00:00:02 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bemcGain table with beginTime 1999-01-01 00:08:00 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bprsPed table with beginTime 2008-03-04 10:30:56 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bprsStatus table with beginTime 2008-12-15 00:00:00 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bprsCalib table with beginTime 1999-01-01 00:10:00 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bprsGain table with beginTime 1999-01-01 00:08:00 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bsmdePed table with beginTime 2009-03-24 15:42:29 and endTime 2009-03-25 11:24:55
StEmcSimulatorMaker:INFO - loaded a new bsmdeStatus table with beginTime 2009-03-24 15:42:29 and endTime 2009-03-25 11:24:55
StEmcSimulatorMaker:INFO - loaded a new bsmdeCalib table with beginTime 2002-11-14 00:01:00 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bsmdeGain table with beginTime 1999-01-01 00:08:00 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bsmdpPed table with beginTime 2009-03-24 15:42:29 and endTime 2009-03-25 11:24:55
StEmcSimulatorMaker:INFO - loaded a new bsmdpStatus table with beginTime 2009-03-24 15:42:29 and endTime 2009-03-25 11:24:55
StEmcSimulatorMaker:INFO - loaded a new bsmdpCalib table with beginTime 2002-11-14 00:01:00 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bsmdpGain table with beginTime 1999-01-01 00:08:00 and endTime 2037-12-31 12:00:00
StEmcSimulatorMaker:INFO - loaded a new bemcTriggerStatus table with beginTime 2009-03-23 07:50:04 and endTime 2009-04-01 18:10:03
StEmcSimulatorMaker:INFO - loaded a new bemcTriggerPed table with beginTime 2009-03-20 04:11:43 and endTime 2009-03-30 20:00:05
StEmcSimulatorMaker:INFO - loaded a new bemcTriggerLUT table with beginTime 2009-03-23 07:50:04 and endTime 2009-04-03 22:08:11

2010.01.13 W test sample QA: Pass 2

http://drupal.star.bnl.gov/STAR/node/16704

QA of the test W-sample from Lidia/Jason.
generated MuDst are from /star/simu/jwebb/01-11-2010-w-test-production/

QA plots for the previous pass can be found here

Two channels being analyzed:

wtest10000 W+ --> e+ nu (shown by black line)
wtest10001 W- --> e- nu (shown by red line)

Cuts: |geant_eta_lepton| < 1

Discussions can be found here:
http://www.star.bnl.gov/HyperNews-star/protected/get/starsimu/440.html

Figure 1: (left) Reconstructed vertex z distribution
(right) reconstructed minus geant z-vertex

Figure 2: E2x2/E_geant energy ratio
Black: positron from W+, mean value= 0.972973;
Red - electron from W- mean value = 0.969773

Figure 3: E1x1/E_geant (highest tower) energy ratio
Black: positron from W+, mean value= 0.815287;
Red - electron from W- mean value = 0.812098

Update on Jan 14, 2010

Figure 4: Lepton E2x2/E_geant energy ratio

Parameter	black: positron from W+	red: electron from W-
gaus-Constant	1.60709e+01 , err=3.08565e+00	1.56834e+01 , err=4.71967e+00
gaus-Mean	9.85514e-01 , err=4.94309e-03	9.86118e-01 , err=5.43577e-03
gaus-Sigma	3.15205e-02 , err=3.73952e-03	2.52009e-02 , err=6.57793e-03
Hist-Mean	0.972973	0.969773

Figure 5: Lepton E3x3/E_geant energy ratio

Parameter	black: positron from W+	red: electron from W-
gaus-Constant	1.48719e+01 , err=2.82186e+00	1.35741e+01 , err=3.36776e+00
gaus-Mean	9.89924e-01 , err=5.72959e-03	9.83056e-01 , err=6.28736e-03
gaus-Sigma	3.37758e-02 , err=4.24983e-03	3.00597e-02 , err=6.06841e-03
Hist-Mean	0.975662	0.974163

2010.01.15 W test sample QA: Pass 3

QA of the test W-sample from Lidia/Jason.
generated MuDst are from /star/data08/users/starreco/recowtest/

QA plots for the previous pass 2 can be found here

QA plots for the previous pass 1 can be found here

Two channels being analyzed:

wtest10000 W+ --> e+ nu (shown by black line)
wtest10001 W- --> e- nu (shown by red line)

Cuts: |geant_eta_lepton| < 1

Discussions can be found here:
http://www.star.bnl.gov/HyperNews-star/protected/get/starsimu/443.html

Figure 1: Reconstructed minus geant z-vertex

Figure 2: Lepton E2x2/E_geant energy ratio

Parameter	black: positron from W+	red: electron from W-
gaus-Constant	6.24023e+01 , err=3.46979e+00	4.73536e+01 , err=3.16029e+00
gaus-Mean	9.79982e-01 , err=1.69854e-03	9.79787e-01 , err=1.68813e-03
gaus-Sigma	3.52892e-02 , err=1.40963e-03	3.15336e-02 , err=1.40759e-03
Hist-Mean	0.972122	0.975073

Figure 3: Lepton E3x3/E_geant energy ratio

Parameter	black: positron from W+	red: electron from W-
gaus-Constant	6.33596e+01 , err=3.58862e+00	4.72335e+01 , err=3.19552e+00
gaus-Mean	9.83287e-01 , err=1.72276e-03	9.83632e-01 , err=1.67661e-03
gaus-Sigma	3.45514e-02 , err=1.48019e-03	3.05224e-02 , err=1.38944e-03
Hist-Mean	0.974372	0.977858

2010.01.18 EEMC response 2006 vs. 2009: typo in TPC volume size bug fix

Click here for previous study before TPC typo fix

Monte-Carlo setup

Throwing one photon per event
Full STAR y2006h/y2009a configurations
Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and energy (5-35 GeV)
Using A2Emaker to get reconstructed Tower/SMD energy (no EEMC SlowSimulator in chain)
Vertex z=0
~50K/per particle type
Non-zero energy: 3 sigma above pedestal

Before and after the fix comparison

Geometry configurations:

y6h:10KeV:old (red) - y2006h with emc_10KeV option
y6h:10KeV (black) - y2006h with emc_10KeV option, after TPC typo fixed

Figure 1: Pre-shower migration (1.2 < eta < 1.9)

y2006h vs. y2009a comparison after a TPC typo fix

Geometry configurations:

y9a:10KeV (red) - y2009a with emc_10KeV option
y6h:10KeV (black) - y2006h with emc_10KeV option

Figure 2: Pre-shower migration (1.2 < eta < 1.9): y2006h vs. y2009a

2010.01.18 W test sample: Cluster ratios and skewed gaussian fits

W test sample from Lidia/Jason. MuDst's from /star/data08/users/starreco/recowtest/

Two channels being analyzed:

wtest10000 W+ -> e+ nu
wtest10001 W- -> e- nu

Figure 1: Lepton yield vs. rapidity (no cuts)

Figure 2: Lepton yield vs. pt and energy
(left) no rapidity cuts
(right) |lepton_eta| < 1

Cluster energy vs. original lepton energy

All plots below with |lepton_eta| < 1

Skewed gaussian fits: [const]*exp(-0.5*((x-[mean])/([sigma]*(1+[skewness]*(x-[mean]))))**2)

Figure 3: Lepton E1x1/E_geant energy ratio

Figure 4: Lepton E2x2/E_geant energy ratio

Figure 5: Lepton E3x3/E_geant energy ratio

2010.01.26 Endcap/Barrel clustering with official W-MC

Simulations: official pp 500GeV pythia W production

Two channels being analyzed:

W+ -> e+ nu (rcf10010*.root)
W- -> e- nu (rcf10011*.root)

Lepton from W in the Endcap: 1.2 < eta_lepton < 1.9

Figure 1: Lepton yield vs. energy

Figure 2: Lepton (left) E1x1/E_thrown and (right) E2x2/E_thrown energy ratio
Skewed gaussian fits: [const]*exp(-0.5*((x-[mean])/([sigma]*(1+[skewness]*(x-[mean]))))**2)

Figure 3: Endcap EMC lepton E3x3/E_thrown energy ratio

Figure 4: Endcap 2x2 sampling fraction (s.f.) vs. thrown lepton (left) energy and (right) eta
S.f. is defined as an average E_2x2/E_thrown for E_2x2/E_thrown>0.8

Figure 5: Endcap 3x3 sampling fraction (s.f.) vs. thrown lepton (left) energy and (right) eta
S.f. is defined as an average E_3x3/E_thrown for E_3x3/E_thrown>0.8

Lepton from W in the Barrel: |eta_lepton| <1

Figure 6: Lepton yield vs. energy

Figure 7: Lepton (left) E1x1/E_thrown and (right) E2x2/E_thrown energy ratio

Figure 8: Barrel EMC lepton E3x3/E_thrown energy ratio

Figure 9: Barrel 2x2 s.f. vs. thrown lepton (left) energy and (right) eta
S.f. is defined as an average E_2x2/E_thrown for E_2x2/E_thrown>0.8

Figure 10: Barrel 3x3 s.f. vs. thrown lepton (left) energy and (right) eta
S.f. is defined as an average E_3x3/E_thrown for E_3x3/E_thrown>0.8

02 Feb

February 2010 posts

2010.02.08 StEemcGammaFilterMaker QA: QCD vs. gamma-jet ccept/reject rates

StEemcGammaFilterMaker QA

Pythia generated processes

Pythia gamma-jet	Pythia QCD 2->2 processes
14 f + fbar -> g + gamma	11 f + f' -> f + f' (QCD)
18 f + fbar -> gamma + gamma	12 f + fbar -> f' + fbar'
29 f + g -> f + gamma	13 f + fbar -> g + g
114 g + g -> gamma + gamma	28 f + g -> f + g
115 g + g -> g + gamma	53 g + g -> f + fbar
	68 g + g -> g + g

Number of generated events per parton pt bin

	Number of generated events (100events/job)
Parton pt range (GeV)	2-3	3-4	4-6	6-9	9-15	15-25
Pythia gamma-jet	50K	50K	50K	50K	50K	50K
Pythia QCD 2->2 processes	100K	100K	50K	50K	50K	50K

Filter configuration

Filter parameter	Value	Notes
mConeRadius	0.24
mSeedThreshold	2.5	Cluster seed energy threshold
mClusterThreshold	3.7	Cluster Et threshold
mEtaLow	0.9	EEMC acceptance
mEtaHigh	2.1	EEMC acceptance
mSmearEnergy	0	Disabled
mThrowTracks	0	Disabled
mCalDepth	279.5	ZDC SMD depth
mMinPartEnergy	1e-05	Disabled by mThrowTracks=0
mHadronScale	0.4	Downscale factor for hadron energy
mFilterMode	0	Accepting all events

stasim/Pythia options

detp geometry y2006h
Calorimeter cut for electromagnetic processes: emc_10keV
call pytune(100): PYTUNE v1.015; CDF Tune A

Figures

Figure 1: Pythia Eemc gamma filter QA:
(left) False rejection, (right) fraction of accepted events

Comments

Overall false rejection is < 0.02% (10^-4)
for both QCD and gamma-jet simulations.
In parton pt range 6-15 GeV the acceptance rate
for the gamma-jet MC is ~ 50%
(I think this is due to rapidity fiducial cut,
otherwise it should be closer to 100%).
For the QCD sample (in the same pt region)
acceptance rate is 8-25%.
Somehow for highest parton pt filter acceptance is the
same for gamma-jet and QCD Monte-Carlo.

2010.02.09 BFC and Pythia QA: relative to trigger

QCD and gamma-jet data samples are described here

Pythia filter configuration

StEemcGammaFilter:: running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events
StEemcGammaFilter:: mConeRadius 0.22 mSeedThreshold 2.1 mClusterThreshold 3.25 mEtaLow 0.95 mEtaHigh 2.1
StEemcGammaFilter:: mCalDepth 279.5 mMinPartEnergy 1e-05 mHadronScale 0.4 mFilterMode 0 mPrintLevel 0

BFC filter configuration

StChain:INFO - Init() : Seed energy threshold = 2.8 GeV
StChain:INFO - Init() : Cluster eT threshold = 4.2 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

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

Figure 1: Fraction of accepted events

Figure 2: False rejection (Y-axis scale is 10^-3)

Accept/Reject relative to the number of triggered events

Figure 3: Fraction of accepted events (relative to triggered events)

Figure 4: False rejection relative to triggered events

2010.02.10 Money plots for W cross section

2010.02.11 BFC and Pythia QA: Gain no-gain-spread, mean=1.05

Click here for discussion and results with spread=0.05/gain=0.95

QCD and gamma-jet data samples are described here

Resultys without gain shift can be found here
(Note: ignore parton pt=25-35GeV for the gamma-jet sample since all jobs failed)

Pythia filter configuration

StEemcGammaFilter:: running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events
StEemcGammaFilter:: mConeRadius 0.22 mSeedThreshold 2.1 mClusterThreshold 3.25 mEtaLow 0.95 mEtaHigh 2.1
StEemcGammaFilter:: mCalDepth 279.5 mMinPartEnergy 1e-05 mHadronScale 0.4 mFilterMode 0 mPrintLevel 0

BFC filter configuration

StChain:INFO - Init() : Seed energy threshold = 2.8 GeV
StChain:INFO - Init() : Cluster eT threshold = 4.2 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

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.0023	0.06264	0.00148
pt=3-4	0.0242285	0.250601	0.0126854
pt=4-6	0.103111	0.427535	0.0571313
pt=6-9	0.16828	0.48368	0.13918
pt=9-15	0.1692	0.50118	0.1619
pt=15-25	0.12708	0.42904	0.11786
pt=25-35	0.05702	0.24854	0.0509
QCD
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	3.003e-05	0.00426426	2.002e-05
pt=3-4	0.0001001	0.0122923	1.001e-05
pt=4-6	0.00078	0.03166	0.00014
pt=6-9	0.00622	0.10538	0.00216
pt=9-15	0.02822	0.27666	0.01022
pt=15-25	0.07568	0.4405	0.03086
pt=25-35	0.0761	0.35556	0.04116

Figure 2: False rejection (Y-axis scale is 10^-3)

False reject: fract. of generated events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	2e-05	0	0
pt=3-4	2.00401e-05	0	0
pt=4-6	8.08081e-05	0	0
pt=6-9	6e-05	4e-05	0
pt=9-15	0.0002	0.00018	0
pt=15-25	0.0001	4e-05	0
pt=25-35	0.00018	2e-05	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	4e-05	6e-05	0
pt=9-15	4e-05	0.00026	0
pt=15-25	0.00016	0.00018	0
pt=25-35	4e-05	4e-05	0

Accept/Reject relative to the number of triggered events

Figure 3: Fraction of accepted events (relative to triggered events)

Figure 4: False rejection relative to triggered events

Accept/Reject relative to the number of Pythia filter accepted events

Figure 5: Fraction of accepted events relative to Pythia filter accepted events

Accept rate: fract. of Pythia-filtered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.0363985	1	0.0220307
pt=3-4	0.0966813	1	0.0501399
pt=4-6	0.241081	1	0.133488
pt=6-9	0.347461	1	0.287587
pt=9-15	0.336286	1	0.322639
pt=15-25	0.29573	1	0.274613
pt=25-35	0.228776	1	0.204635
QCD
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.00704225	1	0.00234742
pt=3-4	0.00488599	1	0.000814332
pt=4-6	0.0214782	1	0.00442198
pt=6-9	0.053141	1	0.0191687
pt=9-15	0.0965806	1	0.0356394
pt=15-25	0.16958	1	0.0695573
pt=25-35	0.213185	1	0.115592

Figure 6: False rejection relative to Pythia filter accepted events

False reject: fract. of Pythia-filtered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.000319285	0	0
pt=3-4	7.9968e-05	0	0
pt=4-6	0.000189009	0	0
pt=6-9	0.000124049	0	0
pt=9-15	0.000399058	0	0
pt=15-25	0.000233079	0	0
pt=25-35	0.000724229	0	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.000379579	0	0
pt=9-15	0.000144582	0	0
pt=15-25	0.000363224	0	0
pt=25-35	0.000112499	0	0

2010.02.11 BFC and Pythia QA: Gain spread=0.05, mean=0.95

QCD and gamma-jet data samples are described here

Resultys without gain spread can be found here
(Note: ignore parton pt=25-35GeV for the gamma-jet sample since all jobs failed)

Gain spread implementation in StEEmcSlowMaker.cxx (private version):

void StEEmcSlowMaker::setTowerGainSpread(Float_t s, Float_t mTowerGainMean)
{
LOG_INFO << "setTowerGainSpread(): gain spread: " << s << "; gain mean value: " << mTowerGainMean << endm;
// initialize tower gain factors to 1
for ( Int_t sec=0;sec<kEEmcNumSectors;sec++ )
    for ( Int_t sub=0;sub<kEEmcNumSubSectors;sub++ )
      for ( Int_t eta=0;eta<kEEmcNumEtas;eta++ )
    {
    //    mTowerGainFact[sec][sub][eta]=1.0;

    Float_t f = -1.0E9;
    while ( f <= -1. || f >= 1.0 )
        f = gRandom->Gaus(0., s);

    mTowerGainFact[sec][sub][eta] = mTowerGainMean + f;

    }
}

Pythia filter configuration

StEemcGammaFilter:: running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events
StEemcGammaFilter:: mConeRadius 0.22 mSeedThreshold 2.1 mClusterThreshold 3.25 mEtaLow 0.95 mEtaHigh 2.1
StEemcGammaFilter:: mCalDepth 279.5 mMinPartEnergy 1e-05 mHadronScale 0.4 mFilterMode 0 mPrintLevel 0

BFC filter configuration

StChain:INFO - Init() : Seed energy threshold = 2.8 GeV
StChain:INFO - Init() : Cluster eT threshold = 4.2 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

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

Figure 1: Fraction of accepted events

Figure 2: False rejection (Y-axis scale is 10^-3)

Accept/Reject relative to the number of triggered events

Figure 3: Fraction of accepted events (relative to triggered events)

Figure 4: False rejection relative to triggered events

2010.02.12 Final Pythia and BFC EEMC-gamma-filter paramter settings

Pythia generated processes

Prompt photons (gamma-jets)		2->2 QCD
id	Process	id	Process
14	f + fbar -> g + gamma	11	f + f' -> f + f' (QCD)
18	f + fbar -> gamma + gamma	12	f + fbar -> f' + fbar'
29	f + g -> f + gamma	13	f + fbar -> g + g
114	g + g -> gamma + gamma	28	f + g -> f + g
115	g + g -> g + gamma	53	g + g -> f + fbar
		68	g + g -> g + g

Number of generated events per parton pt bin

	Number of generated events (100events/job)
Parton pt range (GeV)	2-3	3-4	4-6	6-9	9-15	15-25	25-35
gamma-jets	50K	50K	50K	50K	50K	50K	50K
2->2 QCD processes	100K	100K	50K	50K	50K	50K	50K

Pythia Filter configuration

StRoot/StMCFilter/StEemcGammaFilter.cxx
StRoot/StMCFilter/StEemcGammaFilter.h

Filter parameter	Value	Notes
mConeRadius	0.22
mSeedThreshold	2.6	Cluster seed energy threshold (GeV)
mClusterThreshold	3.6	Cluster Et threshold (GeV)
mEtaLow	0.95	EEMC acceptance
mEtaHigh	2.1	EEMC acceptance
mMaxVertex	120.0	Vertex z cut (cm)
mCalDepth	279.5	EEMC SMD depth (cm)
mMinPartEnergy	1e-05	Ignore track with minim energy (GeV)
mHadronScale	0.4	Downscale factor for hadron energy
mFilterMode	0 / 1	0=Accept all events; 1=reject events

BFC Filter configuration

StRoot/StFilterMaker/StEemcGammaFilterMaker.cxx
StRoot/StFilterMaker/StEemcGammaFilterMaker.h

Filter parameter	Value	Notes
mSeedEnergyThreshold	3.4	Cluter seed energy threshold (GeV)
mClusterEtThreshold	4.5	Cluster Et threshold (GeV)
mEemcSamplingFraction	0.05	Assume 5% sampling fraction for EEMC
mMaxVertex	120.0	Vertex z cut (cm)
mFilterMode	0 / 1	0=Accept all events; 1=reject events

GammaMaker configuration

Filter parameter	Value	Notes
ConeRadius	0.7
ClusterEtThreshold	5.5	(GeV)
SeedEnergyThreshold	4.2	(GeV)
ClusterEnergyThreshold	5.5	(GeV)

EEMC SlowSimulator configuration

(for a moment private version) of StEEmcSlowMaker.cxx
with modified setTowerGainSpread(Float_t s, Float_t mTowerGainMean)

Filter parameter	Value	Notes
mTowerGainMean	1.05	Overall +5% gain shift
GainSpread	0	No gain spread

GSTAR/Pythia options

detp geometry y2006h
Calorimeter cut for electromagnetic processes: emc_10keV
call pytune(100): PYTUNE v1.015; CDF Tune A
or
call pytune(320): PYTUNE Perugia; Perugia 0 tune

2010.02.16 Pythia and BFC filter QA vs. gamma candidate pt and eta

QCD and gamma-jet data samples and filter configurtions are given here

Note: for this study I have used ideal EEMC gains (no gain shift/spread)

Note on trigger effect intepretation:
There is no requirement for the Pythia gamma-jet sample to have direct gamma
headed to the EEMC, only requirement is to have a gamma candidate in the EEMC,
so the away side jet may also contribute.

Figure 1: pt distribution of the gamma candidates
for Pythia/Bfc level filter and triggered events
Event cuts: at least one gamma candidate, |v_z| <120
Left: Pythia gamma-jet MC; (right) 2->2 Pythia QCD
Lower plots: fraction of accepted gamma candidates by filter/trigger
No parton pt weights (= ignore bumps in pt distribution for gamma-jet sample)

Figure 2: Rapidity distribution of the gamma candidates (Same conditions as in Fig. 1)

Figure 3: pt distribution of false rejection for Pythia/Bfc filters
Candidate cuts: at least one gamma candidate, l2gamm-trigger=fired, |v_z| <120
Most of false rejection (~ 1-3% for QCD) is for gamma candidates with pt < 8GeV

03 Mar

March 2010 posts

2010.03.02 Endcap photon-jets simulation request (draft)

Request last updated on Jul 21, 2010

Request summary

Total resources estimate for QCD with 1/pb and prompt-photon with 10/pb suimulations:

CPU: 4.2 CPU years (2.2 weeks of running on a 100 CPUs)
Disk space: 0.15Tb
Numbe of filtered events: 0.74M

partonic pt	QCD	prompt photon
2-3	0	30K
3-4	0	36K
4-6	130K	76K
6-9	240K	40K
9-15	150K	10K
15-35	23K	1K

Latest filter bias/timing test and simulation request spreasheet

Note: These and all other studies are linked from here

Filter code in cvs

Pythia gamma-filter code:
StEemcGammaFilter.cxx
StEemcGammaFilter.h
BFC gamma-filter code:
StEemcGammaFilterMaker.cxx
StEemcGammaFilterMaker.h
eemcGammaFilterMakerParams.idl

Further information related to this request

Lidia added "y2006h" tag (latest Endcap geometry fixes, and Calorimeters with LOW_EM cuts)
http://www.star.bnl.gov/HyperNews-star/protected/get/starsimu/452/1.html
x/y beam offset:
Run 6: x=0.0cm, y=-0.3cm (from /STAR/comp/calib/Beamline/Run6)
Run 9: x=0.3cm, y= 0.0cm (from /STAR/comp/calib/BeamLine/Run9)
Vertex z cut:
+/- 120cm
Vertex z spread:
    Run 6: 55cm (gaus fit to Fig.1 from this post: /STAR/node/13276)
    Run 9: 63cm are taken from Pibero's embedding study:
    www4.rcf.bnl.gov/~pibero/spin/dijets/2009.10.23/embedding.html
Vertex x/y spread set to zero for all runs.
FYI, Run 9 x/y spread is x=0.57, y=0.58
http://www4.rcf.bnl.gov/~pibero/spin/dijets/2009.10.23/XYVertexJetTriggers.png
Vertex option:
    Use option consistent with bfc tags used for data production (VFPPV/Run-6 or VFPPVnoCTB/Run-9):
       Run 6: Leave vertex to be reconstructed vertex, and use VFPPV with beamline
   Run 9: Leave vertex to be reconstructed vertex, and use VFPPVnoCTB with beamline

FYI: bfc options for different years:
http://www.star.bnl.gov/devcgi/dbProdOptionRetrv.pl
Use the latest available "SLXXy" library tag
No sdt option for bfc with Monte-Carlo. See note from Jreome's:
http://www.star.bnl.gov/HyperNews-star/protected/get/starsoft/7905/1/1/2.html
Need to add new bfc tag. Request send to starsimu list:
http://www.star.bnl.gov/HyperNews-star/protected/get/starsimu/453.html
Using Perugia0 tunes (i.e. call pytune(320))
GMT timestamp update
http://www.star.bnl.gov/HyperNews-star/protected/get/phana/481.html

Run 6 200 GeV
sdt20060512.043500 (GMT during run 7132005)
sdt20060513.064000 (GMT during run 7133011)
sdt20060514.090000 (GMT during run 7134015)
sdt20060516.152000 (GMT during run 7136022)
sdt20060518.073700 (GMT during run 7138010)
sdt20060520.142000 (GMT durign run 7140024)
sdt20060521.052000 (GMT during run 7141011)
sdt20060522.124500 (GMT during run 7142029)
sdt20060523.204400 (GMT during run 7143044)
sdt20060525.114000 (GMT during run 7145023)
sdt20060526.114000 (GMT during run 7146020)
sdt20060528.144500 (GMT during run 7148028)
sdt20060602.071500 (GMT during run 7153015)
sdt20060604.191200 (GMT during run 7155043)

Run 9 500 GeV
sdt20090320.014942
sdt20090321.095723
sdt20090324.064934
sdt20090328.040659
sdt20090329.014902
sdt20090404.194055
sdt20090407.030832
sdt20090410.020208
sdt20090411.103512
sdt20090413.021450

Run 9 200 GeV
std20090506.083000 (GMT during run 10126017)
std20090508.152000 (GMT during run 10128053)
std20090514.145500 (GMT during run 10134035)
std20090516.182500 (GMT during run 10135070)
std20090517.214000 (GMT during run 10137052)
std20090518.111600 (GMT during run 10138027)
std20090519.173200 (GMT during run 10139069)
std20090520.100500 (GMT during run 10140011)
std20090522.141000 (GMT during run 10142043)
std20090523.183500 (GMT during run 10143065)
std20090524.112000 (GMT during run 10144035)
std20090525.062000 (GMT during run 10145012)
std20090526.140000 (GMT during run 10146052)
FYI:
simulation request posted to SPIN PWG:
http://www.star.bnl.gov/HyperNews-star/protected/get/starspin/3982.html
Code status:
http://www.star.bnl.gov/HyperNews-star/protected/get/starsoft/8073.html
Code per review (by Pibero and Victor):
http://www.star.bnl.gov/HyperNews-star/protected/get/starsoft/8073/2.html
Note: code being approved

Original disk space estimate (see spreadsheed linked above for the latest estimates):

http://www.star.bnl.gov/HyperNews-star/protected/get/starspin/3982.html

------------------------ REQUEST DRAFT BELOW ----------------------------------------

Endcap photon-jets / QCD 2->2 simulations

Request TypeEvent generator simulation, with filtering

General Information


Request ID
Priority: EC	0
Priority: pwg	High
Status	New
Physics Working Group	Spin
Requested by	Photon group for SPIN PWG
Contact email(s)	ilya.selyuzhenkov@gmail.com, bridgeman@hep.anl.gov
Contact phone(s)
PWG email(s)	starspin-hn@www.star.bnl.gov
Assigned Deputy:	Not assigned
Assigned Helper:	Not assigned

Description

Endcap photon-jets request

Global Simulation Settings


Request type:	Event generator simulation, with filtering
Number of events	See list for each partonic pt bins
Magnetic Field	Run 6: Full-Field Run 9: Reversed Full-Field
Collision Type	Run 6: pp@200GeV Run 9: pp@200GeV Run 9: pp@500GeV
Centrality	---- SELECT CENTRALITY ----
BFC tags	Run 6: trs fss y2006h Idst IAna l0 tpcI fcf ftpc Tree logger ITTF Sti VFPPV bbcSim tofsim tags emcY2 EEfs evout -dstout IdTruth geantout big fzin MiniMcMk eemcDb beamLine clearmem Run 9: tpcrs y2009a MakeEvent ITTF NoSsdIt NoSvtIt Idst BAna l0 Tree logger Sti VFPPVnoCTB tpcDB TpcHitMover TpxClu bbcSim tofsim tags emcY2 EEfs evout IdTruth geantout big fzin McEvOut MiniMcMk eemcDb beamLine clearmem
Production	---- SELECT PRODUCTION TAG ----
Geometry: simu	Run 6: y2006h Run 9: y2009a
Geometry: reco	Run 6: y2006h Run 9: y2009a
Library	use library with approved filter code checked in
Vertex option	Run 6: Leave vertex to be reconstructed vertex, and use VFPPVnoCTB with beamline Run 9: Leave vertex to be reconstructed vertex, and use VFPPVnoCTB with beamline
Pileup option	No
Detector Set	Run 6: TPC, ETOW, BTOW, BSMD, ESMD, BPRS, EPRE1, EPRE2, EPOST, TOF, BBC, SVT, SSD Run 9: TPC, ETOW, BTOW, BSMD, ESMD, BPRS, EPRE1, EPRE2, EPOST, TOF, BBC

Data Sources

MC Event Generator


Event generator	Pythia
Extra options	Additional libraries required for Eemc-gamma Pythia-level filter gexec $ROOTSYS/lib/libCint.so gexec $ROOTSYS/lib/libCore.so gexec $ROOTSYS/lib/libMathCore.so gexec $ROOTSYS/lib/libMatrix.so gexec $ROOTSYS/lib/libPhysics.so gexec .sl53_gcc432/lib/StMCFilter.so // filter library Select prompt photon Pythhia processes: MSUB (14)=1 MSUB (18)=1 MSUB (29)=1 MSUB (114)=1 MSUB (115)=1 Select QCD 2->2 Pythhia processes: MSUB (11) = 1 MSUB (12) = 1 MSUB (13) = 1 MSUB (28) = 1 MSUB (53) = 1 MSUB (68) = 1 Perugia0 Pythia tune: call pytune(320)
Vertex Z, cm	-120 < Vertex < 120
Gaussian sigma in X,Y,Z if applicable	x/y spread use 0 Run 6: 0, 0, 55 200 GeV Run 9: 0, 0, 63 200 GeV Run 9: 0, 0, 42 500 GeV
Vertex offset: x, mm	Run 6: 0.0cm Run 9: 0.3cm (note: values in cm)
Vertex offset: y, mm	Run 6: -0.3cm (note: values in cm) Run 9: 0.0cm
Φ (phi), radian	0 < Φ < 6.29
η (eta)	Default (include Barrel, Endcap, BBC)
Pt bin, GeV	See list above for QCD and g-jet samples
Macro file	Pythia gamma-filter code: StEemcGammaFilter.cxx StEemcGammaFilter.h BFC gamma-filter code: StEemcGammaFilterMaker.cxx StEemcGammaFilterMaker.h eemcGammaFilterMakerParams.idl Private bfc: /star/u/seluzhen/star/spin/MCgammaFilter/scripts/bfc.C

04 Apr

April 2010 posts

2010.04.09 BFC and Pythia QA: 10% gain shift

See this post from Alice for QCD sample rates

QCD and gamma-jet data samples are described here (filter parameters are listed below)

Pythia filter configuration

StEemcGammaFilter:: running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events
StEemcGammaFilter:: mConeRadius 0.22 mSeedThreshold 2.6 mClusterThreshold 3.6 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.4 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

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

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.00646	0.09532	0.00656
pt=3-4	0.03042	0.2401	0.02772
pt=4-6	0.09438	0.42552	0.07568
pt=6-9	0.165984	0.54004	0.149137
pt=9-15	0.167329	0.559137	0.162972
pt=15-25	0.12486	0.45662	0.11744
pt=25-35	0.0562525	0.269499	0.0536273

Figure 2: False rejection (Y-axis scale is 10^-3)

False reject: fract. of generated events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	6e-05	0	0
pt=3-4	0.00016	6e-05	0
pt=4-6	0.0003	8e-05	0
pt=6-9	0.000261044	6.0241e-05	0
pt=9-15	0.000220884	6.0241e-05	0
pt=15-25	0.00028	0.00012	0
pt=25-35	0.000280561	0.000140281	0

Accept/Reject relative to the number of Pythia filter accepted events

Figure 3: Fraction of accepted events relative to Pythia filter accepted events

Accept rate: fract. of Pythia-filtered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.0677717	1	0.0612673
pt=3-4	0.126697	1	0.105623
pt=4-6	0.221799	1	0.16991
pt=6-9	0.307318	1	0.26809
pt=9-15	0.299264	1	0.281882
pt=15-25	0.273444	1	0.244317
pt=25-35	0.20873	1	0.181217

Figure 4: False rejection relative to Pythia filter accepted events

False reject: fract. of Pythia-filtered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.000629459	0	0
pt=3-4	0.000416493	0	0
pt=4-6	0.000517014	0	0
pt=6-9	0.00037183	0	0
pt=9-15	0.000287305	0	0
pt=15-25	0.000350401	0	0
pt=25-35	0.000520524	0	0

2010.04.17 BFC and Pythia QA: 10% gain shift: lowered thresholds

See this post from Alice for QCD sample rates with slightly lower thresholds

QCD and gamma-jet data samples are described here (filter parameters are listed below)

Pythia filter configuration

StEemcGammaFilter:: running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events
StEemcGammaFilter:: mConeRadius 0.22 mSeedThreshold 2.4 mClusterThreshold 3.3 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 = 2.8 GeV
StChain:INFO - Init() : Cluster eT threshold = 3.8 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

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

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.01732	0.1343	0.00628
pt=3-4	0.05818	0.3001	0.0261
pt=4-6	0.1317	0.46492	0.07864
pt=6-9	0.17804	0.56314	0.15092
pt=9-15	0.17226	0.57516	0.15964
pt=15-25	0.13356	0.46894	0.1179
pt=25-35	0.062	0.28546	0.05482

Figure 2: False rejection (Y-axis scale is 10^-3)

False reject: fract. of generated events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	2e-05	2e-05	0
pt=3-4	6e-05	6e-05	0
pt=4-6	6e-05	4e-05	0
pt=6-9	0.0001	6e-05	0
pt=9-15	2e-05	2e-05	0
pt=15-25	2e-05	0	0
pt=25-35	2e-05	2e-05	0

Accept/Reject relative to the number of Pythia filter accepted events

Figure 3: Fraction of accepted events relative to Pythia filter accepted events

Accept rate: fract. of Pythia-filtered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.128816	1	0.0415488
pt=3-4	0.193735	1	0.0808397
pt=4-6	0.283232	1	0.161877
pt=6-9	0.316049	1	0.260291
pt=9-15	0.29943	1	0.26876
pt=15-25	0.284813	1	0.240415
pt=25-35	0.217053	1	0.175156

Figure 4: False rejection relative to Pythia filter accepted events

False reject: fract. of Pythia-filtered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0	0	0
pt=3-4	0	0	0
pt=4-6	4.30182e-05	0	0
pt=6-9	7.10303e-05	0	0
pt=9-15	0	0	0
pt=15-25	4.26494e-05	0	0
pt=25-35	0	0	0

2010.04.17 Pythia/BFC gamma-filter accaptance vs. gamma candidate pt, energy, and eta

Data sample used:
Pythia prompt photon Monte-Carlo (partonic pt bins are combined without weights)

Common event cuts:
reconstruct at least one gamma candidate, |v_z| <120, !=0 l2e-gamma-trigger=fired

Figure 1:

upper panel:
Number of candidates vs. gamma candidate pt for specific filter/trigger fired.
lower panel:
Pythia/bfc filter false rejection vs. gamma candidate pt
left panel:
Pythis (bfc) filter configured with seed=2.6 (3.4) and cluster=3.6 (4.5)
right panel:
Pythis (bfc) filter configured with seed=2.4 (2.8) and cluster=3.3 (3.8)

Figure 2: Same as Fig. 1 vs. gamma candidate energy

Figure 3: Same as Fig. 1 vs. gamma candidate pseudo-rapidity

2010.04.30 BFC and Pythia QA after gamma-maker 3x3 cluser fix

Related inks:

3x3 clusering and gamma maker fix
Alice' study with QCD sample
QCD and gamma-jet data samples are described here
(filter parameters are listed below)

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 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.00288	0.01472	0.00616
pt=3-4	0.01504	0.06192	0.02576
pt=4-6	0.06824	0.22112	0.07548
pt=6-9	0.15848	0.45836	0.15016
pt=9-15	0.1584	0.50416	0.15812
pt=15-25	0.12112	0.42076	0.11916
pt=25-55	0.05356	0.2292	0.0538
QCD
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	2e-05	0.00258	6e-05
pt=3-4	6e-05	0.00854	0.00022
pt=4-6	0.00072	0.03492	0.00076
pt=6-9	0.00564	0.144	0.00496
pt=9-15	0.0242	0.36036	0.0186
pt=15-25	0.07368	0.47592	0.05008
pt=25-55	0.0684553	0.323374	0.0557724

Figure 2: False rejection (Y-axis scale is 10^-3)

False reject: fract. of generated events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.0002	0.00012	0
pt=3-4	0.0002	0.00056	0
pt=4-6	0.00096	0.00132	0
pt=6-9	0.0006	0.00032	0
pt=9-15	0.0002	4e-05	0
pt=15-25	4e-05	0	0
pt=25-55	8e-05	8e-05	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.00012	0	0
pt=9-15	8e-05	8e-05	0
pt=15-25	0.00016	4e-05	0
pt=25-55	0.000203252	0	0

Accept/Reject relative to the number of Pythia filter accepted events

Figure 3: Fraction of accepted events relative to Pythia filter accepted events

Accept rate: fract. of Pythia-filtered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.17663	1	0.171196
pt=3-4	0.215762	1	0.197028
pt=4-6	0.291787	1	0.259407
pt=6-9	0.344096	1	0.313814
pt=9-15	0.314107	1	0.301254
pt=15-25	0.28786	1	0.269798
pt=25-55	0.233333	1	0.211693
QCD
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.00775194	1	0
pt=3-4	0.00702576	1	0.00936768
pt=4-6	0.0183276	1	0.0160367
pt=6-9	0.0386111	1	0.0294444
pt=9-15	0.0667111	1	0.048618
pt=15-25	0.154816	1	0.102706
pt=25-55	0.211691	1	0.163545

Figure 4: False rejection relative to Pythia filter accepted events

False reject: fract. of Pythia-filtered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.0108696	0	0
pt=3-4	0.00258398	0	0
pt=4-6	0.00325615	0	0
pt=6-9	0.00113448	0	0
pt=9-15	0.00031736	0	0
pt=15-25	9.50661e-05	0	0
pt=25-55	0.00017452	0	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.000833333	0	0
pt=9-15	0.000222	0	0
pt=15-25	0.000252143	0	0
pt=25-55	0.000628536	0	0

Figure 5: False rejection relative to trigger accepted events

False reject: fract. of triggered events
GammaJet
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.623377	0.558442	0
pt=3-4	0.517081	0.464286	0
pt=4-6	0.240594	0.192369	0
pt=6-9	0.0348961	0.0167821	0
pt=9-15	0.0111308	0.00151783	0
pt=15-25	0.00872776	0.00134273	0
pt=25-55	0.00966543	0.00148699	0
QCD
pt bin	Bfc Filter	Pythia Filter	L2gamma Trigger
pt=2-3	0.333333	0.333333	0
pt=3-4	0.363636	0.181818	0
pt=4-6	0.631579	0.263158	0
pt=6-9	0.282258	0.0887097	0
pt=9-15	0.202151	0.0301075	0
pt=15-25	0.0686901	0.000798722	0
pt=25-55	0.0291545	0.000728863	0

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:

BFC and Pythia filter biases vs. partonic pt
3x3 clusering and gamma maker fix
Alice' study with QCD sample
QCD and gamma-jet data samples are described here
(filter parameters are listed below)

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:

Filter's configurations and bias/acceptance tests used for this timing study

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:

Results before bug fix (See Fig. 1 upper right)

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:

Run 6 Murad's note:
http://www.star.bnl.gov/HyperNews-star/protected/get/jetfinding/664/1/1.html
Run 9 Pibero's note:
http://www.star.bnl.gov/HyperNews-star/protected/get/jetfinding/1010.html
Jan's note on "tower ADC - pedestal > 3 * RMS":
http://www.star.bnl.gov/HyperNews-star/protected/get/jetfinding/1010/2.html
Ily'as post on EEMC only jets for 5-hit branch:
http://www.star.bnl.gov/HyperNews-star/protected/get/jetfinding/1010/3/1/1/1/2.html
Trigger list for Run 9:
http://www.star.bnl.gov/protected/common/common2009/trigger2009/triggers2009.html
Pibero's jet tree format:
http://www4.rcf.bnl.gov/~pibero/spin/jets/2010.04.20
Disk space estimates for Run 9:
http://www.star.bnl.gov/HyperNews-star/protected/get/jetfinding/1010/1/1/2/1.html
http://www.star.bnl.gov/HyperNews-star/protected/get/jetfinding/1010/3/1/1/1.html
http://www.star.bnl.gov/HyperNews-star/protected/get/jetfinding/1010/3/1/1/1/2/1.html

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

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

06 Jun

June 2010 posts

2010.06.15 First look at data vs. TuneA/Perugia0 filtered MC with latest EEMC geometry

Data samples and colour coding

black: pp2006 data
open green: MC-QCD-TuneA, partonic pt 4-35
solid green: MC-QCD-Perugia0, partonic pt 4-35
(these not shown yet -> still generating data points)
open red MC-prompt-photon-TuneA, partonic pt 3-35
solid red MC-prompt-photon-Perugia0, partonic pt 3-35

Event selection

di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data: L2e-gamma triggered events
No trigger emulation in Monte-Carlo yet
MC scaled to 3.164^pb based on Pythia luminosity (no fudge factors)

Figure 1: Reconstructed photon candidate pt, pt_gamma (no cut on pt_gamma, pt_jet > 5GeV)

Figure 2: Partonic pt distribution (pt_gamma>7GeV, pt_jet > 5GeV)

Figure 3: Estimate of the contribution from low partonic pt,
only QCD-TuneA MC are shown (pt_gamma>7GeV, pt_jet > 5GeV)
Black line: Exponential fit to partonic pt distribution in 4-7GeV range
                   (pt_gamma>7GeV cut for the photon candidate)
Red line: Exponential fit extrapolated to the partonic pt range below 4GeV.
                Ratio of the area under the red line (integral pt=0-4geV)
                to the area under the green line (integral pt=4-35GeV) is 0.0028.

Comments

(based on Fig. 3)
I would propose we drop both of the lowest parton pt bins,
i.e. pt=2-3 and pt=3-4 (Inherited error for pt_gamma>7GeV < 0.3%)
and instead use our CPU time to produce more
statistics in the 4-35 partonic pt range.
(based on Fig. 2)
There is a small difference between CDF-Tune-A and
Perugia0 tunes partonic pt distributions
even for the prompt photon Monte-Carlo.
Comparison with Perugia0 QCD MC is coming.
Hopefully after that we will be able to decide what
Pythia tune is better match the L2e-gamma data.

2010.06.17 Pythia TuneA/Perugia0 filtered MC vs. pp2006 data

Data samples and colour coding

black circles: pp2006 data
open green: MC-QCD-TuneA, partonic pt 4-35
solid green: MC-QCD-Perugia0, partonic pt 4-35
open red MC-prompt-photon-TuneA, partonic pt 3-35
solid red MC-prompt-photon-Perugia0, partonic pt 3-35

Event selection

di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data: L2e-gamma triggered events
No trigger emulation in Monte-Carlo yet
MC scaled to 3.164^pb based on Pythia luminosity (no fudge factors)

Plots before cuts on photon candidate pt

Figure 1: Reconstructed photon candidate pt, pt_gamma (no cut on pt_gamma, pt_jet > 5GeV)

Figure 2: Partonic pt distribution (no cut on pt_gamma, pt_jet > 5GeV)

Plots with pt_gamma>7GeV cut

Figure 3: Partonic pt distribution (pt_gamma>7GeV, pt_jet > 5GeV)

Figure 4: Away side jet pt (pt_gamma>7GeV, pt_jet > 5GeV)

Figure 5: Reconstructed z vertex (pt_gamma>7GeV, pt_jet > 5GeV)

Figure 6: Partonic pt distribution for Pythia CDF-Tune-A QCD simulations (pt_gamma>7GeV, pt_jet > 5GeV)

Estimate of the contribution from low partonic pt:
Black line: Exponential fit to partonic pt distribution in 4-7GeV range
Red line: Exponential fit extrapolated to the partonic pt range below 4GeV.
Ratio of the area under the red line (integral over pt=0-4GeV)
to the area under the green line (integral over pt=4-35GeV) is 0.0028 (<0.3%)

Comments

Simulations with Perugia0 tune has a higher yield than that from CDF-Tune-A simulations
Shapes vs. partonic pt are different for Perugia0 and CDF-TuneA simulations
Shapes vs. reconstructed variables are similar for Perugia0 and CDF-TuneA simulations
(based on Fig. 6) I would propose we drop both of the lowest parton pt bins,
i.e. pt=2-3 and pt=3-4 (Inherited error for pt_gamma>7GeV < 0.3%)
and instead use CPU time to produce more statistics in the 4-35 partonic pt range.
More discussion at phana hyper news:
http://www.star.bnl.gov/HyperNews-star/protected/get/phana/496.html

Additional figures

Figure 7a: Photon candidate yield vs. rapidity (pt_gamma>7GeV, pt_jet > 5GeV)
Left: pt_gamma>7GeV; right: zoom into eta < 1 region

Figure 7b: yield vs. jet1 momentum (pt_gamma>7GeV, pt_jet > 5GeV)
Figure 7c: eta yield without pt_gamma cut
Yields ratio for eta <0.95 to the total yield is ~ 1.7% (1004/58766 = 0.0171)

Figure 8: Photon candidate yield vs. rapidity (pt_gamma>7GeV, pt_jet > 5GeV)

Note: trigger condition is not applied in simulations yet
but at high pt the data to Pythia CDF-Tune-A ratio is about 1.28 (at 9GeV: 3200/2500),
what is consistent with an additional 25% scaling factor
used for CIPANP 2009 presentation (see slide 6)

2010.06.18 L2e-gamma trigger effect: Py-CDF-Tune-A, Py-Perugia0, and pp2006 data comparison

Data samples and colour coding

black circles: pp2006 data
open green: MC-QCD-TuneA, partonic pt 4-35
solid green: MC-QCD-Perugia0, partonic pt 4-35
open red MC-prompt-photon-TuneA, partonic pt 3-35
solid red MC-prompt-photon-Perugia0, partonic pt 3-35

Event selection

di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data: L2e-gamma triggered events
No trigger emulation in Monte-Carlo yet
MC scaled to 3.164^pb based on Pythia luminosity (no fudge factors)

Figure 1: Reconstructed photon candidate pt, pt_gamma (no cut on pt_gamma, pt_jet > 5GeV)

Figure 2: Same as Fig. 1 with L2e-gamma condition simulated in Monte-Carlo

Figure 3: Same as Fig. 1, added distribution for photon pt from geant record (prompt photon MC only)

Figure 4: raw jet pt from jet trees: QCD pt=6-9
upper plot: mit0043 M. Betancourt simulations (MIT Simulation Productions)
bottom plot: new filtered MC

2010.06.28 Tests of L2e-gamma trigger emulation with single photon Monte-Carlo

Monte-Carlo configuration

Single photon in the EEMC (flat in eta, pt, phi)
Narrow vertex distribution with sigma=1cm
10 muons thrown in Barrel (|eta|<0.5) to reconstruct vertex
3 muons thrown in each BBC (|eta|~4) to fire the trigger
Run 6 L2e-gamma-trigger id = 137641
STAR geometry tag: y2006h
Photon cuts:
1.1 < eta < 1.95
3 < pt < 15 GeV
0 < phi < 2pi

Trigger effect vs. thrown photon pt, eta, and energy

Figure 1:
Yields vs. thrown photon pt
left: Yields with (red) and without (black) L2e-gamma trigger condition
right: Yield ratio (with/without trigger)

Figure 2: Same as Fig. 1 vs. thrown eta

Figure 3: Same as Fig. 1 vs. thrown energy

Trigger effect vs. reconstructed energy in the EEMC (high tower, 2x1, 3x3, energy and total E_T)

Figure 4: Same as Fig. 1 vs. total reconstructed transverse energy

Figure 5: Same as Fig. 1 vs. reconstructed high tower energy

Figure 6: Same as Fig. 1 vs. reconstructed energy of the 2x1 tower cluster

Figure 7: Same as Fig. 1 vs. reconstructed energy of the 3x3 tower cluster

2010.06.30 Py-tunes (GEANT+L2e-gamma trigger) vs. Run 6 data

Data samples and colour coding

black circles: pp2006 data
open green: MC-QCD-TuneA, partonic pt 4-35
solid green: MC-QCD-Perugia0, partonic pt 4-35
open red MC-prompt-photon-TuneA, partonic pt 3-35
solid red MC-prompt-photon-Perugia0, partonic pt 3-35

Event selection

di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data : L2e-gamma triggered events
Monte-Carlo: emulated L2e-gamma triggered condition
MC scaled to 3.164^pb based on Pythia luminosity (no fudge factors)

Figure 1: Reconstructed photon candidate pt (no pt_gamma cut, pt_jet > 5GeV)
L2e-gamma condition simulated in Monte-Carlo

Figure 2: Yield ratios (no pt_gamma cut, pt_jet > 5GeV)
Black: data[pp2006] / QCD[Perigia0]
Green: QCD[Perigia0] / QCD[CDF-Tune-A]
Red: g-jet[Perigia0] / g-jet[CDF-Tune-A]

Figure 3: Vertex z distribution (pt_gamma>7GeV, pt_jet > 5GeV)

Figure 4: Simulation yield vs. partonic pt (no pt_gamma cut, pt_jet > 5GeV)

07 Jul

July 2010 posts

2010.07.02 Tests of L2e-gamma trigger emulation with full Pythia+Geant Monte-Carlo

Event selection

di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data : L2e-gamma triggered events
Run 6 L2e-gamma trigger algo: E_T[2x2] > 5.2, E_T[high tower]>3.7
Monte-Carlo: emulated L2e-gamma triggered condition
MC scaled to 3.164^pb based on Pythia luminosity (no fudge factors)

Plots for the ratio of N[passed L2] to N[before trigger]

Figure 1: Trigger effect vs. reconstructed photon candidate pt (3x3 patch) (no pt_gamma cut, pt_jet > 5GeV)
Dashed lines: Pythia Tune A, solid lines Pythia Perugia0 tune

Figure 2: Trigger effect vs. simulated direct photon pt (no pt_gamma cut, pt_jet > 5GeV)

Figure 3: Trigger effect vs. simulated direct photon eta (no pt_gamma cut, pt_jet > 5GeV)

Figure 4: Trigger effect vs. reconstructed vertex z (no pt_gamma cut, pt_jet > 5GeV)

Figure 5: Trigger effect vs. reconstructed vertex z (with additional pt_gamma >7GeV, pt_jet > 5GeV)

Figure 6: Trigger effect vs. 2x2 cluster Et (no pt_gamma cut, pt_jet > 5GeV)

Figure 7: Trigger effect vs. 1x1 cluster (high tower) Et (no pt_gamma cut, pt_jet > 5GeV)

2010.07.09 Table of L0-BBC, L0-EEMC, and L2Egamma triggers biases

Real and simulated trigger decisions:

BBC stands for emulated L0 BBC trigger condition
EEMC stands for emulated Run 6 L0 EEMC (137832) trigger condition
L2 stands for emulated L2E-gamma (137641) trigger condition
Trig event satisfied to all three simulated trigger conditions: BBC+EEMC+L2
data-EEMC stands for real data L0 EEMC (137832) trigger condition
(available only for fast offline data, not filled in yet)
data-L2 stands for real data L2E-gamma (137641) trigger condition

Data samples:

pp2006 data (3.4K events from st_physics production)
Pythia prompt photon simulations (147 events gamma-filtered for partonic pt=3-25GeV)
Pythia QCD 2->2 simulations (45 events gamma-filtered for partonic pt=6-9GeV)

Notations in the table:

XXX=0 - stands for XXX trigger did not fired
XXX=1 - stands for XXX trigger did fired
XXX=0 YYY=1 stands for XXX trigger did not fired, but YYY did fired


sample	Total	BBC=1	EEMC=1	L2=1	L2=1 EEMC=0	L2=1 BBC=0	L2=0 BBC=1	L2=0 EEMC=1

pp2006, st_physics
(counts)	3396	3306	568	549	5	2	2759	24
(%)	1.0000	0.9735	0.1673	0.1617	0.0015	0.0006	0.8124	0.0071

gamma-jets (3-25GeV)
(counts)	147	119	127	122	3	24	21	8
(%)	1	0.80952	0.86395	0.82993	0.020	0.16327	0.14286	0.054

QCD (6-9GeV)
(counts)	45	38	19	19	1	4	23	1
(%)	1	0.84444	0.42222	0.42222	0.022	0.08889	0.51111	0.022


simu vs. real triggers	data-EEMC=1	data-L2=1	L2=1 data-L2=0	L2=0 data-L2=1	Trig=1 data-L2=0	Trig=0 data-L2=1	EEMC=0 data-EEMC=1	EEMC=1 data-EEMC=0
	572	548	7	6	0	6	0	4
	0.1673	0.1617	0.0021	0.0018	0.0000	0.0018	0.0000	0.0012

2010.07.14 Pythia/BFC gamma-filter bias tests (vs. gamma pt, eta, energy, and phi)

Related inks:

Tests of L2e-gamma trigger emulation with full Pythia+Geant Monte-Carlo
Previous study before filter emulator fixes
QCD and gamma-jet data samples are described here
(filter parameters are listed below)

Number of generated events per partnic pt bin (pt binsa are: 2-3, 3-4, 4-6, 6-9, 9-15, 15-35):
gamma-jets (2-4): 25K/bin
gamma-jets (4-35): 12.5K/bin
QCD(2-4): 50K/bin
QCD(4-35): 25K/bin

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 (Fig. 1 only)
BFC:INFO - setTowerGainSpread(): gain spread: 0; gain mean value: 1.1 (Fig. 2 and below)

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

(Includes all recent fixes to trigger emulator configuration/software)
emulated L2E-gamma trigger for Run 2006 [eemc-http-mb-l2gamma:: id 137641]
Trigger conditions:
cluster Et (3x3) = 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 1:
(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 (shown in blue)]
No gain shifts in StEEmcSlowMaker
Note: statistics without gain shifts is smaller because ~10-20% jobs died at RCF

Figure 2: Same as Fig. 1, but with +10% (1.1) gain shifts in StEEmcSlowMaker.

Figure 3: Same as Fig. 2 vs. candidate eta (with +10% (1.1) gain shifts in StEEmcSlowMaker).

Figure 4: Same as Fig. 2 vs. candidate azimuthal angle (with +10% (1.1) gain shifts in StEEmcSlowMaker).

Figure 5: Same as Fig. 2 vs. energy (with +10% (1.1) gain shifts in StEEmcSlowMaker).

2010.07.16 Pythia/BFC gamma-filter bias tests with realistic gain variation

Related inks:

Previous study with ideal gains and gain shift only of +10%
Tests of L2e-gamma trigger emulation with full Pythia+Geant Monte-Carlo
Previous study before filter emulator fixes
QCD and gamma-jet data samples are described here
(filter parameters are listed below)

Number of generated events per partnic pt bin (pt binsa are: 2-3, 3-4, 4-6, 6-9, 9-15, 15-35):
gamma-jets (2-4): 25K/bin
gamma-jets (4-35): 12.5K/bin
QCD(2-4): 50K/bin
QCD(4-35): 25K/bin

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 with realistic gain shift/smearing

Figure 1: Error for gains from MIP study minus ideal
Data digitized from Scott's presentation at 2008 Calibartion workshop

BFC:INFO - setTowerGainSpread(): gain spread: 0.1; gain mean value: 1.05 (Fig. 1,3 only)
BFC:INFO - setTowerGainSpread(): gain spread: 0.1; gain mean value: 0.95 (Fig. 2,4 and below)

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

(Includes all recent fixes to trigger emulator configuration/software)
emulated L2E-gamma trigger for Run 2006 [eemc-http-mb-l2gamma:: id 137641]
Trigger conditions:
cluster Et (3x3) = 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 2:
(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 (shown in blue)]
StEEmcSlowMaker configured with +5% (scale factor=1.05) gain shifts and 0.1 sigma
Previous plots: 125K events per pt-bin, 250K/pt-bin
(figure below combines previous statisitcs + 18K for partonic pt=6-9 and pt=9-15 GeV bins)

Figure 2b: Filter bias per partonic pt bin (QCD simulations only)

Figure 3: Same as Fig. 1 with gain shift=0.95 and sigma=0.1

Figure 4: Same as Fig. 1 vs. candidate eta with gain shift=1.05 and sigma=0.1

Figure 5: Same as Fig. 1 vs. candidate eta with gain shift=0.95, sigma=0.1

2010.07.20 EEMC simulation spreadsheet: prompt photons and QCD (Updated)

Simulation request spreadsheet (QCD@L=2/pb, photons@L=10/pb)

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	Number of Pythia filtered events
g-jets
2-3	0.00870	0.2663	0.00232	1280000	10.0	29659	18.75	4.63	111360
3-4	0.03300	0.3787	0.01250	290000	10.0	36237	15.38	5.90	95700
4-6	0.10920	0.5521	0.06029	126700	10.0	76387	25.66	12.67	138356
6-9	0.22320	0.6688	0.14928	26860	10.0	40096	10.74	6.85	59952
9-15	0.25360	0.6274	0.15911	4636	10.0	7376	2.17	1.28	11757
15-35	0.21360	0.5394	0.11522	347	10.0	399	0.17	0.07	740
totals:				1728543		190154	72.9	31.40	417865
						0.19	0.2	years
QCD
2-3	0.00067	0.0185	0.00001	8089000000	0.0		0	0.00
3-4	0.00298	0.0268	0.00008	1302000000	0.0		0	0.00
4-6	0.01312	0.0240	0.00031	413600000	2.0	260469	1428.35	51.72	10852864
6-9	0.06140	0.0640	0.00393	60620000	2.0	476425	1023.1	102.54	7444136
9-15	0.17692	0.1120	0.01982	7733000	2.0	306459	440.74	64.02	2736245
15-35	0.25480	0.2260	0.05758	404300	2.0	46563	34.53	9.72	206031
totals:				9873357300		1089916	2926.72	228.00	21239276
						1.09	8.02	years

QCD lumi, 1/pb	number of events, x 10e6	CPU years	disk space, Gb			total time with 50 CPU, weeks	total time with 100 CPU, weeks
1	.74	4.2	145.4			4.4	2.2
2	1.28	8.2	259.4			8.6	4.3

Timing tests

Figure 1: Timing tests for BFC and Pythia gamma-filters (in seconds)

2010.07.22 Run 6 EEMC gamma-filtered simulation request

Submitted run-6 photon-jet simulation request for spin physics

Request last updated on Aug 19, 2010

Run 6 EEMC gamma-filtered simulation request summary

Total resources estimate for QCD with 1/pb and prompt-photon with 10/pb suimulations:

CPU: 4.2 CPU years (2.2 weeks of running on a 100 CPUs)
Disk space: 0.15Tb
Numbe of filtered events: 0.74M

partonic pt	QCD		prompt photon
	total	Pythia	total	Pythia
2-3	0	0	30K	110K
3-4	0	0	36K	95K
4-6	130K	5.5M	76K	140K
6-9	240K	3.7M	40K	60K
9-15	150K	1.4M	10K	12K
15-35	23K	100K	1K	3K

Latest filter bias/timing test and simulation request spreasheet

Note: These and all other studies are linked from here

Filter code in cvs

Pythia gamma-filter code:
StEemcGammaFilter.cxx
StEemcGammaFilter.h
BFC gamma-filter code:
StEemcGammaFilterMaker.cxx
StEemcGammaFilterMaker.h
eemcGammaFilterMakerParams.idl

Run 6 GMT timestamps

See this study for more details and plots

Request an equal fraction (10%) for each of the 10 timestamps below:
sdt20060516.152000 (GMT during run 7136022)
sdt20060518.073700 (GMT during run 7138010)
sdt20060520.142000 (GMT during run 7140024)
sdt20060521.052000 (GMT during run 7141011)
sdt20060522.124500 (GMT during run 7142029)
sdt20060523.204400 (GMT during run 7143044)
sdt20060525.114000 (GMT during run 7145023)
sdt20060526.114000 (GMT during run 7146020)
sdt20060528.144500 (GMT during run 7148028)
sdt20060602.071500 (GMT during run 7153015)

------------------------ REQUEST DETAILS BELOW ----------------------------------------

prompt photons and QCD simulations

Request TypeEvent generator simulation, with filtering

General Information


Request ID
Priority: EC	0
Priority: pwg	High
Status	New
Physics Working Group	Spin
Requested by	Photon group for SPIN PWG
Contact email(s)	ilya.selyuzhenkov@gmail.com, bridgeman@hep.anl.gov
Contact phone(s)
PWG email(s)	starspin-hn@www.star.bnl.gov
Assigned Deputy:	Not assigned
Assigned Helper:	Not assigned

Description

Endcap photon-jets request

Global Simulation Settings


Request type:	Event generator simulation, with filtering
Number of events	See list for each partonic pt bins
Magnetic Field	Full-Field
Collision Type	pp@200GeV
Centrality	---- SELECT CENTRALITY ----
BFC tags	trs fss y2006h Idst IAna l0 tpcI fcf ftpc Tree logger ITTF Sti VFPPV bbcSim tofsim tags emcY2 EEfs evout -dstout IdTruth geantout big fzin MiniMcMk eemcDb beamLine clearmem
Production	---- SELECT PRODUCTION TAG ----
Geometry: simu	y2006h
Geometry: reco	y2006h
Library	use library with approved filter code checked in
Vertex option	Leave vertex to be reconstructed vertex, and use VFPPVnoCTB with beamline
Pileup option	No
Detector Set	TPC, ETOW, BTOW, BSMD, ESMD, BPRS, EPRE1, EPRE2, EPOST, TOF, BBC, SVT, SSD

Data Sources

MC Event Generator


Event generator	Pythia
Extra options	Additional libraries required for Eemc-gamma Pythia-level filter gexec $ROOTSYS/lib/libCint.so gexec $ROOTSYS/lib/libCore.so gexec $ROOTSYS/lib/libMathCore.so gexec $ROOTSYS/lib/libMatrix.so gexec $ROOTSYS/lib/libPhysics.so gexec .sl53_gcc432/lib/StMCFilter.so // filter library Prompt photon Pythia processes: MSUB (14)=1 MSUB (18)=1 MSUB (29)=1 MSUB (114)=1 MSUB (115)=1 QCD 2->2 Pythia processes: MSUB (11) = 1 MSUB (12) = 1 MSUB (13) = 1 MSUB (28) = 1 MSUB (53) = 1 MSUB (68) = 1 Pro-pT0 Pythia tune: call pytune(329)
Vertex Z, cm	-120 < Vertex < 120
Gaussian sigma in X,Y,Z if applicable	0, 0, 55 200 GeV
Vertex offset: x, mm	0.0cm
Vertex offset: y, mm	-0.3cm
Φ (phi), radian	0 < Φ < 6.29
η (eta)	Default (include Barrel, Endcap, BBC)
Pt bin, GeV	See list above for QCD and g-jet samples
Macro file	Pythia gamma-filter code: StEemcGammaFilter.cxx StEemcGammaFilter.h BFC gamma-filter code: StEemcGammaFilterMaker.cxx StEemcGammaFilterMaker.h eemcGammaFilterMakerParams.idl Private bfc: /star/u/seluzhen/star/spin/MCgammaFilter/scripts/bfc.C

2010.07.23 PyTune comparison with photon candidates: Perugia0 vs. Pro-PT0

Data samples and colour coding

black Pythia QCD Monte-Carlo with Pro-Pt0 tune (pytune=329), partonic pt 9-15
red Pythia QCD Monte-Carlo with Perugia0 tune (pytune=320), partonic pt 9-15

Event selection

Ran full Pythia+GSTAR simulation and require at least one
reconstucred EEMC photon candidate in the gamma Maker.

Figure 1:
Left: Reconstructed photon candidate transverse momentum (no normalization factor applied)
Right: ratio of Pro-Pt0/Perugia0 simulations (solid Line: "a+b*x" fit to ratio)
Event selections: require at least one reconstucred EEMC photon candidate

Figure 2:
Same as in Fig. 1 with different event selection criteria:
L2E-gamma, |v_z| < 120cm, at least one EEMC gamma candidate

Pytune parameters comparison table

pytune(320) Perugia 0
P. Skands, Perugia MPI workshop October 2008
and T. Sjostrand & P. Skands, hep-ph/0408302
CR by M. Sandhoff & P. Skands, in hep-ph/0604120
LEP parameters tuned by Professor

pytune(329) Pro-pT0
See T. Sjostrand & P. Skands, hep-ph/0408302
and M. Sandhoff & P. Skands, in hep-ph/0604120
LEP/Tevatron parameters tuned by Professor

Red text indicates the parameter which are different between tunes

Parameter	Perugia 0	Pro-pT0	Parameter description
MSTP(51)	7	7	PDF set
MSTP(52)	1	1	PDF set internal (=1) or pdflib (=2)
MSTP(64)	3	2	ISR alphaS type
PARP(64)	1.0000	1.3000	ISR renormalization scale prefactor
MSTP(67)	2	2	ISR coherence option for 1st emission
PARP(67)	1.0000	4.0000	ISR Q2max factor
MSTP(68)	3	3	ISR phase space choice & ME corrections (Note: MSTP(68) is not explicitly (re-)set by PYTUNE)
MSTP(70)	2	2	ISR IR regularization scheme
MSTP(72)	1	0	ISR scheme for FSR off ISR
PARP(71)	2.0000	2.0000	FSR Q2max factor for non-s-channel procs
PARJ(81)	0.2570	0.2570	FSR Lambda_QCD scale
PARJ(82)	0.8000	0.8000	FSR IR cutoff
MSTP(81)	21	21	UE model
PARP(82)	2.0000	1.8500	UE IR cutoff at reference ecm (Note: PARP(82) replaces PARP(62).)
PARP(89)	1800.0000	1800.0000	UE IR cutoff reference ecm
PARP(90)	0.2600	0.2200	UE IR cutoff ecm scaling power
MSTP(82)	5	5	UE hadron transverse mass distribution
PARP(83)	1.7000	1.8000	UE mass distribution parameter
MSTP(88)	0	0	BR composite scheme
MSTP(89)	1	1	BR colour scheme
PARP(79)	2.0000	1.1800	BR composite x enhancement
PARP(80)	0.0500	0.0100	BR breakup suppression
MSTP(91)	1	1	BR primordial kT distribution
PARP(91)	2.0000	2.0000	BR primordial kT width <\|kT\|>
PARP(93)	10.0000	7.0000	BR primordial kT UV cutoff
MSTP(95)	6	6	FSI colour (re-)connection model
PARP(78)	0.3300	0.1700	FSI colour reconnection strength
PARP(77)	0.9000	0.0000	FSI colour reco high-pT dampening streng
MSTJ(11)	5	5	HAD choice of fragmentation function(s)
PARJ(21)	0.3130	0.3130	HAD fragmentation pT
PARJ(41)	0.4900	0.4900	HAD string parameter a
PARJ(42)	1.2000	1.2000	HAD string parameter b
PARJ(46)	1.0000	1.0000	HAD Lund(=0)-Bowler(=1) rQ (rc)
PARJ(47)	1.0000	1.0000	HAD Lund(=0)-Bowler(=1) rb

08 Aug

August 2010 posts

2010.08.09 PyTune comparison with gamma candidates from dijets: Perugia0 vs. Pro-PT0

Tunes compared

CDF Tune A
Perugia0
Pro-pT0

Event selection

di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data : L2e-gamma triggered events
Monte-Carlo: emulated L2e-gamma triggered condition
MC scaled to 3.164^pb based on Pythia luminosity (no fudge factors)

Figure 1a: Reconstructed photon candidate pt (L2e-gamma condition simulated in Monte-Carlo)

Figure 1b: Same as Fig. 1 on a linear scale and zoom into low pt

Figure 2: QCD Monte-Carlo yield to pp2006 data ratio

Figure 3: Prompt photon Monte-Carlo yield to pp2006 data ratio

Figure 4: Simulation yield vs. partonic pt (on a linear scale)

Figure 5: Simulation yield vs. partonic pt (log scale)

2010.08.10 Timestamps study for the simulation request

Summary for the gamma filtered simulation request

Generate Monte-Carlo events for
10 different timestamps and 10% of statistics each:

sdt20060516.152000 (GMT during run 7136022)
sdt20060518.073700 (GMT during run 7138010)
sdt20060520.142000 (GMT durign run 7140024)
sdt20060521.052000 (GMT during run 7141011)
sdt20060522.124500 (GMT during run 7142029)
sdt20060523.204400 (GMT during run 7143044)
sdt20060525.114000 (GMT during run 7145023)
sdt20060526.114000 (GMT during run 7146020)
sdt20060528.144500 (GMT during run 7148028)
sdt20060602.071500 (GMT during run 7153015)

Original list of timestamps

(from http://www.star.bnl.gov/HyperNews-star/protected/get/phana/481.html)

sdt20060512.043500 (GMT during run 7132005)
sdt20060513.064000 (GMT during run 7133011)
sdt20060514.090000 (GMT during run 7134015)
sdt20060516.152000 (GMT during run 7136022)
sdt20060518.073700 (GMT during run 7138010)
sdt20060520.142000 (GMT durign run 7140024)
sdt20060521.052000 (GMT during run 7141011)
sdt20060522.124500 (GMT during run 7142029)
sdt20060523.204400 (GMT during run 7143044)
sdt20060525.114000 (GMT during run 7145023)
sdt20060526.114000 (GMT during run 7146020)
sdt20060528.144500 (GMT during run 7148028)
sdt20060602.071500 (GMT during run 7153015)
sdt20060604.191200 (GMT during run 7155043)

Figure 1: Number of events from Run 6 golden runs
which fired L2e-gamma trigger (trigger id 137641 and 127641)
(using jet trees regenerated in new format by Wayne/Renee)

7132005: 0
7133011: 0
7134015: 10474
7136022: 171217
7138010: 221567
7140024: 62826
7141011: 174207
7142029: 187048
7143044: 189752
7145023: 142799
7146020: 133758
7148028: 181269
7153015: 145129
7155043: 89428

Figure 2: Fraction of events per time stamps

09 Sep

September 2010 posts

2010.09.08 First look at the official EEMC gamma filtered production

Event selection

Official EEMC gamma filtered Monte-Carlo with Pro-pT0 tune
di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data : L2e-gamma triggered events
Monte-Carlo: emulated L2e-gamma triggered condition
MC scaled to 3.164^pb based on Pythia luminosity (no fudge factors)

Figure 1: Reconstructed photon candidate pt (L2e-gamma condition simulated in Monte-Carlo)

Figure 2: Partonic pt

Figure 3: Thrown photon pt (from Geant record, prompt photon sample only)

2010.09.10 Data to MC comparison with official EEMC gamma filtered production

Event selection

Official EEMC gamma filtered Monte-Carlo with Pro-pT0 tune
di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data : L2e-gamma triggered events
Monte-Carlo: emulated L2e-gamma triggered condition
QCD Monte-Carlo scaled to the yield in the data (MC down scaled by a factor of 1.8)
All plots with gamma pt >7GeV cut

Figure 1: Reconstructed photon candidate pt

Figure 2: Reconstructed away side jet pt

Figure 3: z vertex distribution

Figure 4: 3x3 tower cluster energy

Figure 5: 2x1 tower cluster energy

Figure 6: Reconstructed photon candidate (detector) rapidity

Figure 7: Reconstructed away side jet rapidity

2010.09.13 Data vs. official filtered MC: cluster energy ratios

Event selection

Official EEMC gamma filtered Monte-Carlo with Pro-pT0 tune
di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data : L2e-gamma triggered events
Monte-Carlo: emulated L2e-gamma triggered condition
QCD Monte-Carlo scaled to the yield in the data (MC down scaled by a factor of 1.8)
All plots with gamma pt >7GeV and jet pt >5GeV cuts

Figure 1: Cluster energy ratio: 2x1/2x2

Figure 2: Cluster energy ratio: 2x1/3x3

Figure 3: Cluster energy ratio: 2x2/3x3

Figure 4: Cluster energy ratio: 2x1/E(R=0.7)

Figure 5: Cluster energy ratio: 2x2/E(R=0.7)

Figure 6: Cluster energy ratio: 3x3/E(R=0.7)

2010.09.15 Data vs. official filtered MC: energy deposition in various EEMC layers

Event selection

Official EEMC gamma filtered Monte-Carlo with Pro-pT0 tune
di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data : L2e-gamma triggered events
Monte-Carlo: emulated L2e-gamma triggered condition
QCD Monte-Carlo scaled to the yield in the data (MC down scaled by a factor of 1.8)
All plots with gamma pt >7GeV and jet pt >5GeV cuts

Figure 1: 3x3 tower cluster energy

Figure 2: 25 central SMD u-strip energy

Figure 3: 3x3 pre-shower1 cluster energy

Figure 4: 3x3 pre-shower2 cluster energy

Figure 5: 3x3 post-shower cluster energy

2010.09.23 Neutral energy jet shape comparison for various tunes

Event selection

Official EEMC gamma filtered Monte-Carlo with Pro-pT0 tune vs. Ilya's private production for Pro-pT0, Perugia0, and CDF Tune A
di-jets from the cone jet-finder algorithm
photon and jet are opposite in phi:
cos (phi_gamma-phi_jet) < -0.8
pt away side jet > 5GeV
detector eta of the away side jet: |eta_jet_det| < 0.8
data : L2e-gamma triggered events
Monte-Carlo: emulated L2e-gamma triggered condition
QCD Monte-Carlo scaled to the yield in the data (All Monte-Carlo normalized to the total yield in the data)
All plots with gamma pt >7GeV and jet pt >5GeV cuts

Figure 1: 3x3 tower cluster energy to jet R=0.7 energy ratio
Left: Official production (Pro-Pt0), Right: Pro-Pt0 (Ilya private production)

Figure 2: 3x3 tower cluster energy to jet R=0.7 energy ratio - Perugia0 (Ilya private production)

Figure 3: 3x3 tower cluster energy to jet R=0.7 energy ratio - CDF Tune A (Ilya private production)

Figure 4: 3x3 tower cluster with jet thresholds energy to jet R=0.7 energy ratio - Official production (Pro-Pt0)

2010

Year 2010 posts

01 Jan

January 2010 posts

2010.01.04 y2006 vs y2009 EEMC pre-shower migration

Monte-Carlo setup:

2010.01.07 EEMC response to single photons with y2006h vs y2009a geometries

Monte-Carlo setup

2010.01.08 EEMC response 2006 vs. 2009: phi cuts

EEMC migration plots with cuts on TPC sector boundaries

Monte-Carlo setup

2010.01.12 W test sample QA

Data base setup

2010.01.13 W test sample QA: Pass 2

Update on Jan 14, 2010

2010.01.15 W test sample QA: Pass 3

2010.01.18 EEMC response 2006 vs. 2009: typo in TPC volume size bug fix

Monte-Carlo setup

Before and after the fix comparison

y2006h vs. y2009a comparison after a TPC typo fix

2010.01.18 W test sample: Cluster ratios and skewed gaussian fits

Cluster energy vs. original lepton energy

2010.01.26 Endcap/Barrel clustering with official W-MC

Lepton from W in the Endcap: 1.2 < eta_lepton < 1.9

Lepton from W in the Barrel: |eta_lepton| <1

02 Feb

February 2010 posts

2010.02.08 StEemcGammaFilterMaker QA: QCD vs. gamma-jet ccept/reject rates

StEemcGammaFilterMaker QA

Pythia generated processes

Number of generated events per parton pt bin

Filter configuration

stasim/Pythia options

Figures

Comments

2010.02.09 BFC and Pythia QA: relative to trigger

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

Accept/Reject relative to the number of triggered events

2010.02.10 Money plots for W cross section

2010.02.11 BFC and Pythia QA: Gain no-gain-spread, mean=1.05

Pythia filter configuration

BFC filter configuration

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

Accept/Reject relative to the number of triggered events

Accept/Reject relative to the number of Pythia filter accepted events

2010.02.11 BFC and Pythia QA: Gain spread=0.05, mean=0.95

Gain spread implementation in StEEmcSlowMaker.cxx (private version):

Pythia filter configuration

BFC filter configuration

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

Accept/Reject relative to the number of triggered events

2010.02.12 Final Pythia and BFC EEMC-gamma-filter paramter settings

Pythia generated processes

Number of generated events per parton pt bin

Pythia Filter configuration

BFC Filter configuration

GammaMaker configuration

EEMC SlowSimulator configuration

GSTAR/Pythia options

2010.02.16 Pythia and BFC filter QA vs. gamma candidate pt and eta

03 Mar

March 2010 posts

2010.03.02 Endcap photon-jets simulation request (draft)

Request summary

Latest filter bias/timing test and simulation request spreasheet

Filter code in cvs

Further information related to this request

Endcap photon-jets / QCD 2->2 simulations

04 Apr

April 2010 posts

2010.04.09 BFC and Pythia QA: 10% gain shift

Pythia filter configuration

BFC filter configuration

StEEmcSlowMaker

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

Accept/Reject relative to the number of Pythia filter accepted events

2010.04.17 BFC and Pythia QA: 10% gain shift: lowered thresholds

Pythia filter configuration

BFC filter configuration

StEEmcSlowMaker