1 / 15

Robert Ciesielski (DESY)

CC e+p Analysis Status. Robert Ciesielski (DESY). Introduction Event selection Control plots Summary and Outlook. High-Q2 Group Meeting 30/05/2007. Intro. Analysis goal: CC e+p 2006-2007 data with L=~138 pb-1 Last report: first look at 2006 e+p RH sample Progress since last report:

cecil
Download Presentation

Robert Ciesielski (DESY)

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. CC e+p Analysis Status Robert Ciesielski (DESY) • Introduction • Event selection • Control plots • Summary and Outlook High-Q2 Group Meeting 30/05/2007

  2. Intro • Analysis goal: • CC e+p 2006-2007 data with L=~138 pb-1 • Last report: • first look at 2006 e+p RH sample • Progress since last report: • Both: RH (L=79.34 pb-1, P= 31.88 %) • LH (L=58.90 pb-1, P = -36.41 %) Final Sample! • Energy Scale correction to data • Reweighted 2006 vertex distribution with stand-alone version routine (thanks Shima!) CC e+p Analysis Status

  3. Event Selection MC sample (for the time being): Djangoh, Ariadne, CTEQ5D Q2> 100 GeV2, 2004 e+p, num05t0.2 No other cleaning cuts: Remaing cosmic and halo muons rejected after Visual Scanning CC e+p Analysis Status

  4. TSUBAME halo muon finder (BCAL) 99% 1% Scanned 640 RH and 340 LH events with Tsu_halo>0 Very high purity! 99% of them with clear halo muon signal Only 6RH/4LH are CC candidates with ~1/2/3 jets aligned along phi (pure geometry, MC should follow) CC e+p Analysis Status

  5. Event Scanning (reject cosmic/halo) CC e+p Analysis Status

  6. Final samples & Control Plots Number of events RH LH Data: 2325 846 MC: 2411 863 Selection efficiency: 69.05% (D-MC)/MC: -3.6% -2.0% No cosmic/halo muon in the sample (rejected after Visual Scan) The remaining background small, will be stimated with MC (PhP, Di-tracks) Final BeamGas rejection cuts under study (Amita) Our worry at the moment: what happens at low gamma (<0.7) for RH data?... CC e+p Analysis Status

  7. Control Plots Zvertex RH: LH: Migrations from high to low gamma for positive Zvtx? CC e+p Analysis Status

  8. Control Plots Kinematic variables RH CC e+p Analysis Status

  9. Control Plots Kinematic variables LH Migrations from high to low gamma for medium Q2? CC e+p Analysis Status

  10. Control Plots Highest Et jet RH CC e+p Analysis Status

  11. Control Plots Highest Et jet LH CC e+p Analysis Status

  12. Amita's control Plots background from beam gas removed with tracking cut: gtrk > 0 AND gtrk > 0.2*(trk_ntrack - 24) CC e+p Analysis Status

  13. Summary && Outlook • Working on the final RH+LH data samples • Energy scale corrections applied • 2006 vertex distribution in MC (rewegthing with stand alone routine) • Trying to understand discrepances at low gamma (RH sample) • Data/MC differences in REGULAR Tracking due to STT? • Data: starting from 2006e+p data STT used in Pattern Recognition • MC: STT not used in Mozart num4.0 or num5.0 versions • Data/MC differences: • - tracks/good tracks? (Amita) • - beam gas events overlayed on CC events?... • Once we understand it -> go for total cross sections • But we definitely need the up-to-date vertex and Monte Carlo! CC e+p Analysis Status

  14. Timing cuts bool reject=false; if(Cal_tr_n >1 && Cal_tr_e>0.6 && Tmath::Abs(Cal_tr) >6) reject = true; if(Cal_tb_n>1 && Cal_tb_e>0.6 && TMath::Abs(Cal_tb)>6) reject = true; if(Cal_tf_n >1 && Cal_tf_e >0.6 && Tmath::Abs(Cal_tf) >6) reject = true; if(Cal_tg_n>1 && Cal_tg_e>2.0 && TMath::Abs(Cal_tg)>6) reject = true; if(Cal_tr_n >1 && Cal_tr_e>10. && Tmath::Abs(Cal_tr) >5) reject = true; if(Cal_tb_n>1 && Cal_tb_e>10. && TMath::Abs(Cal_tb)>5) reject = true; if(Cal_tf_n >1 && Cal_tf_e >10. && Tmath::Abs(Cal_tf) >5) reject = true; if(Cal_tg_n>1 && Cal_tg_e>10. && TMath::Abs(Cal_tg)>5) reject = true; if(Cal_tu_n>1 && Cal_tu_e>2.0 && TMath::Abs(Cal_tu)>5) reject = true; if(Cal_td_n>1 && Cal_td_e>2.0 && TMath::Abs(Cal_td)>5) reject = true; if(Cal_tu_n>1 && Cal_td_n>1 && Cal_tu_e>1. && Cal_td_e>1.){ if(TMath::Abs(Cal_td-Cal_tu)>5) reject = true; } if(Cal_tf_n>1 && Cal_tr_n>1 && Cal_tf_e>1. && Cal_tr_e>1.){ if(TMath::Abs(Cal_tf-Cal_tr)>5) reject = true; } return reject; DATA ONLY CC e+p Analysis Status

  15. Cosmic/Halo Muon Rejection FCALShowerShape (halo): bool muon=false; if(Hietcell_nr > 16384) return muon; if(Cc_fclus_dt < =3 && Cc_fclus_dm < =3){ if(Cc_fclus_et < 10) muon = true; }else{ if(Cc_rclus_asoe > 0.2) muon = true; } return muon; RatioEemcHac (cosmic && halo): float eremc = Remc[3]; float erhac = Rhac[3]; float er = eremc + erhac; float ebemc = Bemc[3]; float ebhac = Bhac[3]; float eb = ebemc + ebhac; float efemc = Femc[3]; float efhac = Fhac[3]; float efhac1 = Fhac[3]-Fhac2[3]; float efhac2 = Fhac2[3]; float ef = efemc + efhac; bool muon=false; if(er>2){ if(erhac/er > 0.9) muon = true;} if(eb>2){ if(ebhac/eb > 0.85) muon = true;} if(ef>2){ if(efhac2/ef > 0.6) muon = true; if(efhac1/ef < 0.1 || efhac1/ef > 0.8) muon = true; if(efemc/ef < 0.1 || efemc/ef > 0.9) muon = true; } return muon; EmaxCell (halo): bool muon=false; if(Hietcell_et/Cc_et > 0.7) muon = true; if(Ptmhietcell/Cc_pt < 0.2) muon = true; return muon; CC e+p Analysis Status

More Related