sboeser/pisa — pull request #2 by tarlen5 (2014-03-26 13:26:57+00:00)
Title: Initial Stage1 Implementation, with oscProbMaps

Hi Sebastian,

This is my implementation of my stage1 analysis: creating the oscillated flux maps. This is also my first crack at writing the code to generate the oscillation probability lookup tables, found in $PISA/resources/oscProbMaps/. So far, until we handle the standards for naming these oscillation probability lookup tables and where to store it/how to access it, I am simply working with 2 oscillation flux map files in this development mode, saving them to a local directory, as you can see in OscillatedFluxMC.py->get_osc_prob_maps().

Implemented so far is the "Smoothing" (or "downsampling" as you have called it) in a very general way, so that I can use logarithmic bins to create the lookup tables and then any linear or log binning to do the actual analysis. And so far I've been able to have enough accuracy with 500x500 bins in the lookup table from energy: 1-80 GeV and coszen from -1 to 1.

Thanks,
Tim