OscillatorSubSampling.h

Go to the documentation of this file.

#ifndef __OSCILLATOR_SUBSAMPLING_BASE_H__
#define __OSCILLATOR_SUBSAMPLING_BASE_H__
 
#include "OscillatorBase.h"
 
class OscillatorSubSampling : public OscillatorBase {
 public:
 
  OscillatorSubSampling(std::string ConfigName_);
 
  OscillatorSubSampling(YAML::Node Config_);
 
  virtual ~OscillatorSubSampling();
 
  // ========================================================================================================================================================================
  // Public functions which are calculation implementation agnostic
 
  const FLOAT_T* ReturnWeightPointer(int InitNuFlav, int FinalNuFlav, FLOAT_T EnergyVal, FLOAT_T CosineZVal=DUMMYVAL) final;
 
  std::vector<FLOAT_T> ReturnBinEdgesForPlotting(bool ReturnEnergy) final;
  
  // ========================================================================================================================================================================
  // Public virtual functions which need calculater specific implementations
 
 protected:
 
  // ========================================================================================================================================================================
  // Protected functions which are calculation implementation agnostic   
 
  // ========================================================================================================================================================================
  // Protected virtual functions which are calculation implementation agnostic
 
  // ========================================================================================================================================================================
  // Basic protected variables required for oscillation probability calculation
 
 private:
 
  void Initialise();
 
  void PostCalculateProbabilities() final;
 
  void SetupOscillatorImplementation() final;
  
  int FindBinIndexFromEdges(FLOAT_T Val, std::vector<FLOAT_T> BinEdges);
 
  std::vector<FLOAT_T> AveragedOscillationProbabilities;         
 
  std::vector< std::vector<const FLOAT_T*> > OscillationProbabilitiesToAverage; 
 
  size_t nCoarseEnergyBins;
 
  size_t nCoarseCosineZBins;
 
  size_t nFineEnergyBins;
 
  size_t nFineCosineZBins;
 
  size_t TotalCoarseBins;
 
  std::vector<NuOscillator::OscillationChannel> OscillationChannels;
 
  size_t nOscillationChannels;
 
  std::vector<int> NeutrinoTypes;
 
  size_t nNeutrinoTypes;
  
  // ========================================================================================================================================================================
  // Basic private variables required for oscillation probability calculation
 
  std::string FileName;
 
  std::string CoarseEnergyAxisHistName;
 
  std::string FineEnergyAxisHistName;
 
  std::string CoarseCosineZAxisHistName;
 
  std::string FineCosineZAxisHistName;
 
  std::vector<FLOAT_T> CoarseEnergyAxisBinEdges;
 
  std::vector<FLOAT_T> FineEnergyAxisBinEdges;
 
  std::vector<FLOAT_T> CoarseCosineZAxisBinEdges;
 
  std::vector<FLOAT_T> FineCosineZAxisBinEdges;
 
  std::vector<FLOAT_T> FineEnergyAxisBinCenters;
  
  std::vector<FLOAT_T> FineCosineZAxisBinCenters;
 
};
 
#endif