90 likes | 219 Views
IFIC meeting. 01/Feb/2012. Agustín Sánchez Losa. BEST LASER vs OFFICIAL OFFSETS. BEST LASER vs OFFICIAL OFFSETS. BEST LASER vs OFFICIAL OFFSETS. BEST LASER vs OFFICIAL OFFSETS. Last Fitted Signal. Means. > nsig0: 0.115321 > nsig1: 0.959034 > nsig3: 2.869490 > nsig5: 4.768660
E N D
IFIC meeting 01/Feb/2012 Agustín Sánchez Losa
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;
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;
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));
Adding beta • Likelihood?: