Changes between Version 232 and Version 233 of DatabaseBasedAnalysis/Spectrum


Ignore:
Timestamp:
12/09/19 18:11:08 (5 years ago)
Author:
tbretz
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DatabaseBasedAnalysis/Spectrum

    v232 v233  
    2020
    2121Note that the exact calculation of the efficiency \(\epsilon(\Delta E)\) depends on prior knowledge of the correct source spectrum \(N_0(E)\). Therefore, it is strictly speaking only correct if the simulated spectrum and the real spectrum are identical. As the real spectrum is unknown, special care has to be taken of the systematic introduced by the assumption of \(N_0(E)\).
    22 
    23 === Zenith Angle Weights and Spectral Weights ===
    24 
    25 While the source spectrum \(N_0\) is of course independent of the zenith angle \(\theta\) of the observation, the efficiency is not, so that
    26 
    27 \[\epsilon(E) = \epsilon(E, \theta)\]
    28 
    29 In addition, the observed differential flux depends on the zenith angle. That means that for a correct calculation of the efficiency, the number of simulated events per zenith angle interval has to match the observation time distribution.
    30 
    31 \[\epsilon(E, \theta) = \epsilon(E, \theta, \Delta T(\theta))\]
    32 
    33 This can be achieved by applying zenith angle dependent weights \(\omega(\theta)\) to each simulated event. Similarly, the discussed requirement of the match between result spectrum and simulated spectrum can be achieved applying energy dependent weights \(\omega(E)\).
    34 
    35 As the simulated energy spectrum is independent of zenith angle, it can be expressed as
    36 
    37 \[\phi(E,\theta) = N_0\cdot \eta(E)\cdot \eta(\theta)\]
    38 
    39 with the differential energy spectrum \(\eta(E)\) and the zenith angle distribution \(\eta(\theta)\) and the normalization
    40 
    41 \[\int_E\int_\theta\phi(E,\theta)dEd\theta = N_0\cdot \int_E\eta(E)dE\cdot \int_\theta\eta(\theta)d\theta\]
    42 
    43 The contents of one energy bin \(\Delta E\) and one zenith angle bin \(\Delta\theta\) can then be written as
    44 
    45 \[N(\Delta E, \Delta\theta) = \sum_{\Delta E}\sum_{\Delta\theta} \omega_{i,j}(E, \theta)= \sum_{\Delta E}\sum_{\Delta\theta} \omega_i(E)\cdot\omega_j(\theta)= \sum_{\Delta E} \omega_i(\theta)\cdot\sum_{\Delta\theta}\omega_i(E) \]
    46 
    47 and the normalization
    48 
    49 \[\sum_E\sum_\theta N(E,\theta) = N_0 \sum_E\sum_\theta \omega_{i,j}(E,\theta) = N_0 \sum_E\sum_\theta \omega_{i}(E)\omega_j(\theta) = N_0\cdot \sum_E\omega_i(\theta)\cdot \sum_\theta\omega_i(E)\]
    50 
    51 
    52 
    5322
    5423=== Excess ===