Beispiel #1
0
import ROOT
import os
from TH1Wrapper import TH2F, TH1F
import re
import mytools2 as mt2
import time
import numpy as np

# dir = "/Users/jeffreyburggraf/PycharmProjects/TwoNeutronCorrelations/Analysis/Other"
dir = "/Volumes/JeffMacEx/2nCorrData/Production/Forward/"

ROOT.gROOT.ProcessLine(".L twontree.h+")

ergs = TH1F(0.4, 12, nbinss=20)
histos = []


def gen_FREYA_dict():
    ROOT.gROOT.ProcessLine("""
    struct coinc_hit {
        float erg[2]; 
        double theta_abs[2];
        double coinc_theta;
    };
    struct neutrons {
        std::vector<coinc_hit> coinc_hits;
        double xs;
    };
    
    struct FF{
        int A,Z;