rand = TRandom()

muonMass = {}

for i in range(10000):
    muonMass[0] = rand.Exp(30)*90
    muonMass[1] = abs(rand.Gaus(91.2,1.3) - muonMass[0])
    muonMass[2] = rand.Exp(30)
    muonMass[3] = rand.Exp(30)
    
    muons = [[i] for i in range(4)]
    for i in range(4):
        muons[i] = muonMass[i]
    shuffle(muons)

    mystruct.f_muonMass1   = muons[0]
    mystruct.f_muonMass2   = muons[1]
    mystruct.f_muonMass3   = muons[2]
    mystruct.f_muonMass4   = muons[3]

    tree.Fill()




############################
######### ANALYSIS #########
############################
allMuonHist = TH1F("allMuonHist","All Muon Mass",1000,0,1000)
diMuonHist12 = TH1F("diMuonHist12","Di Muon Mass",500,0,500)
diMuonHist23 = TH1F("diMuonHist23","Di Muon Mass",500,0,500)