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)