del_R_b_anti_b = []
del_R_mu_other_b = []
del_R_mu_b = []

minimal_dR = []
minimal_pt = []

print "Reading in the data...."
collisions = hep.get_collisions(f)

print len(collisions)

count = 0
for collision in collisions:

    gen_particles, ak5jets, ca8jets = hep.get_truth_particles_from_collision(
        collision)

    #print gen_particles
    #print ak5jets
    #print ca8jets

    #exit()

    #print "--------------------------"
    pdgs = []
    pts = []
    etas = []
    phis = []
    jet_pt = []
    jet_eta = []
    jet_phi = []
del_Theta=[]
del_R_mu_b=[]





print "Reading in the data...."
collisions = hep.get_collisions(f)

print len(collisions)

count = 0
for collision in collisions:

    gen_particles,ak5jets,ca8jets = hep.get_truth_particles_from_collision(collision)


    pdgs=[]
    pts=[]
    etas=[]
    phis=[]
    jet_pt=[]
    jet_eta=[]
    jet_phi=[]

    for p in gen_particles:
        pdg,status,pt,eta,phi = p
        pdgs.append(pdg)
        pts.append(pt)
        etas.append(eta)