SND@LHC Software
Loading...
Searching...
No Matches
makeMuonEM Namespace Reference

Functions

 getMasssq (pid)
 

Variables

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

Function Documentation

◆ getMasssq()

makeMuonEM.getMasssq (   pid)

Definition at line 12 of file makeMuonEM.py.

12def getMasssq(pid):
13 apid = abs(int(pid))
14 if not apid in masssq:
15 masssq[apid] = PDG.GetParticle(apid).Mass()**2
16 return masssq[apid]
17
18# prepare muon input for FairShip/Geant4 processing
19# incoming muon, id:px:py:pz:x:y:z:ox:oy:oz:pythiaid:parentid:ecut:w
20
21# just duplicate muon n times, rather stupid job
22

Variable Documentation

◆ dTree

makeMuonEM.dTree = ROOT.TNtuple("pythia8-Geant4","muons for EM studies","id:px:py:pz:x:y:z:ox:oy:oz:pythiaid:parentid:ecut:w")

Definition at line 24 of file makeMuonEM.py.

◆ E

makeMuonEM.E = ROOT.TMath.Sqrt(getMasssq(pid)+p*p)

Definition at line 37 of file makeMuonEM.py.

◆ fin

makeMuonEM.fin = ROOT.TFile(muonIn)

Definition at line 27 of file makeMuonEM.py.

◆ fout

makeMuonEM.fout = ROOT.TFile('muonEm_'+str(nJob)+'.root','recreate')

Definition at line 23 of file makeMuonEM.py.

◆ m

makeMuonEM.m = array('d',[pid,px,py,pz,E])

Definition at line 44 of file makeMuonEM.py.

◆ masssq

dict makeMuonEM.masssq = {}

Definition at line 10 of file makeMuonEM.py.

◆ mu

makeMuonEM.mu = array('d',[pid,px,py,pz,E,x,y,z,w])

Definition at line 38 of file makeMuonEM.py.

◆ muonIn

str makeMuonEM.muonIn = '/media/Data/HNL/muVetoDIS/muDISVetoCounter.root'

Definition at line 5 of file makeMuonEM.py.

◆ muPart

makeMuonEM.muPart = ROOT.TVectorD(9,mu)

Definition at line 39 of file makeMuonEM.py.

◆ nJob

int makeMuonEM.nJob = 1

Definition at line 3 of file makeMuonEM.py.

◆ nMult

int makeMuonEM.nMult = 1000

Definition at line 4 of file makeMuonEM.py.

◆ nPart

makeMuonEM.nPart = dPart.GetEntries()

Definition at line 47 of file makeMuonEM.py.

◆ p

makeMuonEM.p = ROOT.TMath.Sqrt(px*px+py*py+pz*pz)

Definition at line 36 of file makeMuonEM.py.

◆ part

makeMuonEM.part = ROOT.TVectorD(5,m)

Definition at line 45 of file makeMuonEM.py.

◆ PDG

makeMuonEM.PDG = ROOT.TDatabasePDG.Instance()

Definition at line 9 of file makeMuonEM.py.

◆ pid

makeMuonEM.pid

Definition at line 35 of file makeMuonEM.py.

◆ px

makeMuonEM.px

Definition at line 33 of file makeMuonEM.py.

◆ py

makeMuonEM.py

Definition at line 33 of file makeMuonEM.py.

◆ pz

makeMuonEM.pz

Definition at line 33 of file makeMuonEM.py.

◆ rc

makeMuonEM.rc = sTree.GetEvent(k)

Definition at line 31 of file makeMuonEM.py.

◆ sTree

makeMuonEM.sTree = fin.muons

Definition at line 28 of file makeMuonEM.py.

◆ w

makeMuonEM.w

Definition at line 35 of file makeMuonEM.py.

◆ x

makeMuonEM.x

Definition at line 34 of file makeMuonEM.py.

◆ y

makeMuonEM.y

Definition at line 34 of file makeMuonEM.py.

◆ z

makeMuonEM.z

Definition at line 34 of file makeMuonEM.py.