OscProbCalcer_NuSQUIDSLinear.h

Go to the documentation of this file.

#ifndef __OSCILLATOR_NUSQUIDSLINEAR_H__
#define __OSCILLATOR_NUSQUIDSLINEAR_H__
 
#include "OscProbCalcerBase.h"
 
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wdeprecated-declarations"
#include "nuSQuIDS/nuSQuIDS.h"
#pragma GCC diagnostic pop
#include "examples/Decoherence/nuSQUIDSDecoh.h"
#include "examples/NSI/NSI.h"
#include "examples/LV/LV.h"
 
class OscProbCalcerNuSQUIDSLinear : public OscProbCalcerBase {
 public:
 
  OscProbCalcerNuSQUIDSLinear(YAML::Node Config_);
 
  OscProbCalcerNuSQUIDSLinear(std::string ConfigName_) : OscProbCalcerNuSQUIDSLinear(YAML::LoadFile(ConfigName_)) {}
  
  virtual ~OscProbCalcerNuSQUIDSLinear();
 
  // ========================================================================================================================================================================
  // Functions which need implementation specific code
 
  void SetupPropagator() final;
  
  void CalculateProbabilities() final;
 
  int ReturnWeightArrayIndex(int NuTypeIndex, int OscNuIndex, int EnergyIndex, int CosineZIndex=-1) final;
  
  long DefineWeightArraySize() final;
 
  void SetOscParams(nusquids::nuSQUIDS* base);
 
  void CalculateNuSQUIDS(nusquids::nuSQUIDS* base, nusquids::marray<double,2>& inistate, int& index_counter);
 
  // ========================================================================================================================================================================
  // Functions which help setup implementation specific code
 
  // ========================================================================================================================================================================
  // Variables which are needed for implementation specific code
 
  enum OscParams_PMNS{kTH12, kTH23, kTH13, kDM12, kDM23, kDCP, kPATHL, kDENS, kELECDENS, kNOscParams_PMNS};
 
  enum OscParams_Decoh{kEnergyStrength=kNOscParams_PMNS, kEnergyDep, kEnergyScale, kNOscParams_Decoh};
  
  enum OscParams_LIV{kEMuReal=kNOscParams_PMNS, kEMuImg, kMuTauReal, kMuTauImg, kEnergyPower, kNOscParams_LIV};
 
  int PMNS_StrToInt(const std::string& OscModel);
  
  void SetOscParams(int OscType);
 
  enum OscModels{kSM=0, kDecoherence=1, kLIV=2, kNSI=3};
  
  int fOscModel;
 
  enum NuType{Nu=1,Nubar=-1};
 
  squids::Const units;
  
  nusquids::nuSQUIDS* nus_base;
 
  nusquids::nuSQUIDS* nubars_base;
 
  std::string decoherence_model;
 
  nusquids::nuSQUIDSDecoh::DecoherenceModel nusquids_decoherence_model;
 
  FLOAT_T nsi_mutau_coupling;
 
  double integration_step;
 
  double rel_error;
 
  double abs_error;
 
  int nNeutrinoFlavours;
};
 
#endif