230 likes | 325 Views
PROGRESS ON WATER PROPERTIES ON TRACKS RECONSTRUCTION Harold Yepes-Ramirez 09/11/2011. DATA/MonteCarlo SELECTION: Data 2008 – 2010 data from the official SeaTray production May 2011 (5997 runs). First run: 31051. Last run: 54244. Lifetime: 618.96 days.
E N D
PROGRESS ON WATER PROPERTIES ON TRACKS RECONSTRUCTION Harold Yepes-Ramirez 09/11/2011
DATA/MonteCarlo SELECTION: Data 2008 – 2010 data from the official SeaTray production May 2011 (5997 runs). First run: 31051. Last run: 54244. Lifetime: 618.96 days. MonteCarlo SoS prepared (C. Bogazzi) with the previous runs (5997 data runs). Mupage for muons + Geasim for neutrinos. MonteCarlo statistics (up to now): sc = scattering centers; aa = om angular acceptance; abs = absorption; sca = scattering; eta = fraction of Rayleigh scattering.
SANITY CHECKS WITH PREVIOUS PRODUCTIONS: • Harold: • It is not a run-by-run simulation. • 5997 data runs (2008-2010). • 312 mupage muon runs. • 90 neutrino + 40 anti-neutrino Geasim files. • TE May 2011. • Down-going neutrinos not used in this MC. • Juan Pablo: • Run-by-run simulation. • 5997 data runs (2008-2010). • 5941 mupage muon runs. • 5898 neutrino + 5900 anti-neutrino files. • TE September 2010. • Down-going neutrinos are used in the r-b-r MC. The L cut at -5.2 remove some muons-neutrinos close to the horizon, not enough affected by JP ones due to large statistics. L > -5.4 relax this zone (SEE NEXT PAGE).
> -5.4 relaxes the zone close to the horizon for neutrinos and muons. • Left My distribution for 2008-2010 data. • Right JP distribution ONLY for 2007-2008 data (CM Amsterdam). • Not directly comparable but both distributions have a nice agreement no matter the selected period .
UPDATED PLOTS FOR THE NUMBER OF HITS USED IN THE FIT: DISTRIBUTIONS FOR A LARGE MC SAMPLES AND OVERVIEW TO THE WATER MODELS
NUMBER OF HITS USED IN THE FIT: • Lesson learnt since the CM in Moscow: the extreme scattering models (lsca<22) shows the worst agreement to data: • Lower values (~0.02) of contribution of Rayleigh scattering (eta) over-estimates the data, and higher (~0.17) values under-estimates it. • The ANTARES site seems to have a large scattering length. • Values of labs higher than 55 m could not be an good approach. • The ANTARES standard model, works reasonably well specially at nhit<80 where first photons are expected. • The uncertainty in the labs for a set of common scattering parameters seems to be less than 20%. • The best agreement to data is then expected for large scattering lengths and not enough higher absorption lengths.
UPDATED PLOTS FOR THE ZENITH ANGLE CONTRIBUTION: UNCERTAINITY ESTIMATION OF RECONSTRUCTED TRACKS DUE WATER OPTICAL PARAMETERS
INFLUENCE OF labs ON RECONSTRUCTED TRACKS: Strategy For a couple of water models with different labs but same scattering parameters, estimate the difference on the reconstructed track rate Uncertainty on labs Vs uncertainty on the muon rate. Previous systematic studies in ANTARES (J.A et al / Astroparticle Physics 34, 2010, 179-184, Pag. 182) “The uncertainty of the light absorption length in water is assumed to be ±10% over the whole wavelength spectrum and yields a variation of ±20% on the number of expected events”. CASE 2 CASE 1 Case 1: Uncertainty on labs: ~13 % (8 m). Uncertainty on muon rate: ~25 % (0.14 Hz). Case 2: Uncertainty on labs: ~13 % (8 m). Uncertainty on muon rate: ~13 % (0.07 Hz).
CASE 3 CASE 4 Case 3: Uncertainty on labs: ~13 % (8 m). Uncertainty on muon rate: ~10 % (0.07 Hz). Case 4: Uncertainty on labs: ~13 % (8 m). Uncertainty on muon rate: ~13 % (0.07 Hz). Case 5: Uncertainty on labs: ~13 % (8 m). Uncertainty on muon rate: ~15 % (0.08 Hz). CASE 5
UNCERTAINTY ON labs (~13%) and UNCERTAINTY ON MUON RATES: FUNCTION OF ZENITH ANGLE Case 1: sm_rate ~25 % Case 2: sm_rate ~13 % Case 3: sm_rate ~14 % Case 4: sm_rate ~14 % Case 5: sm_rate ~14 % • UNCERTAINTY ON labs ±13% • MEAN UNCERTAINTY mrate ± 15%
INFLUENCE OF lsca,eff ON RECONSTRUCTED TRACKS: Strategy Two optical parameters fixed (absorption, eta) and one free parameter (scattering length), for both absorption lengths. CASE 2 CASE 1 Case 1: Uncertainty on lsca_eff: ~23 % (51.5 m). Uncertainty on muon rate: ~16 % (0.08 Hz). Uncertainty on lsca_eff: ~46 % (81.5 m). Uncertainty on muon rate: ~17 % (0.1 Hz). Uncertainty on lsca_eff: ~58 % (133 m). Uncertainty on muon rate: ~31 % (0.18 Hz). Case 1: Uncertainty on lsca_eff: ~23 % (51.5 m). Uncertainty on muon rate: ~4 % (0.02 Hz). Uncertainty on lsca_eff: ~46 % (81.5 m). Uncertainty on muon rate: ~14 % (0.09 Hz). Uncertainty on lsca_eff: ~58 % (133 m). Uncertainty on muon rate: ~17 % (0.11 Hz). • UNCERTAINTY ON lsca_eff [23-58]% UNCERTAINTY mrate [4-17]%
INFLUENCE OF h ON RECONSTRUCTED TRACKS: Strategy One optical parameters fixed (absorption) and two free parameters (scattering length and eta), for both absorption lengths. CASE 2 CASE 1 Case 1: Uncertainty on h: ~35 % (0.06). Uncertainty on muon rate: ~13 % (0.06 Hz). Uncertainty on h: ~81 % (0.09). Uncertainty on muon rate: ~2 % (0.01 Hz). Uncertainty on h: ~88 % (0.15). Uncertainty on muon rate: ~15 % (0.07 Hz). Case 1: Uncertainty on h: ~35 % (0.06). Uncertainty on muon rate: ~2 % (0.01 Hz). Uncertainty on h: ~81 % (0.09). Uncertainty on muon rate: ~2 % (0.01 Hz). Uncertainty on h: ~88 % (0.15). Uncertainty on muon rate: ~4 % (0.02 Hz). • UNCERTAINTY ON h [35-88]% UNCERTAINTY mrate [2-4]%
CONCLUSIONS/COMMENTS: For the water models different to the “standard water model” a bug was founded, a factor scale of ~15 % was not applied, now the correct ones are running. So, for the preliminary comparison such water models have to be considered with a 15% less events. The preliminary “comments in boxes” as a kind of conclusions on the slides, are referred to the expected results. Preliminary and raw results point that, for the current physics conditions simulated in the ANTARES KM3 code (whichever they are), the impact of water parameters could be summarized as follow (s = uncertainty): The effective scattering length seems to be the most relevant parameter. These uncertainties could be more lower if we consider that scattering lengths less than 22 m are not possible in ANTARES site. Data/MC agreement could mainly come from other systematics? Once the corrected data sample is going to be ready, the analysis will be done once again and it will prove the preliminary and raw conclusions. A first look on effective areas and angular resolution will be performed. An internal note will be prepared with a dedicated description of the analysis.