Ilya Selyuzhenkov January 30, 2008
jet trees by Murad Sarsour for pp2006 run, runId=7136022 (~60K events, no triggerId cuts yet)
R_EM1 >0.9 and R_EM2 < 0.9
cos(phi1 - phi2) < -0.9
nChargeTracks1 < 2
0 < nEEMCtowers1 < 3
Ilya Selyuzhenkov February 13, 2008
Gamma-jet selection cuts are discussed here. There are 278 candidates found for runId=7136022.
Transverse momentum distribution for the gamma-jet candidates can be found here.
Figure 1: Vertex z distribution.
Red line presents gamma-jet candidates (scaled by x50). Black is for all di-jet events.
Same data on a log scale is here.
EEMC response event by event for all 278 gamma-jet candidate can be found in this pdf file.
Each page shows SMD/Tower energy distribution for a given event:
First row on each page shows SMD response
for the sector which has a maximum energy deposited in the EEMC Tower
(u-plane is on the left, v-plane is on the right).
In the left plot (u-plane energy distribution) numerical values for
pt of the first jet (with maximum EM fraction) and the second jet are given.
In addition, fit results assuming gamma (single Gaussian, red line) or
neutral pion (double Gaussian, blue line ~ red+green) hypotheses are given.
If m_{gamma gamma} value is negative, then the reconstruction procedure has failed
(for example, no uv-strips intersection found, or tower energy and uv-strips intersection point mismatch, etc).
EEMC response for these "bad" events can be found in this pdf file.
If reconstruction procedure succeded, then
m_{gamma gamma} gives reconstructed invariant mass assuming that two gammas hit the calorimeter.
Figure 3: invariant mass distribution (assuming pi0 hypothesis).
Note, that I'm still working on my fitting algorithm (which is not explained here),
and fit results and the invariant mass distribution will be updated.
It is also shown the ratio for each u/v plane
of the integrated single Gaussian fit (red line) to the total energy in the plane
(look for "gamma U/V " values on the right v-plane plot).
Second and third rows on each page show the energy deposition in the
tower, pre-shower1, pre-shower2, and post-shower as a function of eta:phi (etaBin:phiBin).
Trying to isolate the real gammas which hits the calorimeter,
I have sorted events into different subsets based on the following set of cuts:
if (invMass < 0) reject
if (jet2_pt > jet1_pt) reject
if (jet1_pt < 7) reject
if (minFraction < 0.75) reject
(minFraction = gamma U/V - is a fraction of the integrated single Gaussian peak to the total energy in the uv-plane)
Figure 4: Sample gamma-jet candidate EEMC response
(all gamma-jet candidates selected according to these conditions can be found in this pdf file):
if (invMass < 0) reject
if (jet2_pt < jet1_pt) reject
if (jet2_pt < 7) reject
if (minFraction < 0.75) reject
Event by event EEMC response for pi0 (di-jet) candidates
selected according to these conditions can be found in this pdf file.
Ilya Selyuzhenkov February 20, 2008
After processing all available jet-trees for pp2006 (ppProductionLong),
and applying all "gamma-jet" cuts (which are described below):
there are 47K di-jet events selected
for pt1>7GeV there are 5,4K gamma-jet candidates (3,7K with an additional cut of pt1>pt2)
Figure 1: 2,4K events with both pt1, pt2 > 7GeV
721 candidates within a range of pt1>pt2 and both pt1, pt2 > 7 GeV
jet trees by Murad Sarsour for pp2006 run, number of runs processed: 323
4.7M di-jet events found (no triggerId cuts yet)
R_EM1 >0.9 and R_EM2 < 0.9
cos(phi1 - phi2) < -0.9
nChargeTracks1 < 2
0 < nEEMCtowers1 < 3
Ilya Selyuzhenkov February 27, 2008
Gamma-jet isolation cuts:
R_EM1 >0.9 and R_EM2 < 0.9
cos(phi1 - phi2) < -0.8
nCharge1 = 0
for each gamma-jet candidate finding a tower with a maximum energy
associated with a jet1 (jet with a maximum EM fraction).
Calculating energy of the cluster by finding all adjacent towers and adding their energy together.
R_cluster is defined as a ratio of the cluster energy
to the total energy in the calorimeter associated with a jet1.
Note, that with a cut Ncharge1 =0, energy in the calorimeter is equal to the jet energy.
Figure 4: R_cluster vs number of towers fired in EEMC (left) and BEMC (right). No pt cuts.
Figure 5: R_cluster vs number of towers fired in EEMC (left) and BEMC (right). Additional cut: pt1>7GeV
Figure 6: jet1 pseudorapidity vs number of towers fired in EEMC (left) and BEMC (right).
EEMC candidates: nTowerFiredBEMC=0
BEMC candidates: nTowerFiredEEMC=0
Figure 7: Pseudorapidity (left EEMC, right BEMC candidates)
Figure 8: Azimuthal angle (left EEMC, right BEMC candidates)
Figure 9: Transverse momentum (left EEMC, right BEMC candidates)
Figure 10: Transverse momentum (left EEMC, right BEMC candidates): pt1>7GeV
Figure 11: Transverse momentum (left EEMC, right BEMC candidates): pt1>7 and pt2>7
Ilya Selyuzhenkov March 03, 2008
Data set: ppLongitudinal, runId = 7136033.
Some observations/questions:
In general distributions look clean and good
Sectors 7 and 9 for v-plane and sector 7 for u-plane are noise.
Sector 9 has a hot strip (id ~ 120)
In sector 3, strips id=0-5 in v-plane are hot (see figure 2 right, bottom)
Sectors 2 and 8 in u-plane and sectors 3 and 9 in v-plane have missing strips id=283-288?
Strips 288 are always empty?
Figure 1:Average energy E in the strip vs sector and strip number (max < E > = 0.0027)
same figure on a log scale
Figure 2: Average energy E for E>0.02 (max < E > = 0.0682)
Same figure on a log scale
Ilya Selyuzhenkov March 12, 2008
Note: Di-jet transverse momentum distribution for these candidates can be found on figure 11 at this page
Figure 1:Invariant mass distribution for gamma-jet candidates assuming pi0 (2-gammas) hypothesys
Figure 2:Invariant mass distribution for gamma-jet candidates assuming pi0 (2-gammas) hypothesys
with an additional SMD isolation cut: gammaFraction >0.75
GammaFraction is defined as ratio of the integral
other SMD strips for the first peak to the total energy in the sector
pdf file (first 100 events) with event by event EEMC response for the candidates reconstructed into pion mass (gammaFraction >0.75)
pdf file with event by event EEMC response for the candidates not reconstructed into pion mass
(second peak not found), but has a first peak with gammaFraction >0.75.
Ilya Selyuzhenkov March 20, 2008
The procedure to discriminate gamma candidate from pions (and other background)
based on the SMD response is described at Pibero's web page.
Figure 1: Fit integral vs maximum residual for gamma-jet candidates requesting
no energy deposited in the EEMC pre-shower 1 and 2
(within a 3x3 clusters around tower with a maximum energy).
Black line is defined from MC simulations (see Jason's simulation web page, or Pibero's page above).
Figure 2: Fit integral vs maximum residual for gamma-jet candidates requesting requesting
no energy deposited in pre-shower 1 cluster and
no energy deposited in post-shower cluster (this cut is not really essential in demonstrating the main idea)
Figure 3: Fit integral vs maximum residual for gamma-jet candidates requesting requesting
non-zero energy deposited in both clusters of pre-shower 1 and 2
Event by event EEMC response for gamma-jet candidates for the case of
no energy deposited in the EEMC pre-shower 1 and 2 can be found in this pdf file
Figure 4: Fit integral vs maximum residual for gamma-jet candidates requesting
no energy deposited in the EEMC pre-shower 1 and 2
Figure 5: Fit integral vs maximum residual for gamma-jet candidates requesting requesting
no energy deposited in pre-shower 1 cluster and
no energy deposited in post-shower cluster
Figure 6: Fit integral vs maximum residual for gamma-jet candidates requesting requesting
non-zero energy deposited in both clusters of pre-shower 1 and 2
Event Monte Carlo shape allows to distinguish gammas from background by cutting at chi2/ndf < 0.5
(although the distribution looks wider than for the case of Will's shape).
Figure 7: Chi2/ndf for gamma-jet candidates using Monte Carlo shape requesting
no energy deposited in both clusters of pre-shower 1 and 2
Figure 8: Chi2/ndf for gamma-jet candidates using Monte Carlo shape requesting
non-zero energy deposited in both clusters of pre-shower 1 and 2
Less clear where to cut on chi2?
Figure 9: Chi2/ndf for gamma-jet candidates using Monte Carlo shape requesting
no energy deposited in both clusters of pre-shower 1 and 2
Figure 10: Chi2/ndf for gamma-jet candidates using Monte Carlo shape requesting
non-zero energy deposited in both clusters of pre-shower 1 and 2
Ilya Selyuzhenkov March 26, 2008
Definitions:
All results are for combined distributions from u and v planes: ([u]+[v])/2
Gamma-jet isolation cuts described here
Additional quality cuts:
Figure 1: F_peak vs maximum residual
for various cuts on energy deposited in the EEMC pre-shower 1 and 2
(within a 3x3 clusters around tower with a maximum energy).
Figure 2: F_data vs D_tail^max
Note:This plot is fit independend (only the peak position is defined based on the fit)
Figure 3: F_data vs D_tail^max-D_tail^max
Figure 4: Gamma transverse momentum vs jet transverse momentum
Figure 5: F_peak vs maximum residual
for various cuts on energy deposited in the EEMC pre-shower 1 and 2
(within a 3x3 clusters around tower with a maximum energy).
Figure 6: F_data vs D_tail^max
Note:This plot is fit independend (only the peak position is defined based on the fit)
Figure 7: F_data vs D_tail^max-D_tail^max
Figure 8: Gamma transverse momentum vs jet transverse momentum
Figure 9: Gamma pseudorapidity vs jet pseudorapidity
Figure 10: Gamma azimuthal angle vs jet azimuthal angle
Note: for the case of Pre1>1 && Pre2==0 there is an enhancement around phi_gamma = 0?
Figure 11: maximum strip in v-plane vs maximum strip in u-plane
Figure 12:Chi2/ndf for gamma-jet candidates using Monte Carlo shape (combined for [u+v]/2 plane )
Figure 13:Chi2/ndf for gamma-jet candidates (combined for [u+v]/2 plane ) using Will's shape
Ilya Selyuzhenkov March 28, 2008
One interpretation of this can be that in Monte Carlo simulations
the contribution from the material in front of the detector is underestimated
Triple Gaussian fit gives a better description of the data shapes,
compared to the double Gaussian function (compare red and black lines on Figure 1-4)
Figure 1: EEMC SMD shape comparison for various preshower cuts
(black points shows u-plane shape only, v-plane results can be found here)
Figure 2: EEMC SMD shape comparison for various preshower cuts with gamma-jet pt cut of 7GeV
(black points shows u-plane shape only, v-plane results can be found here)
Figure 3: Shapes with an additional cut on 2-gamma candidates within pi0 invariant mass range.
Sample invariant mass distribution using "simple" pi0 finder can be found here
(black points shows u-plane shape only, v-plane results can be found here)
Figure 4: Shapes for the candidates when "simple" pi0 finder failed to find a second peak
(black points shows u-plane shape only, v-plane results can be found here)
Figure 5: Strip by strip SMD energy distribution.
Only 12 strips from the right side of the maximum are shown.
Zero strip (first upper left plot) corresponds to the high strip in the shape
Note, that already at the 3rd strip from a peak,
RMS values are comparable to those for a mean, and for a higher strips numbers RMS starts to be bigger that mean.
(results for u-plane only, v-plane results can be found here)
Results for side residual (together with pt, eta, phi distributions) for gamma-jet candidates can be found at this web page
Red histograms on Figures 6-8 shows chi2 distribution from MC-photons (normalized at chi2=1.4)
Blue histograms on Figures 6-8 shows chi2 distribution from MC-pions (normalized at chi2=1.4)
Figure 6: Chi2/ndf for gamma-jet candidates using Monte Carlo shape
Figure 7: Chi2/ndf for gamma-jet candidates using Will's shape (derived from eta candidates based on Weihong's pi0-finder)
Figure 8: Chi2/ndf for gamma-jet candidates using Pibero's shape (derived from eta candidates)
Ilya Selyuzhenkov April 02, 2008
Figure 1: EEMC SMD shape comparison for various preshower cuts
Note:Only MC gamma-jet shape (open red squares) is the same on all plots
Ilya Selyuzhenkov April 02, 2008
Figure 1: Side residual for various cuts on energy deposited in the EEMC pre-shower 1 and 2
No EEMC SMD based cuts
Figure 2: Side residual for various cuts on energy deposited in the EEMC pre-shower 1 and 2
"Simple" pi0 finder can not find a second peak
Figure 3: Side residual for various cuts on energy deposited in the EEMC pre-shower 1 and 2
"Simple" pi0 finder reconstruct the invarian mass within [0.1,0.18] range
Figure 4: Side residual distribution (Projection for side residual in Figs.1-3 on vertical axis)
Figure 5: Signal (green: m < 0) vs background (black, red) separation
Ilya Selyuzhenkov April 02, 2008
Figure 1: Side residual for various cuts on energy deposited in the EEMC pre-shower 1 and 2
No EEMC SMD based cuts
Ilya Selyuzhenkov April 03, 2008
Hi Ilya,
I think you gave up on the chi-squared method too quickly, and am sorry I missed the phone meeting last week. So, I would like to make a request that will hopefully take a minimal amount of your time to show that all is okay. Then, if there is a delay in getting the sided residual information out and into the beam use request, you can still fall back on the chi-squared method.
In your March 28 posting, Figure 8 at the bottom, I would like to get numerical values for the events per bin for the black curves. I won't use the preshower1>0 and preshower2=0 data, so those you don't need to send. Also, I won't use the red or blue curve information.
I think your problem has been that you normalized your curves at chi-squared/ndf = 1.4 instead of the peak. What I plan to do is to normalize the (pre1=0, pre2=0) to the (pre1=0, pre2>0) data in the peak and subtract. The (pre1=0, pre2=0) set should have some single photons, but also some multiple photons. The (pre1=0, pre2>0) should also have single photons, and more multiple photons, since the chance that one of them will convert is larger. The difference should look roughly like your blue curve, though perhaps not exactly if Pibero's mean shower shape is not perfect (which it isn't). I will do the same thing with taking the difference between (pre1>0, pre2>0) and (pre1=0, pre2=0), and again the difference should look roughly like your blue curve. The (pre1>0, pre2>0) data should have even larger fraction of multiple photons than either of the other two data sets. I would expect the two difference curves to look approximately the same.
Hope this is possible for you to do. Since our reduced chi-squared curve looks so much like the one from CDF, I am pretty confident that we are okay, but this should be checked to convince people that we are not doing anything terribly wrong.
Dear Hal,
I have tried to implement your idea and produce a figure attached.
There are 4 plots in it:
1. Upper left plot shows normalized to unity (at maximum) chi2 distribution (obtained with Pibero shape for gamma-jet candidates) for a different pre1, pre2 conditions
2. Upper right plot shows bin-by-bin difference: a) between normalized chi2 for pre1=0, pre2>0 and pre1=0, pre2=0 (red) and b) between normalized chi2 for pre1>0, pre2>0 and pre1=0, pre2=0 (blue)
3. Bottom left Same as upper right, but normalization were done based on the integral within [-4,4] bins around maximum.
4. Bottom right Same as for upper right, but with a different normalization ([-4,4] bins around maximum)
I have also tried to normalized by the total integral, but the results looks similar.
Figure 1: See description above
Ilya Selyuzhenkov April 09, 2008
Data sample:
Monte-Carlo gamma-jet sample for partonic pt range of 5-7, 7-9, 9-11,11-15, 15-25, 25-35 GeV.
Analysis: Simulated MuDst files were first processed through jet finder algorithm (thanks to Renee Fatemi),
and later analyzed by applying gamma-jet isolation cuts (see this link for details) and studying EEMC SMD response (see below).
To test the algorithm, Geant records were not used in this analysis.
Further studies based on Geant records (yield estimates, etc) are ongoing.
Figure 1:Comparison between shower shape profile for data and MC.
Black circles shows results for MC gamma-jet sample (all partonic pt).
For v-plane results see this figure
Figure 2:Gamma vs jet transverse momentum.
Figure 3:Gamma vs jet azimuthal angle.
Figure 4:Gamma vs jet pseudo-rapidity.
Definitions for F_peak, D_peak, D_tail^max (D_tail^min) can be found here
Figure 5:F_peak vs maximum residual
for various cuts on energy deposited in the EEMC pre-shower 1 and 2
(within a 3x3 clusters around tower with a maximum energy).
Shower shape used to fit data is fixed to the shape from the previous gamma-jet study of real events
(see black point on Fig.1 [upper left plot] at this page)
Figure 6: F_peak vs D_tail^max: click here
Figure 7: F_peak vs D_tail^max-D_tail^min: click here
Figure 8:Logarithmic fraction of energy in post shower (3x3 cluster) to the total energy in SMD u- and v-planes
Figure 8a:
Same as figure 8, but for gamma-jet candidates from the real data (no pt cuts).
Logarithmic fraction of energy in post shower (3x3 cluster) to the total energy in SMD u- and v-planes
Figure 8b:
Comparison between gamma-jet candidates from data with different preshower conditions.
Points are normalized in peak to the case of pre1 > 0, pre2 > 0
Logarithmic fraction of energy in post shower (3x3 cluster) to the total energy in SMD u- and v-planes
Figure 8c:
Comparison between gamma-jet candidates from Monte-Carlo simulations with different preshower conditions.
Points are normalized in peak to the case of pre1 > 0, pre2 > 0
Logarithmic fraction of energy in post shower (3x3 cluster) to the total energy in SMD u- and v-planes
Figure 9: Jet neutral energy fraction
Figure 10: High v-strip vs u-strip
Figure 11: energy post shower (3x3 cluster)
Figure 12: Peak energy SMD-u
Figure 13: Peak energy SMD-v
Figure 14: Gamma phi
Figure 15: Gamma pt
Figure 16: Gamma eta
Figure 17: Delta gamma-jet pt
Figure 18: Delta gamma-jet eta
Figure 19: Delta gamma-jet phi
Figure 20:chi2 distribution using "standard" MC shape
Ilya Selyuzhenkov April 16, 2008
Figure 1: Sided residual for raw MC (partonic pt 9-11)
Figure 2: Sided residual for data-driven MC (partonic pt 9-11)
Figure 3: Sided residual for data (pp Longitudinal 2006)
Figure 4: Number of events which passed various cuts (MC data, partonic pt 9-11)
Ilya Selyuzhenkov April 17, 2008
MC data for different pt weigted according to Michael Betancourt web page:
weight = xSection[ptBin] / xSection[max] / nFiles
Figure 1: Sided residual for raw MC (partonic pt 5-35)
(same plot for partonic pt 9-11)
Figure 2: Sided residual for data-driven MC (partonic pt 5-35)
(same plot for partonic pt 9-11)
Figure 3: Sided residual for data (pp Longitudinal 2006)
Figure 4: Sided residual for data (pp Longitudinal 2006)
Figure 5: Sided residual for data (pp Longitudinal 2006)
Figure 6: pt(gamma) from geant record vs
pt(gamma) from energy in 3x3 tower cluster and position for uv-intersection wrt vertex
(same on a linear scale)
Figure 7: pt(gamma) from geant record vs
pt(jet) as found by the jet-finder
Figure 8: gamma pt distribution:
data-driven MC (red) vs gamma-jet candidates from pp2006 longitudinal run (black).
MC distribution normalized to data at maximum for each preshower condition
Figure 9: Number of events which passed various cuts (MC data, partonic pt 5-35)
Red: cuts applied independent
Black: cuts applied sequential from left to right
Ilya Selyuzhenkov April 23, 2008
MC data for different partonic pt are weigted according to Michael Betancourt web page:
weight = xSection[ptBin] / xSection[max] / nFiles
Figure 1:Sided residual for data-driven gamma-jet MC events (partonic pt 5-35)
Figure 2:Sided residual for data-driven jet-jet MC events (partonic pt 3-55)
Figure 3:Sided residual for data (pp Longitudinal 2006)
Figure 4:pt(gamma) vs pt(jet) for data-driven gamma-jet MC events (partonic pt 5-35)
Figure 5:pt(gamma) vs pt(jet) for data-driven jet-jet MC events (partonic pt 3-55)
Figure 6:pt(gamma) vs pt(jet) for data (pp Longitudinal 2006)
Ilya Selyuzhenkov May 05, 2008
Data samples:
Figure 1:pt distribution. MC data are scaled to the same luminosity as data
(Normalization factor: Luminosity * sigma / N_events).
Figure 2:Integrated gamma yield vs pt.
For each pt bin yield is defined as the integral from this pt up to the maximum available pt.
MC data are scaled to the same luminosity as data.
Figure 3:Signal to background ratio (all results divided by the data)
You can find sided residual 2-D plots here
Figure 4:Maximum sided residual for pt_gamma>7GeV; pt_jet>7GeV
Figure 5:Fitted peak for pt_gamma>7GeV; pt_jet>7GeV
Figure 6:Max data tail for pt_gamma>7GeV; pt_jet>7GeV
Figure 7:Max minus min data tails for pt_gamma>7GeV; pt_jet>7GeV
Ilya Selyuzhenkov May 08, 2008
Figure 1:y:x EEMC position for gamma-jet candidates:
Pythia gamma-jet sample (~170K events). Partonic pt range 5-35 GeV.
Figure 2:y:x EEMC position for gamma-jet candidates:
Pythia QCD bg sample (~4M events). Partonic pt range 3-65 GeV.
Figure 3:y:x EEMC position for gamma-jet candidates:
pp2006 (long) data [eemc-http-mb-l2gamma:137641 trigger]
Figure 3b:y:x EEMC position for gamma-jet candidates:
pp2006 (long) data [eemc-http-mb-l2gamma:137641 trigger]
pt cut of 7 GeV for gamma and 5GeV for the away side jet has been applied.
Figure 4:High v-strip vs high u-strip.
Pythia gamma-jet sample (~170K events). Partonic pt range 5-35 GeV.
Figure 5:High v-strip vs high u-strip:
Pythia QCD bg sample (~4M events). Partonic pt range 3-65 GeV.
Figure 6:High v-strip vs high u-strip:
pp2006 (long) data [eemc-http-mb-l2gamma:137641 trigger]
Figure 6b:High v-strip vs high u-strip:
pp2006 (long) data [eemc-http-mb-l2gamma:137641 trigger]
pt cut of 7 GeV for gamma and 5GeV for the away side jet has been applied.
Ilya Selyuzhenkov May 09, 2008
For a three data samples (pp2006 [long], MC gamma-jet, and MC QCD background events)
the EEMC detector eta cut of 1< eta < 1.4 has been applied.
Although a poor statistics available for MC background QCD sample,
the signal to background ratio (red to green line ratio)
getting closer to 1:3 (expected signal to background ratio from Les study).
Figure 1:Gamma pt distribution. MC data are scaled to the same luminosity as data
(Normalization factor: Luminosity * sigma / N_events).
Figure 2:Gamma yield vs pt. MC data are scaled to the same luminosity as data.
Figure 3:Signal to background ratio (MC results are normalized to the data)
Ilya Selyuzhenkov May 14, 2008
Correlation between gamma-candidate 3x3 cluster energy ratio (R_cluster) and
number of EEMC towers in a jet1 can be found here (Fig. 4).
Gamma pt distribution, yield and signal to background ratio plots
for a cut of R_cluster >0.9 can be found here (Figs. 1-3).
Gamma pt distribution, yield and signal to background ratio plots
for a cut of R_cluster >0.99 are shown below in Figs. 1-3.
One can see that by going from R_cluster>0.9 to R_cluster>0.99
improves signal to background ratio from ~ 1:10 to ~ 1:5 for gamma pt>10 GeV
Figure 1:Gamma pt distribution for R_cluster >0.99.
MC results scaled to the same luminosity as data
(Normalization factor: Luminosity * sigma / N_events).
Figure 2:Integrated gamma yield vs pt for R_cluster >0.99
For each pt bin yield is defined as the integral from this pt up to the maximum available pt.
MC results scaled to the same luminosity as data.
Figure 3:Signal to background ratio for R_cluster >0.99 (all results divided by the data)
Compare this figure with that for R_cluster>0.9 (Fig. 3 at this link)
Figure 4: pt asymmetry between gamma and the away side jet (R_cluster >0.9)
for a three data samples (pp2006[long] data, gamma-jet MC, QCD jets background).
pt cut of 7 GeV for both gamma and jet has been applied.
Figure 5: signal to background ratio (R_cluster >0.9)
as a function of pt asymmetry between gamma and the away side jet
pt cut of 7 GeV for both gamma and jet has been applied.
Figure 6: pt asymmetry between gamma and the away side jet (R_cluster >0.99)
for a three data samples (pp2006[long] data, gamma-jet MC, QCD jets background).
pt cut of 7 GeV for both gamma and jet has been applied.
Figure 7: signal to background ratio
as a functio of pt asymmetry between gamma and the away side jet (R_cluster >0.99)
pt cut of 7 GeV for both gamma and jet has been applied.
Figure 8: pt asymmetry between gamma and the away side jet (R_cluster >0.99)
for a three data samples (pp2006[long] data, gamma-jet MC, QCD jets background).
pt cut of 7 GeV for gamma and 5GeV for the away side jet has been applied.
Figure 9: signal to background ratio
as a function of pt asymmetry between gamma and the away side jet (R_cluster >0.99)
pt cut of 7 GeV for gamma and 5GeV for the away side jet has been applied.
Ilya Selyuzhenkov May 15, 2008
Figure 1:Vertex z distribution for pp2006 (long) data [eemc-http-mb-l2gamma:137641 trigger]
Note: In the upper right plot (pre1=0, pre2>0) one can see
a hole in the acceptance in the range bweeeen z_vertex -10 to 30 cm (probably due to SVT construction)
Figure 1b:Vertex z distribution for pp2006 (same as Fig. 1, but on a linear scale)
Figure 2:Vertex z distribution for three different data samples
MC results scaled to the same luminosity as data
Figure 3:Vertex z distribution for three different data samples
pt cut of 7 GeV for gamma and 5GeV for the away side jet has been applied.
Ilya Selyuzhenkov May 20, 2008
Selecting only di-jet events identified by the STAR jet finder algorithm,
with jets pointing opposite in azimuth:
cos(phi_jet1 - phi_jet2) < -0.8
Data sample:
Note, that all shapes are normalized at peak to unity
Figure 1:Shower shapes for different detector eta bins
Figure 2:Shower shapes for different detector phi bins
Figure 3:Shower shapes for different gamma pt bins
Figure 4:Shower shapes for different z-vertex bins
Ilya Selyuzhenkov May 21, 2008
Data sample:
Subset of 441 eta-meson candidates from Will's analysis.
additional QA info (detector eta, pre1, pre2, etc)
has been added to pi0-tree reader script:
/star/institutions/iucf/wwjacobs/newEtas_fromPi0finder/ReadEtaTree.C
pi0 trees from this RCF directory has been used to regenerate etas NTuple:
/star/institutions/iucf/wwjacobs/newEtas_fromPi0finder/out_23/
Some observations:
eta-meson purity within the invariant mass region [0.5, 0.65] is about 72%
Most of the eta-candidates has detector pseudorapidity less or about 1.4,
what may limits applicability of data-driven shower shapes
derived from these candidates for higher pseudo-rapidity region,
where we have most of the background for the gamma-jet
analysis due to lack of TPC tracking
z-vertex distribution is very asymmetric, and peaked around -50cm.
Only a few candidates has a positive z-vertex values.
Figure 1: Eta-meson invariant mass with signal and background fits and ratio (upper left).
Pseudorapidity [detector and wrt vertex] distributions (right top and bottom plots),
vertex z distributions (bottom left)
Figure 2:2D plots for the eta-meson invariant mass vs
azimuthal angle (upper left), pseudorapidity (upper right),
z-vertex (bottom right), and detector pseudorapidity (bottom right)
Ilya Selyuzhenkov May 27, 2008
Figure 1: Shower shapes and triple Gaussian fits for photons from eta-meson
sorted by different conditions of EEMC 1st and 2nd pre-shower layers.
Note: All shapes have been normalized at peak to unity
Triple Gaussian fit parameters:
Pre1=0 Pre2=0
0.669864*exp(-0.5*sq((x-0.46016)/0.574864))+0.272997*exp(-0.5*sq((x-0.46016)/-1.84608))+0.0585682*exp(-0.5*sq((x-0.46016)/5.49802))
Pre1=0 Pre2>0
0.0694729*exp(-0.5*sq((x-0.493468)/5.65413))+0.615724*exp(-0.5*sq((x-0.493468)/0.590723))+0.314777*exp(-0.5*sq((x-0.493468)/2.00192))
Pre1>0 Pre2>0
0.0955638*exp(-0.5*sq((x-0.481197)/5.59675))+0.558661*exp(-0.5*sq((x-0.481197)/0.567596))+0.345896*exp(-0.5*sq((x-0.481197)/1.9914))
Shower shapes comparison between different data sets:
Some observations:
Shapes for gammas from eta-meson decay
are in a good agreement with those from MC gamma-jet sample
(compare red squares with blue triangle in Fig. 2 and 3).
MC gamma-jet shapes obtained by running a full gamma-jet reconstruction algorithm,
and this agreement indicates that we are able to reconstruct gamma shapes
which we put in with data-driven shower shape library.
MC gamma-jet shapes match pp2006 data shapes
for pre1=0 condition, where we expect to be very efficient in background rejection
(compare red squares with black circles in upper plots of Fig. 2 and 3).
This indicates that we are able to reproduce EEMC SMD of direct photons with data-driven Monte-Carlo.
There is no match between Monte-Carlo QCD background jets and pp2006 data
for the case when both pre-shower layer fired (pre1>0 and pre2>0).
(compare green triangles with black circes in bottom right plots of Fig.2 and 3).
This is the region where we know background dominates our gamma-jet candidates.
This shows that we still do not reproduce SMD response for our background events
in our data-driven Monte-Carlo simulations
(note, that in Monte-Carlo we replace SMD response with real shapes for all background photons
the same way we do it for direct gammas).
Figure 2: Shower shapes comparison between different data sets.
Shapes for gamma-jet candidates obtained with the same gamma-jet reconstruction algorithm
for three different data samples (pp2006, gamma-jet and QCD jets MC).
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 3:Same as Fig. 2, but with no cuts on gamma and jet pt.
All shapes are similar to those in Fig. 2 with an additional pt cuts.
Note, that blue triangles are the same as in Fig. 2.
Ilya Selyuzhenkov May 30, 2008
Three data sets:
Figure 1: Gamma eta distribution.
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 2: Gamma pt distribution.
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 3: Gamma phi distribution.
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 4: Away side jet eta distribution.
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 5: Away side jet pt distribution.
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 6: Gamma-jet delta pt distribution.
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 7: Gamma-jet delta eta distribution.
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 8: Gamma-jet delta phi distribution.
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Ilya Selyuzhenkov June 04, 2008
Three data sets:
Figure 1: Gamma candidate EEMC pre-shower 1 energy (3x3 cluster).
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 2: Gamma candidate EEMC pre-shower 2 energy (3x3 cluster).
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 3: Gamma candidate EEMC tower energy (3x3 cluster).
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 4: Gamma candidate EEMC post-shower energy (3x3 cluster).
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 5: Gamma candidate EEMC SMD u-plane energy [5-strip cluster] (Figure for v-plane)
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 6: Total minus gamma candidate (3x3 cluster) energy in EEMC pre-shower 1 layer
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 7: Total minus gamma candidate (3x3 cluster) energy in EEMC pre-shower 2 layer
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 8: Total minus gamma candidate (3x3 cluster) energy in EEMC tower
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 9: Total minus gamma candidate (3x3 cluster) energy in EEMC post-shower layer
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Figure 10: Total (sector) energy minus gamma candidate (5-strip cluster) energy in EEMC SMD[u-v] layer
pt cuts of 7GeV for the gamma and of 5 GeV for the away side jet have been applied.
Ilya Selyuzhenkov June 09, 2008
Three data sets:
Numerical values for different pt-bins from Fig. 1-2
Figure 1: Gamma pt distribution for R_cluster >0.9.
No energy in both pre-shower layer (left plot), and
No energy in pre-shower1 and non-zero energy in pre-shower2 (right plot)
Same figure for R_cluster>0.99 can be found here
Figure 2: Gamma pt distribution for R_cluster >0.9.
No energy in first EEMC pre-shower1 layer (left plot), and
non-zero energy in pre-shower1 (right plot)
For more details (yield, ratios, all pre12 four conditions, etc) see figures 1-3 here.
Figure 3: Gamma pt distribution for R_cluster >0.99.
For more details (yield, ratios, all pre12 four conditions, etc) see figures 1-3 here.
Ilya Selyuzhenkov June 10, 2008
Note: No background subtraction has been done yet
The case of pre-shower1=0 (left plots) roughly has 1:1 signal to background ratio,
while pre-shower1>0 (right plots) have 1:10 ratio (See MC to data comparison for details).
Data sets:
Figure 1: Gamma-jet candidate A_LL vs gamma pt.
Figures for related epsilon_LL and 1/Lum scaled by a factor 10^7
(see pdf/html links above for epsilon_LL and 1/Lum definitions)
Figure 2: Gamma-jet candidate A_LL vs x_gluon.
Figures for related epsilon_LL and 1/Lum scaled by a factor 10^7
Figure 3: Gamma-jet candidate A_LL vs x_quark.
Figures for related epsilon_LL and 1/Lum scaled by a factor 10^7
Figure 4: Gamma-jet candidate A_LL vs away side jet pt.
Figures for related epsilon_LL and 1/Lum scaled by a factor 10^7
Ilya Selyuzhenkov June 18, 2008
Photon-jet reconstruction with the EEMC detector - Part 1: pdf or odp
Data samples (pp2006, MC gJet, MC QCD bg)
and gamma-jet reconstruction algorithm:
Comparing pp2006 with Monte-Carlo simulations scaled to the same luminosity
(EEMC pre-shower sorting):
EEMC SMD shower shapes from different data samples
(pp2006 and data-driven Monte-Carlo):
Sided residual plots: pp2006 vs data-driven Monte-Carlo
(gammas from eta meson: 3 gaussian fits)
Various cuts study:
Some QA plots:
A_LL reconstruction technique:
Work in progress... To do list:
Ilya Selyuzhenkov July 07, 2008
Data sets:
Figure 1: Correlation between 3x3 cluster energy in pre-shower2 vs. pre-shower1 layers
Figure 1a: Distribution of the 3x3 cluster energy in pre-shower1 layer (zoom in for Epre1<0.03 region)
(pp2006 data vs. MC gamma-jet and QCD events)
Figure 2: Shower shapes after pre-shower1 < 5MeV cut.
Shapes are narrower than those without pre1 cut (see Fig. 2)
Figure 3: Gamma pt distribution with pre-shower1 < 5MeV cut.
Compare with distribution withoud pre-shower1 (Fig. 3)
Figure 4: Fitted peak vs. maximum sided residual (no pre-shower1 cuts)
Only points for pp2006 data are shown.
Figure 5: Fitted peak vs. maximum sided residual (after pre-shower1 < 5MeV cut).
Only points for pp2006 data are shown.
Note that distribution for pre1>0,pre2>0 case are narrower
compared to that in Fig.4 (without pre-shower1 cuts).
Figure 6: Distribution of maximum sided residual with pre-shower1 < 5MeV cut.
Ilya Selyuzhenkov July 16, 2008
Three data sets:
My simple gamma-gamma finder is trying to
find a second peaks (clusters) in each SMD u and v planes,
match u and v plane high strip intersections,
and calculate the invaraint mass from associated tower energies (3x3 cluster)
according to the energy sharing between SMD clusters.
Figure 1: Gamma-gamma invariant mass plot.
Only pp2006 data are shown: black: no pt cuts, red: gamma pt>7GeV and jet pt>5 GeV.
Clear pi0 peak in the [0.1,0.2] invariant mass region.
Same data on the log scale
Figure 2: Gamma pt distribution (no inv mass cuts).
Figure 3: Gamma pt distribution (m_invMass<0.11 or no second peak found).
This cut improves signal to background ratio.
Figure 4: Gamma pt distribution (m_invMass>0.11).
Mostly background events.
Figure 5: Shower shapes (no pre1 and no invMass cuts).
Good match between shapes in case of no energy in pre-shower1 layer (pre1=0 case).
Figure 6: Shower shapes (pre1<5MeV, no invMass cuts).
For pre1&2>0 case shapes getting closer to ech other, but still do not match.
Figure 7: Shower shapes (cuts: pre1<5MeV, invMass<0.11 or no second peak found).
Note, the surprising agreement between eta-meson shapes (blue) and data (black).
Figure 8: Invariant mass distribution (MC vs. pp2006 data): no pre1 cut
Figure 9: Invariant mass distribution (MC vs. pp2006 data): pre1<5MeV
Left side is the same as in Figure 8
Figure 10: Invariant mass distribution (MC vs. pp2006 data): pre1>5MeV
Left side plot is empty, since there is no events with [pre1=0 and pre1>5MeV]
Ilya Selyuzhenkov July 22, 2008
Figure 1: Shower shapes: no energy cuts, only 12 strips from peak (left u-plane, right v-plane).
Figure 1a: Shower shapes: no energy cuts, 150 strips from peak (left u-plane, right v-plane).
Figure 2: Shower shapes Energy>8GeV (left u-plane, right v-plane).
Figure 3: Shower shapes Energy<=8GeV (left u-plane, right v-plane).
Figure 8: Number of library candidates per sector.
Figure 9: Transverse momentum vs. energy.
Figure 10: Distance from center of the detector vs. energy.
Ilya Selyuzhenkov July 29, 2008
Data sets:
Latest data-driven shower shape replacement library:
Figure 1: Shower shapes for u-plane [12 strips]
Shower shapes for the library are for the E>8GeV bin.
Figure 2: Shower shapes for v-plane [12 strips]
Figure 3: Shower shapes for u-plane [expanded to 40 strips]
Figure 4: Shower shapes for v-plane [expanded to 40 strips]
Ilya Selyuzhenkov August 14, 2008
Data sets:
Data-driven maker with bug fixed multi-shape replacement:
Figure 1: Shower shapes for u-plane [12 strips]
Shower shapes for the library are for the E>8GeV bin.
Open squares and triangles represents raw [without dd-maker]
MC gamma-jet and QCD background shower shapes respectively
Figure 2: Shower shapes for v-plane [12 strips]
Figure 3: Shower shapes for u-plane [expanded to 40 strips]
Dashed red and green lines represents raw [without dd-maker]
MC gamma-jet and QCD background shower shapes respectively
Figure 4: Shower shapes for v-plane [expanded to 40 strips]
Ilya Selyuzhenkov August 19, 2008
Data sets:
Note: Due to lack of statistics for 2008 data, no pt cuts applied on gamma-jet candidates for both 2006 and 2008 date.
Figure 1: Shower shapes for u-plane [pp2006 data: eemc-http-mb-l2gamma trigger]
Figure 2: Shower shapes for v-plane [pp2006 data: eemc-http-mb-l2gamma trigger]
Figure 3: Shower shapes for u-plane [pp2008 data: fmsslow trigger]
Figure 4: Shower shapes for v-plane [pp2008 data: fmsslow trigger]
Ilya Selyuzhenkov August 25, 2008
Data sets:
Event selection:
Figure 1: Vertex z distribution (left: pp2008; right: 2006 data)
Figure 2: eta vs. phi distribution for the jet1 (jet with largest NEF) .
Figure 3: eta vs. z-vertex distribution for the jet1 (jet with largest NEF) .
Figure 4: eta vs. z-vertex distribution for the second jet.
Figure 5: Transverse momentum distribution for jet1.
Ilya Selyuzhenkov August 26, 2008
Data sets:
Data-driven library:
Figure 1: Pre-shower1 energy distribution for Pre1>4 MeV:
Eta meson library for E>8GeV bin [left] and data vs. MC results [right].
Figure 2: Shower shapes for v-plane [Pre1<10MeV cut]
Figure 3: Shower shapes for u-plane [Pre1<10MeV cut]
Definitions for side residual plot (F_peak, F_tal, D_tail) can be found here
For a moment same 3-gaussian shape is used to fit SMD response for all pre-shower bins.
Algo needs to be updated with a new shapes sorted by pre-shower bins.
Figure 4: Sided residual plot for pp2006 data only [Pre1<10MeV cut]
Figure 5: Sided residual projection on "Fitted Peak" axis [Pre1<10MeV cut]
Figure 6: Sided residual projection on "tail difference" axis [Pre1<10MeV cut]
Ilya Selyuzhenkov August 27, 2008
Data sets:
Figure 1: High u vs. v strip id distribution for different pre-shower conditions.
Left column: QCD jets, middle column: gamma-jet, right columnt: pp2006 data
Figure 2: x vs. y position of the gamma-candidate within EEMC detector
for different pre-shower conditions.
Left column: QCD jets, middle column: gamma-jet, right columnt: pp2006 data
Figure 3:Reconstructed vs. generated (from geant record) gamma pt for the MC gamma-jet sample.
Pre-shower1<10MeV cut applied.
Figure 4:Reconstructed vs. generated (from geant record) gamma eta for the MC gamma-jet sample.
Pre-shower1<10MeV cut applied.
Figure 5:Reconstructed vs. generated (from geant record) gamma phi for the MC gamma-jet sample.
Pre-shower1<10MeV cut applied.
Ilya Selyuzhenkov September 02, 2008
Data sets:
Shower shape fitting procedure:
Fitting function "[0]*(exp ( -0.5*((x-[1])/[2])**2 )+[3]*exp ( -0.5*((x-[4])/[5])**2 )+[6]*exp ( -0.5*((x-[7])/[8])**2 ))"
Figure 1: MC gamma-jet shower shapes and fits for u-plane
Results from single, double and triple Gaussian fits (using from 9 to 15 strips) are shown.
Figure 2: Same as figure 1. but from v-plane
Figure 3: MC gamma-jet results using triple Gaussian fits within 12 strips from a peak.
Left: u-plane. Right: v-plane
Figure 4: Combined fit results from MC gamma-jet sample
Figure 5: Fitting parameters [see equation for the fit function above].
Note, that parameters 1, 4, and 7 (peak position) has the same value.
Numerical fit results:
Figure 6: Same as Fig. 3, but for gamma-jet candidates from pp2006 data
Figure 7: Same as Fig. 5, but for gamma-jet candidates from pp2006 data
Ilya Selyuzhenkov September 09, 2008
Data sets:
Procedure to calculate maximum sided residual:
For each event fit SMD u and v energy distributions with
triple Gaussian functions from shower shapes analysis:
[0]*(exp(-0.5*((x-[1])/[2])**2)+[3]*exp(-0.5*((x-[1])/[4])**2)+[6]*exp(-0.5*((x-[1])/[5])**2))
Fit parameters sorted by various pre-shower conditions and u and v-planes can be found here
There are only two free parameters in a final fit: overall amplitude [0] and mean value [1]
Fit range is +-2 strips from the high strip (5 strips total).
Integrate energy from a fit within +-2 strips from high strip.
This is our peak energy from fit, F_peak.
Calculate tail energies on left and right sides from the peak for both data, D_tail, and fit, F_tail.
Tails are integrated up to 30 strips excluding 5 highest strips.
Determine maximum difference between D_tail and F_tail:
max(D_tail-F_tail). This is our maximum sided residual.
Plot F_peak vs. max(D_tail-F_tail). This is sided residual plot.
(implementation for this item is in progress)
Based on MC gamma-jet sided residual plot find a line (some polynomial function)
which will serve as a cut to separate signal and background.
Use that cut line to calculate signal to background ratio
and apply it for the real data analysis.
Figure 1: Maximum sided residual plots for different data sets and various pre-shower condition.
Columns [data sets]: 1. MC QCD background; 2. gamma-jet; 3. pp2006 data
Rows [pre-shower bins]: 1. pre1=0 pre2=0; 2. pre1=0, pre2>0; 3. 0<pre1<4MeV; 4. 4<pre1<10MeV
Results from u and v plane are combined as [U+V]/2
Figure 2: max(D_tail-F_tail) distribution (projection on horizontal axis from Fig.1)
Some observations:
Results for pp2006 and MC gamma-jet are consistent for pre1=0 pre2=0 case (upper left plot)
Results for pp2006 and MC QCD background jets are also in agrees for pre1>0 case (lower left and right plots)
Figure 3: F_peak distribution (projection on vertical axis from Fig.1)
Ilya Selyuzhenkov September 16, 2008
These results are obsolete.
Please use this link instead
Data sets:
Notations used in the plots:
Figure 1: D_peak from [U+V]/2.
Figure 2: U/V asymmetry for D_peak: [U-V]/[U+V]
Figure 3: F_peak from [U+V]/2.
Figure 4: U/V asymmetry for F_peak: [U-V]/[U+V]
Figure 5: (D_peak - F_peak)/D_peak asymmetry
Ilya Selyuzhenkov September 23, 2008
Data sets:
Notations used in the plots:
Figure 1: D_peak from [U+V]/2.
Figure 2: (D_peak - F_peak)/D_peak asymmetry
Ilya Selyuzhenkov September 23, 2008
Figure 1: D_peak vs. [right-left] D_tail
Ilya Selyuzhenkov September 23, 2008
Data sets:
Notations used in the plots:
Figure 1: Maximum sided residual plot
Top get more statistics for MC QCD sample plot is redone with a softer R_cluster > 0.9 cut
Figure 2: D_peak (projection on vertical axis for Fig. 1)
Upper left plot (no pre-shower fired case) reveals some difference
between MC gamma-jet and pp2006 data at lower D_peak values.
This difference could be due to background contribution at low energies.
Still needs more statistics for MC QCD jet sample to confirm that statement.
Figure 3: max(D_tail-F_tail) (projection on horisontal axis for Fig. 1)
One can get an idea of signal/background separation (red vs. black) depending on pre-shower condition.
Figure 4: Mean < max(D_tail-F_tail) > vs. D_peak (profile on vertical axis from Fig. 1)
For gamma-jet sample average sided residual is independent on D_peak energy
and has a slight positive shift for all pre-shower>0 conditions.
For large D_peak values (D_peak>0.16) MC gamma-jet and pp2006 data results are getting close to each other.
This corresponds to higher energy gammas, where we have a better signal/background ratio,
and thus more real gammas among gamma-jet candidates from pp2006 data.
(Note: legend's color coding is wrong, colors scheme is the same as in Fig. 3)
Figure 5: Mean < D_peak > vs. max(D_tail-F_tail) (profile on horisontal axis from Fig. 1)
For "no-preshower fired" case MC gamma-jet sample has a large average values than that from pp2006 data.
This reflects the same difference between pp2006 and MC gamma-jet sample at small D_peak values (see Fig. 2, upper left plot).
(Note: legend's color coding is wrong, colors scheme is the same as in Fig. 3)
Figure 7: D_peak vs. gamma 3x3 tower cluster energy
Figure 8: 3x3 cluster tower energy distribution
Figure 9: Gamma pt distribution
The simplest way to get signal/background separation is to draw a straight line
on sided residual plot (Fig. 1) in such a way that
it will contains most of the counts (signal) on the left side,
and use a distance to that line for both MC and pp2006 data samples
as a discriminant for signal/background separation.
To get the distance to the straight line one can rotate sided residual plot
by the angle which corresponds to the slope of this line,
and then project it on "rotated" max(D_tail-F_tail) axis.
Figure 10: Shows "rotated" sided residual plot by "5/6*(pi/2)" angle (this angle has been picked by eye).
One can see that now most of the counts for gamma-jet sample (middle column)
are on the left side from vertical axis.
Figure 11: "Rotated" max(D_tail-F_tail) [projection on horizontal axis for Fig. 10]
Cut on "Rotated" max(D_tail-F_tail) can be used for signal/background separation.
From figure below one can see much better signal/background separation than in Fig. 3
Figure 12: "Rotated" D_peak [projection on vertical axis for Fig. 10]
Ideally, instead of straight line one needs to use
an actual shape of side residual distribution for MC gamma-jet sample.
This shape can be extracted and parametrized by the following procedure:
The distance to that polynomial function can be used to determine our signal/background rejection efficiency.
This work is in progress...
Just last one figure showing shapes for 6 slices from sided plot.
Figure 13: max(D_tail-F_tail) for different slices in D_peak (scaled by the integral for each slice)
Ilya Selyuzhenkov September 30, 2008
Data sets:
Notations used in the plots:
Figure 1: Sided residual plot: D_peak vs. max(D_tail-F_tail)
Red lines show 4th order polynomial functions, a*x^4,
which have 80% of MC gamma-jet counts on the left side.
These lines are obtained independently for each of pre-shower condition
based on fit procedure shown in Fig. 3 below.
Figure 2: max(D_tail-F_tail) distribution
(projection on horizontal axis from sided residual plot, see Fig. 1 above)
Figure 3: max(D_tail-F_tail) [at 80%] vs. D_peak.
For each slice (bin) in D_peak variable, the max(D_tail-F_tail) value
which has 80% of gamma-jet candidates on the left side are plotted.
Lines represent fits to MC gamma-jet points (shown in red) using different fit functions
(linear, 2nd, 4th order polynomials: see legend for color coding).
Note, that in this plot D_peak values are shown on horizontal axis.
Consequently, to get 2nd order polynomial fit on sided residual plot (Fig. 1),
one needs to use sqrt(D_peak) function.
The same apply to 4th order polynomial function.
Figure 4: D_peak vs. horisontal distance from 4th order polinomial function to max(D_tail-F_tail) values.
(compare with Fig. 1: Now 80% of MC gamma-jet counts are on the left side from vertical axis)
Figure 5: Horizontal distance from 4th order polynomial function to max(D_tail-F_tail)
[Projection on horizontal axis from Fig. 4]
Based on this plot one can obtain purity, efficiency, and rejection plots (see Fig. 6 below)
Horizontal distance plotted in Fig. 5 can be used as a cut
separating gamma-jet signal and QCD-jets background,
and for each value of this distance one can define
gamma-jet purity, efficiency, and QCD-background rejection:
gamma-jet purity is defined as the ratio of
the integral on the left for MC gamma-jet data sample, N[g-jet]_left,
to the sum of the integrals on the left for MC gamma-jet and QCD jets, N[QCD]_left, data samples:
Purity[gamma-jet] = N[g-jet]_left/(N[g-jet]_left+N[QCD]_left)
gamma-jet efficiency is defined as the ratio of
the integral on the left side for MC gamma-jet data sample, N[g-jet]_left,
to the total integral for MC gamma-jet data sample, N[g-jet]:
Efficiency[gamma-jet] = N[g-jet]_left/N[g-jet]
QCD background rejection is defined as the ratio of
the integral on the right side for MC QCD jets data sample, N[QCD]_right,
to the total integral for MC QCD jets data sample, N[QCD]:
Rejection[QCD] = N[QCD]_right/N[QCD]
Figure 6: Shows:
purity[g-jet] vs. efficiency[g-jet] (upper left);
rejection[QCD] vs. efficiency[g-jet] (upper right);
purity[g-jet] vs. rejection[QCD] (lower left);
pp2006 to MC ratio, N[pp2006]/(N[g-jet]+N[QCD]), vs. horizontal distance from Fig. 5 (lower right)
Ilya Selyuzhenkov October 13, 2008
I have finished production of jet trees for Michael's gamma filtered events
You can find jet and skim file lists in my directory at IUCF disk (RCF):
Number of jet events is 1284581 (1020 files).
Production size, including archived log files, is 4.0G.
The script to run jet finder:
/star/institutions/iucf/IlyaSelyuzhenkov/simu/JetTrees/20081008_gJet/StRoot/macros/RunJetSimuSkimFinder.C
JetFinder and JetMaker code:
/star/institutions/iucf/IlyaSelyuzhenkov/simu/JetTrees/20081008_gJet/StRoot/StJetFinder
/star/institutions/iucf/IlyaSelyuzhenkov/simu/JetTrees/20081008_gJet/StRoot/StJetMaker
For more details see these threads of discussions:
Ilya Selyuzhenkov October 14, 2008
Data sets:
Figure 1: Horizontal distance from 4th order polynomial function to max(D_tail-F_tail)
See this page for definition and more details
Figure 2:
purity[g-jet] vs. efficiency[g-jet] (upper left);
rejection[QCD] vs. efficiency[g-jet] (upper right);
purity[g-jet] vs. rejection[QCD] (lower left);
pp2006 to MC ratio, N[pp2006]/(N[g-jet]+N[QCD]), vs. horizontal distance (lower right)
Ilya Selyuzhenkov October 15, 2008
Data sets:
Some observations:
Figure 1a: partonic pt for gamma-jet [old] events
after analysis cuts and partonic pt bin weighting
(Note:Arbitrary absolute scale)
Figure 1b: partonic pt for gamma-jet [gamma-filtered] events after analysis cuts.
Michael's StBetaWeightCalculator has been used to caclulate partonic pt weights
Figure 1c: partonic pt for QCD jets [old] events
after analysis cuts and partonic pt bin weighting
(Note:Arbitrary absolute scale)
Figure 1d: partonic pt for QCD jets [gamma-filtered] events after analysis cuts.
Michael's StBetaWeightCalculator has been used to caclulate partonic pt weights
Figure 2: reconstructed gamma pt: old MC vs. pp2006 data (scaled to the same luminosity)
Figure 3: reconstructed gamma pt: gamma-filtered MC vs. pp2006 data (scaled to the same luminosity)
Ilya Selyuzhenkov October 15, 2008
Data sets:
Figure 1: Horizontal distance from sided residual plot: R_cluster > 0.9
(see Figs. 1-5 from this post for horizontal distance definition)
Figure 2: Purity/efficiency/rejection, and data to MC[gamma-jet+QCD] ratio plots: R_cluster > 0.9
(see text above Fig. 6 from this post for purity, efficiency, and background rejection definition)
Figure 3: Reconstructed gamma pt: R_cluster > 0.98
Figure 4: Horizontal distance from sided residual plot: R_cluster > 0.98
Figure 5: Purity/efficiency/rejection, and data to MC[gamma-jet+QCD] ratio plots: R_cluster > 0.98
Ilya Selyuzhenkov October 21, 2008
Data sets:
Some comments:
Fig. 1-3, upper left plots (pre1=0 pre2=0) show that
average energy per strip in data-driven gamma-jet MC (i.e. solid red square in Fig. 3)
is systematically higher than that for pp2006 data (black circles in Fig. 3).
Note, that there is an agreement between SMD shower shapes
for pp2006 data and data-driven gamma-jet simulations
if one scales them to the same peak value
(Compare red vs. black in upper left plot from Fig. 1 at this link)
Fig. 4, upper left plot (pre1=0 pre2=0):
Integrated SMD energy from 25 strips
in raw gamma-jet simulations (red line) match pp2006 data (black line)
in the region where signal to background ratio is high, E_smd(25-strips)>0.1GeV.
This indicates that raw MC does a good job in
reproducing total energy deposited by direct photon.
Fig. 5, upper left plot (pre1=0 pre2=0):
There is mismatch between distributions of energy in 25 strips cluster
from data-driven gamma-jet simulations and pp2006 data.
This probably reflects the way we scale our library shower shapes
in data-driven shower shape replacement procedure.
Currently, the scaling factor for the library shape is calculated based on the ratio
of direct photon energy from Geant record to the energy of the library photon.
Our library is build out of photons from eta-meson decay,
which has been reconstructed by running pi0 finder.
The purity of the library is about 70% (see Fig. 1 at this post for more details).
The improvement of scaling procedure could be to
preserve total SMD energy deposited within 25 strips from raw MC,
and use that energy to scale shower shapes from the library.
Fig. 6, upper left plot (pre1=0 pre2=0):
Mismatch between integrated 5-strip energy for raw MC and pp2006 in Fig. 6
corresponds to "known" difference in shower shapes from raw Monte-Carlo and real data.
Figure 1: SMD shower shapes: data, raw, and data-driven MC (40 strips).
Vertical axis shows average energy per strip (no overall shower shapes scaling)
Figure 2: Shower shapes: data, raw, and data-driven MC (12 strips)
Figure 3: Shower shapes: data, raw, and data-driven MC (5 strips)
Figure 4: 25 strips SMD cluster energy for raw Monte-Carlo
(Note: type in x-axis lables, should be "25 strip peak" instead of 5)
Figure 5: 25 strips SMD cluster energy for data-driven Monte-Carlo
Figure 6: 5 strips SMD peak energy for raw Monte-Carlo
Figure 7: 5 strips SMD peak energy for data-driven Monte-Carlo
Figure 8:Energy from the right tail (up to 30 strips) for raw Monte-Carlo
Figure 9:Energy from the right tail (up to 30 strips) for data-driven Monte-Carlo
Ilya Selyuzhenkov October 27, 2008
Data sets:
Shower shapes scaling options in data-driven maker:
In all figures below (exept for pp2006 data and raw Monte-Carlo)
the SMD based shower shape scaling has been used.
Figure 1: SMD shower shapes: data, raw, and data-driven MC (40 strips).
Vertical axis shows average energy per strip (no overall shower shapes scaling)
Figure 2: Shower shapes: data, raw, and data-driven MC (12 strips)
Figure 3: Shower shapes: data, raw, and data-driven MC (5 strips)
Figure 4: 25 strips SMD cluster energy for data-driven Monte-Carlo
(SMD based shower shape scaling)
Figure 5: 25 strips SMD cluster energy for raw Monte-Carlo
Note, the difference between results in Fig. 4 and 5. for MC gamma-jets (shown in red)
at low energy (Esmd < 0.04) for pre1=0 pre2=0 case.
This effect is due to the "Number of strips fired in 5-strips cluster > 3" cut.
In data-driven Monte-Carlo we may have shower shapes
with small number of strips fired (rejected in raw Monte-Carlo)
to be replaced by library shape with different (bigger) number of strips fired.
This mostly affects photons which starts to shower
later in the detector and only fires few strips (pre1=0 pre2=0 case)
Ilya Selyuzhenkov October 30, 2008
Below are links to drupal pages
with various SMD energy distributions and shower shapes
for the following set of cuts/conditions:
Figure 1: Case A
Figure 2:Case B
Figure 3:Case C
Figure 4:Case D
Figure 1: Case A
Figure 2:Case B
Figure 3:Case C
Figure 4:Case D
Figure 1: Case A
Figure 2: Case B
Figure 3: Case C
Figure 4: Case D
Figure 1: Case A
Figure 2:Case B
Figure 3:Case C
Figure 4:Case D
Figure 1: Case A
Figure 2:Case B
Figure 3:Case C
Figure 4:Case D
Figure 1: Case A
Figure 2:Case B
Figure 3:Case C
Figure 4:Case D
Figure 1: Case A
Figure 2:Case B
Figure 3:Case C
Figure 4:Case D
Figure 1: Case A
Figure 2:Case B
Figure 3:Case C
Figure 4:Case D
Ilya Selyuzhenkov November 06, 2008
Ilya Selyuzhenkov November 11, 2008
Data sets:
Figure 1: Reconstructed gamma pt for di-jet events and
Geant cuts: pt_gamma[Geant] > 7GeV and 1.05 < eta_gamma[Geant] < 2.0
Total integral for the histogram is: N_total = 5284
(after weighting different partonic pt bins and scaled to 3.164pb^-1).
Compare with number from Jim Sowinski study for
Endcap East+West gamma-jet and pt>7 GeV: N_Jim = 5472
( Jim's numbers are scaled to 3.164pb^-1 : [2539+5936]*3.164/4.9)
Figure 2: Reconstructed jet pt for di-jet events and the same cuts as in Fig. 1
List of cuts (sorted by bin number in Figs. 2 and 3):
Figure 3: Number of accepted events vs. various analysis cuts
The starting number of events (shown in first bin of the plots) is
the number of di-jets with reconstructed gamma_pt>7 GeV and jet_pt>5 GeV
upper left: cuts applied independently
upper right: cuts applied sequentially
lower left: ratio of pp2006 data vs. MC sum of gamma-jet and QCD-jets events (cuts applied independently)
lower right:ratio of pp2006 data vs. MC sum of gamma-jet and QCD jets events (cuts applied sequentially)
Figure 4: Number of accepted events vs. various analysis cuts
Data from Fig. 3 (upper plots) scaled to the initial number of events from first bin:
left: cuts applied independently
right: cuts applied sequentially
Ilya Selyuzhenkov November 18, 2008
all gamma-jet candidate selection cuts except 3x3/r=0.7 energy isolation cut
Figure 4b: 3x3/0.7 ratio but only using towers which passed jet finder threshold
Ilya Selyuzhenkov November 21, 2008
all gamma-jet candidate selection cuts except 3x3/r=0.7 energy isolation cut
There are two sets of figures in links below:
Yield is defined as the integral above given energy fraction
up to the maximum value of 1
Gamma candidate detector eta < 1.5
(eta region where we do have most of the TPC tracking):
Gamma candidate detector eta > 1.5:
(smaller tower size)
Figure 1a: 2x1/0.7 energy fraction [number of counts per given fraction]
Figure 2a: 2x2/0.7 energy fraction [number of counts per given fraction]
Figure 3a: 3x3/0.7 energy fraction [number of counts per given fraction]
Figure 4a: 2x1/0.7 energy fraction [yield]
Figure 1a: 2x1/0.7 energy fraction [number of counts per given fraction]
Figure 2a: 2x2/0.7 energy fraction [number of counts per given fraction]
Figure 3a: 3x3/0.7 energy fraction [number of counts per given fraction]
Figure 4a: 2x1/0.7 energy fraction [yield]
Figure 1a: 2x1/3x3 energy fraction [number of counts per given fraction]
Figure 2a: 2x2/3x3 energy fraction [number of counts per given fraction]
Figure 1a: 2x1/3x3 energy fraction [number of counts per given fraction]
Figure 2a: 2x2/3x3 energy fraction [number of counts per given fraction]
Ilya Selyuzhenkov November 25, 2008
Fig.2 [lower right, 5th bin] shows that
charge particle veto also independent from other cuts
List of cuts (sorted according to bin number in Figs. 1-3. [No SMD sided residual cuts]):
Figure 1: Number of accepted events vs. various analysis cuts
The starting number of events (shown in first bin of the plots) is
the number of di-jets with reconstructed gamma_pt>7 GeV and jet_pt>5 GeV
upper left: cuts applied independently
upper right: expept this cut fired
(event passed all other cuts and being rejected by this cut)
lower left: "cuts applied independently" normalized by the total number of events
lower right: "expept this cut fired" normalized by the total number of events
Figure 2: Same as Fig.1 except: no R_cluster cut and photon detector eta < 1.5
(eta region where we do have most of the TPC tracking)
Figure 3: Same as Fig.1 except: no R_cluster cut and photon detector eta > 1.5
Ilya Selyuzhenkov December 08, 2008
Figure 1: EEMC x vs. y position of photon candidate for 2008 data sample
Problem with pre-shower layer in Sector 10 can been seen in the upper left corner
Figure 2: EEMC x vs. y position of photon candidate for 2006 data sample
Figure 3: Average < E_pre1 * E_pre2 > for 3x3 cluster around high tower
vs run number for sectors 9, 10 and 11
Note, zero pre-shower energy for sector 10 (black points) for days 61, 62, 64, and 67.
All di-jet events for pp2008 data are shown (no gamma-jet cuts)
Figure 3a: Same as Fig.3, zoom into day 61
Figure 3b: Same as Fig.3, zoom into day 62
Figure 3c: Same as Fig.3, zoom into day 64
Figure 3d: Same as Fig.3, zoom into day 67
Figure 4: EEMC x vs. y position of photon candidate for 2008 data sample
Same as Fig. 1, but excluding days: 61, 62, 64, and 67
No problem with pp2008 data have been found,
except that for some runs (mostly on days 61, 62, 64, and 67)
EEMC pre-shower layer for sector 10 was off.
Figure 5: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 6: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 5):
Figure 7: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 5, no scaling):
Ilya Selyuzhenkov December 09, 2008
Ilya Selyuzhenkov December 09, 2008
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
Ilya Selyuzhenkov December 09, 2008
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
Ilya Selyuzhenkov December 09, 2008
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
Ilya Selyuzhenkov December 09, 2008
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
Ilya Selyuzhenkov December 09, 2008
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
Ilya Selyuzhenkov December 11, 2008
Presentation in pdf or open office file format
Ilya Selyuzhenkov December 16, 2008
Gamma-jet isolation cuts except 3x3/r=0.7 energy isolation cut
Figure 1: (reproducing old results with dd-maker fix)
transverse momentum and vertex z distributions
before (with ideal gains/pedestals) and
after (with realistic gains/pedestal tables) dd-maker fix are in a good agreement.
For details on "dd-maker problem", read these hyper news threads:
emc2:2905, emc2:2900, and phana:294
Now we can run L2gamma trigger emulation and Eemc SMD ddMaker
in the same analysis chain.
Figure 2: Vertex z distribution with and without trigger condition in simulations
(emulated trigger: eemc-http-mb-L2gamma [id:137641]).
Solid red/green symbols show results with l2gamma condition applied,
while red/green lines show results for the same analysis cuts but without trigger condition.
Note, good agreement between MC QCD jets with trigger condition on (green solid squared)
and pp2006 data (black solid circles) for pre-shower1>0 case.
Figure 3: pt distribution with/without trigger condition in simulations.
Same color coding as in Fig. 2
Figure 4: Same as Fig. 3 just on a log scale
One can clearly see large trigger effect when applied for QCD jet events,
and a little effect for direct gammas.
Figure 5: gamma candidate pt QCD (right) and prompt photon (left) Monte-Carlo:
no (upper) with (lower) L2e-gamma trigger condition
No photon pt and no jet pt cuts
Figure 6: gamma candidate pt for QCD Monte-Carlo: no L2e-gamma trigger condition
No photon pt and no jet pt cuts
Figure 7: gamma candidate pt for QCD Monte-Carlo: L2e-gamma trigger condition applied (id:137641)
No photon pt and no jet pt cuts
Ilya Selyuzhenkov December 19, 2008
Gamma-jet isolation cuts except 3x3/r=0.7 energy isolation cut
Figure 1: Parton pt distibution for gamma-jet candidates from Pythia QCD sample
with various pt and l2gamma trigger conditions
Figure 2: Parton pt distibution for gamma-jet candidates from Pythia prompt photon sample
with various pt and l2gamma trigger conditions
Ilya Selyuzhenkov January 08, 2009
(concentrated on pre-shower1>0 case
which has better statistics for QCD Monte-Carlo):
Figure 1: Vertex z distribution
Figure 4: gamma-jet pt asymmetry: (pt_gamma - pt_jet)/pt_gamma
Figure 5: gamma pt vs. away side jet pt
1st column: triggered pp2006 data
2nd column: gamma-jet MC (l2gamma trigger on)
3rd column: QCD background MC (l2gamma trigger on)
Ilya Selyuzhenkov January 20, 2009
Note:
this is an update with 10x more statitstics for QCD 9-15GeV parton pt bin.
See this post for old results.
Figure 1: Vertex z distribution
Figure 4: gamma-jet pt asymmetry: (pt_gamma - pt_jet)/pt_gamma
Figure 5: gamma pt vs. away side jet pt
1st column: triggered pp2006 data
2nd column: gamma-jet MC (l2gamma trigger on)
3rd column: QCD background MC (l2gamma trigger on)
Figure 6: pp2006 data to Monte -Carlo sum [QCD + gamma-jet] ratio
for pre-shower1>4MeV (most of statistics)
Left: data to MC ratio vs. reconstructed gamma pt.
Solid line shows constant line fit (p0 ~ 1.3)
Right: data to MC ratio vs. reconstructed vertex position
Ilya Selyuzhenkov January 27, 2009
All figures:
Ilya Selyuzhenkov January 27, 2009
All figures:
Figure 1: Vertex z distribution
Figure 4:Photon pt
Same as in Fig.4 on a log scale: no gamma pt cut and pt>7GeV
Figure 5: Away side jet pt
Same as in Fig.5 on a log scale: no gamma pt cut and pt>7GeV
Ilya Selyuzhenkov February 02, 2009
All figures:
Figure 1: Vertex z distribution with pt_jet>7 cut (left) and without pt_jet cut (rigth)
Figure 2: Photon (left) and away side jet (right) pt
Figure 3: Photon detector eta (left) and corrected for vertex eta (right)
Figure 4: Away side jet detector eta (left) and corrected for vertex eta (right)
Ilya Selyuzhenkov February 03, 2009
Each figure has:
Figure 1: Vertex z distribution
Figure 3: Corrected for vetrex photon eta
Figure 4: Away side jet detector eta
Ilya Selyuzhenkov February 06, 2009
Partonic pt range 2-25 GeV.
Each figure has:
Note: There is a "pre-shower sector 10 problem" for pp2008 data,
which results in migration of small fraction of events with pre-shower>0 into
pre-shower=0 bin (first zero bins in Fig.1 and 2. below).
For pre-shower>0 case this only affects overall normalization of pp2008 data,
but not the shape of pre-shower energy distributions.
I'm running jet-finder+my software to get more statistics from pp2008 data set,
and after more QA will produce list of runs with "pre-shower sector 10 problem",
so to exclude them in the next iteration of my plots.
Ilya Selyuzhenkov February 06, 2009
Partonic pt range 2-25 GeV.
Each figure has:
Figure 6: gamma-jet pt balance
Figure 7: Photon neutral energy fraction
Figure 8: Jet neutral energy fraction
Figure 9: cos(phi_gamma-phi_jet)
Figure 13: Number of charged track associated with photon candidate
Figure 14: Number of Barrel towers associated with photon candidate
Figure 15: Number of Endcap towers associated with photon candidate
Ilya Selyuzhenkov February 16, 2009
Partonic pt range 2-25 GeV.
Each figure has:
Slides: download pdf
Link for CIPANP abstract
Title:
"Photon-jet coincidence measurements
in polarized pp collisions at sqrt{s}=200GeV
with the STAR Endcap Calorimeter"
Abstract: download pdf
Previous versions: v1, v2, v3, v4
Conference link: CIPANP 2009
Multilayer perceptron (MLP) is feedforward neural networks
trained with the standard backpropagation algorithm.
They are supervised networks so they require a desired response to be trained.
They learn how to transform input data into a desired response,
so they are widely used for pattern classification.
With one or two hidden layers, they can approximate virtually any input-output map.
They have been shown to approximate the performance of optimal statistical classifiers in difficult problems.
TMultiLayerPerceptron class in ROOT
mlpHiggs.C example
Netwrok structure:
r3x3, (pt_gamma-pt_jet)/pt_gamma, nCharge, bBtow, eTow2x1: 10 hidden layers: one output later
Figure 1:
Figure 2: Input parameters vs. network output
Row: 1: MC QCD, 2: gamma-jet, 3 pp2006 data
Vertical axis: r3x3, (pt_gamma-pt_jet)/pt_gamma, nCharge, bBtow, eTow2x1
Horisontal axis: network output
ROOT implementation for LDA and MLP:
LDA configuration: default
MLP configuration:
Input parameters (same for both LDA and MLP):
Figure 1: Signal efficiency and purity, background rejection (left),
and significance: Sig/sqrt[Sig+Bg] (right) vs. LDA (upper plots) and MLP (lower plots) classifier discriminants
Figure 2:
Figure 3: Data to Monte-Carlo comparison for LDA (upper plots) and MLP (lower plots)
Good (within ~ 10%) match between data nad Monte-Carlo
a) up to 0.8 for LDA discriminant, and b) up to -0.7 for MLP.
Figure 4: Data to Monte-Carlo comparison for input parameters
from left to right
1) pt_gamma 2) pt_jet 3) r3x3 4) gamma-jet pt balance 5) N_ch[gamma] 6) N_eTow[gamma] 7) N_bTow[gamma]
Colour coding: black pp2006 data, red gamma-jet MC, green QCD MC, blue gamma-jet+QCD
Figure 5: Data to Monte-Carlo comparison:
correlations between input variables (in the same order as in Fig. 4)
and LDA classifier discriminant (horizontal axis).
1st raw: QCD MC; 2nd: gamma-jet MC; 3rd: pp2006 data; 4th: QCD+gamma-jet MC
Endcap photon-jet update at the STAR Collaboration meeting
April 2009 posts
The STAR spin program with longitudinally polarized proton beams
ROOT implementation for LDA:
LDA configuration: default
LDA input parameters:
Figure 1: LDA discriminant (no SMD involved in training)
Figure 2: LDA (no SMD): Efficiency, rejection, purity vs. discriminant
Figure 3: SMD energy in 25 central strips (LDA-dsicriminant>0, no pre-shower1 cut)
Figure 4: SMD energy in 25 central strips (LDA-dsicriminant>0, pre-shower1 < 10MeV)
Figure 5: Maximum residual (LDA-dsicriminant>0, no pre-shower1 cut)
Figure 6: Maximum residual (LDA-dsicriminant>0, pre-shower1 < 10MeV)
SMD info added:
a) energy in 5 central srtips
b) maximum sided residual
Figure 7:LDA with SMD: Efficiency, rejection, purity vs. LDA discriminant
Figure 8: LDA discriminant with SMD
Figure 9: Maximum residual (SMD LDA-dsicriminant>0, pre-shower1 < 10MeV)
Figure 10:LDA (no SMD): Efficiency, rejection, purity plots
Figure 11: LDA with SMD: Efficiency, rejection, purity plots
ROOT implementation for LDA:
LDA configuration: default
LDA input parameters (includes SMD inromation of the distance from max sided residual plot):
The number of strips in SMD u or v planes is required to be greater than 3
Figure 1: SMD energy in 25 central strips sorted by pre-shower energy
Right plot for each pre-shower condition shows the ratio of pp2006 data to sum of the Monte-Carlo samples
Colour coding:
black pp2006 data, red gamma-jet MC, green QCD MC, blue gamma-jet+QCD
(combined plot for all pre-shoer bins can be found here)
Figure 2: SMD energy in 5 central strips sorted by pre-shower energy
(combined plot can be found here)
Figure 3: Maximum residual sorted by pre-shower energy
(combined plot can be found here)
Figure 4: LDA discriminant. Note: LDA algo trained for each pre-shower condition independently
Figure 5: LDA: Efficiency, rejection, purity vs. discriminant, sorted by pre-shower energy
Figure 6: LDA: Efficiency, rejection, purity plots sorted by pre-shower energy
For each pre-shower condition each plot has 4 figures:
ROOT implementation for LDA:
LDA configuration: default
LDA input parameters Set0:
The number of strips in SMD u or v planes is required to be greater than 3
Pre-shower sorting (energy in tiles under 3x3 tower patch):
Photon pt and rapidity cuts:
Figure 0: photon pt distribution for pre-shower1<0.01
Colour coding:
black pp2006 data, red gamma-jet MC, green QCD MC, blue gamma-jet+QCD
Figure 1: LDA discriminant with Set0: Data to Monte-Carlo comparison (pt>7GeV cut)
Right plot for each pre-shower condition shows the ratio of pp2006 data to sum of the Monte-Carlo samples
Colour coding:
black pp2006 data, red gamma-jet MC, green QCD MC, blue gamma-jet+QCD
Figure 2: efficiency, purity, rejection vs. LDA discriminant (pt>7GeV cut)
Figure 3: rejection vs. efficiency
Figure 4: purity vs. efficiency
Figure 5: purity vs. rejection
Figure 6: LDA discriminant with Set1: Data to Monte-Carlo comparison
Figure 7: rejection vs. efficiency
Figure 8: purity vs. efficiency
Figure 9: purity vs. rejection (click link to see the figure)
Figure 10: rejection vs. efficiency (click link to see the figure)
Figure 11: purity vs. efficiency
Figure 12: purity vs. rejection (click link to see the figure)
Figure 13: rejection vs. efficiency (click link to see the figure)
Figure 14: purity vs. efficiency
Figure 15: purity vs. rejection (click link to see the figure)
ROOT implementation for LDA:
LDA configuration: default
LDA input parameters Set0:
The number of strips in SMD u or v planes is required to be greater than 3
Pre-shower sorting (energy in tiles under 3x3 tower patch):
Integrated yields per pre-shower bin:
sample | total integral | pre1=0,pre2=0 | pre1=0, pre2>0 | 0 < pre1 < 0.004 | 0.004 < pre1 < 0.01 | pre1 < 0.01 | pre1 >= 0.01 |
photon-jet | 2.5640e+03 | 3.5034e+02 | 5.2041e+02 | 5.6741e+02 | 5.2619e+02 | 1.9644e+03 | 5.9994e+02 |
QCD | 5.6345e+04 | 1.3515e+03 | 4.3010e+03 | 1.2289e+04 | 1.5759e+04 | 3.3701e+04 | 2.2644e+04 |
pp2006 | 6.2811e+04 | 6.8000e+02 | 2.4310e+03 | 1.2195e+04 | 1.6766e+04 | 3.2072e+04 | 3.0739e+04 |
Photon pt and rapidity cuts:
Figure 1: LDA discriminant with Set0: Data to Monte-Carlo comparison (pt>7GeV cut)
Right plot for each pre-shower condition shows the ratio of pp2006 data to sum of the Monte-Carlo samples
Colour coding:
black pp2006 data, red gamma-jet MC, green QCD MC, blue gamma-jet+QCD
Figure 2: rejection vs. efficiency
Figure 3: purity vs. efficiency
Figure 4: purity vs. rejection
Figure 5: Correlation matrix (pt>7GeV cut)
pre1=0, pre2=0
pre1=0, pre2>0
0 < pre1 < 0.004
0.004 < pre1 < 0.01
pre1 < 0.01
pre1 >= 0.01
For this post LDA input parameters Set4 has been used
LDA for various pre-shower bins is trained independetly,
and later results with pre-shower1<0.01 are combined.
There are a set of plots for various photon pt cuts (pt> 7, 8, 9 10 GeV)
and with different selection of cutoff for LDA
(either based on purity or efficiency).
Number in brackets shows the total yield for the sample.
Link to all plots (16 total) as a single pdf file
Figure 1: pt > 7GeV, efficiency@70
Figure 2: pt > 7GeV, purity@35
Figure 3: pt > 7GeV, purity@40
Figure 4: pt > 7GeV, purity@25 (Note: very similar to results with efficiency@70)
Figure 5: pt > 9GeV, efficiency@70
Figure 6: pt > 9GeV, purity@35
Figure 7: pt > 10GeV, efficiency@70
Figure 8: pt > 10GeV, purity@40
(analysis status update for Spin PWG)
Slides in pdf format:
For this post LDA results with Set1 and Set2 has been used
Note, that LDA for various pre-shower bins is trained independetly
pdf-links with results for pre1=0 and pre2=0 (pre-shower bin 1):
Figures below are for 0.004<pre-shower1<0.01 (pre-shower bin 4).
Photon pt cut: pt> 7, pre-shower bin: 0.004 < pre1 < 0.01
LDA cut with efficiency @ 70%
What is added in Set2 compared to Set1:
smaller cluster size information (r2x1, r2x2), post-shower energy
Figure 1: r2x1
before LDA cut
LDA cut for Set1
LDA cut for Set2
Figure 2: r2x2
before LDA cut
LDA cut for Set1
LDA cut for Set2
Figure 3: r3x3
before LDA cut
LDA cut for Set1
LDA cut for Set2
Figure 4: Residual distance
before LDA cut
LDA cut for Set1
LDA cut for Set2
Note: Only plos for LDA cut @70 efficiency for Set2 are shown
Figure : number of charge particles around photon
Figure 5: number of EEMC tower around photon
Figure 6: number of BEMC tower around photon
Figure 7: photon-jet pt balance
Figure 8: SMD energy in 5 centrapl strips
Figure 9: SMD energy in 25 central strips: u and v plane separately (plot for V plane)
Figure 13: tower energy in r=0.7 radius
Figure 14: 3x3 pre-shower1 energy
Figure 15: 3x3 pre-shower2 energy
Figure 16: 3x3 post-shower energy
Title:
"Photon-jet coincidence measurements
in polarized pp collisions at sqrt{s}=200GeV
with the STAR Endcap Calorimeter"
Title:
"Photon-jet coincidence measurements
in polarized pp collisions at sqrt{s}=200 GeV
with the STAR Endcap Calorimeter"
Data set and cuts:
Figure 1: Average ratio: pt_true / (pt_reco/1.3) vs. pt_reco (GeV/c)
Figure 2:
Average momentum difference: pt_true - (pt_reco/1.3) vs. pt_reco (GeV/c)
Figure 3: Average ratio: (pt_true -1.06) (pt_reco/1.3) vs. pt_reco (GeV/c)
Similar to Fig. 1, but with the true photon pt reduced by 1.06 GeV
Resulting true/reco pt ratio is flat in 4-6 GeV range.
Before further pursuing our efforts in tuning the tower energy response in the Monte-Carlo,
needs to address the observed energy loss difference in the fisrt layer of the BEMC/EEMC detector.
See Jason's blog post from 2009.07.16 for more details:
Comparison muon energy deposit in the 1st BEMC/EEMC layers
Monte-Carlo setup:
Some definitions:
Notations used in the plots:
Note: compare "Left" plots with Brians old results
Figure 1a: Et correction factor vs. pt thrown
Figure 1b: Et correction factor vs. eta thrown
Figure 1c: Et correction factor vs. phi thrown
Note: compare "Right" plots with Jason results with EEMC only geometry
Figure 2a: Sampling fraction vs. pt thrown
Figure 2b: Sampling fraction vs. energy thrown
Figure 2c: Sampling fraction vs. eta thrown
Figure 2d: Sampling fraction vs. phi thrown
Figure 3a: SMD energy vs. energy thrown
Figure 3b: SMD energy vs. eta thrown
Monte-Carlo setup is desribed here
Figure 1:Single photon shower shape before (red) and after (black) EEMC cAir bug fixed
pt=7-8GeV, eta=1.2-1.4 (left), eta=1.6-1.8 (right)
Figure 2: Single photon shower shape vs. data
Monte-Carlo: pt=7-10GeV, eta=1.6-1.8
data: no pre-shower1,2; pt_photon>7, pt_jet>5. no eta cuts.
(see Fig. 1 from here for other pre-shower conditions)
Monte-Carlo setup is desribed here
Color coding:
Single particle shower shape before (left) and after (right) EEMC cAir bug fixed
Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
Eta-meson shower shapes (blue) taken from Fig. 1 from here of this post
All shapes are normalized to 1 at peak (central strip).
Figure 1: Pre-shower bin 0: E_pre1=0; E_pre2=0
Figure 2: Pre-shower bin 1: E_pre1=0; E_pre2>0
Figure 3: Pre-shower bin 2: E_pre1>0; E_pre1<0.004
Figure 4: Pre-shower bin 3: E_pre1>0.004; E_pre1<0.01
Results only for corrected EEMC geometry
All shapes are divided by MC single-photon shower shape.
Figure 5a: Pre-shower bin 0: E_pre1=0; E_pre2=0
Figure 5b: Pre-shower bin 1: E_pre1=0; E_pre2>0
Figure 5c: Pre-shower bin 2: E_pre1>0; E_pre1<0.004
Figure 5d: Pre-shower bin 3: E_pre1>0.004; E_pre1<0.01
Figure 6: Single photon to eta-meson shape ratios only (with error bars):
Pre-shower bins 0 (upper-left),1 (upper-right),2 (lower-left), and 3 (lower-right)
Real data, and signal/background Monte-Carlo samples:
pp@200GeV collisions, STAR produnctionLong.
Trigger: eemc-http-mb-L2gamma [id:137641] (L ~ 3.164 pb^1)
Pythia prompt photon (signal) Monte-Carlo sample.
Filtered Prompt Photon p6410EemcGammaFilter.
Partonic pt range 2-25 GeV.
Pythia 2->2 hard QCD processes (background) Monte-Carlo sample.
Filtered QCD p6410EemcGammaFilter.
Partonic pt range 2-25 GeV.
No studies yet
No studies yet
No studies yet
No studies yet
No studies yet
No comparison yet
Monte-Carlo setup:
Color coding:
Pre-shower bins:
Ep1/Ep2 is the energy deposited in the 1st/2nd EEMC pre-shower layer.
For a single particle MC it is a sum over
all pre-shower tiles in the EEMC with energy of 3 sigma above pedestal.
For eta-meson from pp2006 data the sum is over 3x3 tower patch
Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
Eta-meson shower shapes (blue) taken from Fig. 1 from here of this post
All shapes are normalized to 1 at peak (central strip)
Figure 1: Shower shape sorted by pre-shower conditions.
cAir-Fixed EEMC geometry (NO slow simulator, WITH SVT)
Ratio plot
Figure 2: Shower shape sorted by pre-shower conditions.
cAir-Fixed EEMC geometry (NO slow simulator, NO SVT)
Ratio plot
Figure 3: Shower shape sorted by pre-shower conditions.
cAir-Fixed EEMC geometry (WITH slow simulator, WITH SVT)
Ratio plot
Figure 4: Shower shape sorted by pre-shower conditions.
Old cAir-bug EEMC geometry (NO slow simulator, WITH SVT)
Click here to see the plot
Starting with a fixed (50K events) for each type of particle.
Change in number of counts for a given pre-shower bin
with different detector configuration shows pre-shower migration
Figure 5: Pre-shower migration.
cAir-Fixed EEMC geometry (WITH SVT)
Figure 6: Pre-shower migration.
cAir-Fixed EEMC geometry (WITHOUT SVT)
Figure 7: Sampling fraction (0.05 E_reco / E_thrown).
cAir-Fixed EEMC geometry (WITHOUT Slow-simulator)
Figure 8: Sampling fraction (0.05 E_reco / E_thrown).
cAir-Fixed EEMC geometry (WITH Slow-simulator)
Slow simulator introduce some non-linearity in the sampling fraction
Figure 9: Sampling fraction (0.05 E_reco / E_thrown).
cAir-Fixed EEMC geometry (WITHOUT SVT, WITHOUT Slow-simulator)
Click here to see the plot
Figure 10: Sampling fraction (0.05 E_reco / E_thrown).
Old cAir-bug EEMC geometry (NO slow simulator, WITH SVT)
Click here to see the plot
Monte-Carlo setup:
Color coding:
Pre-shower bins:
Note: Ep1/Ep2 is the energy deposited in the 1st/2nd EEMC pre-shower layer.
For a single photon MC it is a sum over
all pre-shower tiles in the EEMC with energy of 3 sigma above pedestal.
For eta-meson/gamma-jet candidates from pp2006 data the sum is over 3x3 tower patch
Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
Eta-meson shower shapes (blue) taken from Fig. 1 from here of this post
All shapes are normalized to 1 at peak (central strip)
Figure 1: Shower shape sorted by pre-shower conditions.
Figure 2: Shower shape ratio. All shapes in Fig. 1 are divided by single photon shape
for "SVT+LOW_EM" configuration (black circles in Fig. 1)
Figure 3: Sampling fraction (0.05 * E_reco/ E_thrown)
Figure 4: Pre-shower1 energy (all tiles)
Figure 5: Pre-shower2 energy (all tiles)
Figure 6: Post-shower energy (all tiles)
Figure 7: Pre-shower bin photon migration
Figure 8a: Energy ratio in 2x1 to 3x3 cluster
For the first 4 pre-shower bins total yield in MC is normalized to that of the data
Blue circles indicate photon-jet candidates [pp2006] (points from this post)
Same data on a linear scale
Figure 8b: Energy ratio in 2x1 to 3x3 cluster: 7 < pt < 8 and 1.2 < eta < 1.4
Figure 8c: Energy ratio in 2x1 to 3x3 cluster: 7 < pt < 8 and 1.6 < eta < 1.8
Figure 9: Average energy ratio in 2x1 to 3x3 cluster vs. thrown energy
Figure 10: Average energy ratio in 2x1 to 3x3 cluster vs. thrown energy
LOW_EM option for the STAR geometry (Low cuts on Electro-Magnetic processes)
is equivalent to the following set of GEANT cuts:
All these values are for kinetic energy in GeV.
Cut meaning and GEANT default values:
Some details can be found at this link and in the GEANT manual
Monte-Carlo setup:
Color coding:
Figure 1: Sampling fraction vs. thrown energy (upper plot)
and vs. azimuthal angle (lower left) and rapidity (lower right)
Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
Eta-meson shower shapes (blue) taken from Fig. 1 from here of this post
All shapes are normalized to 1 at peak (central strip)
Pre-shower bins:
Ep1/Ep2 is the energy deposited in the 1st/2nd EEMC pre-shower layer.
For a single particle MC it is a sum over
all pre-shower tiles in the EEMC with energy of 3 sigma above pedestal.
For eta-meson from pp2006 data the sum is over 3x3 tower patch
The number of nested volumes (nv),
is the total number of parent volumes for the sensitive volume
(sensitive volume is indicated by the HITS in the tree structure below).
For the Jason and CVS files this nv number seems to be the same
(see block tree structures below).
Then why volume ids id in g2t tables has changed?
The answer I found (which seems trivial to me know)
is that in the original (CVS) file ECAL
block has been instantiated (positioned) twice.
The second appearance is the prototype (East) version of the Endcap
(Original ecalgeo.g from CVS)
if (emcg_OnOff==1 | emcg_OnOff==3) then Position ECAL in CAVE z=+center endif if (emcg_OnOff==2 | emcg_OnOff==3) then Position ECAL in CAVE z=-center ThetaZ=180 endif
In Jason version the second appearance has been removed
(what seems natural and it should not have any effect)
(ecalgeo.g Jason edits, g23):
IF (emcg_OnOff>0) THEN Create ECAL ..... IF (emcg_OnOff==2 ) THEN Prin1 ('East Endcap has been removed from the geometry' ) ENDIF EndIF! emcg_OnOff
Unfortunately, this affects the way GEANT counts nested volumes
(effectively the total number was reduced by 1, from 8 to 7)
and this is the reason why the volume numbering scheme
in g2t tables has been affected.
I propose to put back these East Endcap line back,
since in this case it will not require any additional
changes to the EEMC decoder and g2t tables.
blue - added volumes in Jason file
red - G10 volume removed in Jason file
HITS - sensitive volumes
---- Jason file ----
ECAL
EAGA
|EMSS
| -EFLP
| |ECVO
| | |EMOD
| | | |ESEC
| | | | |ERAD
| | | | | -ELED
| | | | |EMGT
| | | | | |EPER
| | | | | | |ETAR
| | | | | | | -EALP
| | | | | | | -ESCI -> HITS
| |ESHM
| | |ESPL
| | | |EXSG
| | | | -EXPS
| | | | -EHMS -> HITS
| | | | -EBLS
| | | | -EFLS
| | |ERSM
| -ESSP
| -ERCM
| -EPSB
|ECGH
| -ECHC
---- CVS file ----
ECAL
EAGA
|EMSS
| -EFLP
| |ECVO
| | |EMOD
| | | |ESEC
| | | | |ERAD
| | | | | -ELED
| | | | |EMGT
| | | | | |EPER
| | | | | | |ETAR
| | | | | | | -EALP
| | | | | | | -ESCI -> HITS
| |ESHM
| | |ESPL
| | | |EXSG
| | | | -EHMS -> HITS
| | | -EXGT
| | -ERSM
| -ESSP
| -ERCM
| -EPSB
|ECGH
| -ECHC
Create ECAL Block ECAL is one EMC EndCap wheel Create and Position EAGA AlphaZ=halfi EndBlock Block EAGA IS HALF OF WHEEL AIR VOLUME FORTHE ENDCAP MODULE Create AND Position EMSS konly='MANY' Create AND Position ECGH alphaz=90 kOnly='ONLY' EndBlock Block EMSS is the steel support of the endcap module Create AND Position EFLP z=zslice-center+zwidth/2 Create AND Position ECVO z=zslice-center+zwidth/2 Create AND Position ESHM z=zslice-center+zwidth/2 kOnly='MANY' Create AND Position ECVO z=zslice-center+zwidth/2 Create AND Position ESSP z=zslice-center+zwidth/2 Create ERCM Create EPSB EndBlock Block ECVO is one of endcap volume with megatiles and radiators Create AND Position EMOD alphaz=d3 ncopy=i_sector EndBlock Block ESHM is the shower maxsection Create and Position ESPL z=currentk Only='MANY' Create ERSM EndBlock Block ECGH is air gap between endcap half wheels Create ECHC EndBlock Block ECHC is steel endcap half cover EndBlock Block ESSP is stainless steelback plate EndBlock Block EPSB IS A PROJECTILE STAINLESS STEEL BAR EndBlock Block ERCM is stainless steel tie rod in calorimeter sections EndBlock Block ERSM is stainless steel tie rod in shower max EndBlock Block EMOD (fsect,lsect) IS ONE MODULEOF THE EM ENDCAP Create AND Position ESEC z=section-curr+secwid/2 EndBlock Block ESEC is a single em section Create AND Position ERAD z=length+(cell)/2+esec_deltaz Create AND Position EMGT z=length+(gap+cell)/2+esec_deltaz Create AND Position ERAD z=length+cell/2+esec_deltaz EndBlock Block EMGT is a 30 degree megatile Create AND Position EPER alphaz=myPhi EndBlock Block EPER is a 5 degree slice of a 30 degree megatile (subsector) Create and Position ETAR x=(rbot+rtop)/2ort=yzx EndBlock Block ETAR is a single calorimeter cell, containing scintillator, fiber router, etc... Create AND Position EALP y=(-megatile+emcs_alincell)/2 Create AND Position ESCI y=(-megatile+g10)/2+emcs_alincell _ EndBlock Block ESCI is the active scintillator (polystyrene) layer EndBlock Block ERAD is the lead radiator with stainless steel cladding Create AND Position ELED EndBlock Block ELED is a lead absorber plate EndBlock Block EFLP is the aluminum (aluminium) front plate of the endcap EndBlock Block EALP is the thin aluminium plate in calorimeter cell EndBlock Block ESPL is the logical volume containing an SMD plane Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' Create and Position EXSG alphaz=d3 ort=x-y-z ncopy=isec kOnly='MANY' Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' Create and Position EXSG alphaz=d3 ort=x-y-z ncopy=isec kOnly='MANY' Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' EndBlock Block EXSG Is another logical volume... this one acutally creates the planes Create and Position EXPS kONLY='MANY' Create and Position EHMS x=xc y=yc alphaz=-45 kOnly='ONLY' Create and Position EBLS x=xc y=yc z=(+esmd_apex/2+esmd_back_layer/2) alphaz=-45 kOnly='ONLY' Create and Position EHMS x=xc y=yc alphaz=-45 ort=x-y-z kOnly='ONLY' Create and Position EFLS x=xc y=yc z=(-esmd_apex/2-esmd_front_layer/2) alphaz=-45 ort=x-y-z kOnly='ONLY' EndBlock Block EHMS defines the triangular SMD strips Endblock! EHMS Block EFLS is the layer of material on the front of the SMD planes EndBlock! EFLS Block EBLS is the layer of material on the back of the SMD planes EndBlock! EFLS Block EXPS is the plastic spacer in the shower maximum section EndBlock
Create ECAL Block ECAL is one EMC EndCap wheel Create and Position EAGA AlphaZ=halfi EndBlock Block EAGA is half of wheel air volume forthe EndCap module Create and Position EMSS konly='MANY' Create and Position ECGH AlphaZ=90 konly='ONLY' EndBlock Block EMSS is steel support of the EndCap module Create and Position EFLP z=zslice-center+slcwid/2 Create and Position ECVO z=zslice-center+slcwid/2 Create and Position ESHM z=zslice-center+slcwid/2 Create and Position ECVO z=zslice-center+slcwid/2 Create and Position ESSP z=zslice-center+slcwid/2 Create ERCM Create EPSB EndBlock Block ECVO is one of EndCap Volume with megatiles and radiators Create and Position EMOD AlphaZ=d3 Ncopy=J_section EndBlock Block ESHM is the SHower Maxsection Create and Position ESPL z=current Create ERSM Endblock Block ECGH is air Gap between endcap Half wheels Create ECHC EndBlock Block ECHC is steel EndCap Half Cover EndBlock Block ESSP is Stainless Steelback Plate endblock Block EPSB is Projectile Stainless steel Bar endblock Block ERCM is stainless steel tie Rod in CaloriMeter sections endblock Block ERSM is stainless steel tie Rod in Shower Max endblock Block EMOD is one moduleof the EM EndCap Create and Position ESEC z=section-curr+secwid/2 endblock Block ESEC is a single EM section Create and Position ERAD z=len + (cell)/2 Create and Position EMGT z=len +(gap+cell)/2 Create and Position ERAD z=len + cell/2 Endblock Block EMGT is a megatile EM section Create and Position EPER AlphaZ=(emcs_Nslices/2-isec+0.5)*dphi Endblock Block EPER is a EM subsection period (super layer) Create and Position ETAR x=(RBot+RTop)/2ORT=YZX EndBlock Block ETAR is one CELL of scintillator, fiber and plastic Create and Position EALP y=(-mgt+emcs_AlinCell)/2 Create and Position ESCI y=(-mgt+G10)/2+emcs_AlinCell _ EndBlock Block ESCI is the active scintillator (polystyren) layer endblock Block ERAD is radiator Create and PositionELED endblock Block ELED is lead absorber Plate endblock Block EFLP is First Aluminium plate endblock Block EALP is ALuminiumPlate in calorimeter cell endblock Block ESPL is one of the Shower maxPLanes Create and position EXSG AlphaZ=d3Ncopy=isec Create and position EXSG AlphaZ=d3Ncopy=isec Create and position EXGT z=msecwd AlphaZ=d3 Create and position EXSG AlphaZ=d3 ORT=X-Y-Z Ncopy=isec Create and position EXGT z=-msecwd AlphaZ=d3 Create and position EXSG AlphaZ=d3Ncopy=isec Create and position EXGT z=msecwd AlphaZ=d3 Create and position EXSG AlphaZ=d3 ORT=X-Y-Z Ncopy=isec Create and position EXGT z=-msecwd AlphaZ=d3 Endblock Block EXSG is the Shower maxGap for scintillator strips Create EHMS endblock Block EHMS is sHower Max Strip Endblock Block EXGT is the G10 layer in the Shower Max EndBlock
Original (ecalgeo.g) file from CVS
****************************************************************************** Module ECALGEO is the EM EndCap Calorimeter GEOmetry Created 26 jan 1996 Author Rashid Mehdiyev * * Version 1.1, W.J. Llope * - changed sensitive medium names... * * Version 2.0, R.R. Mehdiyev 16.04.97 * - Support walls included * - intercell and intermodule gaps width updated * - G10 layers inserted * Version 2.1, R.R. Mehdiyev 23.04.97 * - Shower Max Detector geometry added * - Variable eta grid step size introduced * Version 2.2, R.R. Mehdiyev 03.12.97 * - Eta grid corrected * - Several changes in volume's dimensions * - Material changes in SMD * * Version 3.0, O. Rogachevsky 28.11.99 * - New proposal for calorimeter SN 0401 * * Version 4.1, O.Akio 3 Jan 01 * - Include forward pion detectors * Version 5.0, O. Rogachevsky 20.11.01 * - FPD is eliminated in this version * - More closed to proposal description * of calorimeter and SMD structure * ****************************************************************************** +CDE,AGECOM,GCONST,GCUNIT. * Content EAGA,EALP,ECAL,ECHC,ECVO,ECGH,EFLP,EHMS, ELED,EMGT,EMOD,EPER,EPSB,ERAD,ERCM,ERSM, ESHM,ESEC,ESCI,ESGH,ESPL,ESSP,EMSS, ETAR,EXGT,EXSG * Structure EMCG { Version, int Onoff, int fillMode} Structure EMCS { Type,ZOrig,ZEnd,EtaMin,EtaMax, PhiMin,PhiMax,Offset, Nsupsec,Nsector,Nsection,Nslices, Front,AlinCell,Frplast,Bkplast,PbPlate,LamPlate, BckPlate,Hub,Rmshift,SMShift,GapPlt,GapCel, GapSMD,SMDcentr,TieRod(2),Bckfrnt,GapHalf,Cover} * Structure EETR { Type,Etagr,Phigr,Neta,EtaBin(13)} * Structure ESEC { Isect, FPlmat, Cell, Scint, Nlayer } * Structure EMXG {Version,Sapex,Sbase,Rin,Rout,F4} * Structure EXSE {Jsect,Zshift,Sectype(6)} * Integer I_section,J_section,Ie,is,isec,i_str,Nstr,Type,ii,jj, cut,fsect,lsect,ihalf,filled * Real center,Plate,Cell,G10,diff,halfi, tan_low,tan_upp,Tanf,RBot,Rtop,Deta,etax,sq2,sq3, dup,dd,d2,d3,rshift,dphi,radiator,orgkeep,endkeep * Real maxcnt,msecwd,mxgten,curr,Secwid,Section, curcl,EtaTop,EtaBot,slcwid,zslice,Gap,mgt, xleft,xright,yleft,yright,current, rth,len,p,xc,yc,xx,yy,rbotrad, Rdel,dxy,ddn,ddup Integer N Parameter (N=12) * Tanf(etax) = tan(2*atan(exp(-etax))) * * ---------------------------------------------------------------------------- * * FillMode =1 only 2-5 sectors (in the first half) filled with scintillators * FillMode =2 all sectors filled (still only one half of one side) * FillMode =3 both halves (ie all 12 sectors are filled) Fill EMCG ! EM EndCAp Calorimeter basic data Version = 5.0 ! Geometry version OnOff = 3 ! Configurations 0-no, 1-west 2-east 3-both FillMode = 3 ! sectors fill mode Fill EMCS ! EM Endcap Calorimeter geometry Type = 1 ! =1 endcap, =2 fpd edcap prototype ZOrig = 268.763 ! calorimeter origin in z ZEnd = 310.007 ! Calorimeter end in z EtaMin = 1.086 ! upper feducial eta cut EtaMax = 2.0, ! lower feducial eta cut PhiMin = -90 ! Min phi PhiMax = 90 ! Max phi Offset = 0.0 ! offset in x Nsupsec = 6 ! Number of azimuthal supersectors Nsector = 30 ! Number of azimutal sectors (Phi granularity) Nslices = 5 ! number of phi slices in supersector Nsection = 4 ! Number of readout sections Front = 0.953 ! thickness of the front AL plates AlinCell = 0.02 ! Aluminim plate in cell Frplast = 0.015 ! Front plastic in megatile Bkplast = 0.155 ! Fiber routing guides and back plastic Pbplate = 0.457 ! Lead radiator thickness LamPlate = 0.05 ! Laminated SS plate thickness BckPlate = 3.175 ! Back SS plate thickness Hub = 3.81 ! thickness of EndCap hub Rmshift = 2.121 ! radial shift of module smshift = 0.12 ! radial shift of steel support walls GapPlt = 0.3/2 ! HALF of the inter-plate gap in phi GapCel = 0.03/2 ! HALF of the radial inter-cell gap GapSMD = 3.400 ! space for SMD detector SMDcentr = 279.542 ! SMD position TieRod = {160.,195} ! Radial position of tie rods Bckfrnt = 306.832 ! Backplate front Z GapHalf = 0.4 ! 1/2 Gap between halves of endcap wheel Cover = 0.075 ! Cover of wheel half * Rmshift = 2.121 ! radial shift of module * -------------------------------------------------------------------------- Fill EETR ! Eta and Phi grid values Type = 1 ! =1 endcap, =2 fpd EtaGr = 1.0536 ! eta_top/eta_bot tower granularity PhiGr = 0.0981747 ! Phi granularity (radians) NEta = 12 ! Eta granularity EtaBin = {2.0,1.9008,1.8065,1.7168,1.6317,1.5507,1.4738, 1.4007,1.3312,1.2651,1.2023,1.1427,1.086}! Eta rapidities *--------------------------------------------------------------------------- Fill ESEC ! First EM section ISect = 1 ! Section number Nlayer = 1 ! Number of Sci layers along z Cell = 1.505 ! Cell full width in z Scint = 0.5 ! Sci layer thickness * Fill ESEC ! First EM section ISect = 2 ! Section number Nlayer = 1 ! Number of Sci layers along z Cell = 1.505 ! Cell full width in z Scint = 0.5 ! Sci layer thickness * Fill ESEC ! Second EM section ISect = 3 ! Section number Nlayer = 4 ! Number of Sci layers along z Cell = 1.405 ! Cell full width in z Scint = 0.4 ! Sci layer thickness * Fill ESEC ! Third EM section ISect = 4 ! Section Nlayer = 18 ! Number of layers along z Cell = 1.405 ! Cell full width in z Scint = 0.4 ! Sci layer thickness * Fill ESEC ! 4th EM section ISect = 5 ! Section Nlayer = 1 ! Number of layers along z Cell = 1.505 ! Cell full width in z Scint = 0.5 ! Sci layer thickness *---------------------------------------------------------------------------- Fill EMXG ! EM Endcap SMD basic data Version = 1 ! Geometry version Sapex = 0.7 ! Scintillator strip apex Sbase = 1.0 ! Scintillator strip base Rin = 77.41 ! inner radius of SMD plane Rout = 213.922 ! outer radius of SMD plane F4 = .15 ! F4 thickness *---------------------------------------------------------------------------- Fill EXSE ! First SMD section JSect = 1 ! Section number Zshift = -1.215 ! Section width sectype = {4,1,0,2,1,0} ! 1-V,2-U,3-cutV,4-cutU * Fill EXSE ! Second SMD section JSect = 2 ! Section number Zshift = 0. ! Section width sectype = {0,2,1,0,2,3} ! 1-V,2-U,3-cutV,4-cutU * Fill EXSE ! Third SMD section JSect = 3 ! Section number Zshift = 1.215 ! Section width sectype = {1,0,2,1,0,2} ! 1-V,2-U,3-cutV,4-cutU *---------------------------------------------------------------------------- * Use EMCG * sq3 = sqrt(3.) sq2 = sqrt(2.) prin1 emcg_version ('ECALGEO version ', F4.2) * Endcap USE EMCS type=1 USE EETR type=1 orgkeep = emcs_ZOrig endkeep = emcs_ZEnd if(emcg_OnOff>0) then diff = 0.0 center = (emcs_ZOrig+emcs_ZEnd)/2 Tan_Upp = tanf(emcs_EtaMin) Tan_Low = tanf(emcs_EtaMax) rth = sqrt(1. + Tan_Low*Tan_Low) rshift = emcs_Hub * rth dup=emcs_Rmshift*Tan_Upp dd=emcs_Rmshift*rth d2=rshift + dd radiator = emcs_Pbplate + 2*emcs_LamPlate * d3=emcs_Rmshift-2*emcs_smshift dphi = (emcs_PhiMax-emcs_PhiMin)/emcs_Nsector Create ECAL if (emcg_OnOff==1 | emcg_OnOff==3) then Position ECAL in CAVE z=+center endif if (emcg_OnOff==2 | emcg_OnOff==3) then Position ECAL in CAVE z=-center ThetaZ=180 endif if(section > emcs_Zend) then prin0 section,emcs_Zend (' ECALGEO error: sum of sections exceeds maximum ',2F12.4) endif prin1 section (' EndCap calorimeter total depth ',F12.4) endif prin1 ('ECALGEO finished') * * ---------------------------------------------------------------------------- Block ECAL is one EMC EndCap wheel Material Air Medium standard Attribute ECAL seen=1 colo=7 ! lightblue shape CONE dz=(emcs_Zend-emcs_ZOrig)/2, Rmn1=orgkeep*Tan_Low-d2 Rmn2=endkeep*Tan_Low-d2, Rmx1=orgkeep*Tan_Upp+dup Rmx2=endkeep*Tan_Upp+dup do ihalf=1,2 filled=1 halfi = -105 + (ihalf-1)*180 if (ihalf=2 & emcg_FillMode<3) filled = 0 Create and Position EAGA AlphaZ=halfi enddo * EndBlock * ---------------------------------------------------------------------------- Block EAGA is half of wheel air volume for the EndCap module Attribute EAGA seen=1 colo=1 serial=filled ! black Material Air shape CONS dz=(emcs_Zend-emcs_ZOrig)/2, Rmn1=orgkeep*Tan_Low-d2 Rmn2=endkeep*Tan_Low-d2, Rmx1=orgkeep*Tan_Upp+dup Rmx2=endkeep*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax if (filled=1) then Create and Position EMSS konly='MANY' curr = orgkeep ; curcl = endkeep Create and position ECGH AlphaZ=90 konly='ONLY' endif EndBlock * ---------------------------------------------------------------------------- Block EMSS is steel support of the EndCap module Attribute EMSS seen=1 colo=1 ! black Material Iron shape CONS dz=(emcs_Zend-emcs_ZOrig)/2, Rmn1=orgkeep*Tan_Low-d2 Rmn2=endkeep*Tan_Low-d2, Rmx1=orgkeep*Tan_Upp+dup Rmx2=endkeep*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax zslice = emcs_ZOrig prin1 zslice (' Front Al plane starts at: ',F12.4) slcwid = emcs_Front Create and Position EFLP z=zslice-center+slcwid/2 zslice = zslice + slcwid prin1 zslice (' First calorimeter starts at: ',F12.4) fsect = 1; lsect = 3 slcwid = emcs_SMDcentr - emcs_GapSMD/2 - zslice * Create and Position ECVO z=zslice-center+slcwid/2 slcwid = emcs_GapSMD zslice = emcs_SMDcentr - emcs_GapSMD/2 prin1 section,zslice (' 1st calorimeter ends, SMD starts at: ',2F10.5) Create and Position ESHM z=zslice-center+slcwid/2 zslice = zslice + slcwid prin1 zslice (' SMD ends at: ',F10.5) * slcwid = 0 fsect = 4; lsect = 5 do I_section =fsect,lsect USE ESEC Isect=I_section Slcwid = slcwid + esec_cell*esec_Nlayer enddo slcwid = emcs_bckfrnt - zslice * Create and Position ECVO z = zslice-center+slcwid/2 zslice = emcs_bckfrnt prin1 section,zslice (' 2nd calorimeter ends, Back plate starts at: ',2F10.5) slcwid = emcs_BckPlate * Create and Position ESSP z=zslice-center+slcwid/2 zslice = zslice + slcwid prin1 zslice (' BackPlate ends at: ',F10.5) slcwid = emcs_Zend-emcs_ZOrig Create ERCM do i_str = 1,2 do is = 1,5 xx = emcs_phimin + is*30 yy = xx*degrad xc = cos(yy)*emcs_TieRod(i_str) yc = sin(yy)*emcs_TieRod(i_str) Position ERCM z=0 x=xc y=yc enddo enddo rth = orgkeep*Tan_Upp+dup + 2.5/2 xc = (endkeep - orgkeep)*Tan_Upp len = .5*(endkeep + orgkeep)*Tan_Upp + dup + 2.5/2 yc = emcs_Zend-emcs_ZOrig p = atan(xc/yc)/degrad Create EPSB do is = 1,6 xx = -75 + (is-1)*30 yy = xx*degrad xc = cos(yy)*len yc = sin(yy)*len Position EPSB x=xc y=yc AlphaZ=xx enddo EndBlock * ---------------------------------------------------------------------------- Block ECVO is one of EndCap Volume with megatiles and radiators Material Air Attribute ECVO seen=1 colo=3 ! green shape CONS dz=slcwid/2, Rmn1=zslice*Tan_Low-dd Rmn2=(zslice+slcwid)*Tan_Low-dd, Rmx1=zslice*Tan_Upp+dup Rmx2=(zslice+slcwid)*Tan_Upp+dup do J_section = 1,6 if (1 < J_section < 6 | emcg_FillMode > 1)then filled = 1 else filled = 0 endif d3 = 75 - (J_section-1)*30 Create and Position EMOD AlphaZ=d3 Ncopy=J_section enddo * EndBlock * ---------------------------------------------------------------------------- Block ESHM is the SHower Max section * Material Air Attribute ESHM seen=1 colo=4 ! blue Shape CONS dz=SlcWid/2, rmn1=zslice*Tan_Low-dd, rmn2=(zslice+slcwid)*Tan_Low-dd, rmx1=(zslice)*Tan_Upp+dup, rmx2=(zslice+slcwid)*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax USE EMXG Version=1 maxcnt = emcs_SMDcentr prin1 zslice,section,center (' Z start for SMD,section: ',3F12.4) * do J_section = 1,3 USE EXSE Jsect=J_section * current = exse_Zshift secwid = emxg_Sapex + 2.*emxg_F4 section = maxcnt + exse_zshift prin1 j_section,current,section,secwid (' layer, Z, width : ',i3,3F12.4) rbot=section*Tan_Low rtop=section*Tan_Upp prin1 j_section,rbot,rtop (' layer, rbot,rtop : ',i3,2F12.4) Create and position ESPL z=current * end do Create ERSM do i_str = 1,2 do is = 1,5 xx = emcs_phimin + (is)*30 yy = xx*degrad xc = cos(yy)*emcs_TieRod(i_str) yc = sin(yy)*emcs_TieRod(i_str) Position ERSM z=0 x=xc y=yc enddo enddo Endblock * ---------------------------------------------------------------------------- Block ECGH is air Gap between endcap Half wheels Material Air Medium standard Attribute ECGH seen=0 colo=7 ! lightblue shape TRD1 dz=(emcs_Zend-emcs_ZOrig)/2, dy =(emcs_gaphalf+emcs_cover)/2, dx1=orgkeep*Tan_Upp+dup, dx2=endkeep*Tan_Upp+dup rth = emcs_GapHalf + emcs_cover xx=curr*Tan_Low-d2 xleft = sqrt(xx*xx - rth*rth) yy=curr*Tan_Upp+dup xright = sqrt(yy*yy - rth*rth) secwid = yy - xx xx=curcl*Tan_Low-d2 yleft = sqrt(xx*xx - rth*rth) yy=curcl*Tan_Upp+dup yright = sqrt(yy*yy - rth*rth) slcwid = yy - xx xx=(xleft+xright)/2 yy=(yleft + yright)/2 xc = yy - xx len = (xx+yy)/2 yc = curcl - curr p = atan(xc/yc)/degrad rth = -(emcs_GapHalf + emcs_cover)/2 Create ECHC Position ECHC x=len y=rth Position ECHC x=-len y=rth AlphaZ=180 EndBlock * ---------------------------------------------------------------------------- Block ECHC is steel EndCap Half Cover Attribute ECHC seen=1 colo=1 ! black Material Iron shape TRAP dz=(curcl-curr)/2, thet=p, bl1=secwid/2, tl1=secwid/2, bl2=slcwid/2, tl2=slcwid/2, h1=emcs_cover/2 h2=emcs_cover/2, phi=0 alp1=0 alp2=0 EndBlock * ---------------------------------------------------------------------------- Block ESSP is Stainless Steel back Plate * Material Iron Attribute ESSP seen=1 colo=6 fill=1 shape CONS dz=emcs_BckPlate/2, Rmn1=zslice*Tan_Low-dd Rmn2=(zslice+slcwid)*Tan_Low-dd, Rmx1=zslice*Tan_Upp+dup Rmx2=(zslice+slcwid)*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax endblock * ---------------------------------------------------------------------------- Block EPSB is Projectile Stainless steel Bar * Material Iron Attribute EPSB seen=1 colo=6 fill=1 shape TRAP dz=(emcs_Zend-emcs_ZOrig)/2, thet=p, bl1=2.5/2, tl1=2.5/2, bl2=2.5/2, tl2=2.5/2, h1=2.0/2 h2=2.0/2, phi=0 alp1=0 alp2=0 endblock * ---------------------------------------------------------------------------- Block ERCM is stainless steel tie Rod in CaloriMeter sections * Material Iron Attribute ERSM seen=1 colo=6 fill=1 shape TUBE dz=slcwid/2, rmin=0, rmax=1.0425 ! nobody knows exactly endblock * ---------------------------------------------------------------------------- Block ERSM is stainless steel tie Rod in Shower Max * Material Iron Attribute ERSM seen=1 colo=6 fill=1 shape TUBE dz=slcwid/2, rmin=0, rmax=1.0425 endblock * ---------------------------------------------------------------------------- Block EMOD is one module of the EM EndCap Attribute EMOD seen=1 colo=3 serial=filled ! green Material Air Shape CONS dz=slcwid/2, phi1=emcs_PhiMin/emcs_Nsupsec, phi2=emcs_PhiMax/emcs_Nsupsec, Rmn1=zslice*Tan_Low-dd Rmn2=(zslice+slcwid)*Tan_Low-dd, Rmx1=zslice*Tan_Upp+dup Rmx2=(zslice+slcwid)*Tan_Upp+dup * * Running parameter 'section' contains the position of the current section * It should not be modified in daughters, use 'current' variable instead. * SecWid is used in all 'CONS' daughters to define dimensions. * * section = zslice curr = zslice + slcwid/2 Do I_section =fsect,lsect USE ESEC Isect=I_section * Secwid = esec_cell*esec_Nlayer if (I_section = 3 | I_section = 5) then ! no last radiator Secwid = Secwid - radiator else if (I_section = 4) then ! add one more radiator Secwid = Secwid - esec_cell + radiator endif Create and position ESEC z=section-curr+secwid/2 section = section + secwid * enddo endblock * ---------------------------------------------------------------------------- Block ESEC is a single EM section Attribute ESEC seen=1 colo=1 serial=filled Material Air Medium standard * Shape CONS dz=secwid/2, rmn1=(section-diff)*Tan_Low-dd, rmn2=(section+secwid-diff)*Tan_Low-dd, rmx1=(section-diff)*Tan_Upp+dup, rmx2=(section+secwid-diff)*Tan_Upp+dup * len = -secwid/2 current = section mgt = esec_scint + emcs_AlinCell _ + emcs_FrPlast + emcs_BkPlast gap = esec_cell - radiator - mgt prin2 I_section,section (' ESEC:I_section,section',i3,F12.4) Do is = 1,esec_Nlayer * define actual cell thickness: Cell = esec_cell plate = radiator * if (is=nint(esec_Nlayer) & (I_section = 3 | I_section = 5)) then Cell = mgt + gap Plate=0 else if (I_section = 4 & is = 1) then ! radiator only Cell = radiator endif * prin2 I_section,is,len,cell,current (' ESEC:I_section,is,len,cell,current ',2i3,3F12.4) if (I_section = 4 & is = 1) then ! radiator only cell = radiator + .14 Create and Position ERAD z=len + (cell)/2 len = len + cell current = current + cell else cell = mgt if(filled = 1) then Create and Position EMGT z=len +(gap+cell)/2 xx = current + (gap+cell)/2 prin2 I_section,is,xx (' MEGA I_section,is ',2i3,F10.4) endif len = len + cell + gap current = current + cell + gap if (Plate>0) then cell = radiator Create and Position ERAD z=len + cell/2 len = len + cell current = current + cell end if end if end do Endblock * ---------------------------------------------------------------------------- Block EMGT is a megatile EM section Attribute EMGT seen=1 colo=1 Material Air Medium standard * Shape CONS dz=mgt/2, rmn1=(current-diff)*Tan_Low-dd, rmn2=(current+mgt-diff)*Tan_Low-dd, rmx1=(current-diff)*Tan_Upp+dup, rmx2=(current+mgt-diff)*Tan_Upp+dup if (I_section=1 | I_section=2 | I_section=5) then Call GSTPAR (ag_imed,'CUTGAM',0.00001) Call GSTPAR (ag_imed,'CUTELE',0.00001) else Call GSTPAR (ag_imed,'CUTGAM',0.00008) Call GSTPAR (ag_imed,'CUTELE',0.001) Call GSTPAR (ag_imed,'BCUTE',0.0001) end if * Do isec=1,nint(emcs_Nslices) Create and Position EPER AlphaZ=(emcs_Nslices/2-isec+0.5)*dphi End Do Endblock *--------------------------------------------------------------------------- Block EPER is a EM subsection period (super layer) * Material POLYSTYREN Attribute EPER seen=1 colo=1 Shape CONS dz=mgt/2, phi1=emcs_PhiMin/emcs_Nsector, phi2=+emcs_PhiMax/emcs_Nsector, rmn1=(current-diff)*Tan_Low-dd, rmn2=(current+mgt-diff)*Tan_Low-dd, rmx1=(current-diff)*Tan_Upp+dup, rmx2=(current+mgt-diff)*Tan_Upp+dup * curcl = current+mgt/2 Do ie = 1,nint(eetr_NEta) EtaBot = eetr_EtaBin(ie) EtaTop = eetr_EtaBin(ie+1) RBot=(curcl-diff)*Tanf(EtaBot) * if(Plate > 0) then ! Ordinary Sci layer RTop=min((curcl-diff)*Tanf(EtaTop), _ ((current-diff)*Tan_Upp+dup)) else ! last Sci layer in section RTop=min((curcl-diff)*Tanf(EtaTop), _ ((current-diff)*Tan_Upp+dup)) endif check RBot<RTop * xx=tan(pi*emcs_PhiMax/180.0/emcs_Nsector) yy=cos(pi*emcs_PhiMax/180.0/emcs_Nsector) Create and Position ETAR x=(RBot+RTop)/2 ORT=YZX prin2 ie,EtaTop,EtaBot,rbot,rtop (' EPER : ie,EtaTop,EtaBot,rbot,rtop ',i3,4F12.4) enddo * EndBlock * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Block ETAR is one CELL of scintillator, fiber and plastic * Attribute ETAR seen=1 colo=4 ! blue * local z goes along the radius, y is the thickness Shape TRD1 dy=mgt/2 dz=(RTop-RBot)/2, dx1=RBot*xx-emcs_GapCel/yy, dx2=RTop*xx-emcs_GapCel/yy * Create and Position EALP y=(-mgt+emcs_AlinCell)/2 G10 = esec_scint Create and Position ESCI y=(-mgt+G10)/2+emcs_AlinCell _ +emcs_FrPlast EndBlock * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Block ESCI is the active scintillator (polystyren) layer * Material POLYSTYREN Material Cpolystyren Isvol=1 Attribute ESCI seen=1 colo=7 fill=0 ! lightblue * local z goes along the radius, y is the thickness Shape TRD1 dy=esec_scint/2, dz=(RTop-RBot)/2-emcs_GapCel Call GSTPAR (ag_imed,'CUTGAM',0.00008) Call GSTPAR (ag_imed,'CUTELE',0.001) Call GSTPAR (ag_imed,'BCUTE',0.0001) Call GSTPAR (ag_imed,'CUTNEU',0.001) Call GSTPAR (ag_imed,'CUTHAD',0.001) Call GSTPAR (ag_imed,'CUTMUO',0.001) * define Birks law parameters Call GSTPAR (ag_imed,'BIRK1',1.) Call GSTPAR (ag_imed,'BIRK2',0.013) Call GSTPAR (ag_imed,'BIRK3',9.6E-6) * HITS ESCI Birk:0:(0,10) * xx:16:H(-250,250) yy:16:(-250,250) zz:16:(-350,350), * px:16:(-100,100) py:16:(-100,100) pz:16:(-100,100), * Slen:16:(0,1.e4) Tof:16:(0,1.e-6) Step:16:(0,100), * none:16: endblock * ---------------------------------------------------------------------------- Block ERAD is radiator * Material Iron Attribute ERAD seen=1 colo=6 fill=1 ! violet Shape CONS dz=radiator/2, rmn1=(current)*Tan_Low-dd, rmn2=(current+cell)*Tan_Low-dd, rmx1=(current)*Tan_Upp+dup, rmx2=(current+radiator)*Tan_Upp+dup Create and Position ELED endblock * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Block ELED is lead absorber Plate * Material Lead Attribute ELED seen=1 colo=4 fill=1 Shape TUBS dz=emcs_Pbplate/2, rmin=(current)*Tan_Low, rmax=(current+emcs_Pbplate)*Tan_Upp, Call GSTPAR (ag_imed,'CUTGAM',0.00008) Call GSTPAR (ag_imed,'CUTELE',0.001) Call GSTPAR (ag_imed,'BCUTE',0.0001) Call GSTPAR (ag_imed,'CUTNEU',0.001) Call GSTPAR (ag_imed,'CUTHAD',0.001) Call GSTPAR (ag_imed,'CUTMUO',0.001) endblock * ---------------------------------------------------------------------------- Block EFLP is First Aluminium plate * Material Aluminium Attribute EFLP seen=1 colo=3 fill=1 ! green shape CONS dz=emcs_Front/2, Rmn1=68.813 Rmn2=68.813, Rmx1=(zslice-diff)*Tan_Upp+dup, Rmx2=(zslice + slcwid-diff)*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax endblock * ---------------------------------------------------------------------------- Block EALP is ALuminium Plate in calorimeter cell * Material Aluminium Material StrAluminium isvol=0 Attribute EALP seen=1 colo=1 Shape TRD1 dy=emcs_AlinCell/2 dz=(RTop-RBot)/2 Call GSTPAR (ag_imed,'CUTGAM',0.00001) Call GSTPAR (ag_imed,'CUTELE',0.00001) Call GSTPAR (ag_imed,'LOSS',1.) Call GSTPAR (ag_imed,'STRA',1.) endblock * ---------------------------------------------------------------------------- Block ESPL is one of the Shower max PLanes * Material Air Attribute ESPL seen=1 colo=3 ! blue Shape TUBS dz=SecWid/2, rmin=section*Tan_Low-1.526, rmax=(section-secwid/2)*Tan_Upp+dup, phi1=emcs_PhiMin phi2=emcs_PhiMax USE EMXG Version=1 msecwd = (emxg_Sapex+emxg_F4)/2 do isec=1,6 cut=1 d3 = 75 - (isec-1)*30 if (exse_sectype(isec) = 0 | (emcg_FillMode=1 & (isec=6 | isec=1))) then cut = 0 Create and position EXSG AlphaZ=d3 Ncopy=isec else if(exse_sectype(isec) = 1) then ! V Create and position EXSG AlphaZ=d3 Ncopy=isec Create and position EXGT z=msecwd AlphaZ=d3 else if(exse_sectype(isec) = 2) then ! U Create and position EXSG AlphaZ=d3 ORT=X-Y-Z Ncopy=isec Create and position EXGT z=-msecwd AlphaZ=d3 else if(exse_sectype(isec) = 3) then ! cut V cut=2 Create and position EXSG AlphaZ=d3 Ncopy=isec Create and position EXGT z=msecwd AlphaZ=d3 else if(exse_sectype(isec) = 4) then ! cut U cut=2 Create and position EXSG AlphaZ=d3 ORT=X-Y-Z Ncopy=isec Create and position EXGT z=-msecwd AlphaZ=d3 endif enddo Endblock * ---------------------------------------------------------------------------- Block EXSG is the Shower max Gap for scintillator strips * Attribute EXSG seen=1 colo=7 serial=cut ! black Material Air Shape TUBS dz=SecWid/2, rmin=section*Tan_Low-1.526, rmax=(section-secwid/2)*Tan_Upp+dup, phi1=emcs_PhiMin/emcs_Nsupsec, phi2=emcs_PhiMax/emcs_Nsupsec * Rbot = emxg_Rin Rtop = emxg_Rout if(cut > 0) then if(cut = 1) then Rdel = 3.938 Nstr = 288 else Rdel = -.475 Nstr = 285 endif rth = .53*rdel ! .53 --- tentatavily ddn = sq3*1.713 + Rdel ddup = .5*1.846 + 1.713 prin2 Rbot,Rtop,Nstr (' EXSG: Rbot,Rtop,Nstr',2F12.4,I5) mgt = emxg_Sbase + .01 do i_str = 1,nstr p = .5*(i_str-1)*mgt + 41.3655 * if (p <= (.5*rbot*sq3 + rth)) then dxy = 1.9375*sq2 xleft = .5*sq2*p*(sq3 + 1.) - dxy yleft = .5*sq2*p*(sq3 - 1.) - dxy yright = .5*sq2*(sqrt( rbot*rbot - p*p) - p) xright = sq2*p + yright else if ((.5*rbot*sq3 + rth) < p <= (.5*Rtop + 1.5)) then prin2 i_str,p (' EXSG: 2 - -i_str,p:',i3,F12.4) dxy = 1.9375*sq2 xleft = .5*sq2*p*(sq3 + 1.) - dxy yleft = .5*sq2*p*(sq3 - 1.) - dxy dxy = rdel*sq2/sq3 yright = .5*sq2*p*(1.- 1./sq3) xright = sq2*p - yright - dxy yright = -yright - dxy else if (p > (.5*rtop +1.5)) then prin2 i_str,p (' EXSG: 3 - - i_str,p:',i3,F12.4) yleft = (sqrt(rtop*rtop - p*p) - p)/sq2 xleft = sq2*p + yleft dxy = rdel*sq2/sq3 yright = .5*sq2*p*(1.- 1./sq3) xright = sq2*p - yright - dxy yright = -yright - dxy dxy = 0. if ((.5*sq3*160.- ddn) < p <= (.5*sq3*160.+ ddup) ) then prin2 i_str,p (' EXSG: 4 - - i_str,p:',i3,F12.4) xc = .5*(sq3*160.+1.846) yc = xc - .5*sq3*1.713 if (p > yc) then dxy = .5*sq2*(2/sq3*rdel + .5*sq3*1.846 +_ sqrt(1.713*1.713 - (p-xc)*(p-xc))) else dxy = sq2/sq3*(p - .5*sq3* 160. + ddn) endif else if ((.5*sq3*195.- ddn) < p <= (.5*sq3*195. + ddup) ) then prin2 i_str,p (' EXSG: 5 - - i_str,p:',i3,F12.4) xc = .5*(sq3*195.+1.846) yc = xc - .5*sq3*1.713 if (p > yc) then dxy = .5*sq2*(2/sq3*rdel + .5*sq3*1.846 +_ sqrt(1.713*1.713 - (p-xc)*(p-xc))) else dxy = sq2/sq3*(p - .5*sq3*195. + ddn) endif endif xright = xright + dxy yright = yright + dxy endif dxy = section*Tan_Upp - Rtop xc = .5*(xright+xleft) + dxy yc = .5*(yright+yleft) xx = .5*sq2*(xleft+yleft) yy = .5*sq2*(xright+yright) len = xx-yy prin2 i_str,p,yy,xx,len,xc,yc (' EXSG: i_str,x,y1,y2,len,xc,yc:',i3,6F12.4) * Create EHMS if (mod(i_str,2) != 0 ) then Position EHMS x=xc y=yc AlphaZ=-45 else Position EHMS x=xc y=yc AlphaZ=-45 ORT=X-Y-Z endif end do endif * dcut exsg z 0 0 10 0.1 0.1 * dcut exsg y 0 10 -50 0.7 0.7 endblock * ---------------------------------------------------------------------------- Block EHMS is sHower Max Strip * Material POLYSTYREN Material Cpolystyren Isvol=1 Attribute EHMS seen=1 colo=2 serial=cut ! red Shape TRD1 dx1=0 dx2=emxg_Sbase/2 dy=len/2 dz=emxg_Sapex/2 Call GSTPAR (ag_imed,'CUTGAM',0.00008) Call GSTPAR (ag_imed,'CUTELE',0.001) Call GSTPAR (ag_imed,'BCUTE',0.0001) * define Birks law parameters Call GSTPAR (ag_imed,'BIRK1',1.) Call GSTPAR (ag_imed,'BIRK2',0.0130) Call GSTPAR (ag_imed,'BIRK3',9.6E-6) * HITS EHMS Birk:0:(0,10) * xx:16:SH(-250,250) yy:16:(-250,250) zz:16:(-350,350), * px:16:(-100,100) py:16:(-100,100) pz:16:(-100,100), * Slen:16:(0,1.e4) Tof:16:(0,1.e-6) Step:16:(0,100), * none:16: Eloss:0:(0,10) * Endblock * ---------------------------------------------------------------------------- Block EXGT is the G10 layer in the Shower Max * * G10 is about 60% SiO2 and 40% epoxy Component Si A=28.08 Z=14 W=0.6*1*28./60. Component O A=16 Z=8 W=0.6*2*16./60. Component C A=12 Z=6 W=0.4*8*12./174. Component H A=1 Z=1 W=0.4*14*1./174. Component O A=16 Z=8 W=0.4*4*16./174. Mixture g10 Dens=1.7 Attribute EXGT seen=1 colo=7 Shape TUBS dz=emxg_F4/2, rmin=(section-diff)*Tan_Low-1.526, rmax=(section+msecwd-diff)*Tan_Upp, phi1=emcs_PhiMin/emcs_Nsupsec, phi2=emcs_PhiMax/emcs_Nsupsec Call GSTPAR (ag_imed,'CUTGAM',0.00001) Call GSTPAR (ag_imed,'CUTELE',0.00001) EndBlock * ---------------------------------------------------------------------------- * ECAL nice views: dcut ecvo x 1 10 -5 .5 .1 * draw emdi 105 0 160 2 13 .2 .1 * draw emdi 120 180 150 1 14 .12 .12 * --------------------------------------------------------------------------- end
ecalgeo.g geometry file (Jason Webb edits, g23)
c***************************************************************************** Module ECALGEO is the EM EndCap Calorimeter GEOmetry c-- Created 26 jan 1996 Author Rashid Mehdiyev c-- c Version 1.1, W.J. Llope c - changed sensitive medium names... c c Version 2.0, R.R. Mehdiyev 16.04.97 c - Support walls included c - intercell and intermodule gaps width updated c - G10 layers inserted c Version 2.1, R.R. Mehdiyev 23.04.97 c - Shower Max Detector geometry added c - Variable eta grid step size introduced c Version 2.2, R.R. Mehdiyev 03.12.97 c - Eta grid corrected c - Several changes in volumes dimensions c - Material changes in SMD c c Version 3.0, O. Rogachevsky 28.11.99 c - New proposal for calorimeter SN 0401 c c Version 4.1, O.Akio 3 Jan 01 c - Include forward pion detectors c c Version 5.0, O. Rogachevsky 20.11.01 c - FPD is eliminated in this version c - More closed to proposal description c of calorimeter and SMD structure c c***************************************************************************** +CDE,AGECOM,GCONST,GCUNIT. * Content EAGA,EALP,ECAL,ECHC,ECVO,ECGH,EFLP,EHMS, ELED,EMGT,EMOD,EPER,EPSB,ERAD,ERCM,ERSM, ESHM,ESEC,ESCI,ESGH,ESPL,ESSP,EMSS,ETAR, EXGT,EXSG,EXPS,EFLS,EBLS Structure EMCG { Version, int Onoff, int fillMode} Structure EMCS { Version,Type,zorg,zend,EtaMin,EtaMax, PhiMin,PhiMax,Offset, Nsupsec,Nsector,Nsection,Nslices, Front,AlinCell,Frplast,Bkplast,PbPlate,LamPlate, BckPlate,Hub,Rmshift,SMShift,GapPlt,GapCel, GapSMD,SMDcentr,TieRod(2),Bckfrnt,GapHalf,Cover, Rtie,slop} Structure EETR { Type,Etagr,Phigr,Neta,EtaBin(13)} Structure ESEC { Isect, FPlmat, Cell, Scint, Nlayer, deltaz, Jiggle(18) } Structure EMXG {Version,Sapex,Sbase,Rin,Rout,F4} Structure EXSE {Jsect,Zshift,Sectype(6)} Structure ESMD {Version, front_layer, back_layer, spacer_layer, base, apex } Integer I_section,J_section,Ie,is,isec,istrip,Nstr,Type,ii,jj, cut,fsect,lsect,ihalf,filled,i,j,k,i_sector Real center,Plate,Cell,G10,halfi, tan_low,tan_upp,Tanf,RBot,Rtop,Deta,etax,sq2,sq3, dup,dd,d2,d3,rshift,dphi,radiator Real maxcnt,msecwd,mxgten,curr,Secwid,Section, curcl,EtaTop,EtaBot,zwidth,zslice,Gap,megatile, xleft,xright,yleft,yright,current, rth,length,p,xc,yc,xx,yy,rdel,dxy,ddn,ddup Real myPhi Integer N Parameter (N=12) Tanf(etax) = tan(2*atan(exp(-etax))) c-------------------------------------------------------------------------------- c Data c c FillMode =1 only 2-5 sectors (in the first half) filled with scintillators c FillMode =2 all sectors filled (still only one half of one side) c FillMode =3 both halves (ie all 12 sectors are filled) c c OnOff =0 Do not build geometry c OnOff =1 Build West Endcap c OnOff =2 Build East Endcap (disabled) c OnOff =3 Build Both Endcaps (east disabled) c c Note: Fill EMCG ! EM EndCAp Calorimeter basic data Version = 5.0 ! Geometry version OnOff = 3 ! Configurations 0-no, 1-west 2-east 3-both FillMode = 3 ! sectors fill mode c-- Fill EMCS ! EM Endcap Calorimeter geometry Version = 1 ! Versioning Type = 1 ! =1 endcap, =2 fpd edcap prototype ZOrg = 268.763 ! calorimeter origin in z ZEnd = 310.007 ! Calorimeter end in z EtaMin = 1.086 ! upper feducial eta cut EtaMax = 2.0, ! lower feducial eta cut PhiMin = -90 ! Min phi PhiMax = 90 ! Max phi Offset = 0.0 ! offset in x Nsupsec = 6 ! Number of azimuthal supersectors Nsector = 30 ! Number of azimutal sectors (Phi granularity) Nslices = 5 ! number of phi slices in supersector Nsection = 4 ! Number of readout sections Front = 0.953 ! thickness of the front AL plates AlinCell = 0.02 ! Aluminim plate in cell Frplast = 0.015 ! Front plastic in megatile Bkplast = 0.155 ! Fiber routing guides and back plastic Pbplate = 0.457 ! Lead radiator thickness LamPlate = 0.05 ! Laminated SS plate thickness BckPlate = 3.175 ! Back SS plate thickness Hub = 3.81 ! thickness of EndCap hub Rmshift = 2.121 ! radial shift of module smshift = 0.12 ! radial shift of steel support walls GapPlt = 0.3/2 ! HALF of the inter-plate gap in phi GapCel = 0.03/2 ! HALF of the radial inter-cell gap GapSMD = 3.400 ! space for SMD detector << version 2 -- 3.600 >> SMDcentr = 279.542 ! SMD position TieRod = {160.,195} ! Radial position of tie rods Bckfrnt = 306.832 ! Backplate front Z GapHalf = 0.4 ! 1/2 Gap between halves of endcap wheel Cover = 0.075 ! Cover of wheel half Rtie = 1.0425 ! Radius of tie rod Slop = 0.1400 ! Added to cell containing radiator 6 (formerly hardcoded in geom) c-- Fill EMCS ! EM Endcap Calorimeter geometry Version = 2 ! Versioning Type = 1 ! =1 endcap, =2 fpd edcap prototype ZOrg = 268.763 ! calorimeter origin in z ZEnd = 310.007 ! Calorimeter end in z EtaMin = 1.086 ! upper feducial eta cut EtaMax = 2.0, ! lower feducial eta cut PhiMin = -90 ! Min phi PhiMax = 90 ! Max phi Offset = 0.0 ! offset in x Nsupsec = 6 ! Number of azimuthal supersectors Nsector = 30 ! Number of azimutal sectors (Phi granularity) Nslices = 5 ! number of phi slices in supersector Nsection = 4 ! Number of readout sections Front = 0.953 ! thickness of the front AL plates AlinCell = 0.02 ! Aluminim plate in cell Frplast = 0.015 ! Front plastic in megatile Bkplast = 0.155 ! Fiber routing guides and back plastic Pbplate = 0.457 ! Lead radiator thickness LamPlate = 0.05 ! Laminated SS plate thickness BckPlate = 3.175 ! Back SS plate thickness Hub = 3.81 ! thickness of EndCap hub Rmshift = 2.121 ! radial shift of module smshift = 0.12 ! radial shift of steel support walls GapPlt = 0.3/2 ! HALF of the inter-plate gap in phi GapCel = 0.03/2 ! HALF of the radial inter-cell gap GapSMD = 3.600 ! space for SMD detector (* from master_geom_bmp.xls *) SMDcentr = 279.542 ! SMD position TieRod = {160.,195} ! Radial position of tie rods Bckfrnt = 306.832 ! Backplate front Z GapHalf = 0.4 ! 1/2 Gap between halves of endcap wheel Cover = 0.075 ! Cover of wheel half Rtie = 0.75 ! Radius of tie rod Slop = 0.0000 ! Added to cell containing radiator 6 (formerly hardcoded in geom) c-- c--------------------------------------------------------------------------- c-- c-- Supporting documentation: c-- http://drupal.star.bnl.gov/STAR/system/files/SMD_module_stack.pdf c-- Fill ESMD ! shower maximum detector information Version = 1 ! versioning information front_layer = 0.161 ! thickness of front layer back_layer = 0.210 ! thickness of back layer base = 1.0 ! base of the SMD strip apex = 0.7 ! apex of the SMD strip spacer_layer = 1.2 ! spacer layer c-- Fill EETR ! Eta and Phi grid values Type = 1 ! =1 endcap, =2 fpd EtaGr = 1.0536 ! eta_top/eta_bot tower granularity PhiGr = 0.0981747 ! Phi granularity (radians) NEta = 12 ! Eta granularity EtaBin = {2.0,1.9008,1.8065,1.7168,1.6317,1.5507,1.4738, 1.4007,1.3312,1.2651,1.2023,1.1427,1.086}! Eta rapidities c-- c--------------------------------------------------------------------------- c-- Fill ESEC ! Preshower 1 / Radiator 1 ISect = 1 ! Section number Nlayer = 1 ! Number of Sci layers along z Cell = 1.505 ! Cell full width in z Scint = 0.475 ! Sci layer thickness (4.75mm Bicron) deltaz = -0.014 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- c-- Note: Jiggle allows one to shift each megatile by Jiggle(i) degrees, where c-- i indicates the layer within the section of the calorimeter. This feature c-- has only been crudely tested... i.e. it compiles and creates a reasonable c-- set of pictures, but I have not verified that every scintillator shows up... c-- There could be volume conflicts and this would need to be checked. --JW c-- Fill ESEC ! Preshower 2 / Radiator 2 ISect = 2 ! Section number Nlayer = 1 ! Number of Sci layers along z Cell = 1.505 ! Cell full width in z Scint = 0.475 ! Sci layer thickness (4.75mm Bicron) deltaz = -0.0182 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- Fill ESEC ! Megatiles 3-6 / Radiators 3-5 ISect = 3 ! Section number Nlayer = 4 ! Number of Sci layers along z Cell = 1.405 ! Cell full width in z Scint = 0.4 ! Sci layer thickness deltaz = -0.0145 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- Fill ESEC ! Megatiles 7-23 / Radiators 6-23 ISect = 4 ! Section Nlayer = 18 ! Number of layers along z Cell = 1.405 ! Cell full width in z Scint = 0.4 ! Sci layer thickness deltaz = +0.0336 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- Fill ESEC ! Postshower ISect = 5 ! Section Nlayer = 1 ! Number of layers along z Cell = 1.505 ! Cell full width in z Scint = 0.5 ! Sci layer thickness (5.0mm Kurarary) deltaz = +0.036 ! Amount to shift section in z to align with as-built numbers Jiggle = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} ! Degrees to shift EPER in each layer c-- c---------------------------------------------------------------------------- c-- Fill EMXG ! EM Endcap SMD basic data Version = 1 ! Geometry version Sapex = 0.7 ! Scintillator strip apex Sbase = 1.0 ! Scintillator strip base Rin = 77.41 ! inner radius of SMD plane Rout = 213.922 ! outer radius of SMD plane F4 = .15 ! F4 thickness c-- c---------------------------------------------------------------------------- c-- Fill EXSE ! First SMD section JSect = 1 ! Section number Zshift = -1.215 ! Section width sectype = {4,1,0,2,1,0} ! 1-V,2-U,3-cutV,4-cutU c-- Fill EXSE ! Second SMD section JSect = 2 ! Section number Zshift = 0. ! Section width sectype = {0,2,1,0,2,3} ! 1-V,2-U,3-cutV,4-cutU c-- Fill EXSE ! Third SMD section JSect = 3 ! Section number Zshift = 1.215 ! Section width sectype = {1,0,2,1,0,2} ! 1-V,2-U,3-cutV,4-cutU c-- c---------------------------------------------------------------------------- c-- Materials c-- c-- PVC used in the SMD spacer layers c-- Component H A=1 Z=1 W=3.0*1.0/62.453 Component C A=12 Z=6 W=2.0*12.0/62.453 Component Cl A=35.453 Z=17 W=1.0*35.453/62.453 Mixture PVC_Spacer Dens=1.390*(1.20/1.00) c-- c-- Lead alloy used in the radiators c-- Component Sn A=118.710 Z=50 W=0.014 Component Ca A=40.0780 Z=20 W=0.00075 Component Al A=26.9815 Z=13 W=0.0003 Component Pb A=207.190 Z=82 W=0.98495 Mixture PbAlloy DENS=11.35 c-- c-- Stainless Steel used in various places c-- Component Cr A=51.9960 Z=24 W=0.19 Component Ni A=58.6934 Z=28 W=0.09 Component Fe A=55.8450 Z=26 W=0.72 Mixture Steel DENS=8.03 c-- c-- Aluminized mylar. According to information which I dug up on a google c-- search, this is typically mylar coated with a thin (1000 angstrom) layer c-- of aluminium on each side. c-- c-- http://www.eljentechnology.com/datasheets/EJ590-B10HH%20data%20sheet.pdf c-- Component Mylar A=12.875 Z=6.4580 w=0.999 Component Al A=26.980 Z=13.000 w=0.001 Mixture AlMylar dens=1.390 c-- c-- G10 Epoxy used in various places c-- Component Si A=28.08 Z=14 W=0.6*1*28./60. Component O A=16 Z=8 W=0.6*2*16./60. Component C A=12 Z=6 W=0.4*8*12./174. Component H A=1 Z=1 W=0.4*14*1./174. Component O A=16 Z=8 W=0.4*4*16./174. Mixture G10 Dens=1.7 c-- c-- Fibreglass cloth used in SMD stackup. I googled this one too... a self- c-- described expert quotes typical densities and percent by volume c-- http://en.allexperts.com/q/Composite-Materials-2430/fiberglass-1.htm c-- c-- glass fiber: 2.6 g/cm3 (17.6%) resin: 1.3 g/cm3 (82.4%) c-- c-- Fiberglass density = 1.529 g/cm3 c-- c-- I will assume that G10 epoxy is close enough to the typical resins c-- used, at least in terms of chemical composition. Then c-- Component G10 A=18.017 Z=9.013 W=1.3*0.824/(1.3*0.824+2.6*0.176) Component Si A=28.08 Z=14 W=2.6*0.176/(1.3*0.824+2.6*0.176)*28.08/60.08 Component O A=16 Z=8 W=2.6*0.176/(1.3*0.824+2.6*0.176)*32.00/60.08 Mixture Fiberglass dens=1.53 c-- c-- c---------------------------------------------------------------------------- c-- Select versions of various geometry data c-- Use EMCG Use EMCS Version=2 Use EETR c-- c---------------------------------------------------------------------------- c-- Calculate frequently used quantities c-- sq3 = sqrt(3.) ! 1/tan(30deg) = sq3 sq2 = sqrt(2.) c-- c-- center = (emcs_zorg+emcs_zend)/2 ! center of the calorimeter tan_upp = tanf(emcs_etamin) ! think this is angle pointing to top of calo tan_low = tanf(emcs_etamax) ! think this is angle pointing to bot of calo rth = sqrt(1. + tan_low*tan_low) ! ?? rshift = emcs_hub * rth ! ?? dup = emcs_rmshift*tan_upp ! dd = emcs_rmshift*rth ! d2 = rshift + dd ! radiator = emcs_pbplate + 2*emcs_lamplate ! thickness of radiator assembly dphi = (emcs_phimax-emcs_phimin)/emcs_nsector ! single endcap sector c-- c---------------------------------------------------------------------------- c---------------------------------------------------------------------------- c-- BEGIN Prin1 emcg_version ('ecalgeo version: ',F4.2) c-- IF (emcg_OnOff>0) THEN c-- c-- Build the EEMC geometry for one half wheel c-- Create ECAL c-- c-- Position the two halves. Bottom half installed in 2003, top c-- half in 2004... so we allow logic to allow for the time c-- evolution of the calorimeter c-- c-- c-- West Endcap c-- IF (emcg_OnOff==1 | emcg_OnOff==3) THEN Position ECAL in CAVE z=+center ENDIF IF (section > emcs_zend) THEN Prin1 section, emcs_zend (' ECALGEO error: sum of sections exceeds maximum ',2F12.4) ENDIF IF (emcg_OnOff==2 ) THEN Prin1 ('East Endcap has been removed from the geometry' ) ENDIF c-- EndIF! emcg_OnOff c-- Prin1 ('ECALGEO finished') c-- c-- END c---------------------------------------------------------------------------- c----------------------------------------------------------------- Block ECAL -- c-- Block ECAL is one EMC EndCap wheel c-- c-- The EEMC is built from two 180 degree half-wheels tilted at an angle c-- with respect to zero in the STAR reference frame. This block is serves c-- as a logical volume which creates the two half wheels. c-- c-- Creates: c-- + EAGA c-- Material Air Attribute ECAL seen=0 colo=7 ! lightblue c-- Shape CONE dz=(emcs_zend-emcs_zorg)/2, rmn1=emcs_zorg*tan_low-d2, rmn2=emcs_zend*tan_low-d2, rmx1=emcs_zorg*tan_upp+dup, rmx2=emcs_zend*tan_upp+dup c-- c-- DO ihalf=1,2 c-- filled = 1 halfi = -105 + (ihalf-1)*180 if (ihalf=2 & emcg_FillMode<3) filled = 0 c-- Create and Position EAGA AlphaZ=halfi c-- ENDDO c-- EndBlock c----------------------------------------------------------------- Block EAGA -- c-- Block EAGA IS HALF OF WHEEL AIR VOLUME FOR THE ENDCAP MODULE c-- c-- The eemc is divided into two halves. one half installed for 2003 run, c-- second half added for 2004 and beyond. the eaga block represents one c-- of these half-wheels. it is an air volume which will be filled in c-- with additional detector components. c-- c-- Creates: c-- + EMSS -- steel support block c-- + ECGH -- air gap between the two halves c-- C-- Material AIR Attribute EAGA seen=0 colo=1 serial=FILLED ! BLACK C-- Shape CONS dz=(emcs_zend-emcs_zorg)/2, rmn1=emcs_zorg*tan_low-d2 rmn2=emcs_zend*tan_low-d2, rmx1=emcs_zorg*tan_upp+dup rmx2=emcs_zend*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- c-- IF ( FILLED .EQ. 1 ) THEN c-- Create AND Position EMSS konly='MANY' c-- curr = emcs_zorg curcl = emcs_zend c-- Create AND Position ECGH alphaz=90 kOnly='ONLY' c-- ENDIF c-- EndBlock c----------------------------------------------------------------- Block EMSS -- c-- Block EMSS is the steel support of the endcap module c-- c-- Creates: c-- + EFLP -- ALUMINIUM FRONT PLATE c-- + ECVO -- VOLUMES TO CONTAIN RADIATORS AND MEGATILES c-- + ESHM -- SHOWER MAX DETECTOR VOLUME c-- + ESSP -- STAINLESS STEEL BACKPLATE c-- + ERCM -- STAINLESS STEEL TIE-RODS PENETRATING ECVO c-- c-- Material Steel c-- Attribute EMSS seen=1 colo=1 ! BLACK Shape CONS dz=(emcs_zend-emcs_zorg)/2, rmn1=emcs_zorg*tan_low-d2 rmn2=emcs_zend*tan_low-d2, rmx1=emcs_zorg*tan_upp+dup rmx2=emcs_zend*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- c-- Aluminium front plate C-- zslice = emcs_zorg zwidth = emcs_front c-- Prin1 zslice+zwidth/2 (' Front Al plate centered at: ', F12.4 ) c-- Create AND Position EFLP z=zslice-center+zwidth/2 zslice = zslice + zwidth C-- Prin1 zslice (' FIRST CALORIMETER STARTS AT: ',F12.4) c-- c-- Preshower 1, preshower 2, and calorimeter tiles up to c-- megatile number six. c-- fsect = 1 ! first section lsect = 3 ! last section c-- zwidth = emcs_smdcentr - emcs_gapsmd/2 - zslice ! width of current slice c-- Prin1 zslice+zwidth/2 ('Sections 1-3 positioned at: ', F12.4 ) c-- Create AND Position ECVO z=zslice-center+zwidth/2 c-- zwidth = emcs_gapsmd zslice = emcs_smdcentr - emcs_gapsmd/2 c-- Prin1 section, zslice (' 1st calorimeter ends, smd starts at: ',2f10.5) Prin1 zwidth (' smd width = ',f10.5 ) c-- Prin1 zslice+zwidth/2 ('SMD section centered at: ', F12.4 ) c-- Do not kill neighbors Create AND Position ESHM z=zslice-center+zwidth/2 kOnly='MANY' zslice = zslice + zwidth c-- Prin1 zslice (' SMD ends at: ',f10.5) c-- c-- fsect = 4 ! first section lsect = 5 ! last section c-- c-- Calculate the width of the last two calorimeter sections c-- zwidth = 0 DO i_section = fsect,lsect c-- USE ESEC isect=i_section zwidth = zwidth + esec_cell*esec_nlayer c-- ENDDO c-- c-- ============================================================= c-- c-- Total width will be between the back plate and the current c-- position... this effectively turns the geometry into an c-- accordian... whatever was defined earlier will compress c-- / expand this section. so correcting the smd gap will c-- result in some small, sub-mm shifts of radiators and c-- megatiles... one would like to actually place these c-- into their absolute positions. c-- c-- ============================================================== c-- zwidth = emcs_bckfrnt - zslice c-- Prin1 zslice+zwidth/2 ('Sections 4-5 positioned at: ', F12.4 ) c-- Create AND Position ECVO z=zslice-center+zwidth/2 c-- zslice = emcs_bckfrnt c-- Prin1 section,zslice (' 2nd calorimeter ends, back plate starts at: ',2f10.5) c-- zwidth = emcs_bckplate c-- Create AND Position ESSP z=zslice-center+zwidth/2 c-- zslice = zslice + zwidth c-- Prin1 zslice ('EEMC Al backplate ends at: ',F12.4 ) c-- c-- Done with the calorimeter stackup. now go back and cut through the c-- calorimeter stack with the tie rods c-- c-- slice width will be full calorimeter depth zwidth = emcs_zend-emcs_zorg c-- Create ERCM c-- DO i = 1,2 ! two tie rods along DO j = 1,5 ! each gap between sectors (5 gaps) xx = emcs_phimin + j*30 yy = xx*degrad xc = cos(yy)*emcs_tierod(i) yc = sin(yy)*emcs_tierod(i) Position ERCM z=0 x=xc y=yc ENDDO ENDDO c-- c-- Now add in projective steel bars which form part of the support c-- structure of the eemc c-- rth = emcs_zorg*tan_upp+dup + 2.5/2 xc = (emcs_zend - emcs_zorg)*tan_upp length = .5*(emcs_zend + emcs_zorg)*tan_upp + dup + 2.5/2 yc = emcs_zend-emcs_zorg p = atan(xc/yc)/degrad c-- Create EPSB DO i = 1,6 c-- xx = -75 + (i-1)*30 yy = xx*degrad xc = cos(yy)*length yc = sin(yy)*length c-- Position EPSB X=XC Y=YC ALPHAZ=XX c-- ENDDO c-- EndBlock c----------------------------------------------------------------- Block ECVO -- c-- Block ECVO is one of endcap volume with megatiles and radiators c-- c-- CreateS: c-- + EMOD -- Responsible for creating esec which, in a glorious example c-- of spaghetti code, turns around and creates esec, which is c-- responsible for creating the radiators before and after the c-- smd layers. C-- Material AIR Attribute ECVO seen=1 colo=3 ! GREEN Shape CONS dz=zwidth/2, rmn1=zslice*tan_low-dd, rmn2=(zslice+zwidth)*tan_low-dd, rmx1=zslice*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup c-- c-- Loop over the SIX SECTORS in the current half-wheel. determine c-- whether the sector is filled or not, and create the "module". c-- By "module", we really mean endcap sector. (Lots of code in the c-- EEMC borrows from the barrel, and so barrel modlues get mapped c-- to EEMC sectors). c-- DO i_sector = 1,6 c-- IF (1 < I_SECTOR < 6 | EMCG_FILLMODE > 1) THEN filled = 1 ELSE filled = 0 ENDIF c-- d3 = 75 - (i_sector-1)*30 Create AND Position EMOD alphaz=d3 ncopy=i_sector c-- ENDDO c-- EndBlock c----------------------------------------------------------------- Block ESHM -- c-- Block ESHM is the shower max section c-- c-- CreateS: c-- + ESPL -- SHOWER MAXIMUM DETECTOR PLANES c-- + ERSM -- TIE RODS W/IN THE SHOWER MAXIMUM DETECTOR c-- Material AIR Attribute ESHM seen=1 colo=4 ! BLUE c-- Shape CONS dz=zwidth/2, rmn1=(zslice*tan_low)-dd, rmn2=(zslice+zwidth)*tan_low-dd, rmx1=(zslice)*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- USE EMXG c-- maxcnt = emcs_smdcentr Prin1 zslice, section, center (' === z start for smd,section: ',3f12.4) c-- c-- Loop over the three possible locations for the smd planes and c-- create them. note that code w/in espl will decide which of c-- 5 types of smd planes are created... u, v, cutu,cutv or spacer. c-- DO j_section = 1,3 c-- USE EXSE jsect=j_section c-- current = exse_zshift secwid = emxg_sapex + 2.*emxg_f4 section = maxcnt + exse_zshift c-- Prin1 j_section,current,section,secwid (' layer, z, width : ',i3,3f12.4) c-- rbot=section*tan_low rtop=section*tan_upp c-- Prin1 j_section,rbot,rtop (' layer, rbot,rtop : ',i3,2f12.4) c-- Prin1 j_section, center+current (' smd layer=',I1,' z=',F12.4 ) c-- Do not kill neighbors Create and Position ESPL z=current kOnly='MANY' c-- ENDDO c-- c-- Add in the tie rods which penetrate the SMD layers c-- Create ERSM c-- DO i = 1,2 DO j = 1,5 xx = emcs_phimin + j*30 yy = xx*degrad xc = cos(yy)*emcs_tierod(i) yc = sin(yy)*emcs_tierod(i) Position ERSM Z=0 X=XC Y=YC END DO END DO C-- EndBlock c----------------------------------------------------------------- Block ECGH -- c-- Block ECGH is air gap between endcap half wheels c-- c-- Creates: c-- + ECHC -- THE STAINLESS STEEL COVER FOR 1/2 OF THE EEMC. c-- Material AIR Medium standard Attribute ECGH seen=0 colo=7 ! LIGHTBLUE Shape TRD1 dz=(emcs_zend-emcs_zorg)/2, dy =(emcs_gaphalf+emcs_cover)/2, dx1=emcs_zorg*tan_upp+dup, dx2=emcs_zend*tan_upp+dup c-- c-- rth = emcs_gaphalf + emcs_cover xx=curr*tan_low-d2 xleft = sqrt(xx*xx - rth*rth) yy=curr*tan_upp+dup xright = sqrt(yy*yy - rth*rth) secwid = yy - xx xx=curcl*tan_low-d2 yleft = sqrt(xx*xx - rth*rth) yy=curcl*tan_upp+dup yright = sqrt(yy*yy - rth*rth) zwidth = yy - xx xx=(xleft+xright)/2 yy=(yleft + yright)/2 xc = yy - xx length = (xx+yy)/2 yc = curcl - curr p = atan(xc/yc)/degrad rth = -(emcs_gaphalf + emcs_cover)/2 c-- Create ECHC c-- Position ECHC X=+LENGTH Y=RTH Position ECHC X=-LENGTH Y=RTH ALPHAZ=180 c-- EndBlock c----------------------------------------------------------------- Block ECHC -- c-- Block ECHC is steel endcap half cover c-- Material steel Attribute ECHC seen=1 colo=1 ! BLACK c-- Shape TRAP dz=(curcl-curr)/2, thet=p, bl1=secwid/2, tl1=secwid/2, bl2=zwidth/2, tl2=zwidth/2, h1=emcs_cover/2, h2=emcs_cover/2, phi=0, alp1=0, alp2=0 c-- EndBlock c----------------------------------------------------------------- Block ESSP -- c-- Block ESSP is stainless steel back plate c-- Material steel Attribute ESSP seen=1 colo=6 fill=1 Shape CONS dz=emcs_bckplate/2, rmn1=zslice*tan_low-dd, rmn2=(zslice+zwidth)*tan_low-dd, rmx1=zslice*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup, phi1=emcs_phimin, phi2=emcs_phimax c-- EndBlock c----------------------------------------------------------------- Block EPSB -- c-- Block EPSB IS A PROJECTILE STAINLESS STEEL BAR C-- Material Steel Attribute EPSB seen=1 colo=6 FILL=1 Shape TRAP dz=(emcs_zend-emcs_zorg)/2, thet=p, bl1=2.5/2, tl1=2.5/2, bl2=2.5/2, tl2=2.5/2, h1=2.0/2, h2=2.0/2, phi=0, alp1=0, alp2=0 c-- c-- EndBlock c----------------------------------------------------------------- Block ERCM -- c-- Block ERCM is stainless steel tie rod in calorimeter sections c-- Material Steel Attribute ERSM seen=1 colo=6 FILL=1 c-- Shape TUBE dz=zwidth/2, rmin=0, rmax=emcs_rtie c-- c-- Looks like the tie rods are meant to engage the 1.525 cm diameter holes c-- piercing the ears of the smd spacer... 1.5 cm may be a better approximation c-- here. c-- c-- http://drupal.star.bnl.gov/star/system/files/smd_spacer_drawings.pdf c-- EndBlock c----------------------------------------------------------------- Block ERSM -- c-- Block ERSM is stainless steel tie rod in shower max c-- Material Steel Attribute ERSM seen=1 colo=6 FILL=1 c-- Shape TUBE dz=zwidth/2, rmin=0, rmax=emcs_rtie c-- c-- see comments above c-- EndBlock c----------------------------------------------------------------- Block EMOD -- c-- Block EMOD (fsect,lsect) IS ONE MODULE OF THE EM ENDCAP c-- c-- Arguements: (do be defined prior to the creation of this block) c-- c-- fsect -- first section to create c-- lsect -- last section to create c-- Attribute EMOD seen=1 colo=3 serial=FILLED ! GREEN Material Air Shape CONS dz=zwidth/2, phi1=emcs_phimin/emcs_nsupsec, phi2=emcs_phimax/emcs_nsupsec, rmn1=zslice*tan_low-dd, rmn2=(zslice+zwidth)*tan_low-dd, rmx1=zslice*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup c-- c-- Running parameter 'section' contains the position of the current section c-- it should not be modified in daughters, use 'current' variable instead. c-- secwid is used in all 'cons' daughters to define dimensions. c-- section = zslice curr = zslice + zwidth/2 c-- c-- DO i_section = fsect, lsect USE ESEC isect=i_section c-- secwid = esec_cell*esec_nlayer c-- c-- Section 3 precedes the smd. section 5 is the post shower. in c-- both cases these sections end with a scintillator layer and no c-- radiator. c-- IF (I_SECTION = 3 | I_SECTION = 5) THEN secwid = secwid - radiator ELSE IF (I_SECTION = 4) THEN ! add one more radiator secwid = secwid - esec_cell + radiator ENDIF c-- Prin1 i_section, section-curr+secwid/2 ('+ ECVO isection=',I1,' zcenter=', F12.4) c-- Create AND Position ESEC z=section-curr+secwid/2 c-- section = section + secwid c-- ENDDO! Loop over sections c-- EndBlock c----------------------------------------------------------------- Block ESEC -- c-- Block ESEC is a single em section Material AIR Medium standard Attribute ESEC seen=1 colo=1 serial=filled lsty=2 c-- Shape CONS dz=secwid/2, rmn1=(section)*tan_low-dd, rmn2=(section+secwid)*tan_low-dd, rmx1=(section)*tan_upp+dup, rmx2=(section+secwid)*tan_upp+dup c-- length = -secwid/2 current = section c-- megatile = esec_scint+emcs_alincell+emcs_frplast+emcs_bkplast c-- gap = esec_cell - radiator - megatile Prin2 i_section,section (' ESEC:i_section,section',i3,f12.4) c-- c-- Loop over all layers in this section c-- DO is = 1,esec_nlayer c-- c-- Define actual cell thickness: cell = esec_cell plate = radiator c-- IF (is=nint(esec_nlayer) & (i_section = 3 | i_section = 5)) THEN c-- cell = megatile + gap plate=0 c-- ELSE IF (i_section = 4 & is = 1) THEN ! RADIATOR ONLY c-- cell = radiator c-- ENDIF c-- Prin2 i_section,is,length,cell,current (' esec:i_section,is,length,cell,current ',2i3,3f12.4) C-- C-- This handles the special case in the section after the smd. c-- this section begins with a lead radiator. the previous section c-- ended with a plastic scintillator c-- IF (i_section = 4 & is = 1) THEN ! radiator only c-- c$$$ cell = radiator + .14 cell = radiator + emcs_slop ! ^^^^ probably the fiber router layer... but is this needed here? c-- Prin1 is, current + cell/2+esec_deltaz ( ' + ESEC radiator ilayer=',I2,' z=',F12.4 ) Create AND Position ERAD z=length+(cell)/2+esec_deltaz c-- length = length + cell current = current + cell c-- c-- All other cases are standard radiator followed by scintillator c-- ELSE c-- cell = megatile IF (FILLED = 1) THEN c-- Create AND Position EMGT z=length+(gap+cell)/2+esec_deltaz c-- xx = current + (gap+cell)/2+esec_deltaz prin2 i_section,is,xx (' mega i_section,is ',2i3,f10.4) Prin1 is, xx (' + ESEC megatile ilayer=',I2,' z=',F12.4) c-- ENDIF c-- length = length + cell + gap current = current + cell + gap c-- IF (PLATE>0) THEN c-- cell = radiator Prin1 is, current + cell/2+esec_deltaz ( ' + ESEC radiator ilayer=',I2,' z=',F12.4 ) Create AND Position ERAD z=length+cell/2+esec_deltaz c-- length = length + cell current = current + cell c-- ENDIF c-- ENDIF c-- ENDDO c-- c-- EndBlock c----------------------------------------------------------------- Block EMGT -- c-- Block EMGT is a 30 degree megatile c-- Material Air Medium Standard Attribute EMGT seen=1 colo=1 lsty=2 c-- Shape CONS dz=megatile/2, rmn1=(current)*tan_low-dd, rmn2=(current+megatile)*tan_low-dd, rmx1=(current)*tan_upp+dup, rmx2=(current+megatile)*tan_upp+dup c-- c-- DO isec=1,nint(emcs_nslices) c-- myPhi = (emcs_nslices/2-isec+0.5)*dphi + esec_jiggle(is) c-- Create AND Position EPER alphaz=myPhi c-- END DO c-- EndBlock c----------------------------------------------------------------- Block EPER -- c-- Block EPER is a 5 degree slice of a 30 degree megatile (subsector) c-- c-- Creates: c-- + ETAR -- The pseudo-rapidity divivisions in the megatiles c-- Material Polystyren Attribute EPER seen=1 colo=1 lsty=1 c-- c-- c-- Shape CONS dz=megatile/2, phi1=emcs_phimin/emcs_nsector, phi2=emcs_phimax/emcs_nsector, rmn1=(current)*tan_low-dd, rmn2=(current+megatile)*tan_low-dd, rmx1=(current)*tan_upp+dup, rmx2=(current+megatile)*tan_upp+dup c-- curcl = current+megatile/2 DO ie = 1, nint(eetr_neta) c-- etabot = eetr_etabin(ie) etatop = eetr_etabin(ie+1) rbot=(curcl)*tanf(etabot) rtop=min((curcl)*tanf(etatop), ((current)*tan_upp+dup)) c-- check rbot<rtop c-- xx=tan(pi*emcs_phimax/180.0/emcs_nsector) yy=cos(pi*emcs_phimax/180.0/emcs_nsector) Create and Position ETAR x=(rbot+rtop)/2 ort=yzx prin2 ie,etatop,etabot,rbot,rtop (' EPER : ie,etatop,etabot,rbot,rtop ',i3,4f12.4) c-- ENDDO c-- EndBlock c----------------------------------------------------------------- Block ETAR -- c-- c-- ETAR is a single cell of scintillator, including fiber router, plastic, c-- etc... c-- c-- local z is radially outward in star c-- local y is the thickness of the layer c-- Block ETAR is a single calorimeter cell, containing scintillator, fiber router, etc... c-- Material POLYSTYREN Attribute ETAR seen=1 colo=4 lsty=1 ! BLUE c-- Shape TRD1 dy=megatile/2 dz=(rtop-rbot)/2, dx1=rbot*xx-emcs_gapcel/yy, dx2=rtop*xx-emcs_gapcel/yy c-- Create AND Position EALP y=(-megatile+emcs_alincell)/2 g10 = esec_scint Create AND Position ESCI y=(-megatile+g10)/2+emcs_alincell _ +emcs_frplast c-- EndBlock c----------------------------------------------------------------- Block ESCI -- c-- Block ESCI is the active scintillator (polystyrene) layer c-- c-- Obtain the definition of polystyrene on this line, next line clones Material Polystyren Material Ecal_scint isvol=1 Medium Ecal_active isvol=1 c-- Attribute ESCI seen=1 colo=7 fill=0 lsty=1 ! LIGHTBLUE c-- local z goes along the radius, y is the thickness Shape TRD1 dy=esec_scint/2, dz=(rtop-rbot)/2-emcs_gapcel c-- c-- Call ecal_set_cuts( ag_imed, 'detector' ) c-- c-- HITS ESCI BIRK:0:(0,10) c-- c-- EndBlock c----------------------------------------------------------------- Block ERAD -- c-- Block ERAD is the lead radiator with stainless steel cladding c-- c-- Creates: c-- + ELED -- the business end of the calorimeter... c-- Material STEEL c-- Attribute ERAD seen=1 colo=6 fill=1 lsty=1 ! VIOLET Shape CONS dz=radiator/2, rmn1=(current)*tan_low-dd, rmn2=(current+cell)*tan_low-dd, rmx1=(current)*tan_upp+dup, rmx2=(current+radiator)*tan_upp+dup c-- Create AND Position ELED c-- EndBlock c------------------------------------------------------------------------- c----------------------------------------------------------------- Block ELED -- c-- Block ELED is a lead absorber plate c-- c-- Material PbAlloy Medium Ecal_lead Attribute ELED seen=1 colo=4 fill=1 lsty=1 c-- Shape TUBS dz=emcs_pbplate/2, rmin=(current)*tan_low, rmax=(current+emcs_pbplate)*tan_upp, c-- Call ecal_set_cuts( ag_imed, 'radiator' ) c-- EndBlock c-- c----------------------------------------------------------------------------- c----------------------------------------------------------------- Block EFLP -- c-- Block EFLP is the aluminum (aluminium) front plate of the endcap c-- Material ALUMINIUM Attribute EFLP seen=1 colo=3 fill=1 lsty=1 ! GREEN Shape CONS dz=emcs_front/2, rmn1=68.813 rmn2=68.813, rmx1=(zslice)*tan_upp+dup, rmx2=(zslice+zwidth)*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- EndBlock c----------------------------------------------------------------------------- c----------------------------------------------------------------- Block EALP -- c-- Block EALP is the thin aluminium plate in calorimeter cell c-- c-- Material Aluminium Attribute EALP seen=1 colo=1 lsty=1 c-- c-- Shape TRD1 dy=emcs_alincell/2 dz=(rtop-rbot)/2 c-- c-- Thin aluminium plate in each calorimeter cell. The energy-loss c-- fluctuations are restricted in this thin material. c-- CALL GsTPar (AG_IMED,'CUTGAM',0.00001) CALL GsTPar (AG_IMED,'CUTELE',0.00001) CALL GsTPar (AG_IMED,'LOSS',1.) CALL GsTPar (AG_IMED,'STRA',1.) c-- EndBlock c----------------------------------------------------------------- Block ESPL -- c-- Block ESPL is the logical volume containing an SMD plane c-- Material Air Attribute ESPL seen=1 colo=4 lsty=4 Shape TUBS dz=emcs_gapsmd/3/2, rmin=section*tan_low-1.526, rmax=(section-secwid/2)*tan_upp+dup, phi1=emcs_phimin phi2=emcs_phimax c-- USE EMXG version=1 msecwd = (emxg_sapex+emxg_f4)/2 c-- ^^^^^^ what is this used for? --jw c-- looks like the g10 layer which we are retiring c-- c-- loop over the six sectors in an endcap half wheel c-- DO isec=1,6 cut=1 d3 = 75 - (isec-1)*30 c-- IF (exse_sectype(isec)=0|(emcg_fillmode=1&(isec=6|isec=1))) THEN cut = 0 c -- come back and build spacers -- ElseIF (exse_sectype(isec) = 1) then ! v c-- Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' c-- ElseIF (exse_sectype(isec) = 2) then ! u c-- Create and Position EXSG alphaz=d3 ort=x-y-z ncopy=isec kOnly='MANY' c-- ElseIF (exse_sectype(isec) = 3) then ! cut v c-- cut=2 Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' c-- ElseIF (exse_sectype(isec) = 4) then ! cut u c-- cut=2 Create and Position EXSG alphaz=d3 ort=x-y-z ncopy=isec kOnly='MANY' c-- EndIF c-- EndDO! loop over six sectors in eemc half wheel c-- c-- repeat the loop and add in the spacer layers c-- DO isec=1,6 d3=75 - (isec-1)*30 IF (exse_sectype(isec)=0|(emcg_fillmode=1&(isec=6|isec=1))) then cut = 0 c-- Do not kill neighbors Create and Position EXSG alphaz=d3 ncopy=isec kOnly='MANY' c ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ c potential side effect... may screw up the mapping c of the smd strips into the tables? c EndIF EndDO c-- EndBlock c----------------------------------------------------------------- Block EXSG -- c-- Block EXSG Is another logical volume... this one acutally creates the planes c-- c-- Creates: c-- + EHMS -- shower max strips c-- + EFLS -- front cover for SMD planes c-- + EBLS -- back cover for SMD planes c-- Attribute EXSG seen=1 colo=7 serial=cut lsty=3 ! MEH Material Air c$$$ Medium TMED_EXSG stemax=0.01 Shape TUBS dz=emcs_gapsmd/3/2, rmin=section*tan_low-1.526, rmax=(section-secwid/2)*tan_upp+dup, phi1=emcs_phimin/emcs_nsupsec-5, phi2=emcs_phimax/emcs_nsupsec+5 c-- rbot = emxg_rin rtop = emxg_rout c-- c-- Code to handle smd spacers c-- IF ( cut .eq. 0 ) THEN Create and Position EXPS kONLY='MANY' ENDIF c-- c-- Code to handle smd planes c-- IF (cut > 0) THEN c-- c-- setup which plane we are utilizing c-- IF (cut = 1) THEN nstr = 288 ELSE nstr = 285 ENDIF c-- c-- loop over all smd strips and place them w/in this smd plane c-- DO istrip = 1,nstr c-- Call ecal_get_strip( section, cut, istrip, xc, yc, length ) c-- IF (mod(istrip,2) != 0 ) THEN Create and Position EHMS x=xc y=yc alphaz=-45 kOnly='ONLY' Create and Position EBLS x=xc y=yc z=(+esmd_apex/2+esmd_back_layer/2) alphaz=-45 kOnly='ONLY' ELSE Create and Position EHMS x=xc y=yc alphaz=-45 ort=x-y-z kOnly='ONLY' Create and Position EFLS x=xc y=yc z=(-esmd_apex/2-esmd_front_layer/2) alphaz=-45 ort=x-y-z kOnly='ONLY' ENDIF c-- Prin1 istrip, xc, yc, length ( 'SMD Plane: strip=',I3,' xc=',F5.1,' yc=,'F5.1,' length=',F5.1 ) c-- ENDDO c-- ENDIF c-- c-- * dcut exsg z 0 0 10 0.1 0.1 * dcut exsg y 0 10 -50 0.7 0.7 c-- EndBlock c-- c-- c----------------------------------------------------------------------------- c----------------------------------------------------------------- Block EHMS -- c-- Block EHMS defines the triangular SMD strips c-- Material Ecal_scint Medium Ecal_active isvol=1 Attribute EHMS seen=1 colo=2 serial=cut lsty=1 ! red c-- Shape TRD1 dx1=0 dx2=emxg_Sbase/2 dy=length/2 dz=emxg_Sapex/2 c-- HITS EHMS Birk:0:(0,10) c-- Endblock! EHMS c----------------------------------------------------------------------------- c--- c-- Several thin layers of material are applied to the front and back of the c-- SMD planes to provide structural support. We combine these layers into c-- a single effective volume, which is affixed to the base of the SMD c-- strips. As with the SMD strips, z along the depth, y is length c-- c-- http://drupal.star.bnl.gov/STAR/system/files/SMD_module_stack.pdf c-- c-- 1.19 mm G10 c-- 0.25 mm Fiberglass and epoxy c-- 0.17 mm Aluminized mylar c-- c-- Weight in mixture by mass = (depth)*(Area) c-- c-- Weighted density is given by sum (density)_i * (depth)_i / sum (depth)_i c-- c----------------------------------------------------------------- Block EFLS -- c-- Block EFLS is the layer of material on the front of the SMD planes c-- c-- Component G10 A=18.017 Z=9.013 w=1.19*1.700/(1.19*1.700+0.25*1.530+0.17*1.390) Component Fiberglass A=19.103 Z=9.549 w=0.25*1.530/(1.19*1.700+0.25*1.530+0.17*1.390) Component AlMylar A=12.889 Z=6.465 w=0.17*1.390/(1.19*1.700+0.25*1.530+0.17*1.390) Mixture EFLS dens=(1.19*1.7+0.25*1.53+0.17*1.39)/(1.19+0.25+0.17) Attribute EFLS seen=1 colo=22 lsty=1 Shape TRD1 dz=esmd_front_layer/2 dy=length/2 dx1=esmd_base/2 dx2=esmd_base/2 c-- EndBlock! EFLS c-- c-- see link above for documentation c-- c-- 0.10 mm aluminized mylar c-- 0.25 mm fiberglass and epoxy c-- 1.50 mm WLS fiber router layer (polystyrene) c-- 0.25 mm aluminum c-- c----------------------------------------------------------------- Block EBLS -- c-- Block EBLS is the layer of material on the back of the SMD planes c-- Component AlMylar A=12.889 Z=6.465 w=0.10*1.390/(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699) Component Fiberglass A=19.103 Z=9.549 w=0.25*1.530/(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699) Component Polystyren A=11.154 Z=5.615 w=1.50*1.032/(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699) Component Al A=28.08 Z=14.00 w=0.25*2.699/(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699) Mixture EBLS dens=(0.10*1.390+0.25*1.530+1.50*1.032+0.25*2.699)/(0.10+0.25+1.50+0.25) c-- Attribute EFLS seen=1 colo=22 lsty=1 Shape TRD1 dz=esmd_back_layer/2 dy=length/2 dx1=esmd_base/2 dx2=esmd_base/2 c-- EndBlock! EFLS c----------------------------------------------------------------- Block EXPS -- c-- Block EXPS is the plastic spacer in the shower maximum section c-- c-- Simple implementation of the spacer in the shwoer maximum detector. c-- This implmentation neglects the ears and the source tube. c-- c-- n.b. There may be a side effect in the way this gets created... c-- it could overwrite SMD strips which extend into this plane. c-- Probably need to go with a different approach here. c-- c-- Scanned Drawings: c-- + http://drupal.star.bnl.gov/STAR/system/files/SMD_spacer_drawings.pdf c-- c-- thickness is 1.2 cm, as given by detail B and C... but I do not want c-- to do alot of complicated recoding of the geometry. So I am limiting c-- it to be the same width as a normal SMD volume. c-- Material PVC_Spacer Attribute EXPS seen=1 colo=6 lsty=1 lwid=2 c-- c-- Spacer layers are extended by +/- 5 degrees into the adjacent sectors. c-- The kONLY='Many' option at creation time should mean that conflicts c-- in volume will be resolved in favor of the SMD strips. c-- Shape TUBS dz=esmd_apex/2, rmin=(section)*Tan_Low-1.526, rmax=(section+msecwd)*Tan_Upp, phi1=emcs_PhiMin/emcs_Nsupsec, phi2=emcs_PhiMax/emcs_Nsupsec c-- EndBlock c-- END c----------------------------------------------------------------- End Module -- c------------------------------------------------------------------------------ c-- Helper subroutines and functions c------------------------------------------------------------------------------ c-- c-- Subroutine ecal_set_cuts(id, medium) c-- c-- id -- integer ID idetifying the current tracking medium c-- medium -- character switch selecting the type of cuts to be c-- used in this tracking volumne c-- c------------------------------------------------------------------------------ Subroutine ecal_set_cuts(id,medium) c-- Implicit NONE Integer id Character medium*(*) c-- Integer radiator, megatile, detector Save radiator, megatile, detector c-- IF ( medium == 'print' ) THEN c-- Write (*,400) radiator Write (*,401) megatile Write (*,402) detector c-- Call GpTMed( +radiator ) Call GpTMed( -megatile ) Call GpTMed( -detector ) c-- Return c-- ENDIF c-- 400 Format('radiator cuts set for ag_imed=',I3) 401 Format('megatile cuts set for ag_imed=',I3) 402 Format('detector cuts set for ag_imed=',I3) c-- c-- c-- Setup common cuts for neutrons, hadrons and muons c-- Call GsTPar (id,'CUTNEU',0.001) Call GsTPar (id,'CUTHAD',0.001) Call GsTPar (id,'CUTMUO',0.001) c-- IF ( medium == 'radiator' ) THEN Call GsTPar (id,'CUTGAM',0.00008) Call GsTPar (id,'CUTELE',0.001) Call GsTPar (id,'BCUTE' ,0.0001) radiator = id C-- c-- ELSEIF ( medium == 'megatile' ) THEN Call GsTPar (id,'CUTGAM',0.00008) Call GsTPar (id,'CUTELE',0.001) Call GsTPar (id,'BCUTE' ,0.0001) megatile = id c-- c-- ELSEIF ( medium == 'detector' ) THEN Call GsTPar (id,'CUTGAM',0.00008) Call GsTPar (id,'CUTELE',0.001) Call GsTPar (id,'BCUTE' ,0.0001) c-- Call GsTPar (id,'BIRK1',1.) Call GsTPar (id,'BIRK2',0.0130) Call GsTPar (id,'BIRK3',9.6E-6) detector = id c-- c-- ELSE Call GsTPar (id,'CUTGAM',0.00008) Call GsTPar (id,'CUTELE',0.001) Call GsTPar (id,'BCUTE' ,0.0001) Write(*,300) 300 Format('Warning: unknown medium[',A20,'] in ecal_set_cuts') c-- c-- ENDIF c-- Return End c----------------------------------------------------------------------- c----------------------------------------------------------------------- c-- c-- Subroutine ecal_get_strip( section, cut, istrip, xcenter, ycenter, length ) c-- in in in out out out Implicit NONE c-- Real section Integer cut ! 0=no plane 1=normal plane 2=cut plane Integer istrip ! strip index Real xcenter ! output Real ycenter ! output Real length ! output c-- Integer nstrips Real rdel ! shift in radius (?) Real rth Real ddn, ddup Real megatile, p c-- Real xleft, yleft, xright, yright Real dxy, xx, yy Real sqrt2, sqrt3 c-- c-- SMD data copied from data structures above c-- Real base, apex Data base, apex / 1.0, 0.7/ !cm c-- Real Rbot, Rtop Data Rbot, Rtop / 77.41, 213.922 / c-- Real EtaMin, EtaMax Data EtaMin, EtaMax / 1.086, 2.000 / c-- Real tan_theta_min, tan_theta_max c-- Real tanf, eta tanf(eta) = tan(2*atan(exp(-eta))) c-- tan_theta_min = tanf( EtaMax ) tan_theta_max = tanf( EtaMin ) c-- IF (cut = 1) THEN rdel = 3.938 nstrips = 288 ELSE rdel = -.475 nstrips = 285 ENDIF c-- xcenter=0. ycenter=0. length=0. c-- IF ( cut = 0 ) THEN RETURN ENDIF c-- sqrt2 = sqrt(2.0) sqrt3 = sqrt(3.0) c-- rth = .53*rdel ! .53 --- tentatavily jcw-- wtf? ddn = sqrt(3.0)*1.713 + rdel ddup = .5*1.846 + 1.713 megatile = base + .01 c-- p = .5*(istrip-1)*megatile + 41.3655 IF (p <= (.5*rbot*sqrt3 + rth)) THEN dxy = 1.9375*sqrt2 xleft = .5*sqrt2*p*(sqrt3 + 1.) - dxy yleft = .5*sqrt2*p*(sqrt3 - 1.) - dxy yright = .5*sqrt2*(sqrt( rbot*rbot - p*p) - p) xright = sqrt2*p + yright ELSEIF ((.5*rbot*sqrt3 + rth) < p <= (.5*rtop + 1.5)) THEN dxy = 1.9375*sqrt2 xleft = .5*sqrt2*p*(sqrt3 + 1.) - dxy yleft = .5*sqrt2*p*(sqrt3 - 1.) - dxy dxy = rdel*sqrt2/sqrt3 yright = .5*sqrt2*p*(1.- 1./sqrt3) xright = sqrt2*p - yright - dxy yright = -yright - dxy ELSEIF (p > (.5*rtop +1.5)) THEN yleft = (sqrt(rtop*rtop - p*p) - p)/sqrt2 xleft = sqrt2*p + yleft dxy = rdel*sqrt2/sqrt3 yright = .5*sqrt2*p*(1.- 1./sqrt3) xright = sqrt2*p - yright - dxy yright = -yright - dxy dxy = 0. c-- IF ((.5*sqrt3*160.- ddn) < p <= (.5*sqrt3*160.+ ddup) ) THEN xcenter = .5*(sqrt3*160.+1.846) ycenter = xcenter - .5*sqrt3*1.713 IF (p > ycenter) THEN dxy = .5*sqrt2*(2/sqrt3*rdel + .5*sqrt3*1.846 +_ sqrt(1.713*1.713 - (p-xcenter)*(p-xcenter))) ELSE dxy = sqrt2/sqrt3*(p - .5*sqrt3* 160. + ddn) ENDIF ELSEIF ((.5*sqrt3*195.- ddn) < p <= (.5*sqrt3*195. + ddup) ) THEN xcenter = .5*(sqrt3*195.+1.846) ycenter = xcenter - .5*sqrt3*1.713 IF (p > ycenter) THEN dxy = .5*sqrt2*(2/sqrt3*rdel + .5*sqrt3*1.846 +_ sqrt(1.713*1.713 - (p-xcenter)*(p-xcenter))) ELSE dxy = sqrt2/sqrt3*(p - .5*sqrt3*195. + ddn) ENDIF ENDIF xright = xright + dxy yright = yright + dxy ENDIF dxy = section*tan_theta_max - rtop xcenter = .5*(xright+xleft) + dxy ycenter = .5*(yright+yleft) xx = .5*sqrt2*(xleft+yleft) yy = .5*sqrt2*(xright+yright) length = xx-yy c-- c-- Return c-- End! Subroutine smd_strip c-- * ---------------------------------------------------------------------------- * ECAL nice views: dcut ecvo x 1 10 -5 .5 .1 * draw emdi 105 0 160 2 13 .2 .1 * draw emdi 120 180 150 1 14 .12 .12 * --------------------------------------------------------------------------- c-- examples of HITS * HITS EHMS Birk:0:(0,10) * xx:16:SH(-250,250) yy:16:(-250,250) zz:16:(-350,350), * px:16:(-100,100) py:16:(-100,100) pz:16:(-100,100), * Slen:16:(0,1.e4) Tof:16:(0,1.e-6) Step:16:(0,100), * none:16: Eloss:0:(0,10) *
Monte-Carlo setup:
Cuts for shower shapes:
Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
All shapes are normalized to 1 at peak (central strip)
Added layer definition from Jason file:
Some comments:
Figure 1: Sampling fraction vs. thrown energy
Figure 3: Sampling fraction vs. thrown energy (left), 2x1/3x3 energy ratio (right)
See legend for details
Figure 4: Shower shapes. See legend for details
Pre-shower bins:
Ep1/Ep2 is the energy deposited in the 1st/2nd EEMC pre-shower layer.
For a single particle MC it is a sum over
all pre-shower tiles in the EEMC with energy of 3 sigma above pedestal.
For eta-meson from pp2006 data the sum is over 3x3 tower patch
Monte-Carlo setup:
Added layer definition from Jason file:
Geometry configurations and notations (shown in the center of the plot):
cross section of 1st SMD plane labeled with "SUV" ordering
Note: u-v ordering scheme can be found here (Fig. 9-11)
Figure 1: Average number of SMD u-strip fired vs. thrown photon's (x,y)
Figure 2:Average number of SMD v-strip fired vs. thrown photon's (x,y)
Monte-Carlo setup:
Geometry configurations and notations (shown in the center of the plot):
Figure 2: Total energy distribution
Figure 3: Shower shapes (left) and shape ratios (right) for 0 < pre-shower1 < 4MeV
For the previous study click here
Monte-Carlo setup:
Geometry configurations and notations (shown in the center of the plot):
Figure 1: Average energy in SMD-u plane vs. position of the thrown photon
SMD v (left) and u (right) sampling fraction (E_smd/E_thrown) vs. E_thrown
Figure 2: Sampling fraction (E_tower^total/E_thrown) vs. position of the thrown photon
Sampling fraction (E_tower^total/E_thrown) vs. E_thrown
Figure 3: Number of towers above threshold vs. position of the thrown photon
Number of towers above threshold vs. E_thrown
Figure 4: (left) Pre-shower1 and (right) Pre-shower2 sampling fraction vs. E_thrown
Figure 5: (left) High tower sampling fraction and (right) residual energy, [E_tot-E_3x3]/E_thrown, vs. E_thrown
Monte-Carlo setup:
Geometry configurations and notations (shown in the center of the plot):
Figure 1: number of post-shower tiles
Figure 2: number of pre-1-shower tiles
Figure 3: number of pre-2-shower tiles
Figure 4: number of towers
Figure 5: Average pre-shower1 energy
Figure 6: Average pre-shower2 energy
Figure 7: Average number of SMD-u strips
Figure 8: Average number of SMD-v strips
Figure 9: Average post-shower energy
Figure 10: Sampling fraction 1x1 vs. thrown energy
Figure 11: Sampling fraction 2x1 vs. thrown energy
Figure 12: Sampling fraction 3x3 vs. thrown energy
Figure 13: Sampling fraction (total energy) vs. thrown energy
Figure 14: Sampling fraction 1x1
Figure 15: Sampling fraction 2x1
Figure 16: Sampling fraction 3x3
Monte-Carlo setup:
Geometry configurations and notations (shown in the center of the plot):
Figure 1: Average number of SMD-u strips
Figure 2: Average number of SMD-v strips
Figure 3: distribution of 1x1 sampling fraction
Figure 4: distribution of 2x1 sampling fraction
Figure 5: distribution of 3x3 sampling fraction
Figure 6: 1x1 sampling fraction vs. thrown energy
Monte-Carlo setup:
Geometry configurations and notations (shown in the center of the plot):
Figure 1: Sampling fraction 1x1
Figure 2: Sampling fraction 2x1
Figure 3: Sampling fraction 3x3
Figure 4: Sampling fraction total energy
Figure 5: Sampling fraction pre1-shower
Figure 6: Sampling fraction pre2-shower
Figure 7: Sampling fraction smd-u
Figure 8: Sampling fraction smd-v
Figure 9: Sampling fraction post-shower
Figure 10: Sampling fraction 1x1 vs. thrown energy
Figure 11: Sampling fraction 2x1 vs. thrown energy
Figure 12: Sampling fraction 3x3 vs. thrown energy
Figure 13: Sampling fraction (tatal energy) vs. thrown energy
FYI: Alice blog on ELED block study
Monte-Carlo setup:
Geometry configurations and notations (shown in the center of the plot):
Figure 1: Sampling fraction 1x1 (up-left), 2x1 (up-right), 3x3 (low-left), total energy (low-right)
Figure 2: Sampling fraction pre1 (up-left), pre2 (up-right), SMD-u (low-left), post (low-right)
Figure 3: Shower shapes (left) and shower shape ratio (right)
Monte-Carlo setup:
Geometry configurations and notations (shown in the center of the plot):
Figure 1: Distribution of the sampling fraction (total energy in EEMC)
Figure 2: Sampling fraction (total energy in EEMC) vs. thrown energy
Figure 3: Sampling fraction (total energy in EEMC) vs. position of the thrown photon
Monte-Carlo setup:
Geometry configurations and notations:
data base settings (same settings in bfc.C (Jan's trick) and in my MuDst reader):
dbMk->SetFlavor("sim","bemcPed");
dbMk->SetFlavor("Wbose","bemcCalib");
dbMk->SetFlavor("sim","bemcGain");
dbMk->SetFlavor("sim","bemcStatus");
dbMk->SetFlavor("sim","bprsPed");
dbMk->SetFlavor("Wbose","bprsCalib");
dbMk->SetFlavor("sim","bprsGain");
dbMk->SetFlavor("sim","bprsStatus");
dbMk->SetFlavor("sim","bsmdePed");
dbMk->SetFlavor("Wbose","bsmdeCalib");
dbMk->SetFlavor("sim","bsmdeGain");
dbMk->SetFlavor("sim","bsmdeStatus");
dbMk->SetFlavor("sim","bsmdpPed");
dbMk->SetFlavor("Wbose","bsmdpCalib");
dbMk->SetFlavor("sim","bsmdpGain");
dbMk->SetFlavor("sim","bsmdpStatus");
Note: for BEMC ideal pedSigma set to 0, so effectively
there is no effect when I apply 3-sigma threshold above pedestal.
Figure 1: E_reco/E_thrown distribution.
E_reco is the total energy in the BEMC towers from mMuDstMaker->muDst()->muEmcCollection()
E_thrown energy of the thrown photon from tne GEant record
No cut (yet) applied to exclude otliers in the average
Outliers in E_reco/E_thrown
Figure 2: Average E_reco/E_thrown vs. thrown photon eta (left) and phi (right)
Average is taken over a slice in eta or phi (no gaussian fits)
Figure 3: Average E_reco/E_thrown vs. thrown position (eta and phi)
Left: without LOW_EM option; right: with LOW_EM option
No cut applied to exclude otliers
Monte-Carlo setup:
Figure 1: Sampling fraction (total energy in EEMC)
Figure 2: Sampling fraction (total energy in EEMC)
black: same black as in Fig. 1, upper plots
red: EEMC geometry with Material PbAlloy isvol=0
(modification suggested by Jason in this post)
Monte-Carlo setup:
Geometry configurations
Figure 1: Tower sampling fraction distribution
Figure 2: Tower sampling fraction vs. thrown energy
Figure 3: Tower sampling fraction vs. position of the thrown photon
Figure 4: Pre1, pre2, post and SMD sampling fraction distribution
Figure 5: Pre1, pre2, post and SMD sampling fraction vs. thrown energy
Figure 6: SMD-v shower shapes
Figure 7: SMD-v shower shape ratios
Figure 8: Number of SMD-u strips
Figure 9: Number of SMD-v strips
Figure 10: Energy ractio of 2x1 to 3x3 cluster vs. gamma-jet data
Figure 11: Pre-shower1 energy
Figure 12: Pre-shower2 energy
Figure 13: Post-shower energy
Figure 14: SMD-v energy
Figure 15: Number of towers
Figure 16: Tower Sampling fraction: LOW_EM option and pre-shower migration
Monte-Carlo setup:
Geometry configurations and notations:
data base settings (same settings in bfc.C (Jan's trick) and in my MuDst reader):
dbMk->SetFlavor("sim","bemcPed");
dbMk->SetFlavor("Wbose","bemcCalib");
dbMk->SetFlavor("sim","bemcGain");
dbMk->SetFlavor("sim","bemcStatus");
Note: for BEMC ideal pedSigma set to 0, so effectively
there is no effect when I apply 3-sigma threshold above pedestal.
Figure 1: Sampling fraction (0.07*E_reco/E_thrown) distribution: average vs. gaussian fit
E_reco is the total energy in the BEMC towers from mMuDstMaker->muDst()->muEmcCollection()
E_thrown energy of the thrown photon from tne GEant record
The difference between fit and using average values is < 0.7%
Figure 2: Otliers vs. eta and phi: (left) no energy reconstrycted, (right) s.f. < 55%
Most outlier are at eta = 0, -1, +1
Figure 3: Sampling fraction (0.07*E_reco/E_thrown) distribution
Effect of LOW_EM cuts
Figure 4: Sampling fraction vs. thrown photon eta (left) and phi (right)
Average is taken over a slice in eta or phi with cut on outliers (events with s.f. < 5.5% rejected)
Figure 5: Sampling fraction vs. thrown position (eta and phi)
Average is taken over a slice in eta or phi with cut on outliers (events with s.f. < 5.5% rejected)
Figure 6: (left) Single tower sampling fraction
and (right) energy ratio of 1x1 cluster to the total BEMC energy
Not much of the effect from LOW_EM cuts on the 1x1 clustering. Need to look at other (2x1, 2x2 clusters)
Monte-Carlo setup:
Geometry configurations
Figure 1: Sampling fraction of various EEMC layers vs. thrown photon energy:
(a) tower s.f.; (b) tower s.f. distribution; (c) pre-shower1; (d) pre-shower2; (e) SMD, (f) post-shower
Figure 2: (left) Shower shapes and (right) shower shape ratios
Monte-Carlo setup:
Geometry configurations
Monte-Carlo setup:
data base settings (same settings in bfc.C (Jan's trick) and in my MuDst reader):
dbMk->SetFlavor("sim","bemcPed");
dbMk->SetFlavor("Wbose","bemcCalib");
dbMk->SetFlavor("sim","bemcGain");
dbMk->SetFlavor("sim","bemcStatus");
Note: for BEMC ideal pedSigma set to 0, so effectively
there is no effect when I apply 3-sigma threshold above pedestal.
Figure 1: Rapidity cuts study (no eta cuts, no cuts on otliers in this figure)
Figure 2: Sampling fraction (0.07*E_reco/E_thrown) distribution
E_reco is the total energy in the BEMC towers from mMuDstMaker->muDst()->muEmcCollection()
E_thrown energy of the thrown photon from tne Geant record
Cuts: |eta| < 0.97 && |eta|>0.01 && s.f. > 0.055
s.f. distribution on the log scale
Monte-Carlo setup:
Geometry configurations
Figure 1: (left) time/event distribution, (right) average time for the particle type
Conclusion: for single particle Monte-Carlo required time in starsim
with LOW_EM option is ~ 2.6 times higher.
Monte-Carlo setup:
Geometry configurations
Figure 1: EEMC sampling fraction (left) distribution (right) vs. thrown photon energy (1.2 < eta < 1.9; no pt cuts)
Figure 2: EEMC sampling fraction (left) distribution (right) vs. thrown photon energy (1.2 < eta < 1.9; pt > 7GeV cut)
Figure 5: Shower shape ratios (v plane)
Figure 6: Shower shape ratios (u plane)
Figure 7: Pre-shower migration (1.2 < eta < 1.9; no pt cuts)
Monte-Carlo setup:
data base settings (same settings in bfc.C (Jan's trick) and in my MuDst reader):
dbMk->SetFlavor("sim","bemcPed");
dbMk->SetFlavor("Wbose","bemcCalib");
dbMk->SetFlavor("sim","bemcGain");
dbMk->SetFlavor("sim","bemcStatus");
Note: for BEMC ideal pedSigma set to 0, so effectively
there is no effect when I apply 3-sigma threshold above pedestal.
Figure 1: Energy sampling of various cluster in the Barrel EMC
E_reco is the total energy in the BEMC towers from mMuDstMaker->muDst()->muEmcCollection()
eta_thrown - rapidity of the thrown photon from the Geant record
Cuts: |eta| < 0.97 && |eta|>0.01 && total energy s.f. > 0.055
Figure 2: Various cluster energy ratios
Conclusions/dicsussion at the emc2 hypernew
http://www.star.bnl.gov/HyperNews-star/get/emc2/3369.html
http://www.star.bnl.gov/HyperNews-star/get/emc2/3375.html
NoCuts: Default STAR geometry EM cuts
Endcap EMC setup is quite non-uniform
(all cuts are set via "Call GSTPAR (ag_imed,'CutName', Value)":
CUTGAM = 0.00001
CUTELE = 0.00001
CUTGAM = 0.00008
CUTELE = 0.001
BCUTE = 0.0001
CUTNEU = 0.001
CUTHAD = 0.001
CUTMUO = 0.001
c-- Define Birks law parameters
BIRK1 = 1.
BIRK2 = 0.013
BIRK3 = 9.6E-6
CUTGAM = 0.00008
CUTELE = 0.001
BCUTE = 0.0001
CUTNEU = 0.001
CUTHAD = 0.001
CUTMUO = 0.001
CUTGAM = 0.00001
CUTELE = 0.00001
LOSS = 1.
STRA = 1.
CUTGAM = 0.00008
CUTELE = 0.001
BCUTE = 0.0001
c-- Define Birks law parameters
BIRK1 = 1.
BIRK2 = 0.0130
BIRK3 = 9.6E-6
CUTGAM = 0.0001
CUTELE = 0.0001
BCUTE = 0.0001
BCUTM = 0.0001
DCUTE = 0.0001
DCUTM = 0.0001
CUTGAM = 0.00001
CUTELE = 0.00001
BCUTE = 0.00001
BCUTM = 0.00001
DCUTE = 0.0001
DCUTM = 0.00001
CUTGAM = 0.00001
CUTELE = 0.00001
BCUTE = 0.00001
BCUTM = 0.00001
DCUTE = 0.00001
DCUTM = 0.00001
Figure 1: Endcap EMC sampling fraction for different cluster sizes:
1x1, 2x1, 3x3, and total energy in the EEMC
Lower right plot shows total s.f. vs. photon thrown energy
Conclusions/dicsussion at the emc2 hypernew:
http://www.star.bnl.gov/HyperNews-star/get/emc2/3374.html
List of LOW_EM cuts and defaults
Pythia QCD Monte-Carlo:
Figure 1: QCD Total (GEANT/GSTAR+bfc) timing (seconds/event)
Figure 2: QCD GEANT/GSTAR timing (seconds/event)
Figure 3: QCD bfc.C timing (seconds/event)
EEMC single photons Monte-Carlo
Figure 4: EEMC single photon Total (GEANT/GSTAR+bfc) timing (seconds/event)
Figure 5: EEMC single photon GEANT/GSTAR timing (seconds/event)
Figure 1: (left) Endcap EMC sampling fraction (total calorimeter energy), (right) SMD-u sampling fraction
Red: (previous) ecalgeo-v6.1 with global LOW_EM option
(Note: same points as in this post, Fig. 1 lower left, label y6:LOW_EM)
Black: (new) ecalgeo-v6.2 (embedded LOW_EM cuts), no global LOW_EM option
Figure 2: Pre-shower migrations
There is only a few events with pre1>4MeV with new simulations: potential problem with TPC geometry?
Results: Update for the previous tests of EMC v6.2 geometry after fixing TPC/EEMC overlap
Figure 1: Endcap EMC sampling fraction: total calorimeter energy, pre1-, pre2-, post- shower layers, and SMD-u energy
Red: (previous) ecalgeo-v6.1 with global LOW_EM option
(Note: same points as in this post, Fig. 1 lower left, label y6:LOW_EM)
Black: (new) ecalgeo-v6.2 (embedded LOW_EM cuts), no global LOW_EM option
Figure 2: Pre-shower migrations
Change in TPC geometry seems to introduce a reasonable (small) change in pre-shower migration
* Update version of the previous simulation
request from December 18, 2008 (see Ref. [1])
Understanding effects of trigger, material budget differences,
and throughout comparison between 2006 and 2009 data
requires to have dedicated Monte-Carlo
data samples with both y2006 and y2009 geometries.
We request to produce the following set of
Monte-Carlo samples for the photon-jet analysis:
Dedicated (gamma filtered, Refs. [2-5]) data sample
for Pythia pp@200GeV prompt photon processes
with y2006 STAR geometry configuration
and partonic pt range 2-25GeV.
Simulations configured with:
LOW_EM option in starsim (Ref. [6]).
Low cuts on electromagnetic processes in GSTAR
y2006h geometry tag, which includes
latest Endcap EMC (v6.1) and TPC (v4) geometry fixes.
Pythia 6.4 CDF Tune A or Perugia tunes (6.4.22)?
Dedicated (gamma filtered, Refs. [2-5]) data sample
for Pythia pp@200GeV hard QCD processes
with y2006 STAR geometry configuration
and partonic pt range 2-25GeV.
Same simulation setup as for the sample S1.
Pythia pp@200GeV prompt photon and hard QCD
processes with y2009 STAR geometry configuration
and partonic pt range 2-25GeV.
Same simulation setup as for the sample S1
but with y2009a geometry tag.
S4: 3rd priority
Pythia pp@500GeV prompt photon and hard QCD
processes with y2009 STAR geometry configuration
and partonic pt range 2-25GeV.
Same simulation setup as for the sample S1
but with y2009a geometry tag.
Below I provide some estimates of CPU and disk space
which are required to produce the data samples listed above.
These estimates are based on the previous (private)
production of the MC gamma-filtered events with y2006
geometry which was done at MIT computer cluster
by Michael Betancourt (Ref. [2,4-5]):
Pythia pp@200GeV prompt photon simulations
with ~7 pb^-1 luminosity:
~60 days running time on a single CPU
~17Gb of disk space to store MuDst/geant files
Number of (filtered) events:
~ 30K for pt range 6-9GeV
~ 15K for pt range 9-15GeV
Pythia pp@200 QCD hard process simulations
with (at least) 1 pb^-1 luminosity:
~ 620 days running on a single CPU
(less than a week on a cluster with 100 CPUs)
~ 150Gb of disk space to store MuDst/geant files
Number of (filtered) events:
~ 650K for pt range 6-9GeV
~ 300K for pt range 9-15GeV
Notes on the estimates:
Enabling LOW_EM option in GSTAR increases
the time estimates by ~40% (Ref. [7]).
Additional production of jet trees will
require a disk space on the order of < 2%
of the total size of the MuDst/geant files.
Additional production of gamma trees will also
require a disk space on the order of a few percents
of the total size of the MuDst/geant files.
Previous simulation request (Date: 2008, Dec 18):
http://www.star.bnl.gov/HyperNews-star/protected/get/starspin/3596.html
Michael's document on
"Targeted MC procedure for the gamma-jet program at STAR":
http://drupal.star.bnl.gov/STAR/system/files/20080729_gammaFilter_by_MichaelBetancourt.pdf
simulations with filtering readiness:
http://www.star.bnl.gov/HyperNews-star/protected/get/starsimu/387/1/1/2/1/1/2/3/1.html
Filtered photon production with y2006 geometry:
http://www.star.bnl.gov/HyperNews-star/protected/get/phana/256.html
More details on statistics needed and disk space estimates:
http://www.star.bnl.gov/HyperNews-star/protected/get/phana/297.html
LOW_EM option in GSTAR:
http://www.star.bnl.gov/HyperNews-star/protected/get/phana/371.html
Time estimates with and without LOW_EM option:
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
Geometry configurations:
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 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
Click here for results before phi cuts
Geometry configurations:
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
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
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
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:
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
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 |
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:
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 |
Click here for previous study before TPC typo fix
Geometry configurations:
Figure 1: Pre-shower migration (1.2 < eta < 1.9)
Geometry configurations:
Figure 2: Pre-shower migration (1.2 < eta < 1.9): y2006h vs. y2009a
W test sample from Lidia/Jason. MuDst's from /star/data08/users/starreco/recowtest/
Two channels being analyzed:
Figure 1: Lepton yield vs. rapidity (no cuts)
Figure 2: Lepton yield vs. pt and energy
(left) no rapidity cuts
(right) |lepton_eta| < 1
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
Simulations: official pp 500GeV pythia W production
Two channels being analyzed:
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
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
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 (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 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 |
Figure 1: Pythia Eemc gamma filter QA:
(left) False rejection, (right) fraction of accepted events
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
Figure 1: Fraction of accepted events
Figure 2: False rejection (Y-axis scale is 10^-3)
Figure 3: Fraction of accepted events (relative to triggered events)
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)
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
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
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 |