OscProbCalcer_OscProb.h

Go to the documentation of this file.

#ifndef __OSCILLATOR_OSCPROB_H__
#define __OSCILLATOR_OSCPROB_H__
 
#include "OscProbCalcerBase.h"
 
#include "inc/PMNS_Base.h"
#include "inc/PremModel.h"
 
class OscProbCalcerOscProb : public OscProbCalcerBase {
 public:
 
  OscProbCalcerOscProb(YAML::Node Config_);
 
  OscProbCalcerOscProb(std::string ConfigName_) : OscProbCalcerOscProb(YAML::LoadFile(ConfigName_)) {}
 
  virtual ~OscProbCalcerOscProb();
 
  // ==========================================================================
  // Functions which need implementation specific code
 
  void SetupPropagator() final;
 
  void CalculateProbabilities() final;
 
  void CalcProbPMNS();
 
  void SetPath(int iCosineZ);
 
  int ReturnWeightArrayIndex(int NuTypeIndex, int OscNuIndex, int EnergyIndex, int CosineZIndex=-1) final;
 
  long DefineWeightArraySize() final;
 
  // ==========================================================================
  // Functions which help setup implementation specific code
 
  int PMNS_StrToInt(const std::string& PMNSType);
 
  void SetOscParams();
 
  OscProb::PMNS_Base* GetPMNSObj();
 
  void SetPMNSParams();
 
  void SetSterileParams();
  void SetDecayParams();
  void SetDecoParams();
  void SetNSIParams();
  void SetIterParams();
  void SetNUNMParams();
  void SetLIVParams();
  void SetSiderealLIVParams();
  void SetOQSParams();
 
  int GetNCosineZ();
 
  // ==========================================================================
  // Variables which are needed for implementation specific code
 
  enum OscParams{kTH12, kTH23, kTH13, kDM12, kDM23, kDCP, kNOscParams};
 
  enum OscParams_Sterile1{kTH14 = kNOscParams, kTH24, kTH34, kDM14,
                          kDelta14, kDelta24, kNOscParams_Sterile1};
  enum OscParams_Sterile2{kTH15 = kNOscParams_Sterile1, kTH25, kTH35, kTH45, kDM15,
                          kDelta15, kDelta25, kDelta35, kNOscParams_Sterile2};
  enum OscParams_Sterile3{kTH16 = kNOscParams_Sterile2, kTH26, kTH36, kTH46, kTH56, kDM16,
                          kDelta16, kDelta26, kDelta36, kDelta46, kNOscParams_Sterile3};
 
  enum OscParams_Decay{kAlpha2 = kNOscParams, kAlpha3, kNOscParams_Decay};
 
  enum OscParams_Deco{kGamma21 = kNOscParams, kGamma31, kDecoAngle, kPower, kNOscParams_Deco};
 
  enum OscParams_NSI{kEps_ee = kNOscParams, kEps_emu, kEps_etau, kEps_mumu,
                     kEps_mutau, kEps_tautau, kDelta_emu, kDelta_etau,
                     kDelta_mutau, kElecCoup, kUpCoup, kDownCoup,
                     kNOscParams_NSI};
 
  enum OscParams_SNSI{kLightMass = kNOscParams_NSI, kNOscParams_SNSI};
 
  enum OscParams_Iter{kPrec = kNOscParams, kNOscParams_Iter};
 
  enum OscParams_NUNM{kAlpha11 = kNOscParams, kAlpha21, kAlpha31, kAlpha22, kAlpha32, kAlpha33,
                      kPhi21, kPhi31, kPhi32, kFracVnc, kNOscParams_NUNM};
 
