FairShip
Loading...
Searching...
No Matches
makeMuonEM.py File Reference

Go to the source code of this file.

Namespaces

namespace  makeMuonEM
 

Functions

def makeMuonEM.getMasssq (pid)
 

Variables

int makeMuonEM.nJob = 1
 
int makeMuonEM.nMult = 1000
 
str makeMuonEM.muonIn = "/media/Data/HNL/muVetoDIS/muDISVetoCounter.root"
 
ROOT makeMuonEM.PDG = ROOT.TDatabasePDG.Instance()
 
dict makeMuonEM.masssq = {}
 
ROOT makeMuonEM.fout = ROOT.TFile("muonEm_" + str(nJob) + ".root", "recreate")
 
ROOT makeMuonEM.dTree = ROOT.TTree("pythia8-Geant4", "muons for EM studies")
 
ROOT makeMuonEM.iMuon = ROOT.TClonesArray("TVectorD")
 
ROOT makeMuonEM.iMuonBranch = dTree.Branch("InMuon", iMuon, 32000, -1)
 
ROOT makeMuonEM.dPart = ROOT.TClonesArray("TVectorD")
 
ROOT makeMuonEM.dPartBranch = dTree.Branch("Particles", dPart, 32000, -1)
 
ROOT makeMuonEM.fin = ROOT.TFile(muonIn)
 
ROOT makeMuonEM.sTree = fin.muons
 
ROOT makeMuonEM.rc = sTree.GetEvent(k)
 
 makeMuonEM.px
 
 makeMuonEM.py
 
 makeMuonEM.pz
 
 makeMuonEM.x
 
 makeMuonEM.y
 
 makeMuonEM.z
 
 makeMuonEM.pid
 
 makeMuonEM.w
 
ROOT makeMuonEM.p = ROOT.TMath.Sqrt(px * px + py * py + pz * pz)
 
ROOT makeMuonEM.E = ROOT.TMath.Sqrt(getMasssq(pid) + p * p)
 
array makeMuonEM.mu = array("d", [pid, px, py, pz, E, x, y, z, w])
 
ROOT makeMuonEM.muPart = ROOT.TVectorD(9, mu)
 
ROOT makeMuonEM.tca_vec = iMuon.ConstructedAt(0)
 
array makeMuonEM.m = array("d", [pid, px, py, pz, E])
 
ROOT makeMuonEM.part = ROOT.TVectorD(5, m)
 
len makeMuonEM.nPart = len(dPart)