Download
ific meeting n.
Skip this Video
Loading SlideShow in 5 Seconds..
IFIC meeting PowerPoint Presentation
Download Presentation
IFIC meeting

IFIC meeting

109 Vues Download Presentation
Télécharger la présentation

IFIC meeting

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

  1. IFIC meeting 01/Feb/2012 Agustín Sánchez Losa

  2. BEST LASER vs OFFICIAL OFFSETS

  3. BEST LASER vs OFFICIAL OFFSETS

  4. BEST LASER vs OFFICIAL OFFSETS

  5. BEST LASER vs OFFICIAL OFFSETS

  6. LastFittedSignal Means > nsig0: 0.115321 > nsig1: 0.959034 > nsig3: 2.869490 > nsig5: 4.768660 > nsig7: 6.777240 > nsig9: 8.693610 Declination: 0°; Fixed Search; Systematics ON; Original code and ingredients;

  7. Last Q stat Thresholds: Mine: > th(3s): 4.65 > th(5s): 12.05 Internal note: > th(3s): 3.53 > th(5s): 15.24 Declination: 0°; Fixed Search; Systematics ON; Original code and ingredients;

  8. PEX structure • Simulating experiments: • add_backgroundfrom DATA histograms: • Bgr& bgr.f_per_sindecl_vs_sindecl->GetRandom() (background rate for sin(δ))) • TH1* nhits_atm_histogram(nhits in data) • add_signal_psfhistogram3d from MC histograms: • TH3& psf_3d_histogram (3d histogram of log(α), sin(δ)and nhits for E-2 neutrinos) • Computing likelihood: • lik_sb_with_hits with: • Psf& psf.per_sr_vs_a(alpha) (psf for sin(δ)) • Bgr& bgr.per_sr_vs_decl(decl) (background rate for sin(δ))) • TH1* hist_E2 (nhitsin MC) • TH1* hist_atm(nhits in data) • Likelihood: r = SUM[log(ls + lb)]- nsig • ls = nsig * psf.per_sr_vs_a(dist(ev.decl, ev.ra, decl, ra))* prob_hit_E2; • lb = bgr.per_sr_vs_decl(ev.decl) * prob_hit_atm; • prob_hit_E2 = hist_E2->GetBinContent(hist_E2->FindBin(ev.nhits)); • prob_hit_atm= hist_atm->GetBinContent(hist_atm->FindBin(ev.nhits));

  9. Adding beta • Likelihood?: