|
| ROOT | makeCascade.timer = ROOT.TStopwatch() |
| |
| argparse | makeCascade.ap = argparse.ArgumentParser(description="Run SHiP makeCascade: generate ccbar or bbar") |
| |
| | makeCascade.type |
| |
| | makeCascade.default |
| |
| | makeCascade.help |
| |
| | makeCascade.int |
| |
| | makeCascade.float |
| |
| | makeCascade.choices |
| |
| | makeCascade.action |
| |
| argparse | makeCascade.args = ap.parse_args() |
| |
| | makeCascade.Fntuple |
| |
| list | makeCascade.idbeam = [2212, 211, 2112, 321, 130, 310] |
| |
| list | makeCascade.target = ["p+", "n0"] |
| |
| float | makeCascade.fracp = 0.40 |
| |
| | makeCascade.else : |
| |
| float | makeCascade.pbeaml = 34.0 |
| |
| list | makeCascade.idsig = [411, 421, 431, 4122, 4132, 4232, 4332, 4412, 4414, 4422, 4424, 4432, 4434, 4444] |
| |
| ROOT | makeCascade.PDG = ROOT.TDatabasePDG.Instance() |
| |
| ROOT | makeCascade.myPythia = ROOT.TPythia6() |
| |
| ROOT | makeCascade.tp = ROOT.tPythia6Generator() |
| |
| dict | makeCascade.h = {} |
| |
| int | makeCascade.id = 0 |
| |
| int | makeCascade.nb = 400 |
| |
| time | makeCascade.t0 = time.time() |
| |
| len | makeCascade.idhist = len(idbeam) * 4 * [0] |
| |
| list | makeCascade.idw = idbeam[idp] * idpm |
| |
| ROOT | makeCascade.name = PDG.GetParticle(idw).GetName() |
| |
| int | makeCascade.idb = id * 10 + idnp * 4 |
| |
| idpn | makeCascade.idadd = idpn * 4 |
| |
| ipbeam | makeCascade.pbw = ipbeam * (args.pbeamh - pbeaml) / (args.nrpoints - 1) + pbeaml |
| |
| dict | makeCascade.ibin = h[str(id * 10 + 1 + idadd)].FindBin(pbw, 0.0, 0.0) |
| |
| float | makeCascade.chimx = 0.0 |
| |
| dict | makeCascade.ibh = h[str(idw)].FindBin(args.pbeamh) |
| |
| ROOT | makeCascade.ftup = ROOT.TFile.Open(args.Fntuple, "RECREATE") |
| |
| ROOT | makeCascade.Ntup |
| |
| ROOT | makeCascade.kc = myPythia.Pycomp(kf) |
| |
| int | makeCascade.stack = 1000 * [0] |
| |
| int | makeCascade.nstack = 0 |
| |
| ROOT | makeCascade.ptot = ROOT.TMath.Sqrt(stack[nstack][1] ** 2 + stack[nstack][2] ** 2 + stack[nstack][3] ** 2) |
| |
| float | makeCascade.prbsig = 0.0 |
| |
| int | makeCascade.idpn = 0 |
| |
| dict | makeCascade.ib = h[str(idw)].FindBin(ptot, 0.0, 0.0) |
| |
| list | makeCascade.charmFound = [] |
| |
| ROOT | makeCascade.idabs = ROOT.TMath.Abs(myPythia.GetK(itrk, 2)) |
| |
| array | makeCascade.vl = array("f") |
| |
| int | makeCascade.nsub = stack[nstack][4] - 1 |
| |
| ROOT | makeCascade.pt2 = myPythia.GetP(itrk, 1) ** 2 + myPythia.GetP(itrk, 2) ** 2 |
| |
| argparse | makeCascade.beta = args.pbeamh / (myPythia.GetP(1, 4) + myPythia.GetP(2, 5)) |
| |
| tuple | makeCascade.gamma = (1 - beta**2) ** (-0.5) |
| |
| tuple | makeCascade.pbcm = -gamma * beta * myPythia.GetP(1, 4) + gamma * myPythia.GetP(1, 3) |
| |
| tuple | makeCascade.pDcm = -gamma * beta * myPythia.GetP(itrk, 4) + gamma * myPythia.GetP(itrk, 3) |
| |
| tuple | makeCascade.xf = pDcm / pbcm |
| |
| int | makeCascade.icas = stack[nstack][4] + 1 |
| |
| copy | makeCascade.anclist = copy.copy(stack[nstack][5]) |
| |
| copy | makeCascade.sublist = copy.copy(stack[nstack][6]) |
| |
| copy | makeCascade.tmp = copy.copy(anclist) |
| |
| copy | makeCascade.stmp = copy.copy(sublist) |
| |
| ROOT | makeCascade.rtime = timer.RealTime() |
| |
| ROOT | makeCascade.ctime = timer.CpuTime() |
| |
1.9.6