Ejemplo n.º 1
0
def style(nu):
    """Return the line style to use for a given neutrino flavour."""
    if oscillations.isNeutrino(nu):
        return 1
    elif oscillations.isAntiNeutrino(nu):
        return 7
    else:
        raise ValueError("Invalid neutrino/anti-neutrino enum.")
Ejemplo n.º 2
0
def style(nu):
	"""Return the line style to use for a given neutrino flavour."""
	if oscillations.isNeutrino(nu):
		return 1
	elif oscillations.isAntiNeutrino(nu):
		return 7
	else:
		raise ValueError("Invalid neutrino/anti-neutrino enum.")
Ejemplo n.º 3
0
	def setUp(self):
		# Example experiment oscillations
		self.osc = oscillations.Oscillations()
		self.osc.setL(295.0 * oscillations.units.km )
		self.osc.setE(0.6   * oscillations.units.GeV)
		self.osc.setDeltaM21( 7.50e-5 * oscillations.units.eV2 )
		self.osc.setDeltaM32( 2.32e-3 * oscillations.units.eV2 )
		self.osc.setTheta12( 33.9 * oscillations.units.degrees )
		self.osc.setTheta13(  9.1 * oscillations.units.degrees )
		self.osc.setTheta23( 45.0 * oscillations.units.degrees )
		self.osc.setDeltaCP(  0.0 * oscillations.units.degrees )
		# get an index which does not represent a neutrino or anti-neutrino
		non_nu = 1
		while( oscillations.isNeutrino(non_nu) or oscillations.isAntiNeutrino(non_nu) ):
			non_nu += 1
		self.non_neutrino = non_nu
Ejemplo n.º 4
0
 def setUp(self):
     # Example experiment oscillations
     self.osc = oscillations.Oscillations()
     self.osc.setL(295.0 * oscillations.units.km)
     self.osc.setE(0.6 * oscillations.units.GeV)
     self.osc.setDeltaM21(7.50e-5 * oscillations.units.eV2)
     self.osc.setDeltaM32(2.32e-3 * oscillations.units.eV2)
     self.osc.setTheta12(33.9 * oscillations.units.degrees)
     self.osc.setTheta13(9.1 * oscillations.units.degrees)
     self.osc.setTheta23(45.0 * oscillations.units.degrees)
     self.osc.setDeltaCP(0.0 * oscillations.units.degrees)
     # get an index which does not represent a neutrino or anti-neutrino
     non_nu = 1
     while (oscillations.isNeutrino(non_nu)
            or oscillations.isAntiNeutrino(non_nu)):
         non_nu += 1
     self.non_neutrino = non_nu
Ejemplo n.º 5
0
 def test_nuEBarIsAntiNeutrino(self):
     self.assertTrue(oscillations.isAntiNeutrino(oscillations.nu_e_bar))
Ejemplo n.º 6
0
 def test_nuTauIsNotAntiNeutrino(self):
     self.assertFalse(oscillations.isAntiNeutrino(oscillations.nu_tau))
Ejemplo n.º 7
0
	def test_nuTauBarIsAntiNeutrino(self):
		self.assertTrue( oscillations.isAntiNeutrino( oscillations.nu_tau_bar ) )
Ejemplo n.º 8
0
	def test_nuTauIsNotAntiNeutrino(self):
		self.assertFalse( oscillations.isAntiNeutrino( oscillations.nu_tau ) )