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:

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