rtang/ANASEN_analysis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

modified: .gitignore

Browse Source 
new file:   ELoss/Eloss.py made a Eloss calculator using the pyCatima library
	new file:   ELoss/alpha_loss_250torr.dat
	new file:   ELoss/proton_loss_250torr.dat
	modified:   MakeVertex.C
This commit is contained in:
Vignesh Sitaraman 2026-04-22 14:46:48 -04:00
parent 6e0100253b
commit 50a8e03eef
5 changed files with 585 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
.gitignore vendored
View File

@ -30,4 +30,5 @@ qqq_gains_det3.dat
qqq_relative_gains.dat qqq_relative_gains.dat
Armory/CorrelateQQQ.h Armory/CorrelateQQQ.h
QQQStage2.C QQQStage2.C
anasen_fem/scalars.dat.names anasen_fem/scalars.dat.names
myenv/

76
ELoss/Eloss.py Normal file
View File

@ -0,0 +1,76 @@
#VSitaraman 2024-06-20
# This script calculates the stopping power and range for protons and alpha particles in a gas mixture of 96% Helium and 4% CO2 at 250 Torr using the Catima library.
import numpy as np
# --- 1. Experimental & Physical Constants ---
P_TORR = 250
TEMP_K = 293.15
R = 8.3144
MEV2U = 1.0 / 931.494 # Conversion factor for mass
# Pressure to Pascal and Molar Density
p_pa = P_TORR * 133.322
molar_density = p_pa / (R * TEMP_K)
# Nuclear Data (Masses in u)
m_he = 4.002603
m_c = 12.00000
m_o = 15.99491
m_co2 = m_c + (2 * m_o)
# Effective mass of mixture (96% He, 4% CO2)
m_mix_avg = (0.96 * m_he) + (0.04 * m_co2)
# Density in g/cm3
rho = (molar_density * m_mix_avg) / 1e6
# --- 2. Create the Material ---
# We treat 0.96 and 0.04 as the stoichiometry/abundance in the mixture
material_def = [
(m_he, 2, 0.96), # Helium
(m_c, 6, 0.04), # Carbon (from CO2)
(m_o, 8, 0.08) # Oxygen (2 * 0.04 from CO2)
]
gas_mix = catima.Material(material_def)
gas_mix.density(rho)
# --- 3. Generation Function ---
def save_loss_file(z, mass_u, label):
filename = f"{label}_loss_250torr.dat"
# Energy range from 0.05 MeV/u to 50 MeV/u
energies_u = np.logspace(-1.3, 1.7, 250)
# Initialize Projectile: (Mass_u, Z)
projectile = catima.Projectile(mass_u, z)
header = (f"{label.capitalize()} Energy Loss | 96% He + 4% CO2 @ {P_TORR} Torr\n"
f"Density: {rho:.6e} g/cm3\n"
f"E_total(MeV) \tdE/dx(MeV/cm) \tRange(cm)")
output_data = []
for e_u in energies_u:
# Set kinetic energy in MeV/u
projectile.T(e_u)
# Calculate stopping power (MeV / (g/cm2))
# Multiply by density to get MeV/cm
stop_power = catima.dedx(projectile, gas_mix)
dedx_mev_cm = stop_power * rho
# Calculate range (g/cm2) and convert to cm
range_g_cm2 = catima.range(projectile, gas_mix)
range_cm = range_g_cm2 / rho
e_total = e_u * mass_u
output_data.append([e_total, dedx_mev_cm, range_cm])
np.savetxt(filename, output_data, fmt='%.6e', delimiter='\t', header=header)
print(f"File created: {filename}")
# --- 4. Execute ---
# Proton: Z=1, Mass ~1.007 u
save_loss_file(1, 1.0078, "proton")
# Alpha: Z=2, Mass ~4.0026 u
save_loss_file(2, 4.0026, "alpha")
print(f"\nCompleted. Calculated gas density: {rho:.6e} g/cm3")

253
ELoss/alpha_loss_250torr.dat Normal file
View File

