|
| | makeCascade.timer = ROOT.TStopwatch() |
| |
| str | makeCascade.R = '' |
| |
| int | makeCascade.mselcb = 4 |
| |
| int | makeCascade.pbeamh = 400. |
| |
| bool | makeCascade.storePrimaries = False |
| |
| int | makeCascade.nevgen = 100000 |
| |
| str | makeCascade.Fntuple = 'Cascade100k-parp16-MSTP82-1-MSEL'+str(mselcb)+'-ntuple.root' |
| |
| | makeCascade.opts |
| |
| | makeCascade.args |
| |
| int | makeCascade.nev = 5000 |
| |
| int | makeCascade.nrpoints = 20 |
| |
| list | makeCascade.idbeam = [2212,211,2112,321,130,310] |
| |
| list | makeCascade.target = ['p+','n0'] |
| |
| float | makeCascade.fracp = 0.43 |
| |
| int | makeCascade.pbeaml = 34. |
| |
| list | makeCascade.idsig = [411, 421, 431,4122,4132,4232,4332,4412,4414,4422,4424,4432,4434,4444] |
| |
| | makeCascade.PDG = ROOT.TDatabasePDG.Instance() |
| |
| | makeCascade.myPythia = ROOT.TPythia6() |
| |
| | makeCascade.tp = ROOT.tPythia6Generator() |
| |
| dict | makeCascade.h = {} |
| |
| int | makeCascade.id = 0 |
| |
| int | makeCascade.nb = 400 |
| |
| | makeCascade.t0 = time.time() |
| |
| int | makeCascade.idhist = len(idbeam)*4*[0] |
| |
| list | makeCascade.idw = idbeam[idp]*idpm |
| |
| | makeCascade.name = PDG.GetParticle(idw).GetName() |
| |
| int | makeCascade.idb = id*10+idnp*4 |
| |
| int | makeCascade.idadd = idpn*4 |
| |
| int | makeCascade.pbw = ipbeam*(pbeamh-pbeaml)/(nrpoints-1)+pbeaml |
| |
| dict | makeCascade.ibin = h[str(id*10+1+idadd)].FindBin(pbw,0.,0.) |
| |
| int | makeCascade.chimx = 0. |
| |
| dict | makeCascade.ibh = h[str(idw)].FindBin(pbeamh) |
| |
| | makeCascade.ftup = ROOT.TFile.Open(Fntuple, 'RECREATE') |
| |
| | makeCascade.Ntup |
| |
| | makeCascade.kc = myPythia.Pycomp(kf) |
| |
| int | makeCascade.stack = 1000*[0] |
| |
| int | makeCascade.nstack = 0 |
| |
| | makeCascade.ptot = ROOT.TMath.Sqrt(stack[nstack][1]**2+stack[nstack][2]**2+stack[nstack][3]**2) |
| |
| int | makeCascade.prbsig = 0. |
| |
| int | makeCascade.idpn = 0 |
| |
| dict | makeCascade.ib = h[str(idw)].FindBin(ptot,0.,0.) |
| |
| list | makeCascade.charmFound = [] |
| |
| | makeCascade.idabs = ROOT.TMath.Abs(myPythia.GetK(itrk,2)) |
| |
| | makeCascade.vl = array('f') |
| |
| int | makeCascade.nsub = stack[nstack][4]-1 |
| |
| int | makeCascade.pt2 = myPythia.GetP(itrk,1)**2+myPythia.GetP(itrk,2)**2 |
| |
| int | makeCascade.beta = 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 |
| |
| | makeCascade.anclist = copy.copy(stack[nstack][5]) |
| |
| | makeCascade.sublist = copy.copy(stack[nstack][6]) |
| |
| | makeCascade.tmp = copy.copy(anclist) |
| |
| | makeCascade.stmp = copy.copy(sublist) |
| |
| | makeCascade.rtime = timer.RealTime() |
| |
| | makeCascade.ctime = timer.CpuTime() |
| |
1.9.8