makeDecay.py File Reference

Go to the source code of this file.

Namespaces
namespace	makeDecay

Variables
str	makeDecay.fname = "/eos/experiment/ship/data/Charm/Cascade-parp16-MSTP82-1-MSEL4-76Mpot_1"
float	makeDecay.nrpotspill = 5.0e13
float	makeDecay.chicc = 1.7e-3
float	makeDecay.chibb = 1.6e-7
bool	makeDecay.setByHand = False
	makeDecay.try :
	makeDecay.opts
	makeDecay.args
str	makeDecay.FIN = fname + ".root"
os	makeDecay.tmp = os.path.abspath(FIN).split("/")
str	makeDecay.FOUT = "Decay-" + tmp[len(tmp) - 1]
ROOT	makeDecay.fin = ROOT.TFile(FIN)
	makeDecay.else :
ROOT	makeDecay.sTree = fin.FindObjectAny("pythia6")
ROOT	makeDecay.nEvents = sTree.GetEntries()
dict	makeDecay.hc = {}
ROOT	makeDecay.fhin = ROOT.TFile(FIN.replace("ntuple", "hists"))
dict	makeDecay.nrcpot = hc["2"].GetBinContent(1) / 2.0
ROOT	makeDecay.P8gen = ROOT.TPythia8()
ROOT	makeDecay.P8 = P8gen.Pythia8()
int	makeDecay.n = 1
ROOT	makeDecay.p = P8.particleData.particleDataEntryPtr(n)
str	makeDecay.command = str(n) + ":mayDecay = false"
ROOT	makeDecay.ftup = ROOT.TFile.Open(FOUT, "RECREATE")
ROOT	makeDecay.Ntup = ROOT.TNtuple("Decay", "pythia8 heavy flavour decays", "id:px:py:pz:E:M:weight:mid:mpx:mpy:mpz:mE:pot:ptGM:pzGM")
dict	makeDecay.h = {}
ROOT	makeDecay.PDG = ROOT.TDatabasePDG.Instance()
int	makeDecay.idhnu = 1000 + idnu
idnu	makeDecay.idw = idnu
ROOT	makeDecay.name = PDG.GetParticle(idw).GetName()
float	makeDecay.pot = 0.0
int	makeDecay.nDsprim = 0
int	makeDecay.ntotprim = 0
ROOT	makeDecay.rc = sTree.GetEvent(n)
float	makeDecay.wspill = nrpotspill * chicc / nrcpot
ROOT	makeDecay.pt = ROOT.TMath.Sqrt(sTree.mpx**2 + sTree.mpy**2)
int	makeDecay.idabs = int(abs(sTree.id))
ROOT	makeDecay.par = P8.event[n]
ROOT	makeDecay.ptGM = ROOT.TMath.Sqrt(sTree.mpx * sTree.mpx + sTree.mpy * sTree.mpy)
ROOT	makeDecay.pt2 = par.px() ** 2 + par.py() ** 2
ROOT	makeDecay.ptot = ROOT.TMath.Sqrt(pt2 + par.pz() ** 2)
min	makeDecay.l10ptot = min(max(ROOT.TMath.Log10(ptot), -0.3), 1.69999)
min	makeDecay.l10pt = min(max(ROOT.TMath.Log10(ROOT.TMath.Sqrt(pt2)), -2.0), 0.4999)
float	makeDecay.fDsP6 = 1.0 * nDsprim / ntotprim
float	makeDecay.fDs = 0.077