SND@LHC Software
Loading...
Searching...
No Matches
makeMuonEM.py File Reference

Go to the source code of this file.

Namespaces

namespace  makeMuonEM
 

Functions

 makeMuonEM.getMasssq (pid)
 

Variables

int makeMuonEM.nJob = 1
 
int makeMuonEM.nMult = 1000
 
str makeMuonEM.muonIn = '/media/Data/HNL/muVetoDIS/muDISVetoCounter.root'
 
 makeMuonEM.PDG = ROOT.TDatabasePDG.Instance()
 
dict makeMuonEM.masssq = {}
 
 makeMuonEM.fout = ROOT.TFile('muonEm_'+str(nJob)+'.root','recreate')
 
 makeMuonEM.dTree = ROOT.TNtuple("pythia8-Geant4","muons for EM studies","id:px:py:pz:x:y:z:ox:oy:oz:pythiaid:parentid:ecut:w")
 
 makeMuonEM.fin = ROOT.TFile(muonIn)
 
 makeMuonEM.sTree = fin.muons
 
 makeMuonEM.rc = sTree.GetEvent(k)
 
 makeMuonEM.px
 
 makeMuonEM.py
 
 makeMuonEM.pz
 
 makeMuonEM.x
 
 makeMuonEM.y
 
 makeMuonEM.z
 
 makeMuonEM.pid
 
 makeMuonEM.w
 
 makeMuonEM.p = ROOT.TMath.Sqrt(px*px+py*py+pz*pz)
 
 makeMuonEM.E = ROOT.TMath.Sqrt(getMasssq(pid)+p*p)
 
 makeMuonEM.mu = array('d',[pid,px,py,pz,E,x,y,z,w])
 
 makeMuonEM.muPart = ROOT.TVectorD(9,mu)
 
 makeMuonEM.m = array('d',[pid,px,py,pz,E])
 
 makeMuonEM.part = ROOT.TVectorD(5,m)
 
 makeMuonEM.nPart = dPart.GetEntries()