makeMuonDIS.py File Reference

Go to the source code of this file.

Namespaces
namespace	makeMuonDIS

Functions
def	makeMuonDIS.getMasssq (pid)
def	makeMuonDIS.rotate (ctheta, stheta, cphi, sphi, px, py, pz)

Variables
int	makeMuonDIS.nJob = 2
int	makeMuonDIS.nMult = 10
str	makeMuonDIS.muonIn = "$SHIPSOFT/data/muConcrete.root"
int	makeMuonDIS.nPerJob = 20000
ROOT	makeMuonDIS.myPythia = ROOT.TPythia6()
ROOT	makeMuonDIS.kc = myPythia.Pycomp(kf)
dict	makeMuonDIS.masssq = {}
int	makeMuonDIS.R = int(time.time() % 900000000)
dict	makeMuonDIS.mutype = {-13: "gamma/mu+", 13: "gamma/mu-"}
ROOT	makeMuonDIS.fout = ROOT.TFile("muonDis_" + str(nJob) + ".root", "recreate")
ROOT	makeMuonDIS.dTree = ROOT.TTree("DIS", "muon DIS")
ROOT	makeMuonDIS.iMuon = ROOT.TClonesArray("TVectorD")
ROOT	makeMuonDIS.iMuonBranch = dTree.Branch("InMuon", iMuon, 32000, -1)
ROOT	makeMuonDIS.dPart = ROOT.TClonesArray("TVectorD")
ROOT	makeMuonDIS.dPartBranch = dTree.Branch("Particles", dPart, 32000, -1)
ROOT	makeMuonDIS.fin = ROOT.TFile(muonIn)
ROOT	makeMuonDIS.sTree = fin.muons
ROOT	makeMuonDIS.nTOT = sTree.GetEntries()
int	makeMuonDIS.nStart = nPerJob * nJob
min	makeMuonDIS.nEnd = min(nTOT, nStart + nPerJob)
int	makeMuonDIS.nMade = 0
ROOT	makeMuonDIS.rc = sTree.GetEvent(k)
	makeMuonDIS.px
	makeMuonDIS.py
	makeMuonDIS.pz
	makeMuonDIS.x
	makeMuonDIS.y
	makeMuonDIS.z
	makeMuonDIS.pid
	makeMuonDIS.w
ROOT	makeMuonDIS.p = ROOT.TMath.Sqrt(px * px + py * py + pz * pz)
ROOT	makeMuonDIS.E = ROOT.TMath.Sqrt(getMasssq(pid) + p * p)
ROOT	makeMuonDIS.theta = ROOT.TMath.ACos(pz / p)
ROOT	makeMuonDIS.phi = ROOT.TMath.ATan2(py, px)
	makeMuonDIS.ctheta
	makeMuonDIS.stheta
	makeMuonDIS.cphi
	makeMuonDIS.sphi
array	makeMuonDIS.mu = array("d", [pid, px, py, pz, E, x, y, z, w])
ROOT	makeMuonDIS.tca_vec = iMuon.ConstructedAt(0)
ROOT	makeMuonDIS.muPart = ROOT.TVectorD(9, mu)
ROOT	makeMuonDIS.did = myPythia.GetK(itrk, 2)
	makeMuonDIS.dpx
	makeMuonDIS.dpy
	makeMuonDIS.dpz
dpx	makeMuonDIS.psq = dpx**2 + dpy**2 + dpz**2
array	makeMuonDIS.m = array("d", [did, dpx, dpy, dpz, E])
ROOT	makeMuonDIS.part = ROOT.TVectorD(5, m)
len	makeMuonDIS.nPart = len(dPart)