Size: 5274
Comment:
|
Size: 11791
Comment:
|
Deletions are marked like this. | Additions are marked like this. |
Line 4: | Line 4: |
''' [[Projectos/AUGER/MARTA/SIM_B2015|Info about available simulations]]''' | ''' [[Projectos/AUGER/MARTA/SIM_B2015|Info about available simulations|target="_blank]]''' |
Line 8: | Line 8: |
=== List of plots === ==== 3.1 Performance of individual stations ==== "The following information should be provided. As applicable to the item, quantitative values should be based on the analysis of individual stations using the reconstructed datasets at fixed energies (F) : 1019 eV and 1019.8 eV and for p, He, N and Fe primaries. The target for these metrics is the March meeting, with items (iii) and (v) having highest priority." |
=== 3.1 Performance of individual stations === "The following information should be provided. As applicable to the item, quantitative values should be based on the analysis of individual stations using the reconstructed datasets at fixed energies (F) : 1019 eV and 1019.8 eV and for p, He, N and Fe primaries. The target for these metrics is the '''March meeting, with items (iii) and (v) having highest priority'''." |
Line 12: | Line 11: |
I. Strategy for calibrating the new detectors. Expected precision of the relative and absolute calibration to be achieved should be provided. Means for validating the calibration should be discussed. | i. '''Strategy for calibrating the new detectors. Expected precision of the relative and absolute calibration to be achieved should be provided. Means for validating the calibration should be discussed.''' |
Line 14: | Line 13: |
I. The time resolution of the signal trace(s) and response time of the detector. | i. '''The time resolution of the signal trace(s) and response time of the detector.''' |
Line 16: | Line 15: |
I. The resolution and the bias in reconstructing the muonic signal (or electromagnetic signal, as applicable to the upgrade) in individual detector stations. The bias is defined as Delta(Smu)=(Smu,rec - Smu,MC)/Smu,MC and the resolution is given by its variance sigma(Delta(Smu)). These quantities should be shown as a function of the distance, r, of the tank to the shower core for each of the different zenith angles of the simulation sets (F). | i. '''The resolution and the bias in reconstructing the muonic signal (or electromagnetic signal, as applicable to the upgrade) in individual detector stations. The bias is defined as Delta(Smu)=(Smu,rec - Smu,MC)/Smu,MC and the resolution is given by its variance sigma(Delta(Smu)). These quantities should be shown as a function of the distance, r, of the tank to the shower core for each of the different zenith angles of the simulation sets (F).''' |
Line 18: | Line 17: |
I. Lateral distance range and corresponding typical range of signal amplitudes within which the detectors are expected to perform optimally. At a minimum this should include determination of the saturation distance, Rsat, defined to be the radius above which the reconstructed signal has less than a 20% departure from linearity in at least 90% of the events, for showers of 1019.8 eV, for the different angles as given in (F). | ==== All plots for 38 deg ==== |
Line 20: | Line 19: |
I. --(For a detector which is planned to trigger independently of the existing water-Cherenkov detectors (WCD), or is a modification of the existing WCD, the trigger probability as function of the lateral distance of the detector should be given. In particular, above what distance and below what typical detector signal does the trigger probability fall below 95% and 50%?)-- | [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iii_bias_en19.000_th38.000_fiducial.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iii_bias_en19.000_th38.000_fiducial.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''E=10^19^ eV''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iii_bias_en19.800_th38.000_fiducial.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iii_bias_en19.800_th38.000_fiducial.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''E=10^19.8^ eV''' |
Line 22: | Line 22: |
==== 3.2 Event-based observables ==== "The performance of the event reconstruction is of great interest even though such performance measures are only approximate for the time being, because some reconstruction codes are clearly not as mature as others. With that being recognized, the following information should be provided by June. It would be helpful if the different proponents would undertake a first evaluation of items (b), (c), (d) and (e) earlier." |
|
Line 25: | Line 23: |
. '''(a) Average lateral distribution of measured signals and any other pertinent quantities, for different energy, zenith angles and composition (F).''' '''(b) The resolution and bias of the reconstructed Smu (or applicable relevant quantity such as Nm) ''' |
[[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/plots/section.3.1/SimChallenge_v1/resolution_en19.000_th38.000_fiducial.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/plots/section.3.1/SimChallenge_v1/resolution_en19.000_th38.000_fiducial.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''E=10^19^ eV''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/plots/section.3.1/SimChallenge_v1/resolution_en19.800_th38.000_fiducial.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/plots/section.3.1/SimChallenge_v1/resolution_en19.800_th38.000_fiducial.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''E=10^19.8^ eV''' |
Line 28: | Line 26: |
'''of the proposed detector, at an “optimum radius” ropt: ''' | i.#4 '''Lateral distance range and corresponding typical range of signal amplitudes within which the detectors are expected to perform optimally. At a minimum this should include determination of the saturation distance, Rsat, defined to be the radius above which the reconstructed signal has less than a 20% departure from linearity in at least 90% of the events, for showers of 1019.8 eV, for the different angles as given in (F).''' |
Line 30: | Line 28: |
'''The optimum radius1 for a given detection method and detector spacing is the radius giving the smallest variance of DSm, or alternatively the best separation power for different primary masses. As ropt depends on the shower properties, it should be estimated for the fixed energy library (F). In general, ropt will not be 1000 m; in the event that ropt has not yet been determined for the given upgrade proposal, the resolution and bias can be reported taking ropt =1000m or other specified distance.''' | ==== Plots for 10^{19.8}, 38 deg, Epos ==== |
Line 32: | Line 30: |
{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/plots/SIgma_FixedQGSJet.gif|caption|width=200, align="bottom"}}Gráfico 1 [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/plots/SIgma_FixedQGSJet.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/plots/SIgma_FixedQGSJet.gif|caption|width=200, align="bottom", target="_blank"}}|target="_blank"]]Gráfico 1 | [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iv_a_nonlinearity_Proton_Epos_en19.800_th38.000_fiducial.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iv_a_nonlinearity_Proton_Epos_en19.800_th38.000_fiducial.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]]'''Proton''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iv_a_nonlinearity_Iron_Epos_en19.800_th38.000_fiducial.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iv_a_nonlinearity_Iron_Epos_en19.800_th38.000_fiducial.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]]'''Iron''' ==== Plots for p,He,N,Fe, 38 deg., QGSjet, Epos ==== [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iv_b_nonlinearity_en19.000_th38.000_fiducial.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iv_b_nonlinearity_en19.000_th38.000_fiducial.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]]'''E=10^19^eV''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iv_b_nonlinearity_en19.800_th38.000_fiducial.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.1/plots/SimChallenge_3_1_iv_b_nonlinearity_en19.800_th38.000_fiducial.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]]'''E=10^19.8^''' i. --(For a detector which is planned to trigger independently of the existing water-Cherenkov detectors (WCD), or is a modification of the existing WCD, the trigger probability as function of the lateral distance of the detector should be given. In particular, above what distance and below what typical detector signal does the trigger probability fall below 95% and 50%?)-- === 3.2 Event-based observables === "The performance of the event reconstruction is of great interest even though such performance measures are only approximate for the time being, because some reconstruction codes are clearly not as mature as others. With that being recognized, the following information '''should be provided by June'''. It would be helpful if the different proponents would '''undertake a first evaluation of items (b), (c), (d) and (e) earlier.'''" a. ''' Average lateral distribution of measured signals and any other pertinent quantities, for different energy, zenith angles and composition (F)''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/LDF_Mix_Composition.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/LDF_Mix_Composition.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/MeanLDF_Fit.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/MeanLDF_Fit.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/0Final_Xmax_beta_LDFm_CH.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/0Final_Xmax_beta_LDFm_CH.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/0Final_MuonRaul_beta_LDFm_CH.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/0Final_MuonRaul_beta_LDFm_CH.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] a. '''The resolution and bias of the reconstructed Smu (or applicable relevant quantity such as Nmu) of the proposed detector, at an “optimum radius” ropt (see [[attachment:DatachallengeV7.pdf|B2015 Simulation Challenge document|target="_blank"]] for more details). As ropt depends on the shower properties, it should be estimated for the fixed energy library (F). In general, ropt will not be 1000 m; in the event that ropt has not yet been determined for the given upgrade proposal, the resolution and bias can be reported taking ropt =1000m or other specified distance.''' ==== Bias; All plots for 38 deg. ==== [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedQGSjetII04BiasDrawing.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedQGSjetII04BiasDrawing.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''QGSJet04''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedEposBiasDrawing.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedEposBiasDrawing.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''Epos''' ==== Resolution; All plots for 38 deg. ==== [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedQGSjetII04ResolutionsDrawing.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedQGSjetII04ResolutionsDrawing.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''QGSJet04''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedEposResolutionsDrawing.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedEposResolutionsDrawing.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''Epos''' a.#3 ''' The resolution and bias of the reconstructed muon production depth distribution (MPD), for different energies, zenith angles and primary composition (F). The reference is the MPD one would obtain from the muons for a given detection energy threshold, arriving at ground as predicted by Monte Carlo simulation.''' a. ''' Scatter plots for each library in (F) to illustrate the correlation and separation power of the different composition sensitive parameters that can be derived with the upgraded array (i.e. Sm, Nm or Sem vs. S(1000), with each shower being one point in the plot).''' a. ''' Merit factor for separating primary mass groups, to be derived from simulation sets (F). The primary combinations p-Fe, N-Fe and p-He should be specifically considered. See [[attachment:DatachallengeV7.pdf|B2015 Simulation Challenge document|target="_blank"]] for more details. In addition it will be important to consider also the limited energy resolution by deriving the same merit factors including a realistic energy reconstruction using simulation sets (C).''' ==== Distributions for rho^MARTA_1000 @ 10^19^ and 10^19.8^, 38^o^ ==== [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedQGSjetII04DistributionDrawing.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedQGSjetII04DistributionDrawing.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''QGSJet04''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedEposDistributionDrawing.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedEposDistributionDrawing.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''Epos''' |
Line 36: | Line 77: |
==== Merit factor using rho^MARTA^_1000_ alone. ==== [[Projectos/AUGER/MARTA/B2015SimulationChallenge/AdditionalPlots|click here for the same plots including the ideal case|target="_blank]] [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedQGSjetII04MeritDrawing.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedQGSjetII04MeritDrawing.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''QGSJet04''' [[http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedEposMeritDrawing.gif|{{http://www.lip.pt/~auger/MARTA/B2015SimulationChallenge/section.3.2/plots/PlotsFixedEposMeritDrawing.gif|caption|width=400, align="bottom", target="_blank"}}|target="_blank"]] '''Epos''' |
|
Line 38: | Line 84: |
. '''(c) The resolution and bias of the reconstructed muon production depth distribution (MPD), for ''' | |
Line 40: | Line 85: |
'''different energies, zenith angles and primary composition (F). The reference is the MPD one would obtain from the muons for a given detection energy threshold, arriving at ground as predicted by Monte Carlo simulation.''' . '''(d) Scatter plots for each library in (F) to illustrate the correlation and separation power of the ''' '''different composition sensitive parameters that can be derived with the upgraded array (i.e. Sm, Nm or Sem vs. S(1000), with each shower being one point in the plot).''' . '''(e) Merit factor, fI,J, for separating primary mass groups, to be derived from simulation sets (F). ''' '''The primary combinations p-Fe, N-Fe and p-He should be specifically considered. This might be, e.g., fp,Fe = jhSFe m i U00100000 hSp mij q s2Fe m + s2pm , or a more general multi-variate discrimination method can be used. In addition it will be important to consider also the limited energy resolution by deriving the same merit factors including a realistic energy reconstruction using simulation sets (C).''' . '''(f) The event reconstruction efficiency should be evaluated (F,C). Are there event geometries or ''' '''other conditions that lead to an improved or much poorer than average event reconstruction?''' |
a.#6 ''' The event reconstruction efficiency should be evaluated (F,C). Are there event geometries or other conditions that lead to an improved or much poorer than average event reconstruction? ''' |
Beyond 2015 Simulation Challenge
B2015 Simulation Challenge document (V7)
Info about available simulations
3.1 Performance of individual stations
"The following information should be provided. As applicable to the item, quantitative values should be based on the analysis of individual stations using the reconstructed datasets at fixed energies (F) : 1019 eV and 1019.8 eV and for p, He, N and Fe primaries. The target for these metrics is the March meeting, with items (iii) and (v) having highest priority."
Strategy for calibrating the new detectors. Expected precision of the relative and absolute calibration to be achieved should be provided. Means for validating the calibration should be discussed.
The time resolution of the signal trace(s) and response time of the detector.
The resolution and the bias in reconstructing the muonic signal (or electromagnetic signal, as applicable to the upgrade) in individual detector stations. The bias is defined as Delta(Smu)=(Smu,rec - Smu,MC)/Smu,MC and the resolution is given by its variance sigma(Delta(Smu)). These quantities should be shown as a function of the distance, r, of the tank to the shower core for each of the different zenith angles of the simulation sets (F).
All plots for 38 deg
Lateral distance range and corresponding typical range of signal amplitudes within which the detectors are expected to perform optimally. At a minimum this should include determination of the saturation distance, Rsat, defined to be the radius above which the reconstructed signal has less than a 20% departure from linearity in at least 90% of the events, for showers of 1019.8 eV, for the different angles as given in (F).
Plots for 10^{19.8}, 38 deg, Epos
Plots for p,He,N,Fe, 38 deg., QGSjet, Epos
For a detector which is planned to trigger independently of the existing water-Cherenkov detectors (WCD), or is a modification of the existing WCD, the trigger probability as function of the lateral distance of the detector should be given. In particular, above what distance and below what typical detector signal does the trigger probability fall below 95% and 50%?
3.2 Event-based observables
"The performance of the event reconstruction is of great interest even though such performance measures are only approximate for the time being, because some reconstruction codes are clearly not as mature as others. With that being recognized, the following information should be provided by June. It would be helpful if the different proponents would undertake a first evaluation of items (b), (c), (d) and (e) earlier."
Average lateral distribution of measured signals and any other pertinent quantities, for different energy, zenith angles and composition (F)
The resolution and bias of the reconstructed Smu (or applicable relevant quantity such as Nmu) of the proposed detector, at an “optimum radius” ropt (see B2015 Simulation Challenge document for more details). As ropt depends on the shower properties, it should be estimated for the fixed energy library (F). In general, ropt will not be 1000 m; in the event that ropt has not yet been determined for the given upgrade proposal, the resolution and bias can be reported taking ropt =1000m or other specified distance.
Bias; All plots for 38 deg.
Resolution; All plots for 38 deg.
The resolution and bias of the reconstructed muon production depth distribution (MPD), for different energies, zenith angles and primary composition (F). The reference is the MPD one would obtain from the muons for a given detection energy threshold, arriving at ground as predicted by Monte Carlo simulation.
Scatter plots for each library in (F) to illustrate the correlation and separation power of the different composition sensitive parameters that can be derived with the upgraded array (i.e. Sm, Nm or Sem vs. S(1000), with each shower being one point in the plot).
Merit factor for separating primary mass groups, to be derived from simulation sets (F). The primary combinations p-Fe, N-Fe and p-He should be specifically considered. See B2015 Simulation Challenge document for more details. In addition it will be important to consider also the limited energy resolution by deriving the same merit factors including a realistic energy reconstruction using simulation sets (C).
Distributions for rho^MARTA_1000 @ 10^19^ and 10^19.8^, 38^o^
Merit factor using rho^MARTA^_1000_ alone.
click here for the same plots including the ideal case
The event reconstruction efficiency should be evaluated (F,C). Are there event geometries or other conditions that lead to an improved or much poorer than average event reconstruction?