osc Namespace Reference

Classes
class	EarthModel
class	IOscCalculator
class	IOscCalculatorAdjustable
class	NoOscillations
class	OscCalculator
class	OscCalculatorGeneral
class	OscCalculatorPMNS
class	OscCalculatorPMNSOpt
class	PMNS
class	PMNS_NSI
class	PMNSOpt

Typedefs
typedef std::complex< double >	val_t
typedef ublas::bounded_array< val_t, kNumFlavours >	alloc_t
typedef ublas::bounded_matrix< val_t, kNumFlavours, kNumFlavours >	ComplexMat
typedef ublas::c_vector< val_t, kNumFlavours >	ComplexVec
typedef ublas::unit_vector< val_t, alloc_t >	UnitVec

Functions
const ublas::zero_matrix< val_t, alloc_t >	kZeroMat (kNumFlavours, kNumFlavours)
const ublas::identity_matrix< val_t, alloc_t >	kIdentity (kNumFlavours)
ComplexMat	GetPMNS (OscCalculatorGeneral::Priv *d)
ComplexMat	VacuumHamiltonian (const ComplexMat &U, std::vector< double > mSq, double E)
ComplexMat	MatterHamiltonianComponent (double Ne, double emutau)
ComplexMat	MatrixExp (const ComplexMat &m2)
ComplexVec	EvolveState (ComplexVec A, const ComplexMat &H, double L)
void	conjugate_elements (ComplexMat &m)

Variables
const unsigned int	kNumFlavours = 3

Typedef Documentation

typedef ublas::bounded_array<val_t, kNumFlavours> osc::alloc_t

Definition at line 22 of file OscCalculatorGeneral.cxx.

typedef ublas::bounded_matrix<val_t, kNumFlavours, kNumFlavours> osc::ComplexMat

Definition at line 24 of file OscCalculatorGeneral.cxx.

typedef ublas::c_vector<val_t, kNumFlavours> osc::ComplexVec

Definition at line 25 of file OscCalculatorGeneral.cxx.

typedef ublas::unit_vector<val_t, alloc_t> osc::UnitVec

Definition at line 26 of file OscCalculatorGeneral.cxx.

typedef std::complex<double> osc::val_t

Definition at line 21 of file OscCalculatorGeneral.cxx.

Function Documentation

void osc::conjugate_elements

(

ComplexMat &

m

)

Definition at line 193 of file OscCalculatorGeneral.cxx.

  {
    for(unsigned int i = 0; i < kNumFlavours; ++i)
      for(unsigned int j = 0; j < kNumFlavours; ++j)
        m(i, j) = std::conj(m(i, j));
  }

ComplexVec osc::EvolveState	(	ComplexVec	A,
		const ComplexMat &	H,
		double	L
	)

Definition at line 186 of file OscCalculatorGeneral.cxx.

  {
    const std::complex<double> i = std::complex<double>(0, 1);

    return ublas::prod(MatrixExp(-H*L*i), A);
  }

ComplexMat osc::GetPMNS

(

OscCalculatorGeneral::Priv *

d

)

Definition at line 106 of file OscCalculatorGeneral.cxx.

  {
    if(d->dirty){
      ublas::noalias(d->pmns) = ublas::prod(d->atmos, ublas::prod<ComplexMat>(d->react, d->solar));
      d->dirty = false;
    }
    return d->pmns;
  }

const ublas::identity_matrix<val_t, alloc_t> osc::kIdentity

(

kNumFlavours

)

const ublas::zero_matrix<val_t, alloc_t> osc::kZeroMat	(	kNumFlavours	,
		kNumFlavours
	)

ComplexMat osc::MatrixExp

(

const ComplexMat &

m2

)

Definition at line 170 of file OscCalculatorGeneral.cxx.

  {
    // We use the identity e^(a*b) = (e^a)^b
    // First: ensure the exponent is really small, 65536 = 2^16
    ComplexMat m = m2/65536;

    // To first order e^x = 1+x
    ComplexMat ret = kIdentity+m;

    // Raise the result to the power 65536, by squaring it 16 times
    for(int n = 0; n < 16; ++n) ret = ublas::prod(ret, ret);

    return ret;
  }

ComplexMat osc::MatterHamiltonianComponent	(	double	Ne,
		double	emutau
	)

Definition at line 144 of file OscCalculatorGeneral.cxx.

  {
    // Need to convert avogadro's constant so that the total term comes out in
    // units of inverse distance. Note that Ne will be specified in g/cm^-3
    // I put the following into Wolfram Alpha:
    // (fermi coupling constant)*((avogadro number)/cm^3)*(reduced planck constant)^2*(speed of light)^2
    // And then multiplied by 1000 because we specify L in km not m.
    const double GF = 1.368e-4;

    ComplexMat H = kZeroMat;

    H(0, 0) = sqrt(2)*GF*Ne;

    // Ignoring conjugates here because we assume e_mutau is real
    H(1, 2) = H(2, 1) = emutau*sqrt(2)*GF*Ne;

    return H;
  }

ComplexMat osc::VacuumHamiltonian	(	const ComplexMat &	U,
		std::vector< double >	mSq,
		double	E
	)

Definition at line 119 of file OscCalculatorGeneral.cxx.

  {
    // Conversion factor for units
    E /= 4*1.267;

    assert(mSq.size() == kNumFlavours);

    ComplexMat H = kZeroMat;

    // Loop over rows and columns of the output
    for(unsigned int a = 0; a < kNumFlavours; ++a){
      for(unsigned int b = 0; b < kNumFlavours; ++b){
        // Form sum over mass states
        for(unsigned int i = 0; i < kNumFlavours; ++i){
          H(a, b) += U(a, i)*mSq[i]/(2*E)*std::conj(U(b, i));
        }
      }
    }

    return H;
  }

Variable Documentation

const unsigned int osc::kNumFlavours = 3

Definition at line 18 of file OscCalculatorGeneral.cxx.