makeMuonDIS.py File Reference

Go to the source code of this file.

Namespaces
namespace	makeMuonDIS

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