Skip to main content

Trying to input counts that are float values into histogram with set bins

I have this set of data here that I am trying to create a histogram for:

| cos_low | cos_up | E_low | E_up | diff_cross |
|-1       |-0.9    |105    |162.9 |21.7861     |
|-1       |-0.9    |162.9  |220.8 |41.0181     |
|-1       |-0.9    |220.8  |278.7 |38.4258     |
|-1       |-0.9    |278.7  |336.6 |23.33       |
|-1       |-0.9    |336.6  |394.5 |7.47169     |
|-1       |-0.9    |394.5  |452.4 |1.41047     |
|-1       |-0.9    |452.4  |510.3 |0.190604    |
|-1       |-0.9    |510.3  |568.2 |0.0381209   |
|-1       |-0.9    |568.2  |626.1 |0           |
|-1       |-0.9    |626.1  |684   |0           |

(There is no header in the original text file)

This is just a portion of the full set of data I have but the format is the same. It contains the bin edges for the angle and energy. The energy bin goes from 105 to 3000 with 50 bins. I am trying to put the diff_cross values as the counts for each energy bin into a histogram but I can't get it to plot. I've tried the code for data that's already been binned and counted from matplotlib's histogram page, specifically:

counts, bins = np.histogram(x)
plt.stairs(bins, counts)

Perhaps I'm not using it correctly. I've also tried to look around to see if I could find a similar issue but I couldn't really find one. Here is the code for what I've tried:

import numpy as np
import sympy as smp
import matplotlib.pyplot as plt
import scipy as spy
import pandas as pd

#Reads text files containing data for each flux
flux_0 = pd.read_csv("2d_bins_flux_0.txt", sep=" ", header=None)

#Naming the columns
flux_0.columns = ['cos_low', ' ', 'cos_up', ' ', 'E_low', ' ', 'E_up', ' ', 
'diff_cross']

#Separating out data according to angles for each flux
flux_0_1 = flux_0.loc[lambda df: df['cos_low'] == -1, :] #Angle -1.0 < x < -0.9
flux_0_2 = flux_0.loc[lambda df: df['cos_low'] == -.9, :] #Angle -0.9 < x < -0.8

#This part is where the trouble starts
#I've tried some different things here, none working
counts = flux_0_1.diff_cross
bins = flux_0_1.E_low
plt.hist(bins[:-1], bins, weights=counts)

The histogram I'm trying to get is very similar to this one. I might be missing some detail that's preventing me from plotting this so I'll continue playing around. Much thanks in advance!



source https://stackoverflow.com/questions/73840325/trying-to-input-counts-that-are-float-values-into-histogram-with-set-bins

Labels

Labels:

Comments

Post a Comment

Popular posts from this blog

ValueError: X has 10 features, but LinearRegression is expecting 1 features as input

So, I am trying to predict the model but its throwing error like it has 10 features but it expacts only 1. So I am confused can anyone help me with it? more importantly its not working for me when my friend runs it. It works perfectly fine dose anyone know the reason about it? cv = KFold(n_splits = 10) all_loss = [] for i in range(9): # 1st for loop over polynomial orders poly_order = i X_train = make_polynomial(x, poly_order) loss_at_order = [] # initiate a set to collect loss for CV for train_index, test_index in cv.split(X_train): print('TRAIN:', train_index, 'TEST:', test_index) X_train_cv, X_test_cv = X_train[train_index], X_test[test_index] t_train_cv, t_test_cv = t[train_index], t[test_index] reg.fit(X_train_cv, t_train_cv) loss_at_order.append(np.mean((t_test_cv - reg.predict(X_test_cv))**2)) # collect loss at fold all_loss.append(np.mean(loss_at_order)) # collect loss at order plt.plot(np.log(al...
Post a Comment
Read more

Sorting large arrays of big numeric stings

I was solving bigSorting() problem from hackerrank: Consider an array of numeric strings where each string is a positive number with anywhere from to digits. Sort the array's elements in non-decreasing, or ascending order of their integer values and return the sorted array. I know it works as follows: def bigSorting(unsorted): return sorted(unsorted, key=int) But I didnt guess this approach earlier. Initially I tried below: def bigSorting(unsorted): int_unsorted = [int(i) for i in unsorted] int_sorted = sorted(int_unsorted) return [str(i) for i in int_sorted] However, for some of the test cases, it was showing time limit exceeded. Why is it so? PS: I dont know exactly what those test cases were as hacker rank does not reveal all test cases. source https://stackoverflow.com/questions/73007397/sorting-large-arrays-of-big-numeric-stings
Post a Comment
Read more

How to load Javascript with imported modules?

I am trying to import modules from tensorflowjs, and below is my code. test.html <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8"> <title>Document</title </head> <body> <script src="https://cdn.jsdelivr.net/npm/@tensorflow/tfjs@2.0.0/dist/tf.min.js"></script> <script type="module" src="./test.js"></script> </body> </html> test.js import * as tf from "./node_modules/@tensorflow/tfjs"; import {loadGraphModel} from "./node_modules/@tensorflow/tfjs-converter"; const MODEL_URL = './model.json'; const model = await loadGraphModel(MODEL_URL); const cat = document.getElementById('cat'); model.execute(tf.browser.fromPixels(cat)); Besides, I run the server using python -m http.server in my command prompt(Windows 10), and this is the error prompt in the console log of my browser: Failed to loa...
Post a Comment
Read more