Skip to main content

code to do a minimised chi squared fit on a function by varying two variables in python for a function with a maximum point not working correctly

I have written a piece of code intended to perform a minimised chi squared fit by varying two parameters, the problem section of code is below:

conversion_factor = 0.3894 * 10**9

def cross_section(parameters, E, Γee):

    mZ, ΓZ = parameters

    return conversion_factor * (12 * np.pi / (mZ ** 2)) * (
                (E ** 2) / ((E ** 2 - mZ ** 2) ** 2) + (mZ ** 2) * (
                    ΓZ ** 2)) * (Γee ** 2)

def chi_squared(parameters, E, sigma, sigma_error, mZ_start, ΓZ_start, Γee):

    mZ, ΓZ = parameters

    model = cross_section(parameters, E,Γee)

    return np.sum(((model - sigma) / sigma_error) ** 2)

E = combined_data['% centre-of-mass energy (GeV)']
sigma = combined_data['cross section (nb)']
sigma_error = combined_data['uncertainty (nb)']

mZ_start, ΓZ_start = [90, 3]

Γee = 83.91

result = fmin(chi_squared, (mZ_start, ΓZ_start), 
              args=(E, sigma, sigma_error, mZ_start, ΓZ_start, Γee),
              full_output=True, disp=False)

print(result[0][0])
print(result[0][1])

The code only returns my initial values of mZ_start and ΓZ_start, not the maximum point of my chi squared fit. Can anybody help me understand why?



source https://stackoverflow.com/questions/77578836/code-to-do-a-minimised-chi-squared-fit-on-a-function-by-varying-two-variables-in
Labels:

Comments

Post a Comment

Popular posts from this blog

Prop `className` did not match in next js app

I have written a sample code ( Github Link here ). this is a simple next js app, but giving me error when I refresh the page. This seems to be the common problem and I tried the fix provided in the internet but does not seem to fix my issue. The error is Warning: Prop className did not match. Server: "MuiBox-root MuiBox-root-1" Client: "MuiBox-root MuiBox-root-2". Did changes for _document.js, modified _app.js as mentioned in official website and solutions in stackoverflow. but nothing seems to work. Could someone take a look and help me whats wrong with the code? Via Active questions tagged javascript - Stack Overflow https://ift.tt/2FdjaAW
Post a Comment
Read more

How to show number of registered users in Laravel based on usertype?

i'm trying to display data from the database in the admin dashboard i used this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count(); echo $users; ?> and i have successfully get the correct data from the database but what if i want to display a specific data for example in this user table there is "usertype" that specify if the user is normal user or admin i want to user the same code above but to display a specific usertype i tried this: <?php use Illuminate\Support\Facades\DB; $users = DB::table('users')->count()->WHERE usertype =admin; echo $users; ?> but it didn't work, what am i doing wrong? source https://stackoverflow.com/questions/68199726/how-to-show-number-of-registered-users-in-laravel-based-on-usertype
Post a Comment
Read more

Why is my reports service not connecting?

I am trying to pull some data from a Postgres database using Node.js and node-postures but I can't figure out why my service isn't connecting. my routes/index.js file: const express = require('express'); const router = express.Router(); const ordersCountController = require('../controllers/ordersCountController'); const ordersController = require('../controllers/ordersController'); const weeklyReportsController = require('../controllers/weeklyReportsController'); router.get('/orders_count', ordersCountController); router.get('/orders', ordersController); router.get('/weekly_reports', weeklyReportsController); module.exports = router; My controllers/weeklyReportsController.js file: const weeklyReportsService = require('../services/weeklyReportsService'); const weeklyReportsController = async (req, res) => { try { const data = await weeklyReportsService; res.json({data}) console
Post a Comment
Read more