diff --git a/data/SciBooNE/flux_fhc.txt b/data/SciBooNE/flux_fhc.txt new file mode 100644 index 0000000..1b6c8f0 --- /dev/null +++ b/data/SciBooNE/flux_fhc.txt @@ -0,0 +1,202 @@ +# From https://www-sciboone.fnal.gov/data_release/flux/flux_numode.txt +# E_low(GeV) E_high(GeV) total nue anti-nue numu anti-numu +0 0.025 2.8822e-11 1.31509e-13 5.37544e-14 1.39903e-11 1.46464e-11 +0.025 0.05 1.72977e-10 5.908e-13 1.76322e-13 8.60701e-11 8.61401e-11 +0.05 0.075 2.83006e-10 1.25253e-12 2.78521e-13 1.56765e-10 1.24709e-10 +0.075 0.1 3.24191e-10 1.44687e-12 2.73281e-13 2.01136e-10 1.21334e-10 +0.1 0.125 2.99657e-10 2.08327e-12 3.48294e-13 2.09584e-10 8.76419e-11 +0.125 0.15 2.94326e-10 1.9994e-12 2.17575e-13 2.26275e-10 6.58333e-11 +0.15 0.175 3.03569e-10 2.29971e-12 2.16336e-13 2.48995e-10 5.20581e-11 +0.175 0.2 3.18462e-10 2.286e-12 2.10615e-13 2.71575e-10 4.43898e-11 +0.2 0.225 3.45154e-10 2.34662e-12 2.17263e-13 3.02614e-10 3.99758e-11 +0.225 0.25 4.03232e-10 2.43333e-12 2.10372e-13 3.6268e-10 3.79087e-11 +0.25 0.275 4.01383e-10 2.49172e-12 1.72814e-13 3.64971e-10 3.37479e-11 +0.275 0.3 4.16553e-10 2.44537e-12 2.06293e-13 3.81392e-10 3.25088e-11 +0.3 0.325 4.36424e-10 2.64459e-12 1.73228e-13 4.023e-10 3.13068e-11 +0.325 0.35 4.52456e-10 2.51904e-12 1.83636e-13 4.20197e-10 2.95559e-11 +0.35 0.375 4.70429e-10 2.5833e-12 1.91672e-13 4.39311e-10 2.83425e-11 +0.375 0.4 4.88394e-10 2.49031e-12 2.19945e-13 4.58628e-10 2.70554e-11 +0.4 0.425 5.08162e-10 2.45046e-12 1.81591e-13 4.79551e-10 2.5978e-11 +0.425 0.45 5.21222e-10 2.39272e-12 1.64697e-13 4.93657e-10 2.50067e-11 +0.45 0.475 5.32337e-10 2.46304e-12 1.92387e-13 5.06058e-10 2.36233e-11 +0.475 0.5 5.42836e-10 2.27778e-12 1.73471e-13 5.18029e-10 2.23558e-11 +0.5 0.525 5.48144e-10 2.29219e-12 2.01434e-13 5.24078e-10 2.15726e-11 +0.525 0.55 5.50852e-10 2.25053e-12 2.06163e-13 5.27468e-10 2.09267e-11 +0.55 0.575 5.53515e-10 2.14778e-12 1.41057e-13 5.30971e-10 2.02548e-11 +0.575 0.6 5.52044e-10 2.16317e-12 1.99858e-13 5.30576e-10 1.91054e-11 +0.6 0.625 5.45932e-10 2.07367e-12 1.73836e-13 5.25508e-10 1.81764e-11 +0.625 0.65 5.39653e-10 2.03829e-12 1.73756e-13 5.19473e-10 1.79679e-11 +0.65 0.675 5.33142e-10 1.98214e-12 1.27021e-13 5.13726e-10 1.73063e-11 +0.675 0.7 5.22956e-10 1.89735e-12 1.96783e-13 5.04391e-10 1.64708e-11 +0.7 0.725 5.11204e-10 1.88347e-12 1.37892e-13 4.93976e-10 1.52058e-11 +0.725 0.75 4.98931e-10 1.76558e-12 1.50334e-13 4.81887e-10 1.51278e-11 +0.75 0.775 4.83511e-10 1.78574e-12 1.25966e-13 4.67202e-10 1.43968e-11 +0.775 0.8 4.6901e-10 1.62353e-12 1.60205e-13 4.53401e-10 1.38256e-11 +0.8 0.825 4.52677e-10 1.71987e-12 1.4512e-13 4.37737e-10 1.30742e-11 +0.825 0.85 4.35423e-10 1.57004e-12 1.40293e-13 4.21029e-10 1.26838e-11 +0.85 0.875 4.18229e-10 1.53015e-12 1.22709e-13 4.04355e-10 1.22214e-11 +0.875 0.9 3.99259e-10 1.49486e-12 9.90342e-14 3.85976e-10 1.16885e-11 +0.9 0.925 3.80367e-10 1.44038e-12 1.33772e-13 3.67788e-10 1.10049e-11 +0.925 0.95 3.62855e-10 1.45615e-12 1.28733e-13 3.50638e-10 1.06321e-11 +0.95 0.975 3.4383e-10 1.37374e-12 1.19223e-13 3.32258e-10 1.00792e-11 +0.975 1 3.24656e-10 1.27294e-12 1.61534e-13 3.13454e-10 9.76789e-12 +1 1.025 3.06879e-10 1.23161e-12 1.17673e-13 2.96532e-10 8.99813e-12 +1.025 1.05 2.88276e-10 1.21993e-12 1.12232e-13 2.78042e-10 8.90264e-12 +1.05 1.075 2.70177e-10 1.17833e-12 9.50004e-14 2.60634e-10 8.26926e-12 +1.075 1.1 2.53206e-10 1.15611e-12 1.20086e-13 2.44073e-10 7.85637e-12 +1.1 1.125 2.35756e-10 1.06533e-12 1.0562e-13 2.27027e-10 7.55748e-12 +1.125 1.15 2.2043e-10 1.0816e-12 8.86255e-14 2.122e-10 7.05972e-12 +1.15 1.175 2.05495e-10 1.08103e-12 1.06474e-13 1.97443e-10 6.86463e-12 +1.175 1.2 1.90534e-10 1.02108e-12 8.25464e-14 1.83198e-10 6.23265e-12 +1.2 1.225 1.76781e-10 9.74136e-13 1.12742e-13 1.69623e-10 6.07123e-12 +1.225 1.25 1.63433e-10 9.73232e-13 9.08882e-14 1.56811e-10 5.55704e-12 +1.25 1.275 1.51423e-10 9.1933e-13 1.01861e-13 1.45301e-10 5.10072e-12 +1.275 1.3 1.39884e-10 8.53899e-13 9.12507e-14 1.34027e-10 4.9112e-12 +1.3 1.325 1.29783e-10 8.85115e-13 8.02695e-14 1.24054e-10 4.7634e-12 +1.325 1.35 1.19712e-10 8.5994e-13 8.4219e-14 1.14314e-10 4.45385e-12 +1.35 1.375 1.10388e-10 8.0388e-13 8.57882e-14 1.05394e-10 4.10471e-12 +1.375 1.4 1.02243e-10 7.9281e-13 7.97369e-14 9.73576e-11 4.01327e-12 +1.4 1.425 9.42337e-11 7.33139e-13 7.13119e-14 8.97386e-11 3.69062e-12 +1.425 1.45 8.70625e-11 7.34457e-13 6.56455e-14 8.27238e-11 3.53862e-12 +1.45 1.475 8.01288e-11 6.95377e-13 6.56668e-14 7.61248e-11 3.24296e-12 +1.475 1.5 7.37059e-11 6.96046e-13 7.02433e-14 6.99214e-11 3.01819e-12 +1.5 1.525 6.84942e-11 7.45743e-13 7.43809e-14 6.4742e-11 2.93211e-12 +1.525 1.55 6.30213e-11 6.39512e-13 5.9934e-14 5.95523e-11 2.76952e-12 +1.55 1.575 5.80933e-11 6.3328e-13 6.92622e-14 5.48378e-11 2.55291e-12 +1.575 1.6 5.34512e-11 6.22272e-13 6.89485e-14 5.03704e-11 2.38953e-12 +1.6 1.625 4.92993e-11 5.90767e-13 5.72316e-14 4.64117e-11 2.2396e-12 +1.625 1.65 4.56293e-11 5.64915e-13 6.00168e-14 4.29206e-11 2.08371e-12 +1.65 1.675 4.22024e-11 5.31521e-13 5.60754e-14 3.9672e-11 1.94281e-12 +1.675 1.7 3.89931e-11 5.46799e-13 6.45607e-14 3.65417e-11 1.83997e-12 +1.7 1.725 3.61471e-11 4.97904e-13 5.43773e-14 3.38388e-11 1.75604e-12 +1.725 1.75 3.33796e-11 5.2927e-13 6.17084e-14 3.12227e-11 1.56586e-12 +1.75 1.775 3.08986e-11 5.13869e-13 5.35416e-14 2.8858e-11 1.47316e-12 +1.775 1.8 2.85133e-11 4.66955e-13 6.56469e-14 2.66125e-11 1.36822e-12 +1.8 1.825 2.67321e-11 4.91123e-13 4.76625e-14 2.48534e-11 1.3399e-12 +1.825 1.85 2.47177e-11 4.36559e-13 4.96381e-14 2.29634e-11 1.26806e-12 +1.85 1.875 2.30567e-11 4.60941e-13 5.13715e-14 2.14096e-11 1.13474e-12 +1.875 1.9 2.12442e-11 4.17003e-13 4.84688e-14 1.97384e-11 1.04034e-12 +1.9 1.925 1.96996e-11 4.11708e-13 4.43265e-14 1.83042e-11 9.39407e-13 +1.925 1.95 1.85507e-11 4.16464e-13 5.09792e-14 1.72243e-11 8.58949e-13 +1.95 1.975 1.72437e-11 3.94176e-13 4.12696e-14 1.59616e-11 8.46697e-13 +1.975 2 1.62087e-11 3.65697e-13 5.64389e-14 1.49702e-11 8.16429e-13 +2 2.025 1.49526e-11 3.59407e-13 3.93851e-14 1.38237e-11 7.30162e-13 +2.025 2.05 1.40632e-11 3.69166e-13 4.14981e-14 1.29782e-11 6.74365e-13 +2.05 2.075 1.3105e-11 3.28769e-13 4.0204e-14 1.21114e-11 6.2469e-13 +2.075 2.1 1.23491e-11 3.13878e-13 3.81462e-14 1.14046e-11 5.92437e-13 +2.1 2.125 1.15635e-11 2.98277e-13 3.75313e-14 1.06785e-11 5.49111e-13 +2.125 2.15 1.08505e-11 2.85466e-13 4.06669e-14 1.00494e-11 4.74889e-13 +2.15 2.175 1.02219e-11 2.85579e-13 3.91155e-14 9.4246e-12 4.72657e-13 +2.175 2.2 9.78664e-12 3.04789e-13 3.46606e-14 8.99942e-12 4.47775e-13 +2.2 2.225 9.1449e-12 2.81255e-13 4.92334e-14 8.41247e-12 4.01946e-13 +2.225 2.25 8.831e-12 2.71582e-13 3.02791e-14 8.13299e-12 3.96144e-13 +2.25 2.275 8.21362e-12 2.70383e-13 3.88703e-14 7.56284e-12 3.41536e-13 +2.275 2.3 7.83474e-12 2.49526e-13 3.811e-14 7.23769e-12 3.09413e-13 +2.3 2.325 7.44709e-12 2.35226e-13 2.91798e-14 6.90651e-12 2.76181e-13 +2.325 2.35 7.06565e-12 2.24695e-13 2.67039e-14 6.54643e-12 2.67816e-13 +2.35 2.375 6.65761e-12 2.21456e-13 2.85896e-14 6.16849e-12 2.39072e-13 +2.375 2.4 6.47984e-12 2.17267e-13 2.7226e-14 6.01798e-12 2.17372e-13 +2.4 2.425 6.20644e-12 2.01579e-13 2.14441e-14 5.77533e-12 2.08092e-13 +2.425 2.45 6.01351e-12 2.08772e-13 2.27141e-14 5.57007e-12 2.11945e-13 +2.45 2.475 5.72589e-12 1.90328e-13 2.64737e-14 5.32265e-12 1.8644e-13 +2.475 2.5 5.47782e-12 1.89915e-13 2.80894e-14 5.09896e-12 1.6086e-13 +2.5 2.525 5.33271e-12 1.7631e-13 3.08683e-14 4.97255e-12 1.52976e-13 +2.525 2.55 5.09848e-12 1.80883e-13 2.17407e-14 4.75541e-12 1.40445e-13 +2.55 2.575 5.00265e-12 1.66733e-13 1.89777e-14 4.68129e-12 1.35647e-13 +2.575 2.6 4.78551e-12 1.5252e-13 2.03882e-14 4.50246e-12 1.1014e-13 +2.6 2.625 4.66066e-12 1.62801e-13 2.02896e-14 4.37397e-12 1.03601e-13 +2.625 2.65 4.52361e-12 1.50746e-13 1.95936e-14 4.25341e-12 9.98601e-14 +2.65 2.675 4.41852e-12 1.43951e-13 1.81542e-14 4.16787e-12 8.85514e-14 +2.675 2.7 4.25497e-12 1.47109e-13 1.74095e-14 4.00295e-12 8.75063e-14 +2.7 2.725 4.14718e-12 1.30507e-13 1.76574e-14 3.91074e-12 8.82827e-14 +2.725 2.75 4.09619e-12 1.21832e-13 1.85008e-14 3.87073e-12 8.5131e-14 +2.75 2.775 4.03924e-12 1.29242e-13 2.13082e-14 3.82877e-12 5.99152e-14 +2.775 2.8 3.89925e-12 1.22346e-13 1.56767e-14 3.69991e-12 6.13133e-14 +2.8 2.825 3.76529e-12 1.09802e-13 1.55563e-14 3.57979e-12 6.01432e-14 +2.825 2.85 3.68704e-12 1.102e-13 1.50583e-14 3.51433e-12 4.74533e-14 +2.85 2.875 3.5984e-12 9.8085e-14 1.58721e-14 3.4387e-12 4.57439e-14 +2.875 2.9 3.55416e-12 1.15661e-13 1.25879e-14 3.37874e-12 4.71765e-14 +2.9 2.925 3.38886e-12 9.3374e-14 1.28854e-14 3.24135e-12 4.12496e-14 +2.925 2.95 3.29228e-12 9.52858e-14 1.2556e-14 3.1524e-12 3.20435e-14 +2.95 2.975 3.31778e-12 8.75579e-14 1.29469e-14 3.18686e-12 3.04143e-14 +2.975 3 3.19942e-12 8.25905e-14 1.04178e-14 3.07622e-12 3.01886e-14 +3 3.025 3.14372e-12 9.00918e-14 9.8719e-15 3.02364e-12 2.0117e-14 +3.025 3.05 3.17015e-12 7.80918e-14 1.93541e-14 3.05303e-12 1.96748e-14 +3.05 3.075 3.11726e-12 7.57707e-14 1.11443e-14 2.99831e-12 3.20366e-14 +3.075 3.1 2.99728e-12 7.47496e-14 1.12012e-14 2.89239e-12 1.89376e-14 +3.1 3.125 2.96392e-12 7.06854e-14 8.44872e-15 2.86223e-12 2.25584e-14 +3.125 3.15 2.90715e-12 6.81669e-14 8.51335e-15 2.81342e-12 1.70544e-14 +3.15 3.175 2.87042e-12 6.49559e-14 9.7821e-15 2.78293e-12 1.27529e-14 +3.175 3.2 2.81936e-12 6.20918e-14 7.7778e-15 2.73417e-12 1.53287e-14 +3.2 3.225 2.79602e-12 6.40616e-14 9.02572e-15 2.70388e-12 1.9059e-14 +3.225 3.25 2.70177e-12 5.29808e-14 7.10035e-15 2.62698e-12 1.47121e-14 +3.25 3.275 2.64771e-12 5.55716e-14 7.05534e-15 2.5712e-12 1.38869e-14 +3.275 3.3 2.61137e-12 5.43971e-14 6.86005e-15 2.54195e-12 8.16257e-15 +3.3 3.325 2.51269e-12 4.55791e-14 6.02029e-15 2.44834e-12 1.27487e-14 +3.325 3.35 2.54876e-12 6.1362e-14 6.13845e-15 2.46931e-12 1.19452e-14 +3.35 3.375 2.48534e-12 4.64747e-14 6.54846e-15 2.41185e-12 2.04663e-14 +3.375 3.4 2.4191e-12 4.39196e-14 5.63623e-15 2.35827e-12 1.12803e-14 +3.4 3.425 2.35355e-12 4.88507e-14 5.34503e-15 2.29327e-12 6.07872e-15 +3.425 3.45 2.29398e-12 4.14293e-14 5.52927e-15 2.23621e-12 1.08132e-14 +3.45 3.475 2.26581e-12 3.80315e-14 5.21859e-15 2.21725e-12 5.31235e-15 +3.475 3.5 2.26604e-12 3.57265e-14 8.26067e-15 2.21251e-12 9.53858e-15 +3.5 3.525 2.17865e-12 3.26542e-14 4.58895e-15 2.1319e-12 9.50783e-15 +3.525 3.55 2.16805e-12 3.39181e-14 5.37686e-15 2.11498e-12 1.37757e-14 +3.55 3.575 2.08992e-12 3.25442e-14 4.68693e-15 2.04751e-12 5.18182e-15 +3.575 3.6 2.00691e-12 2.96052e-14 3.6164e-15 1.96537e-12 8.31469e-15 +3.6 3.625 1.97338e-12 3.35164e-14 3.33579e-15 1.92811e-12 8.4237e-15 +3.625 3.65 1.92172e-12 2.59029e-14 3.97392e-15 1.88858e-12 3.26144e-15 +3.65 3.675 1.90998e-12 2.56568e-14 4.22816e-15 1.87215e-12 7.94533e-15 +3.675 3.7 1.79792e-12 2.53015e-14 3.31157e-15 1.7665e-12 2.80411e-15 +3.7 3.725 1.79792e-12 2.41838e-14 3.4701e-15 1.76772e-12 2.54333e-15 +3.725 3.75 1.7585e-12 2.42643e-14 2.76085e-15 1.7286e-12 2.87469e-15 +3.75 3.775 1.72697e-12 2.28589e-14 3.26209e-15 1.6938e-12 7.04923e-15 +3.775 3.8 1.69222e-12 1.96765e-14 2.60586e-15 1.66746e-12 2.4855e-15 +3.8 3.825 1.66658e-12 2.11844e-14 2.73895e-15 1.64017e-12 2.48236e-15 +3.825 3.85 1.59576e-12 1.86674e-14 2.12227e-15 1.56976e-12 5.20091e-15 +3.85 3.875 1.5729e-12 1.77223e-14 2.53373e-15 1.5458e-12 6.84657e-15 +3.875 3.9 1.5538e-12 1.64324e-14 2.35975e-15 1.53289e-12 2.11842e-15 +3.9 3.925 1.50724e-12 1.6059e-14 2.59838e-15 1.48682e-12 1.76614e-15 +3.925 3.95 1.45152e-12 1.56431e-14 2.27806e-15 1.43146e-12 2.13706e-15 +3.95 3.975 1.416e-12 1.55525e-14 2.22121e-15 1.39191e-12 6.31465e-15 +3.975 4 1.40021e-12 1.3182e-14 2.11597e-15 1.37842e-12 6.4926e-15 +4 4.025 1.34972e-12 1.52378e-14 1.84729e-15 1.33139e-12 1.24003e-15 +4.025 4.05 1.27945e-12 1.20633e-14 2.11946e-15 1.26377e-12 1.49851e-15 +4.05 4.075 1.25184e-12 1.22853e-14 1.70156e-15 1.23644e-12 1.41784e-15 +4.075 4.1 1.19845e-12 1.15967e-14 1.63862e-15 1.18384e-12 1.37024e-15 +4.1 4.125 1.16432e-12 1.13562e-14 1.44119e-15 1.15018e-12 1.34172e-15 +4.125 4.15 1.16352e-12 9.66668e-15 1.34104e-15 1.15123e-12 1.28884e-15 +4.15 4.175 1.11981e-12 9.74513e-15 1.29304e-15 1.10793e-12 8.45802e-16 +4.175 4.2 1.11637e-12 9.03034e-15 1.1307e-15 1.10525e-12 9.59063e-16 +4.2 4.225 1.05792e-12 9.4861e-15 1.1511e-15 1.04637e-12 9.04003e-16 +4.225 4.25 1.02823e-12 8.67947e-15 1.33376e-15 1.01742e-12 8.03884e-16 +4.25 4.275 1.00444e-12 8.73385e-15 9.88021e-16 9.94008e-13 7.14724e-16 +4.275 4.3 9.88215e-13 8.11028e-15 1.20453e-15 9.77892e-13 1.00863e-15 +4.3 4.325 9.45791e-13 7.77913e-15 9.31072e-16 9.36486e-13 5.95165e-16 +4.325 4.35 9.22066e-13 6.98414e-15 1.06215e-15 9.13073e-13 9.46357e-16 +4.35 4.375 8.68839e-13 7.6479e-15 9.5925e-16 8.59371e-13 8.60399e-16 +4.375 4.4 8.52245e-13 6.86665e-15 7.46546e-16 8.43877e-13 7.53864e-16 +4.4 4.425 8.23759e-13 6.62824e-15 8.38834e-16 8.15854e-13 4.37502e-16 +4.425 4.45 8.15444e-13 5.71784e-15 6.90826e-16 8.03871e-13 5.16403e-15 +4.45 4.475 7.68451e-13 5.52603e-15 9.19705e-16 7.61284e-13 7.21315e-16 +4.475 4.5 7.6842e-13 4.80607e-15 7.73047e-16 7.48002e-13 1.48389e-14 +4.5 4.525 7.52038e-13 4.93254e-15 6.38427e-16 7.41285e-13 5.18183e-15 +4.525 4.55 7.20485e-13 4.67466e-15 6.88256e-16 7.14395e-13 7.27115e-16 +4.55 4.575 6.91817e-13 4.50743e-15 6.49754e-16 6.86451e-13 2.08807e-16 +4.575 4.6 6.82552e-13 3.76591e-15 6.41842e-16 6.73022e-13 5.12224e-15 +4.6 4.625 6.48685e-13 4.00868e-15 5.59075e-16 6.43739e-13 3.77919e-16 +4.625 4.65 6.38958e-13 3.68227e-15 4.92241e-16 6.25133e-13 9.65034e-15 +4.65 4.675 6.02433e-13 3.65448e-15 4.32836e-16 5.98144e-13 2.01085e-16 +4.675 4.7 5.82049e-13 2.94739e-15 4.88186e-16 5.78298e-13 3.14827e-16 +4.7 4.725 5.6225e-13 2.82263e-15 4.62162e-16 5.58814e-13 1.51314e-16 +4.725 4.75 5.57914e-13 3.28537e-15 4.38147e-16 5.54039e-13 1.51536e-16 +4.75 4.775 5.26916e-13 3.03382e-15 4.85317e-16 5.23276e-13 1.21156e-16 +4.775 4.8 5.09983e-13 2.75081e-15 2.79262e-16 5.06836e-13 1.16892e-16 +4.8 4.825 4.97122e-13 2.42494e-15 2.37454e-16 4.89633e-13 4.8262e-15 +4.825 4.85 4.73253e-13 1.95223e-15 2.90448e-16 4.66072e-13 4.93813e-15 +4.85 4.875 4.58272e-13 1.88281e-15 3.25333e-16 4.5596e-13 1.04015e-16 +4.875 4.9 4.44385e-13 2.08297e-15 3.68548e-16 4.4179e-13 1.43158e-16 +4.9 4.925 4.21759e-13 1.77438e-15 1.98987e-16 4.19689e-13 9.67065e-17 +4.925 4.95 4.19446e-13 1.7138e-15 2.79072e-16 4.17376e-13 7.68308e-17 +4.95 4.975 3.95305e-13 1.63623e-15 2.93307e-16 3.93275e-13 1.00819e-16 +4.975 5 3.83982e-13 1.55561e-15 2.2816e-16 3.82136e-13 6.21645e-17 diff --git a/data/SciBooNE/makefluxroot.C b/data/SciBooNE/makefluxroot.C new file mode 100644 index 0000000..2565534 --- /dev/null +++ b/data/SciBooNE/makefluxroot.C @@ -0,0 +1,47 @@ +#include "TFile.h" +#include "TH1D.h" + +#include +#include +#include + +void makefluxroot() { + + std::ifstream input("flux_fhc.txt"); + + std::string line; + + double low, up, totval, nueval, nuebarval, numuval, numubarval; + std::vector bins; + std::vector vals; + while (std::getline(input, line)) { + if (line.front() == '#') { + continue; + } + std::istringstream ss(line); + ss >> low; + ss >> up; + ss >> totval; + ss >> nueval; + ss >> nuebarval; + ss >> numuval; + ss >> numubarval; + + if (!bins.size()) { + bins.push_back(low); + } + bins.push_back(up); + vals.push_back(numuval); + } + + TFile f("sciboone_fhc.root", "RECREATE"); + + TH1D *flux = new TH1D("flux_numu", ";E_{#nu} (GeV); #Phi_{#nu};", + bins.size() - 1, bins.data()); + for (size_t i = 0; i < vals.size(); ++i) { + flux->SetBinContent(i + 1, vals[i]); + } + + f.Write(); + f.Close(); +} diff --git a/data/flux/sciboone_fhc.root b/data/flux/sciboone_fhc.root new file mode 100644 index 0000000..ecae166 Binary files /dev/null and b/data/flux/sciboone_fhc.root differ