jet analysis Dec.22/2008 ~ 1 visible , without smear, |eta|<1, [0.1,0.1] cells, pp5500GeV JET PROPERTIES 2 Celljet (the Jet-Finding Method) 1. divide η-φ space in [0.1, 0.1] cells 2. calculate transverse energy (eTcell) in each cell 3. select candidates of jet-seed by eTcell > ”eTseed” 4. calculate sum of eTcell in the cone which center positioned at jet-seed (eTsum=ΣeTcell) 5. requirement : eTsum>”Min-eT” 6. define the survivors as found jets Parameters to input are “eTseed”, “coneRadius”, “Min-eT” 3 Distance btw primary parton and celljet center assumption: – – • • • dr<1にprimary partonからのjetが入っている dr<1に他のjetがない proton dr<1の中でcelljet.eTが最大のjetをprimary parton起源のjetとした celljet.coneradius = 1. distance b/w primary parton & celljet center proton parton found jet • dr quark gluon photon integral(0~r) 80% 4 Distance btw celljet center and each tracks • assumption: – – • • • dr<1にprimary partonからのjetが入っている dr<1に他のjetがない dr<1の中でcelljet.eTが最大のjetをprimary parton起源のjetとした celljet.coneradius = 1. distance b/w celljet center & each tracks (eT weighted) quark gluon photon integral(0~r) 80% 5 jet energy balance • dr(distance btw parton & celljet) • coneradius • qq/gg/γγ : jet2 energy/jet1 energy(jet2 energy<jet1 energy) • qγ/gγ : q/g jet energy/photon jet energy dr={0.45(quark),0.55(gluon),0.3(photon)}, coneradius=1.0 dr = 0.4, coneRadius=0.25 qq gg γγ qγ gγ 6 visible/charged/neutral , with smear, |eta|<3, [0.1,0.1] cells, pp5500GeV, quark jet, dr<{0.45/0.55/0.30}, coneRadius=1.0, SINGLE JET ANALYSIS 7 charged/neutral • MineTcharged = MineTvisible*2/3の確認 • energy fraction of each particle seed against with energy sum of generated particles MinBias pTHatMin=200GeV だいたいOK? 8 influence by smear • assumption: – • ALICE-EMC,TPC resolution setting : – – – – dr<0.45(quark)の中でcelljet.eTが最大のjet dr<0.55(gluon)の中でcelljet.eTが最大のjet dr<0.30(photon)の中でcelljet.eTが最大のjet coneRadius=1.0 charged : TPC, Neutral : EMC particle energy(momentum) is smeared assume pion mass EMC TPC+ITS P 0.068 P 0.608[%] P E 8.4 2 2.4 [%] E E 2 9 energy resolution 1. 2. 3. 4. 見つけたcelljet-eT とjet-eT(parton-eT)の相関を作る 直線Fitから補正係数を決定 見つけたcelljet-eTを補正したもの(E’)とjet-eT(parton-eT)(E)から(E’-E)/Eを計算 resolution = (E’-E)/E分布のRMS eT-eT correlation (E’-E)/E Profile resolution 10 jet energy resolution • setting: – – – – visible : MineT = 20 GeV charged : MineT = 20*2/3 GeV neutral : MineT = 20*1/3 GeV eTseed = 0.3 * MineT quark gluon jet resolution がquarkに比べよい? =>low-pT粒子でsmearingの桁落ち? visible resolutionがあまりによい =>計算違い? gluon なんだか怪しいが、計算ミスは見当たらない… photon 11 parton energy resolution • setting: – – – – visible : MineT = 20 GeV charged : MineT = 20*2/3 GeV neutral : MineT = 20*1/3 GeV eTseed = 0.3 * MineT quark visibleで10~20%程度 ほとんどkinematicsで決まっている gluon 前pageと同じ計算なので少し自信がない… photon 12 visible , |eta|<2, [0.1,0.1] cells, embedded jet : jet center eta < 1. coneradius = 1.0 quark jet (pp5500GeV からのB-to-B jet)+(HYDJET HD+J event), JET FINDING METHOD 13 parameter selection(threshold/coneRadius) • BKG subtraction eTcell(sum of eT in cell) < threshold を集めてくる それらのcellを使いmean-eTを求める(mean-eT) mean-eTを中心に、v2に応じたBKGを計算(mean-eT’) 全てのeTcellからmean-eT’を引く pythia-celljetと同じ作業 1. 2. 3. 4. 5. threshold/coneRadiusはどの程度に設定すべきか? (Ejet’-Ejet)/Ejetの幅が最小となる値を見つけたい coneRadius(0.05~0.5) • • jet energy resolution (0.1~0.75) threshold[GeV](0~22) 14 50GeV 100GeV threshold : 15GeV coneRadius : centralityの関数 150GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 16 60-75% parameter selection(Min-eT/eTseed) • S/(S+B) ≡ Nreal/Nfound • efficiency ≡ Nreal/Nembed Nreal : 見つかったpythiaからのB-to-B jetの数 Nfound : celljetで見つけたB-to-B jetの数 Nembed : 埋め込んだB-to-B jetの数 eTseed/Min-eT(0.1~0.9) • S/(S+B), efficiencyの値からMin-eT, eTseedの最適化 • {S/(S+B)}*{efficiency} 隣のcellまで許容 pythia B-to-B jet : dphi-PI 1.5sigma S/(S+B)(0~1) efficiency(0~1) {S/(S+B)}*{efficiency}(0~0.9) Min-eT/pTHat(0.2~1.0) energy resolution B-to-B jet の決め方/答え合わせ方法 |dφ-π|<0.284 dr(pythia-jet/found jet) < 0.11 1.-|Efound jet/Epythia jet| < 0.45 17 single quark jet 50~100GeV {S/(S+B)}*eff efficiency S/(S+B) 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 21 60-75% BtoB quark jet 50~100GeV {S/(S+B)}*eff efficiency S/(S+B) 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 22 60-75% plan • • • • • gluon jet, photon jet でも同様にconeRadius等の最適化 gamma-jetのS/(S+B),efficiencyの計算 gamma-jetのenergy balanceの計算 (method dependence) jet quenchingした上で、parton energy lossの計算(再現性能) 27 visible , |eta|<2, [0.1,0.1] cells, embedded jet : jet center eta < 1. coneradius = 1.0 gluon jet (pp5500GeV からのB-to-B jet)+(HYDJET HD+J event), JET FINDING METHOD 28 50GeV 100GeV threshold : 15GeV coneRadius : centralityの関数 150GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 29 60-75% parameter selection(Min-eT/eTseed) • S/(S+B) ≡ Nreal/Nfound • efficiency ≡ Nreal/Nembed Nreal : 見つかったpythiaからのB-to-B jetの数 Nfound : celljetで見つけたB-to-B jetの数 Nembed : 埋め込んだB-to-B jetの数 eTseed/Min-eT(0.1~0.9) • S/(S+B), efficiencyの値からMin-eT, eTseedの最適化 • {S/(S+B)}*{efficiency} 隣のcellまで許容 pythia B-to-B jet : dphi-PI 1.5sigma S/(S+B)(0~1) efficiency(0~1) {S/(S+B)}*{efficiency}(0~0.9) Min-eT/pTHat(0.2~1.0) B-to-B jet の決め方/答え合わせ方法 |dφ-π|<0.409 dr(pythia-jet/found jet) < 0.11 1.-|Efound jet/Epythia jet| < 0.45 30 single gluon jet 50~100GeV {S/(S+B)}*eff efficiency S/(S+B) 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 31 60-75% BtoB gluon jet 50~100GeV {S/(S+B)}*eff efficiency S/(S+B) 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 32 60-75% jet(parton) energy resolution • • • • coneRadius : function of centrality threshold : 15GeV Min-eT : pTHat * 0.5 eTseed : 0.3 * Min-eT central 上 : (pythia jet eT) vs (celljet eT(corrected)) 下 : (pythia parton eT) vs (celljet eT(corrected)) peripheral 33 jet(parton) energy resolution jet energy resolution parton energy resolution 100GeV jet energy resolution = 25%程度(central) 100GeV parton energy resolution = 25%程度(central) 34 S/(S+B), efficiency, eff*S/(S+B) • S/(S+B) ≡ Nreal/Nfound • efficiency ≡ Nreal/Nembed Nreal : 見つかったpythiaからのB-to-B jetの数 Nfound : celljetで見つけたB-to-B jetの数 Nembed : 埋め込んだB-to-B jetの数 • Min-eT : pTHat * 0.5 • eTseed : 0.3 * Min-eT S/(S+B) efficiency {S/(S+B)}*eff 35 visible , |eta|<2, [0.1,0.1] cells, embedded jet : jet center eta < 1. coneradius = 1.0 photon (pp5500GeV からのB-to-B jet)+(HYDJET HD+J event), JET FINDING METHOD 36 50GeV 100GeV threshold : 15GeV coneRadius : 0.1 150GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 37 60-75% parameter selection(Min-eT/eTseed) • S/(S+B) ≡ Nreal/Nfound • efficiency ≡ Nreal/Nembed Nreal : 見つかったpythiaからのB-to-B jetの数 Nfound : celljetで見つけたB-to-B jetの数 Nembed : 埋め込んだB-to-B jetの数 eTseed/Min-eT(0.1~0.9) • S/(S+B), efficiencyの値からMin-eT, eTseedの最適化 • {S/(S+B)}*{efficiency} 隣のcellまで許容 pythia B-to-B jet : dphi-PI 1.5sigma S/(S+B)(0~1) efficiency(0~1) {S/(S+B)}*{efficiency}(0~0.9) Min-eT/pTHat(0.2~1.0) B-to-B jet の決め方/答え合わせ方法 |dφ-π|<0.169 dr(pythia-jet/found jet) < 0.11 1.-|Efound jet/Epythia jet| < 0.45 38 single photon 50~100GeV {S/(S+B)}*eff efficiency S/(S+B) 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 39 60-75% BtoB photon 50~100GeV {S/(S+B)}*eff efficiency S/(S+B) 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 40 60-75% jet(parton) energy resolution • • • • coneRadius : function of centrality threshold : 15GeV Min-eT : pTHat * 0.5 eTseed : 0.3 * Min-eT central 上 : (pythia jet eT) vs (celljet eT(corrected)) 下 : (pythia parton eT) vs (celljet eT(corrected)) peripheral 41 jet(parton) energy resolution jet energy resolution parton energy resolution 100GeV jet energy resolution = 25%程度(central) 100GeV parton energy resolution = 25%程度(central) 42 S/(S+B), efficiency, eff*S/(S+B) • S/(S+B) ≡ Nreal/Nfound • efficiency ≡ Nreal/Nembed Nreal : 見つかったpythiaからのB-to-B jetの数 Nfound : celljetで見つけたB-to-B jetの数 Nembed : 埋め込んだB-to-B jetの数 • Min-eT : pTHat * 0.5 • eTseed : 0.3 * Min-eT S/(S+B) efficiency {S/(S+B)}*eff 43 BACK UP SLIDES 44 PYTHIA jet resolution correlation of jet eT correlation of parton/jet eT 45 quark gluon photon jet-Not smeared visible charged neutral 46 quark gluon photon jet-smeared visible charged neutral 47 quark gluon photon parton-Not smeared visible charged neutral 48 quark gluon photon parton- smeared visible charged neutral 49 |dφ-π|<0.38 dr(pythia-jet/found jet) < 0.11 1.-|Efound jet/Epythia jet| < 0.45 50 50~100GeV S/(S+B) 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 51 60-75% 50~100GeV efficiency 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 52 60-75% 50~100GeV {S/(S+B)}*eff 100~150GeV Min-eT = 0.5*pTHat eTseed = 0.3*Min-eT 実際にはひとつの値しか使えない 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 53 60-75% |dφ-π|<0.3 dr(pythia-jet/found jet) < 0.11 1.-|Efound jet/Epythia jet| < 0.2 54 50~100GeV S/(S+B) 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 55 60-75% 50~100GeV efficiency 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 56 60-75% 50~100GeV {S/(S+B)}*eff 100~150GeV 150~200GeV 0-10% 10-20% 20-30% 30-40% 40-50% 50-60% 57 60-75% PYTHIA-jetの作り直し • cellJet configuration : R=0.7 • MineT=20GeV, eTseed=6GeV -> jet candidate • primary partonからの距離dr<0.3内で最大のeTの jet(|eta|<1)を選ぶ • BtoBの決め方 – – – – 2つのprimary partonからのjet(jet1,jet2) |eTjet1-eTpr1|/(eTjet1+eTpr1)<0.1 |eTjet2-eTpr2|/(eTjet2+eTpr2)<0.1 |veTjet1+veTjet2|/(eTjet1+eTjet2)<0.2 • celljet, primary parton, final particlesの情報を記録 PYTHIA-jetの作り直し • parameters – – – – – R=0.7 : 充分広く MineT =20GeV : efficiency ,eTseed=6GeV : MineT*0.3 dr=0.3 : (prparton-celljet)のdphi,detaのRMS~0.2 energy fraction 0.1 : (2~3)*RMS BtoB variable 0.2 : (2~3)*RMS quark gluon 50GeV 100GeV 150GeV gluon 0.68 0.63 0.58 quark 0.32 0.37 0.42 その前に... • quark/gluon jet separationは難しそう – quark(gluon) enrich sample • dijet, gamma-jet dijet はなんとか識別methodを作った方がよい gamma-jetはquark-rich 0.73-1.01e-3x 0.27+1.02e-3x BKGの選択 • mean-eT = <eTcell>|(eTcell<threshold)を求め る • mean-eTを中心とした、dphi分布(eTweighted)の幅を持たせてた分布からBKGを 設定 jet energy resolution • celljet energyによってRを変えるのは変なバイアスがかか りそうで気持ち悪い • quark/gluon jetは同じconeRadiusを選ぶ Min-eT, eTseed • single/di-jet のS/(S+B) * efficiencyの最大値 • 本物の決め方(1jet) – PYTHIA-jetとcelljetの距離<0.15^2 – |eTjet-eTcelljet|/eTjet<0.45 • BtoB の決め方 – |dphi-PI|<0.3 – エネルギーの高いjetから選んでいき、一度使ったものは使わない • 本物の決め方(dijet) – (1jet)+(1jet) • 実際には一種類しか使えない – 50GeV conf.(Min-eT=25GeV, eTseed=7.5GeV)で一度jetを見つけ、その jet-eTに対しseed-cell-eT > 0.3*jet-eTを要請 yield • 前頁のdijetのグラフよりそれぞれのcentrality でのjet-yieldを計算 – Ncollで重み付け • 前に求めたefficiencyをかけることで reconstruct可能なjet数が求まる • acceptanceの効果を加えることでannual (reconstucted) yieldが求まる
© Copyright 2025 ExpyDoc