Changes between Version 247 and Version 248 of DatabaseBasedAnalysis/Spectrum


Ignore:
Timestamp:
12/10/19 14:20:22 (5 years ago)
Author:
tbretz
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DatabaseBasedAnalysis/Spectrum

    v247 v248  
    22
    33= Spectrum Analysis =
     4
     5== Theory ==
    46
    57The differential flux \(\phi(E)\) per area, time and energy interval is defined as
     
    108110\[\rightarrow\quad\sigma^2(\phi) = \phi^2 \cdot\left[\left(\frac{\sigma(M_\textrm{exc})}{M_\textrm{exc} }\right)^2 + \left(\frac{\sigma(N'_\textrm{exc})}{N'_\textrm{exc}}\right)^2 + \left(\frac{\sigma(N'_\textrm{src})}{N'_\textrm{src}}\right)^2\right]\]
    109111
    110 
    111 == Define Binnings ==
    112 
    113 == Get Data File List ==
     112== Conceptual Example ==
     113
     114=== Define Binnings ===
     115
     116=== Get Data File List ===
    114117
    115118A list with file IDs containing the events to be analyed is required a.t.m. The following query retrieves such a list and fills a temporary table (`DataFiles`) with the IDs.
     
    146149}}}
    147150
    148 == Get Observation Time ==
     151=== Get Observation Time ===
    149152
    150153The following query bins the effective observation time of the runs listed above in zenith angle bins and stores the result in a temporary table (`ObservationTime`). Note that the result contains only those bins which have entries.
     
    177180`%0:bins` is a placeholder for the bin boundaries, e.g. `5, 10, 15, 20, 25, 30` (five bins between 5° and 30° plus underflow and overflow).
    178181
    179 == Get Monte Carlo File List ==
     182=== Get Monte Carlo File List ===
    180183
    181184The next query obtains all Monte Carlo runs which have their !ThetaMin or !ThetaMax within one of the bins obtained in the previous query (so all MC runs that correspond to bins in which data is available). Strictly speaking, this step is not necessary, but it accelerats further processing. In addition (here as an example) only runs with even FileIDs are obtained as test-runs (assuming that odd runs were used for training). The resulting FileIDs are stored in a temporary table (`MonteCarloFiles`).
     
    210213
    211214
    212 == Get Zenith Angle Histogram ==
     215=== Get Zenith Angle Histogram ===
    213216
    214217The following table creates a temporaray table (`EventCount`) internally which bins the !MonetCarlo files from the file list in !MonteCarloFiles in zenith angle bins. This temporary table ois then joined with table containing the binning for the data files (!EventCount) and for each bin, the ratio (!ZdWeight) and the corresponding error (!ErrZdWeight) is calculated (assuming an error on the on-time of 1s per 5min).
     
    262265
    263266
    264 == Analyze Data ==
    265 
    266 == Analyze Monte Carlo Data ==
    267 
    268 == Summarize Corsika Production ==
     267=== Analyze Data ===
     268
     269=== Analyze Monte Carlo Data ===
     270
     271=== Summarize Corsika Production ===
    269272
    270273{{{#!Spoiler
     
    300303
    301304
    302 == Result (Spectrum) ==
    303 
    304 == Result (Threshold) ==
    305 
    306 == Result (Migration) ==
     305=== Result (Spectrum) ===
     306
     307=== Result (Threshold) ===
     308
     309=== Result (Migration) ===