Changes between Version 247 and Version 248 of DatabaseBasedAnalysis/Spectrum
- Timestamp:
- 12/10/19 14:20:22 (5 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
DatabaseBasedAnalysis/Spectrum
v247 v248 2 2 3 3 = Spectrum Analysis = 4 5 == Theory == 4 6 5 7 The differential flux \(\phi(E)\) per area, time and energy interval is defined as … … 108 110 \[\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]\] 109 111 110 111 == Define Binnings == 112 113 == Get Data File List == 112 == Conceptual Example == 113 114 === Define Binnings === 115 116 === Get Data File List === 114 117 115 118 A 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. … … 146 149 }}} 147 150 148 == Get Observation Time==151 === Get Observation Time === 149 152 150 153 The 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. … … 177 180 `%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). 178 181 179 == Get Monte Carlo File List==182 === Get Monte Carlo File List === 180 183 181 184 The 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`). … … 210 213 211 214 212 == Get Zenith Angle Histogram==215 === Get Zenith Angle Histogram === 213 216 214 217 The 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). … … 262 265 263 266 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 === 269 272 270 273 {{{#!Spoiler … … 300 303 301 304 302 == Result (Spectrum)==303 304 == Result (Threshold)==305 306 == Result (Migration)==305 === Result (Spectrum) === 306 307 === Result (Threshold) === 308 309 === Result (Migration) ===