Esempio n. 1
0
    def test_fix_param(self):
        m = Minuit(func4, printlevel=-1)
        m.migrad()
        m.minos()
        val = m.values
        self.assertAlmostEqual(val['x'], 2.)
        self.assertAlmostEqual(val['y'], 5.)
        self.assertAlmostEqual(val['z'], 7.)
        err = m.errors  #second derivative
        # self.assertAlmostEqual(err['x'],5.)
        # self.assertAlmostEqual(err['y'],10.)
        # self.assertAlmostEqual(err['z'],4.)

        #now fix z = 10
        m = Minuit(func4, printlevel=-1, y=10., fix_y=True)
        m.migrad()
        val = m.values
        self.assertAlmostEqual(val['x'], 2.)
        self.assertAlmostEqual(val['y'], 10.)
        self.assertAlmostEqual(val['z'], 7.)
        self.assertAlmostEqual(m.fmin(), 10. + 2.5)

        self.assertIn('x', m.free_param)
        self.assertNotIn('x', m.fix_param)
        self.assertIn('y', m.fix_param)
        self.assertNotIn('y', m.free_param)
        self.assertNotIn('z', m.fix_param)
Esempio n. 2
0
 def functesthelper(self, f):
     m = Minuit(f, printlevel=-1)
     m.migrad()
     val = m.values
     self.assertAlmostEqual(val['x'], 2.)
     self.assertAlmostEqual(val['y'], 5.)
     self.assertAlmostEqual(m.fmin(), 10.)