290 likes | 480 Views
NCC Task Force. Richard Seto NCC Task Force Meeting Dec 4, 2007 BNL. Agenda. Reports Its Dec 4–asses whether Dec 20 is a reasonable goal for completion. Report from Dave/Stefan – pi0, dir photon Carla/Sky - chic Ken – delta G
E N D
NCC Task Force Richard Seto NCC Task Force Meeting Dec 4, 2007 BNL
Agenda • Reports • Its Dec 4–asses whether Dec 20 is a reasonable goal for completion. Report from • Dave/Stefan – pi0, dir photon • Carla/Sky - chic • Ken – delta G • Edouard/Ondrej – general state of reconstruction/calibrations • The question is can we start the following: • Some folks will start to write: • Goal? Doc to collaboration, mid -Jan • Doc to DOE – end of Jan, before QM • I will have to give an update for the Wed core week meeting so please provide me with plots by next Tues meeting
Goals • Started with ncc bullets • The goals in simplest terms for the NCC task force is to produce the following figure of merit (FOM) plots from realistic physics simulations. Realistic means a full simulation using PISA with detector response, reconstruction, and data analysis. • a) R_AA for pi0 as a function of pt for various centralities in AuAu and dAu collisions • b) Direct photons cross sections, and likewise R_AA. • c) R_AA for the Chi_C via its decay to J/psi photon as a function of pt for various centralities. Very likely we will limit this to rapidity 1-1.5 • d) A_LL(pt) for direct photons in 500 GeV polarized pp collisions • d') Because we already have jet reconstruction going in pp collisions we should easily be able to use this information to get A_LL(x) giving a greater sensitivity to low-x • e) a demonstration that we can we can extract pi0's at various angles in various centrality collisions with a known efficiency - This should answer the projective Geometry question.
ngammas, from edouards new code, a new chi2 from Ondrej richard seto Dec 4
Ngammas black – mc ngammas red=ngammas where I cut chi2<10 black- from MC red – from NCC edouards new code black – mc ngammas red=ntracks cut E>1000 MeV cut E>300 MeV
Ngammas +Ondrej new chi2 black – mc ngammas red=ngammas where I cut chi2<10 black- from MC red – from NCC edouards new code ondrej’s new code conclusion chi2<10 good Low energy need recalibration black – mc ngammas red=ntracks cut E>1000 MeV cut E>300 MeV
chi2=Lchi2+Tchi2 I think Y=1.5-2.0 Y=1.0-1.5 from this I think maybe I should cut at ln(chi2)=4 or chi2=50 instead of 10? ln(chi2) ln(chi2) black – all tracks from pp min bias file colored lines are from pi0 only files red photons from pt_pi0=1-1.5 blue photons from pt_pi0=5-5.5 magenta photons from pi0=9-9.5 except y=2-2.5 where its 7.5-8 Y=2.0-2.5 ln(chi2)
chi2=Lchi2+Tchi2 I think Y=1.5-2.0 Y=1.0-1.5 Ondrej – new chi2 from this I think maybe I should cut at ln(chi2)=4 or chi2=50 instead of 10? ANSWER:NO now chi2 kill bkg more effectively ln(chi2) ln(chi2) black – all tracks from pp min bias file colored lines are from pi0 only files red photons from pt_pi0=1-1.5 blue photons from pt_pi0=5-5.5 magenta photons from pi0=9-9.5 except y=2-2.5 where its 7.5-8 Y=2.0-2.5 ln(chi2)
Lchi2 Y=1.5-2.0 Y=1.0-1.5 ln(chi2) ln(chi2) black – all tracks from pp min bias file colored lines are from pi0 only files red photons from pt_pi0=1-1.5 blue photons from pt_pi0=5-5.5 magenta photons from pi0=9-9.5 except y=2-2.5 where its 7.5-8 Y=2.0-2.5 ln(chi2)
Lchi2 Y=1.5-2.0 Y=1.0-1.5 Ondrej – new chi2 kills bkg more effectively ln(chi2) ln(chi2) black – all tracks from pp min bias file colored lines are from pi0 only files red photons from pt_pi0=1-1.5 blue photons from pt_pi0=5-5.5 magenta photons from pi0=9-9.5 except y=2-2.5 where its 7.5-8 Y=2.0-2.5 ln(chi2)
Tchi2 Y=1.5-2.0 Y=1.0-1.5 ln(chi2) ln(chi2) black – all tracks from pp min bias file colored lines are from pi0 only files red photons from pt_pi0=1-1.5 blue photons from pt_pi0=5-5.5 magenta photons from pi0=9-9.5 except y=2-2.5 where its 7.5-8 Y=2.0-2.5 ln(chi2)
Tchi2 Y=1.5-2.0 Y=1.0-1.5 Ondrej – new chi2 Start thinking about cutting separately on Lchi2 and Tchi2 ln(chi2) ln(chi2) black – all tracks from pp min bias file colored lines are from pi0 only files red photons from pt_pi0=1-1.5 blue photons from pt_pi0=5-5.5 magenta photons from pi0=9-9.5 except y=2-2.5 where its 7.5-8 Y=2.0-2.5 ln(chi2)
2 gamma pi0After Edouards new code Dec 2BEFORE Ondrej chi2 Cuts: 0 – chisq<10 1 – eta<1.5 chisq<10 2- pt>200 chisq<10 3- eta<1.5 pt>200 chisq<10 good only for y=1-1.25
ppmb y=1-1.5 embedded events subtracted fit signal on full bkg pt=1-1.5 cut3 pt=3-3.5 cut3 pt==2-2.5 cut3 black signal embedded in bkg red – mixed events blue – bkg from events with no embedded pi0 cut is so effective we need to run much more bkg
ppmb higher y embedded events subtracted fit signal on full bkg y=1.5-2 pt=1-1.5 cut2 y=2-2.5 pt=1-1.5 cut2 y=1.5-2 pt=2-2.5 cut2 black signal embedded in bkg red – mixed events blue – bkg from events with no embedded pi0 as high in pt as we can go for 2 gamma at y=1.5-2 as high in pt as we can go for 2 gamma at y=2-2.5
ppmb :inv yield with error bars, pi0 mean/width S/B S/B pi0 width pi0 mass pi0 width pi0 mass y=1-1.5 y=2-2.5 y=1-1.5 Central arm pp pT/GeV S/B ----------- 1 1.5 2 7
ppmb eff all pt bins correct pt bin higher pt bin lower pt bin y=2.0-2.5 y=1-1.5 y=1.5-2.0
dAu mb y=1-1.5 embedded events subtracted fit signal on full bkg pt=1-1.5 cut3 pt=3-3.5 cut3 pt==2-2.5 cut3 black signal embedded in bkg red – mixed events blue – bkg from events with no embedded pi0
dAu mb higher y subtracted subtracted embedded events embedded events signal on full bkg fit signal on full bkg fit y=1.5-2 pt=1-1.5 cut2 y=1.4-2 pt=1.5-3 cut2 black signal embedded in bkg red – mixed events blue – bkg from events with no embedded pi0
dAu mb y=1-1.5 eff all pt bins S/B eff correct pt bin pi0 width pi0 mass efff higher pt bin eff lower pt bin
dAu mb y=1.5-2 eff all pt bins S/B eff correct pt bin pi0 width pi0 mass efff higher pt bin eff lower pt bin
dAu cent y=1.0-1.5 subtracted subtracted embedded events embedded events signal on full bkg fit signal on full bkg fit pt=3.5-4 pt=1-1.5
dAu cent higher y subtracted embedded events signal on full bkg fit y=1.5-2 pt=1-1.5
dAu cent y=1.0-1.5 eff all pt bins S/B eff correct pt bin pi0 width pi0 mass efff higher pt bin eff lower pt bin y=1-1.5
dAu cent 1.5-2.0 all pt bins S/B correct pt bin pi0 width pi0 mass higher pt bin lower pt bin y=1-1.5
AuAu ~50% subtracted subtracted embedded events embedded events signal on full bkg signal on full bkg fit fit
AuAu 50% higher y subtracted embedded events signal on full bkg fit y=1-1.5 pt=1-1.5
AuAu 50% y=1-1.5 all pt bins S/B correct pt bin pi0 width pi0 mass higher pt bin lower pt bin y=1-1.5
plans • 2 track pi0’s - limited pt range • Central ran at ISU last nite – low bkg statistics (need to make plots) • use bigger stats for pp bkg • get 500 GeV S/B etc • run everything with final cuts • Get final R_AA plots out • Big step – merge with single track data