@ -0,0 +1,253 @@
# Alpha Energy Loss | 96% He + 4% CO2 @ 250 Torr
# Density: 7.660749e-05 g/cm3
# E_total(MeV) dE/dx(MeV/cm) Range(cm)
2.006052e-01 1.383778e-01 2.186560e+00
2.062483e-01 1.399877e-01 2.227104e+00
2.120501e-01 1.416105e-01 2.268310e+00
2.180152e-01 1.432457e-01 2.310191e+00
2.241480e-01 1.448927e-01 2.352759e+00
2.304534e-01 1.465509e-01 2.396029e+00
2.369362e-01 1.482198e-01 2.440013e+00
2.436013e-01 1.498984e-01 2.484728e+00
2.504539e-01 1.515859e-01 2.530186e+00
2.574992e-01 1.532815e-01 2.576405e+00
2.647428e-01 1.549839e-01 2.623400e+00
2.721901e-01 1.566921e-01 2.671188e+00
2.798469e-01 1.584046e-01 2.719787e+00
2.877191e-01 1.601202e-01 2.769215e+00
2.958128e-01 1.618373e-01 2.819493e+00
3.041341e-01 1.635541e-01 2.870638e+00
3.126895e-01 1.652689e-01 2.922674e+00
3.214856e-01 1.669797e-01 2.975622e+00
3.305291e-01 1.686844e-01 3.029505e+00
3.398270e-01 1.703807e-01 3.084348e+00
3.493864e-01 1.720663e-01 3.140178e+00
3.592148e-01 1.737385e-01 3.197020e+00
3.693197e-01 1.753948e-01 3.254904e+00
3.797088e-01 1.770322e-01 3.313860e+00
3.903901e-01 1.786477e-01 3.373920e+00
4.013719e-01 1.802381e-01 3.435118e+00
4.126627e-01 1.818002e-01 3.497489e+00
4.242710e-01 1.833306e-01 3.561072e+00
4.362060e-01 1.848257e-01 3.625906e+00
4.484766e-01 1.862817e-01 3.692034e+00
4.610924e-01 1.876951e-01 3.759500e+00
4.740631e-01 1.890617e-01 3.828353e+00
4.873987e-01 1.903778e-01 3.898641e+00
5.011094e-01 1.916393e-01 3.970419e+00
5.152058e-01 1.928421e-01 4.043743e+00
5.296987e-01 1.939822e-01 4.118673e+00
5.445994e-01 1.950555e-01 4.195273e+00
5.599192e-01 1.960580e-01 4.273609e+00
5.756699e-01 1.969855e-01 4.353753e+00
5.918637e-01 1.978343e-01 4.435781e+00
6.085131e-01 1.986005e-01 4.519773e+00
6.256308e-01 1.992803e-01 4.605813e+00
6.432300e-01 1.998703e-01 4.693993e+00
6.613243e-01 2.003670e-01 4.784407e+00
6.799276e-01 2.007674e-01 4.877156e+00
6.990542e-01 2.010685e-01 4.972348e+00
7.187189e-01 2.012678e-01 5.070096e+00
7.389367e-01 2.013629e-01 5.170520e+00
7.597233e-01 2.013518e-01 5.273748e+00
7.810946e-01 2.012329e-01 5.379913e+00
8.030671e-01 2.010048e-01 5.489160e+00
8.256577e-01 2.006667e-01 5.601638e+00
8.488837e-01 2.002180e-01 5.717508e+00
8.727632e-01 1.996587e-01 5.836937e+00
8.973143e-01 1.989890e-01 5.960105e+00
9.225561e-01 1.982096e-01 6.087199e+00
9.485080e-01 1.973219e-01 6.218420e+00
9.751899e-01 1.963272e-01 6.353978e+00
1.002622e+00 1.952275e-01 6.494094e+00
1.030826e+00 1.940253e-01 6.639005e+00
1.059824e+00 1.927232e-01 6.788957e+00
1.089637e+00 1.913244e-01 6.944212e+00
1.120289e+00 1.898324e-01 7.105045e+00
1.151803e+00 1.882507e-01 7.271748e+00
1.184204e+00 1.865836e-01 7.444627e+00
1.217516e+00 1.848353e-01 7.624003e+00
1.251765e+00 1.830102e-01 7.810217e+00
1.286978e+00 1.811131e-01 8.003628e+00
1.323181e+00 1.791489e-01 8.204612e+00
1.360403e+00 1.771224e-01 8.413565e+00
1.398671e+00 1.750388e-01 8.630903e+00
1.438017e+00 1.729032e-01 8.857067e+00
1.478468e+00 1.707205e-01 9.092515e+00
1.520058e+00 1.684960e-01 9.337732e+00
1.562818e+00 1.662347e-01 9.593228e+00
1.606781e+00 1.639414e-01 9.859535e+00
1.651980e+00 1.616211e-01 1.013721e+01
1.698451e+00 1.592785e-01 1.042686e+01
1.746229e+00 1.569182e-01 1.072907e+01
1.795351e+00 1.545444e-01 1.104452e+01
1.845855e+00 1.521615e-01 1.137387e+01
1.897780e+00 1.497734e-01 1.171783e+01
1.951165e+00 1.473839e-01 1.207716e+01
2.006052e+00 1.449966e-01 1.245263e+01
2.062483e+00 1.426148e-01 1.284507e+01
2.120501e+00 1.402416e-01 1.325533e+01
2.180152e+00 1.378799e-01 1.368431e+01
2.241480e+00 1.355323e-01 1.413296e+01
2.304534e+00 1.332012e-01 1.460226e+01
2.369362e+00 1.308888e-01 1.509325e+01
2.436013e+00 1.285970e-01 1.560701e+01
2.504539e+00 1.263277e-01 1.614467e+01
2.574992e+00 1.240824e-01 1.670742e+01
2.647428e+00 1.218624e-01 1.729651e+01
2.721901e+00 1.196688e-01 1.791324e+01
2.798469e+00 1.175028e-01 1.855898e+01
2.877191e+00 1.153651e-01 1.923515e+01
2.958128e+00 1.132564e-01 1.994325e+01
3.041341e+00 1.111773e-01 2.068486e+01
3.126895e+00 1.091281e-01 2.146162e+01
3.214856e+00 1.071092e-01 2.227526e+01
3.305291e+00 1.051208e-01 2.312758e+01
3.398270e+00 1.031628e-01 2.402047e+01
3.493864e+00 1.012354e-01 2.495594e+01
3.592148e+00 9.933846e-02 2.593606e+01
3.693197e+00 9.747178e-02 2.696303e+01
3.797088e+00 9.563517e-02 2.803913e+01
3.903901e+00 9.382837e-02 2.916678e+01
4.013719e+00 9.205106e-02 3.034851e+01
4.126627e+00 9.030289e-02 3.158697e+01
4.242710e+00 8.858349e-02 3.288495e+01
4.362060e+00 8.689244e-02 3.424537e+01
4.484766e+00 8.522932e-02 3.567133e+01
4.610924e+00 8.359369e-02 3.716605e+01
4.740631e+00 8.198509e-02 3.873292e+01
4.873987e+00 8.040308e-02 4.037552e+01
5.011094e+00 7.884717e-02 4.209761e+01
5.152058e+00 7.731693e-02 4.390313e+01
5.296987e+00 7.581190e-02 4.579623e+01
5.445994e+00 7.433163e-02 4.778129e+01
5.599192e+00 7.287567e-02 4.986290e+01
5.756699e+00 7.144362e-02 5.204590e+01
5.918637e+00 7.003504e-02 5.433536e+01
6.085131e+00 6.864953e-02 5.673666e+01
6.256308e+00 6.728672e-02 5.925542e+01
6.432300e+00 6.594622e-02 6.189758e+01
6.613243e+00 6.462768e-02 6.466938e+01
6.799276e+00 6.333075e-02 6.757741e+01
6.990542e+00 6.205510e-02 7.062859e+01
7.187189e+00 6.080043e-02 7.383021e+01
7.389367e+00 5.956643e-02 7.718995e+01
7.597233e+00 5.835282e-02 8.071590e+01
7.810946e+00 5.715932e-02 8.441659e+01
8.030671e+00 5.598567e-02 8.830098e+01
8.256577e+00 5.483164e-02 9.237853e+01
8.488837e+00 5.369698e-02 9.665919e+01
8.727632e+00 5.258146e-02 1.011535e+02
8.973143e+00 5.148486e-02 1.058724e+02
9.225561e+00 5.040698e-02 1.108276e+02
9.485080e+00 4.934762e-02 1.160313e+02
9.751899e+00 4.830657e-02 1.214965e+02
1.002622e+01 4.728365e-02 1.272368e+02
1.030826e+01 4.627866e-02 1.332665e+02
1.059824e+01 4.529142e-02 1.396007e+02
1.089637e+01 4.432176e-02 1.462552e+02
1.120289e+01 4.336949e-02 1.532470e+02
1.151803e+01 4.243444e-02 1.605935e+02
1.184204e+01 4.151642e-02 1.683134e+02
1.217516e+01 4.061526e-02 1.764263e+02
1.251765e+01 3.973079e-02 1.849528e+02
1.286978e+01 3.886282e-02 1.939146e+02
1.323181e+01 3.801117e-02 2.033347e+02
1.360403e+01 3.717566e-02 2.132370e+02
1.398671e+01 3.635611e-02 2.236471e+02
1.438017e+01 3.555233e-02 2.345916e+02
1.478468e+01 3.476413e-02 2.460988e+02
1.520058e+01 3.399133e-02 2.581983e+02
1.562818e+01 3.323372e-02 2.709214e+02
1.606781e+01 3.249112e-02 2.843010e+02
1.651980e+01 3.176332e-02 2.983717e+02
1.698451e+01 3.105012e-02 3.131702e+02
1.746229e+01 3.035133e-02 3.287348e+02
1.795351e+01 2.966673e-02 3.451061e+02
1.845855e+01 2.899612e-02 3.623268e+02
1.897780e+01 2.833929e-02 3.804419e+02
1.951165e+01 2.769603e-02 3.994986e+02
2.006052e+01 2.706613e-02 4.195470e+02
2.062483e+01 2.644937e-02 4.406395e+02
2.120501e+01 2.584555e-02 4.628315e+02
2.180152e+01 2.525445e-02 4.861813e+02
2.241480e+01 2.467585e-02 5.107504e+02
2.304534e+01 2.410954e-02 5.366035e+02
2.369362e+01 2.355531e-02 5.638086e+02
2.436013e+01 2.301294e-02 5.924378e+02
2.504539e+01 2.248222e-02 6.225666e+02
2.574992e+01 2.196293e-02 6.542747e+02
2.647428e+01 2.145488e-02 6.876462e+02
2.721901e+01 2.095784e-02 7.227695e+02
2.798469e+01 2.047161e-02 7.597379e+02
2.877191e+01 1.999599e-02 7.986497e+02
2.958128e+01 1.953076e-02 8.396084e+02
3.041341e+01 1.907573e-02 8.827230e+02
3.126895e+01 1.863070e-02 9.281086e+02
3.214856e+01 1.819548e-02 9.758863e+02
3.305291e+01 1.776986e-02 1.026184e+03
3.398270e+01 1.735365e-02 1.079136e+03
3.493864e+01 1.694666e-02 1.134884e+03
3.592148e+01 1.654872e-02 1.193577e+03
3.693197e+01 1.615963e-02 1.255374e+03
3.797088e+01 1.577921e-02 1.320440e+03
3.903901e+01 1.540728e-02 1.388950e+03
4.013719e+01 1.504319e-02 1.461088e+03
4.126627e+01 1.468318e-02 1.537064e+03
4.242710e+01 1.433174e-02 1.617093e+03
4.362060e+01 1.398869e-02 1.701391e+03
4.484766e+01 1.365387e-02 1.790185e+03
4.610924e+01 1.332712e-02 1.883716e+03
4.740631e+01 1.300827e-02 1.982235e+03
4.873987e+01 1.269716e-02 2.086008e+03
5.011094e+01 1.239366e-02 2.195314e+03
5.152058e+01 1.209759e-02 2.310445e+03
5.296987e+01 1.180881e-02 2.431711e+03
5.445994e+01 1.152718e-02 2.559436e+03
5.599192e+01 1.125255e-02 2.693961e+03
5.756699e+01 1.098479e-02 2.835643e+03
5.918637e+01 1.072375e-02 2.984859e+03
6.085131e+01 1.046930e-02 3.142004e+03
6.256308e+01 1.022132e-02 3.307494e+03
6.432300e+01 9.979656e-03 3.481762e+03
6.613243e+01 9.744203e-03 3.665267e+03
6.799276e+01 9.514839e-03 3.858487e+03
6.990542e+01 9.291427e-03 4.061925e+03
7.187189e+01 9.073864e-03 4.276108e+03
7.389367e+01 8.862036e-03 4.501589e+03
7.597233e+01 8.655815e-03 4.738944e+03
7.810946e+01 8.455107e-03 4.988779e+03
8.030671e+01 8.259805e-03 5.251728e+03
8.256577e+01 8.069788e-03 5.528454e+03
8.488837e+01 7.884972e-03 5.819647e+03
8.727632e+01 7.705254e-03 6.126031e+03
8.973143e+01 7.530523e-03 6.448363e+03
9.225561e+01 7.360702e-03 6.787430e+03
9.485080e+01 7.195693e-03 7.144053e+03
9.751899e+01 7.035395e-03 7.519090e+03
1.002622e+02 6.879736e-03 7.913431e+03
1.030826e+02 6.728624e-03 8.328004e+03
1.059824e+02 6.581968e-03 8.763774e+03
1.089637e+02 6.439700e-03 9.221742e+03
1.120289e+02 6.301734e-03 9.702954e+03
1.151803e+02 6.167990e-03 1.020847e+04
1.184204e+02 6.038403e-03 1.073942e+04
1.217516e+02 5.907706e-03 1.129715e+04
1.251765e+02 5.775150e-03 1.188348e+04
1.286978e+02 5.645666e-03 1.250022e+04
1.323181e+02 5.519178e-03 1.314883e+04
1.360403e+02 5.395621e-03 1.383096e+04
1.398671e+02 5.274938e-03 1.454834e+04
1.438017e+02 5.157056e-03 1.530276e+04
1.478468e+02 5.041919e-03 1.609613e+04
1.520058e+02 4.929469e-03 1.693043e+04
1.562818e+02 4.819639e-03 1.780776e+04
1.606781e+02 4.712376e-03 1.873030e+04
1.651980e+02 4.607627e-03 1.970038e+04
1.698451e+02 4.505327e-03 2.072041e+04
1.746229e+02 4.405427e-03 2.179293e+04
1.795351e+02 4.307877e-03 2.292060e+04
1.845855e+02 4.212616e-03 2.410623e+04
1.897780e+02 4.119599e-03 2.535276e+04
1.951165e+02 4.028778e-03 2.666327e+04
2.006052e+02 3.940096e-03 2.804100e+04