  enum OscParams_LIV{kaT_ee_3 = kNOscParams, kaT_emu_3, kaT_etau_3, kaT_mumu_3, kaT_mutau_3,
                     kaT_tautau_3, kDelta_emu_3, kDelta_etau_3, kDelta_mutau_3,
                     kcT_ee_4, kcT_emu_4, kcT_etau_4, kcT_mumu_4, kcT_mutau_4,
                     kcT_tautau_4, kDelta_emu_4, kDelta_etau_4, kDelta_mutau_4,
                     kaT_ee_5, kaT_emu_5, kaT_etau_5, kaT_mumu_5, kaT_mutau_5,
                     kaT_tautau_5, kDelta_emu_5, kDelta_etau_5, kDelta_mutau_5,
                     kcT_ee_6, kcT_emu_6, kcT_etau_6, kcT_mumu_6, kcT_mutau_6,
                     kcT_tautau_6, kDelta_emu_6, kDelta_etau_6, kDelta_mutau_6,
                     kaT_ee_7, kaT_emu_7, kaT_etau_7, kaT_mumu_7, kaT_mutau_7,
                     kaT_tautau_7, kDelta_emu_7, kDelta_etau_7, kDelta_mutau_7,
                     kcT_ee_8, kcT_emu_8, kcT_etau_8, kcT_mumu_8, kcT_mutau_8,
                     kcT_tautau_8, kDelta_emu_8, kDelta_etau_8, kDelta_mutau_8,
                     kNOscParams_LIV};
 
 
  enum OscParams_SiderealLIV {
    // A coefficients: X,Y,Z components
    aT_ee_X = kNOscParams, aT_ee_Y, aT_ee_Z,
    aT_emu_X, aT_emu_Y, aT_emu_Z,
    aT_etau_X, aT_etau_Y, aT_etau_Z,
 
    aT_mumu_X, aT_mumu_Y, aT_mumu_Z,
    aT_mutau_X, aT_mutau_Y, aT_mutau_Z,
    aT_tautau_X, aT_tautau_Y, aT_tautau_Z,
 
 
    // C coefficients: symmetric spatial tensor XX,XY,XZ,YY,YZ,ZZ
    cT_ee_XX, cT_ee_XY, cT_ee_XZ, cT_ee_YY, cT_ee_YZ, cT_ee_ZZ,
    cT_emu_XX, cT_emu_XY, cT_emu_XZ, cT_emu_YY, cT_emu_YZ, cT_emu_ZZ,
    cT_etau_XX, cT_etau_XY, cT_etau_XZ, cT_etau_YY, cT_etau_YZ, cT_etau_ZZ,
    cT_mumu_XX, cT_mumu_XY, cT_mumu_XZ, cT_mumu_YY, cT_mumu_YZ, cT_mumu_ZZ,
    cT_mutau_XX, cT_mutau_XY, cT_mutau_XZ, cT_mutau_YY, cT_mutau_YZ, cT_mutau_ZZ,
    cT_tautau_XX, cT_tautau_XY, cT_tautau_XZ, cT_tautau_YY, cT_tautau_YZ, cT_tautau_ZZ,
 
 
    // Detector / sidereal parameters
    kLatitude, kLongitude, kAltitude,
    kZenith, kAzimuth,
    kSiderealTime,
    kOscParams_SiderealLIV
  };
 
  enum OscParams_OQS {
    // Decoherence magnitudes (Gell-Mann basis)
    kA1 = kNOscParams, kA2, kA3, kA4, kA5, kA6, kA7, kA8,
    // Angles
    kTheta12, kTheta13, kTheta14, kTheta15, kTheta16, kTheta17, kTheta18,
    kTheta23, kTheta24, kTheta25, kTheta26, kTheta27, kTheta28,
    kTheta34, kTheta35, kTheta36, kTheta37, kTheta38,
    kTheta45, kTheta46, kTheta47, kTheta48,
    kTheta56, kTheta57, kTheta58,
    kTheta67, kTheta68,
    kTheta78,
    kPower_OQS,
    kNOscParams_OQS
  };
 
  enum NuType{Nu=1,Nubar=-1};
 
  enum PMNSMatrix{kFast, kPMNSSterile1, kPMNSSterile2, kPMNSSterile3,
                  kDecay, kDeco, kNSI, kSNSI, kIter, kNUNM, kLIV, kOQS, kSiderealLIV};
 
  int fOscType;
 
private:
 
  OscProb::PremModel fPremModel;
 
  std::string fPremFile;
 
  double fDetDepth;
 
  int fMaxGenFlavour;
 
  int fMaxDetFlavour;
 
  OscProb::PMNS_Base* fPMNSObj;
 
};
 
#endif