410 likes | 593 Views
Angular resolution study of isolated gamma with GLD detector simulation. 2007/Feb/ ACFA ILC Workshop M1 ICEPP, Tokyo Hitoshi HANO. collaborated with Acfa-Sim-J group. Contents. Introduction Angular Resolution Study position resolution of ECAL cluster
E N D
Angular resolution study of isolated gamma with GLD detector simulation 2007/Feb/ ACFA ILC Workshop M1 ICEPP, Tokyo Hitoshi HANO collaborated with Acfa-Sim-J group
Contents • Introduction • Angular Resolution Study • position resolution of ECAL cluster • direction of reconstructed gamma • Calorimeter Component Dependence • Cell size dependence • Material dependence • Summary
Measurement of the direction of non-pointing photon is important for GMSB (gauge mediated supersymmetry breaking) scenarios. IP Motivation and PFA analysis decay length : [m] • To identify a non-pointing photon, we need to know angular resolution of the detector (EM Calorimeter). We have studied angular resolution using full-simulator (Jupiter) ECAL In this study, we have used single-gamma shot from IP to evaluate angular resolution.
Jupiter for GLD detector GLD detector has large-radius and fine-segmented Calorimeter. Calorimeter cell size and absorber material can be changed. It’s important to optimize cost vs. physics performance. ECAL geometry in Jupiter : R [m] Z [m] Structure W/Scinti./gap3/2/1(mm) x 33 layers cell size 1x1(cm2) barrel endcap 2.1-2.3 0.4-2.3 0-2.8 2.8-3.0 ECAL
Method of reconstructed gamma • Clustering • Find an energy-weighted central point of each layer • Fit each point with least-square method • Evaluate an angle between gamma-line and reconstructed gamma Calorimeter γ IP (generated point) reconstructed gamma
Angular Resolution Study position resolution of cluster
θ(φ)resolution [rad] = θ(φ)meas – θ(φ)MC Method (position resolution study of aaaaaaa each hit cluster) • Shoot single-gamma from IP with random direction • Clustering (more details in next page) • Search energy-weighted central point of cluster • Evaluate θ, φ of a central point • Compare with MC truth ECAL central point (θ,φ) IP (generated point) γ clustering
Clustering method • Find the highest energy deposit cell • Make the cone with a focus on it • Define cells which are inside of the cone as one cluster (around all layers) • Find energy-weighted central point highest energy deposit cell central point clustering angle = 10° IP (generated point) γ@10GeV
Position resolution of cluster (cell : 1 cm) barrel endcap barrel endcap 1 GeV 2 GeV 5 GeV 10 GeV σ [mrad] σ [mrad] |cos(θ)| |cos(θ)| θ resolution is better for larger cos(θ) φ resolution is worse for larger cos(θ) ECAL geometrical effect Position resolution : ~0.3 [cm] IP (generated point)
EnergyDependentResult of position resolution 1GeV 2GeV σ [mrad] σ [mrad] 5GeV 10GeV 1/√E 1/√E θ barrel: [mrad] φ barrel: [mrad] θ endcap: [mrad] φ endcap: [mrad]
Angular Resolution Study direction of reconstructed gamma
Method (Angular resolution study of reconstructed gamma) • Clustering • Find an energy-weighted central point of each layer • Fit each point with least-square method • Evaluate an angle between gamma-line and reconstructed gamma Calorimeter γ IP (shot point) reconstructed gamma
r d Histogram and angular resolution γ central point of cluster IP d r reconstructed gamma angle [rad] = r/d • r histogram F(r) fitting function σ = 48.3 ± 0.3 [mrad] gamma@10GeV
[mrad] Energy dependence (1,2,5,10,50GeV) 1GeV 2GeV σ [mrad] Average over full acceptance 10GeV 5GeV 50GeV 1/√E
IP ECAL Shoot from another point gamma@10GeV • Shoot from x=y=20, z=0 reconstructed gamma • Shoot from IP IP ECAL σ= 48.3±0.3[mrad] σ= 48.6±0.3[mrad] If gamma has been shot from another position, we could not observed significant difference.
Structure (cell size dependence) gamma : E = 10GeV How about cell size dependence?
Cell size dependence 1 [cm] : 48.3 ± 0.3 [mrad] 0.5 [cm] : 46.4 ± 0.3 [mrad] <5% gamma @10GeV We could not observed significant improvement from 1cm to 0.5cm
Structure (energy dependence) gamma : E = 1~10GeV How about energy dependence between 1cm and 0.5cm?
Energy dependence (1,2,5,10GeV) 1GeV 2GeV 5GeV 10GeV No significant difference has been observed between 1cm and 0.5cm around all of energy.
Structure (Absorber dependence) gamma : E = 10GeV How about absorber dependence?
Absorber dependence (Tungsten, Lead) Tungsten[3mm] Lead[4.8mm] Same total radiation length Lead[3mm] Tungsten [3mm] : 48.3 ± 0.3 [mrad] Lead [4.8mm] : 45.5 ± 0.3 [mrad] @1x1 [cm] Angular resolution with Lead is better than Tungsten
Cause reconstructed gamma gamma MC gamma MC reconstructed gamma depth depth Angular resolution is better than Tungsten, since Lead has deeper distribution.
Summary • Angular resolution of default-GLD Calorimeter (W:1cm) • The angular resolution is estimated to be 125mrad/√(E/GeV) • Dependence on cell size granularity and material dependence (W, Pb) has been studied • No significant difference has been observed between 1cm and 0.5cm • Lead is better than Tungsten for isolated gamma • Energy resolution is same • How about energy resolution for jet ? Next talk
θ(φ)resolution [rad] = θ(φ)MC – θ(φ)meas θ, φ resolution study of cluster • Shoot single-gamma from IP with random direction • Clustering - use hit data from ECAL(,HCAL) • Search central point of cluster • Find θ, φ of a central point • Compare with MC truth γ central point IP
Judging by inner product IP clustering angle = 10° IP γ@10GeV Clustering • Find the highest energy deposit cell • Make a cone with centering around it • Define cells which are inside of a cone as one cluster • Find a central point by energy weighted mean
IP 1cell ( 1cm x 1cm ) θ, φ z=0 ( |cos(θ)|=0 ) max θ1cell z=280 ( |cos(θ)|=0.8 ) min θ1cell 1cm 1cm 210cm IP θ(φ)1cell = 1/210 ≒ 4.7 [mrad] θ1cell≒ 1.71 [mrad]
Result (θ, φ resolution) gamma@10GeV • θ ,φ resolution θ1cell ≒ 1.71~4.70 [mrad] θbarrel : 0.430±0.004 [mrad] θendcap: 0.282±0.006 [mrad]φbarrel: 0.423±0.004 [mrad] φendcap: 0.699±0.014 [mrad] Angular resolution is good as well as cell size (1x1cm)
Angular resolution study ofreconstructed line • Clustering • Find a central point of each layer by energy weighted mean • Fit each point with least-square method • Find an angle between IP and reconstructed line γ reconstructed line IP
Fitting method Find a central point of each layer by energy weighted mean Fitting 2-dimentions (x-y) y y’ x Fitting new 2-dimentions (y’-z) weighted by energy deposit z y’ Distance[cm]
2-dimension normal distribution • r histogram F(r)
central point of cluster γ r reconstructed line r d d IP Distance (d) and angle angle [rad] = d/r fitting function
2x2 cm 1x1 cm [mrad] Linearity(1,2,5,10,50 GeV) Linearity is kept below 10GeV.
Angular resolution Absorber (Tungsten, Lead) Tungsten[3mm] Lead[4.8mm] Lead[3mm] Tungsten : 48.26 ± 0.29 [mrad] Lead : 45.51 ± 0.28 [mrad] @1x1 [cm]
Average (10000 events) gamma @10GeV
Hitting distribution and Average Fitting of Lead makes successfully than Tungsten, because Lead has deep distribution. gamma @10GeV
Fitting 0.068*x-14 0.031*x+6.8
Cause reconstructed gamma gamma MC gamma MC reconstructed gamma depth (layer) depth (layer) Since Lead has deeper distribution, angular resolution is better than Tungsten. Average of gamma @10GeV