253
ELoss/proton_loss_250torr.dat Normal file
View File

@ -0,0 +1,253 @@
# Proton Energy Loss | 96% He + 4% CO2 @ 250 Torr
# Density: 7.660749e-05 g/cm3
# E_total(MeV) dE/dx(MeV/cm) Range(cm)
5.050965e-02 5.849786e-02 1.113037e+00
5.193051e-02 5.892181e-02 1.137238e+00
5.339133e-02 5.934106e-02 1.161942e+00
5.489325e-02 5.975510e-02 1.187164e+00
5.643742e-02 6.016342e-02 1.212917e+00
5.802502e-02 6.056548e-02 1.239217e+00
5.965729e-02 6.096073e-02 1.266079e+00
6.133547e-02 6.134858e-02 1.293520e+00
6.306086e-02 6.172843e-02 1.321557e+00
6.483479e-02 6.209967e-02 1.350207e+00
6.665861e-02 6.246165e-02 1.379491e+00
6.853375e-02 6.281371e-02 1.409426e+00
7.046163e-02 6.315518e-02 1.440034e+00
7.244374e-02 6.348535e-02 1.471337e+00
7.448161e-02 6.380351e-02 1.503355e+00
7.657681e-02 6.410892e-02 1.536114e+00
7.873094e-02 6.440086e-02 1.569638e+00
8.094567e-02 6.467854e-02 1.603953e+00
8.322270e-02 6.494121e-02 1.639086e+00
8.556379e-02 6.518808e-02 1.675066e+00
8.797073e-02 6.541837e-02 1.711923e+00
9.044538e-02 6.563127e-02 1.749689e+00
9.298965e-02 6.582601e-02 1.788396e+00
9.560548e-02 6.600177e-02 1.828080e+00
9.829490e-02 6.615778e-02 1.868779e+00
1.010600e-01 6.629324e-02 1.910530e+00
1.039028e-01 6.640738e-02 1.953375e+00
1.068257e-01 6.649945e-02 1.997356e+00
1.098307e-01 6.656871e-02 2.042520e+00
1.129203e-01 6.661443e-02 2.088915e+00
1.160968e-01 6.663593e-02 2.136590e+00
1.193626e-01 6.663255e-02 2.185600e+00
1.227203e-01 6.660368e-02 2.236001e+00
1.261725e-01 6.654872e-02 2.287853e+00
1.297218e-01 6.646715e-02 2.341217e+00
1.333709e-01 6.635850e-02 2.396162e+00
1.371227e-01 6.622232e-02 2.452756e+00
1.409800e-01 6.605826e-02 2.511074e+00
1.449458e-01 6.586601e-02 2.571195e+00
1.490232e-01 6.564536e-02 2.633202e+00
1.532153e-01 6.539613e-02 2.697181e+00
1.575253e-01 6.511827e-02 2.763226e+00
1.619565e-01 6.481176e-02 2.831434e+00
1.665124e-01 6.447669e-02 2.901909e+00
1.711965e-01 6.411324e-02 2.974760e+00
1.760123e-01 6.372167e-02 3.050103e+00
1.809636e-01 6.330234e-02 3.128060e+00
1.860542e-01 6.285567e-02 3.208760e+00
1.912879e-01 6.238220e-02 3.292341e+00
1.966690e-01 6.188256e-02 3.378945e+00
2.022013e-01 6.135745e-02 3.468726e+00
2.078893e-01 6.080766e-02 3.561846e+00
2.137373e-01 6.023407e-02 3.658473e+00
2.197498e-01 5.963763e-02 3.758789e+00
2.259315e-01 5.901936e-02 3.862983e+00
2.322870e-01 5.838038e-02 3.971256e+00
2.388214e-01 5.772182e-02 4.083818e+00
2.455395e-01 5.704492e-02 4.200894e+00
2.524466e-01 5.635092e-02 4.322720e+00
2.595480e-01 5.564115e-02 4.449542e+00
2.668492e-01 5.491694e-02 4.581624e+00
2.743558e-01 5.417965e-02 4.719241e+00
2.820735e-01 5.343068e-02 4.862684e+00
2.900083e-01 5.267143e-02 5.012259e+00
2.981664e-01 5.190328e-02 5.168289e+00
3.065539e-01 5.112763e-02 5.331111e+00
3.151774e-01 5.034586e-02 5.501084e+00
3.240435e-01 4.955933e-02 5.678583e+00
3.331589e-01 4.876936e-02 5.864001e+00
3.425308e-01 4.797724e-02 6.057753e+00
3.521663e-01 4.718423e-02 6.260274e+00
3.620729e-01 4.639152e-02 6.472022e+00
3.722582e-01 4.560028e-02 6.693475e+00
3.827299e-01 4.481159e-02 6.925137e+00
3.934963e-01 4.402648e-02 7.167537e+00
4.045655e-01 4.324593e-02 7.421227e+00
4.159461e-01 4.247085e-02 7.686787e+00
4.276468e-01 4.170205e-02 7.964827e+00
4.396766e-01 4.094033e-02 8.255982e+00
4.520449e-01 4.018637e-02 8.560922e+00
4.647611e-01 3.944081e-02 8.880344e+00
4.778350e-01 3.870421e-02 9.214981e+00
4.912767e-01 3.797707e-02 9.565600e+00
5.050965e-01 3.725982e-02 9.933004e+00
5.193051e-01 3.655283e-02 1.031803e+01
5.339133e-01 3.585643e-02 1.072156e+01
5.489325e-01 3.517085e-02 1.114452e+01
5.643742e-01 3.449630e-02 1.158786e+01
5.802502e-01 3.383294e-02 1.205260e+01
5.965729e-01 3.318087e-02 1.253980e+01
6.133547e-01 3.254015e-02 1.305055e+01
6.306086e-01 3.191080e-02 1.358602e+01
6.483479e-01 3.129281e-02 1.414742e+01
6.665861e-01 3.068613e-02 1.473601e+01
6.853375e-01 3.009068e-02 1.535313e+01
7.046163e-01 2.950638e-02 1.600018e+01
7.244374e-01 2.893308e-02 1.667860e+01
7.448161e-01 2.837065e-02 1.738993e+01
7.657681e-01 2.781893e-02 1.813577e+01
7.873094e-01 2.727774e-02 1.891780e+01
8.094567e-01 2.674691e-02 1.973779e+01
8.322270e-01 2.622623e-02 2.059757e+01
8.556379e-01 2.571552e-02 2.149910e+01
8.797073e-01 2.521455e-02 2.244440e+01
9.044538e-01 2.472314e-02 2.343560e+01
9.298965e-01 2.424107e-02 2.447495e+01
9.560548e-01 2.376812e-02 2.556479e+01
9.829490e-01 2.330410e-02 2.670759e+01
1.010600e+00 2.284879e-02 2.790595e+01
1.039028e+00 2.240200e-02 2.916257e+01
1.068257e+00 2.196353e-02 3.048033e+01
1.098307e+00 2.153318e-02 3.186221e+01
1.129203e+00 2.111077e-02 3.331138e+01
1.160968e+00 2.069611e-02 3.483115e+01
1.193626e+00 2.028902e-02 3.642500e+01
1.227203e+00 1.988934e-02 3.809659e+01
1.261725e+00 1.949689e-02 3.984976e+01
1.297218e+00 1.911151e-02 4.168857e+01
1.333709e+00 1.873307e-02 4.361726e+01
1.371227e+00 1.836139e-02 4.564030e+01
1.409800e+00 1.799635e-02 4.776240e+01
1.449458e+00 1.763781e-02 4.998850e+01
1.490232e+00 1.728563e-02 5.232380e+01
1.532153e+00 1.693969e-02 5.477377e+01
1.575253e+00 1.659987e-02 5.734416e+01
1.619565e+00 1.626606e-02 6.004101e+01
1.665124e+00 1.593814e-02 6.287071e+01
1.711965e+00 1.561600e-02 6.583995e+01
1.760123e+00 1.529956e-02 6.895577e+01
1.809636e+00 1.498869e-02 7.222559e+01
1.860542e+00 1.468332e-02 7.565721e+01
1.912879e+00 1.438335e-02 7.925884e+01
1.966690e+00 1.408869e-02 8.303912e+01
2.022013e+00 1.379925e-02 8.700715e+01
2.078893e+00 1.351496e-02 9.117250e+01
2.137373e+00 1.323573e-02 9.554523e+01
2.197498e+00 1.296148e-02 1.001360e+02
2.259315e+00 1.269215e-02 1.049558e+02
2.322870e+00 1.242765e-02 1.100166e+02
2.388214e+00 1.216792e-02 1.153307e+02
2.455395e+00 1.191288e-02 1.209110e+02
2.524466e+00 1.166246e-02 1.267713e+02
2.595480e+00 1.141661e-02 1.329260e+02
2.668492e+00 1.117525e-02 1.393904e+02
2.743558e+00 1.093831e-02 1.461803e+02
2.820735e+00 1.070574e-02 1.533127e+02
2.900083e+00 1.047747e-02 1.608052e+02
2.981664e+00 1.025345e-02 1.686766e+02
3.065539e+00 1.003360e-02 1.769465e+02
3.151774e+00 9.817869e-03 1.856355e+02
3.240435e+00 9.606199e-03 1.947656e+02
3.331589e+00 9.398531e-03 2.043597e+02
3.425308e+00 9.194805e-03 2.144418e+02
3.521663e+00 8.994964e-03 2.250376e+02
3.620729e+00 8.798952e-03 2.361738e+02
3.722582e+00 8.606711e-03 2.478787e+02
3.827299e+00 8.418186e-03 2.601820e+02
3.934963e+00 8.233322e-03 2.731150e+02
4.045655e+00 8.052062e-03 2.867108e+02
4.159461e+00 7.874353e-03 3.010041e+02
4.276468e+00 7.700139e-03 3.160316e+02
4.396766e+00 7.529368e-03 3.318317e+02
4.520449e+00 7.361986e-03 3.484453e+02
4.647611e+00 7.197939e-03 3.659150e+02
4.778350e+00 7.037176e-03 3.842859e+02
4.912767e+00 6.879643e-03 4.036057e+02
5.050965e+00 6.725289e-03 4.239243e+02
5.193051e+00 6.574063e-03 4.452944e+02
5.339133e+00 6.425913e-03 4.677718e+02
5.489325e+00 6.280789e-03 4.914147e+02
5.643742e+00 6.138641e-03 5.162851e+02
5.802502e+00 5.999419e-03 5.424478e+02
5.965729e+00 5.863074e-03 5.699714e+02
6.133547e+00 5.729557e-03 5.989279e+02
6.306086e+00 5.598819e-03 6.293935e+02
6.483479e+00 5.470812e-03 6.614483e+02
6.665861e+00 5.345490e-03 6.951767e+02
6.853375e+00 5.222804e-03 7.306676e+02
7.046163e+00 5.102709e-03 7.680150e+02
7.244374e+00 4.985158e-03 8.073174e+02
7.448161e+00 4.870106e-03 8.486793e+02
7.657681e+00 4.757507e-03 8.922102e+02
7.873094e+00 4.647318e-03 9.380258e+02
8.094567e+00 4.539494e-03 9.862481e+02
8.322270e+00 4.433991e-03 1.037006e+03
8.556379e+00 4.330766e-03 1.090434e+03
8.797073e+00 4.229778e-03 1.146675e+03
9.044538e+00 4.130983e-03 1.205881e+03
9.298965e+00 4.034341e-03 1.268208e+03
9.560548e+00 3.939811e-03 1.333826e+03
9.829490e+00 3.847351e-03 1.402909e+03
1.010600e+01 3.756808e-03 1.475644e+03
1.039028e+01 3.667293e-03 1.552240e+03
1.068257e+01 3.579867e-03 1.632913e+03
1.098307e+01 3.494490e-03 1.717882e+03
1.129203e+01 3.411128e-03 1.807376e+03
1.160968e+01 3.329742e-03 1.901636e+03
1.193626e+01 3.250297e-03 2.000917e+03
1.227203e+01 3.172757e-03 2.105485e+03
1.261725e+01 3.097086e-03 2.215622e+03
1.297218e+01 3.023249e-03 2.331623e+03
1.333709e+01 2.951212e-03 2.453799e+03
1.371227e+01 2.880940e-03 2.582477e+03
1.409800e+01 2.812400e-03 2.718001e+03
1.449458e+01 2.745559e-03 2.860730e+03
1.490232e+01 2.680385e-03 3.011046e+03
1.532153e+01 2.616845e-03 3.169345e+03
1.575253e+01 2.554909e-03 3.336045e+03
1.619565e+01 2.494545e-03 3.511586e+03
1.665124e+01 2.435724e-03 3.696428e+03
1.711965e+01 2.378415e-03 3.891054e+03
1.760123e+01 2.322589e-03 4.095970e+03
1.809636e+01 2.268219e-03 4.311708e+03
1.860542e+01 2.215276e-03 4.538824e+03
1.912879e+01 2.163732e-03 4.777900e+03
1.966690e+01 2.113561e-03 5.029546e+03
2.022013e+01 2.064737e-03 5.294401e+03
2.078893e+01 2.017233e-03 5.573133e+03
2.137373e+01 1.971026e-03 5.866438e+03
2.197498e+01 1.926090e-03 6.175047e+03
2.259315e+01 1.882401e-03 6.499721e+03
2.322870e+01 1.839936e-03 6.841253e+03
2.388214e+01 1.798672e-03 7.200473e+03
2.455395e+01 1.758586e-03 7.578244e+03
2.524466e+01 1.719658e-03 7.975464e+03
2.595480e+01 1.681865e-03 8.393069e+03
2.668492e+01 1.645187e-03 8.832031e+03
2.743558e+01 1.609603e-03 9.293361e+03
2.820735e+01 1.575095e-03 9.778112e+03
2.900083e+01 1.541642e-03 1.028735e+04
2.981664e+01 1.509227e-03 1.082222e+04
3.065539e+01 1.476542e-03 1.138408e+04
3.151774e+01 1.443401e-03 1.197476e+04
3.240435e+01 1.411027e-03 1.259607e+04
3.331589e+01 1.379403e-03 1.324951e+04
3.425308e+01 1.348513e-03 1.393671e+04
3.521663e+01 1.318340e-03 1.465942e+04
3.620729e+01 1.288869e-03 1.541947e+04
3.722582e+01 1.260084e-03 1.621875e+04
3.827299e+01 1.231971e-03 1.705928e+04
3.934963e+01 1.204514e-03 1.794316e+04
4.045655e+01 1.177698e-03 1.887261e+04
4.159461e+01 1.151511e-03 1.984995e+04
4.276468e+01 1.125937e-03 2.087762e+04
4.396766e+01 1.100963e-03 2.195819e+04
4.520449e+01 1.076576e-03 2.309433e+04
4.647611e+01 1.052762e-03 2.428887e+04
4.778350e+01 1.029509e-03 2.554478e+04
4.912767e+01 1.006805e-03 2.686516e+04
5.050965e+01 9.846366e-04 2.825327e+04

2
MakeVertex.C
View File

@ -156,7 +156,7 @@ void MakeVertex::Begin(TTree * /*tree*/)
fflush(stdout); fflush(stdout);
//usleep(4e5); //usleep(4e5);
// Load PC Calibrations // Load PC Calibrations
std::ifstream inputFile("slope_intercept_results_"+dataset+".txt"); std::ifstream inputFile("slope_intercept_results_"+dataset+".dat");
if (inputFile.is_open()) if (inputFile.is_open())
{ {
std::string line; std::string line;