Functions
def	yBeam (float Mproton=0.938272081, float pbeam=400.0)

None	na50 (bool online=True)

None	muflux ()

None	debugging (g)

Variables
int	theSeed = 0

dict	h = {}

ArgumentParser	parser = ArgumentParser()

	dest

	type

	int

	help

	default

	float

	action

	str

ArgumentParser	options = parser.parse_args()

ROOT	X = ROOT.FixedTargetGenerator()

dict	generators = {"p": ROOT.Pythia8.Pythia(), "n": ROOT.Pythia8.Pythia()}

	else :

dict	rc = generators["p"].next()

dict	processes = generators["p"].info.codesHard()

tuple	hname

ROOT	f = ROOT.TFile("ntuple-" + hname + ".root", "RECREATE")

ROOT	signal = ROOT.TNtuple("ntuple", "ntuple", "M:P:Pt:y:p1x:p1y:p1z:p2x:p2y:p2z:cosCS")

ROOT	timer = ROOT.TStopwatch()

dict	ntagged = {"p": 0, "n": 0}

def	ybeam = yBeam(pbeam=options.fMom)

dict	py = generators[g]

dict	nmu = {}

list	ks = list(nmu)

dict	Zstar = py.event[nmu[ks[0]]]

dict	ii = nmu[ks[0]]

dict	d0 = py.event.daughterList(ii)[0]

dict	d1 = py.event.daughterList(ii)[1]

dict	nlep = py.event[d0]

dict	nantilep = py.event[d1]

dict	P1pl = nlep.e() + nlep.pz()

dict	P2pl = nantilep.e() + nantilep.pz()

dict	P1mi = nlep.e() - nlep.pz()

dict	P2mi = nantilep.e() - nantilep.pz()

dict	A = P1pl * P2mi - P2pl * P1mi

tuple	cosCS

dict	M = {}

int	k = 0

dict	G = M[0] + M[1]

ROOT	rtime = timer.RealTime()

ROOT	ctime = timer.CpuTime()

dict	sigma = generators[g].info.sigmaGen(processes[0])

Function Documentation

◆ debugging()

None runPythia8.debugging ( g )

Definition at line 287 of file runPythia8.py.

287def debugging(g) -> None:

288 generators[g].settings.listAll()

◆ muflux()

None runPythia8.muflux ( )

Definition at line 258 of file runPythia8.py.

def muflux() -> None:
    Z_Mo = 96.0
    P_Mo = 42
    fraction = {}
    for g in generators:
        processes = generators[g].info.codesHard()
        name = generators[g].info.nameProc(processes[0])
        sigma = generators[g].info.sigmaGen(processes[0])
        yax = h["M_" + g].GetYaxis()
        xax = h["M_" + g].GetXaxis()
        Mmin = xax.FindBin(0.0)
        Mmax = xax.FindBin(5.0)
        Ymin = yax.FindBin(0.3)
        Ymax = yax.FindBin(3.0)
        h["MA"] = h["M_" + g].ProjectionX("MA")
        h["M"] = h["M_" + g].ProjectionX("M", Ymin, Ymax)
        print(
            g,
            name,
            sigma,
            float(h["MA"].GetEntries()) / options.NPoT,
            h["MA"].Integral(Mmin, Mmax) / float(h["MA"].GetEntries()),
            h["M"].GetEntries() / float(h["MA"].GetEntries()),
        )
        fraction[g] = h["M"].Integral(Mmin, Mmax) / options.NPoT
    meanFraction = (fraction["p"] * P_Mo + fraction["n"] * (Z_Mo - P_Mo)) / Z_Mo * sigma
    print("cross section a la muflux: %5.2F pb" % (0.5 * meanFraction * 1e9))
 
 

◆ na50()

None runPythia8.na50 ( bool online = True )

Definition at line 192 of file runPythia8.py.

def na50(online: bool = True) -> None:
    for g in generators:
        if online:
            processes = generators[g].info.codesHard()
            name = generators[g].info.nameProc(processes[0])
            sigma = generators[g].info.sigmaGen(processes[0])
        else:
            name = ""
            sigma = h["xsec_" + g].GetBinContent(1)
        yax = h["M_" + g].GetYaxis()
        xax = h["M_" + g].GetXaxis()
        Mmin = xax.FindBin(2.9)
        Mmax = xax.FindBin(4.5)
        Ymin = yax.FindBin(-0.425)
        Ymax = yax.FindBin(0.575)
        h["MA"] = h["M_" + g].ProjectionX("MA")
        h["M"] = h["M_" + g].ProjectionX("M", Ymin, Ymax)
        print("generator     sigma   mumu-ratio  in-mass-range  in-y-range")
        print(
            "%s %s %6.2F nbarn, %5.2F, %5.2G, %5.2F    "
            % (
                g,
                name,
                sigma * 1e6,
                float(h["MA"].GetEntries()) / options.NPoT,
                h["MA"].Integral(Mmin, Mmax) / float(h["MA"].GetEntries()),
                h["M"].GetEntries() / float(h["MA"].GetEntries()),
            )
        )
        fraction = h["M"].Integral(Mmin, Mmax) / options.NPoT
        # multiply with 0.5 assuming no polarization -0.5 < cosCS < 0.5
        print(f"cross section a la NA50 for : {g} {0.5 * fraction * sigma * 1e9:5.2F} pb")
    # via cosCS
    for g in generators:
        if online:
            processes = generators[g].info.codesHard()
            name = generators[g].info.nameProc(processes[0])
            sigma = generators[g].info.sigmaGen(processes[0])
        else:
            name = ""
            sigma = h["xsec_" + g].GetBinContent(1)
        yax = h["cosCSJpsi_" + g].GetYaxis()
        xax = h["cosCSJpsi_" + g].GetXaxis()
        Mmin = xax.FindBin(-0.5)
        Mmax = xax.FindBin(0.5)
        Ymin = yax.FindBin(-0.425)
        Ymax = yax.FindBin(0.575)
        h["MA"] = h["cosCSJpsi_" + g].ProjectionX("MA")
        h["M"] = h["cosCSJpsi_" + g].ProjectionX("M", Ymin, Ymax)
        print("generator     sigma   mumu-in-mass-range% cosCS in-y-range")
        print(
            "%s %s %6.2F nbarn, %5.2F, %5.2F, %5.2F   "
            % (
                g,
                name,
                sigma * 1e6,
                float(h["MA"].GetEntries()) / options.NPoT * 100.0,
                h["MA"].Integral(Mmin, Mmax) / float(h["MA"].GetEntries()),
                h["M"].GetEntries() / float(h["MA"].GetEntries()),
            )
        )
        fraction = h["M"].Integral(Mmin, Mmax) / options.NPoT
        # taking polarization into account.
        print("cross section a la NA50, -0.5<cosCS<0.5: %5.2F pb" % (fraction * sigma * 1e9))
 
 

◆ yBeam()

def runPythia8.yBeam	(	float	Mproton = `0.938272081`,
		float	pbeam = `400.0`
	)

Definition at line 43 of file runPythia8.py.

def yBeam(Mproton: float = 0.938272081, pbeam: float = 400.0):
    Ebeam = ROOT.TMath.Sqrt(pbeam**2 + Mproton**2)
    betaCM = pbeam / (Ebeam + Mproton)
    y_beam = 0.5 * ROOT.TMath.Log((1 + betaCM) / (1 - betaCM))  # https://arxiv.org/pdf/1604.02651.pdf
    return y_beam
 
 

Variable Documentation

◆ A

dict runPythia8.A = P1pl * P2mi - P2pl * P1mi

Definition at line 146 of file runPythia8.py.

◆ action

runPythia8.action

Definition at line 28 of file runPythia8.py.

◆ cosCS

tuple runPythia8.cosCS

Initial value:

=  (
                    Zstar.pz() / abs(Zstar.pz()) * 1.0 / Zstar.m() / ROOT.TMath.Sqrt(Zstar.m2() + Zstar.pT() ** 2) * A
                )

Definition at line 147 of file runPythia8.py.

◆ ctime

ROOT runPythia8.ctime = timer.CpuTime()

Definition at line 183 of file runPythia8.py.

◆ d0

dict runPythia8.d0 = py.event.daughterList(ii)[0]

Definition at line 134 of file runPythia8.py.

◆ d1

dict runPythia8.d1 = py.event.daughterList(ii)[1]

Definition at line 135 of file runPythia8.py.

◆ default

runPythia8.default

Definition at line 23 of file runPythia8.py.

◆ dest

runPythia8.dest

Definition at line 23 of file runPythia8.py.

◆ else

runPythia8.else :

Definition at line 87 of file runPythia8.py.

◆ f

ROOT runPythia8.f = ROOT.TFile("ntuple-" + hname + ".root", "RECREATE")

Definition at line 105 of file runPythia8.py.

◆ float

runPythia8.float

Definition at line 26 of file runPythia8.py.

◆ G

dict runPythia8.G = M[0] + M[1]

Definition at line 172 of file runPythia8.py.

◆ generators

dict runPythia8.generators = {"p": ROOT.Pythia8.Pythia(), "n": ROOT.Pythia8.Pythia()}

Definition at line 50 of file runPythia8.py.

◆ h

dict runPythia8.h = {}

Definition at line 11 of file runPythia8.py.

◆ help

runPythia8.help

Definition at line 23 of file runPythia8.py.

◆ hname

tuple runPythia8.hname

Initial value:

=  (
    "pythia8_PDFpset"
    + options.PDFpSet
    + "_Emin"
    + str(options.Emin)
    + "_"
    + generators["p"].info.nameProc(processes[0])
)

Definition at line 93 of file runPythia8.py.

◆ ii

dict runPythia8.ii = nmu[ks[0]]

Definition at line 133 of file runPythia8.py.

◆ int

runPythia8.int

Definition at line 23 of file runPythia8.py.

◆ k

int runPythia8.k = 0

Definition at line 168 of file runPythia8.py.

◆ ks

list runPythia8.ks = list(nmu)

Definition at line 128 of file runPythia8.py.

◆ M

dict runPythia8.M = {}

Definition at line 167 of file runPythia8.py.

◆ nantilep

dict runPythia8.nantilep = py.event[d1]

Definition at line 138 of file runPythia8.py.

◆ nlep

dict runPythia8.nlep = py.event[d0]

Definition at line 137 of file runPythia8.py.

◆ nmu

dict runPythia8.nmu = {}

Definition at line 117 of file runPythia8.py.

◆ ntagged

dict runPythia8.ntagged = {"p": 0, "n": 0}

Definition at line 111 of file runPythia8.py.

◆ options

ArgumentParser runPythia8.options = parser.parse_args()

Definition at line 37 of file runPythia8.py.

◆ P1mi

dict runPythia8.P1mi = nlep.e() - nlep.pz()

Definition at line 144 of file runPythia8.py.

◆ P1pl

dict runPythia8.P1pl = nlep.e() + nlep.pz()

Definition at line 142 of file runPythia8.py.

◆ P2mi

dict runPythia8.P2mi = nantilep.e() - nantilep.pz()

Definition at line 145 of file runPythia8.py.

◆ P2pl

dict runPythia8.P2pl = nantilep.e() + nantilep.pz()

Definition at line 143 of file runPythia8.py.

◆ parser

ArgumentParser runPythia8.parser = ArgumentParser()

Definition at line 22 of file runPythia8.py.

◆ processes

dict runPythia8.processes = generators["p"].info.codesHard()

Definition at line 92 of file runPythia8.py.

◆ py

dict runPythia8.py = generators[g]

Definition at line 115 of file runPythia8.py.

◆ rc

dict runPythia8.rc = generators["p"].next()

Definition at line 91 of file runPythia8.py.

◆ rtime

ROOT runPythia8.rtime = timer.RealTime()

Definition at line 182 of file runPythia8.py.

◆ sigma

dict runPythia8.sigma = generators[g].info.sigmaGen(processes[0])

Definition at line 187 of file runPythia8.py.

◆ signal

ROOT runPythia8.signal = ROOT.TNtuple("ntuple", "ntuple", "M:P:Pt:y:p1x:p1y:p1z:p2x:p2y:p2z:cosCS")

Definition at line 106 of file runPythia8.py.

◆ str

runPythia8.str

Definition at line 32 of file runPythia8.py.

◆ theSeed

int runPythia8.theSeed = 0

Definition at line 10 of file runPythia8.py.

◆ timer

ROOT runPythia8.timer = ROOT.TStopwatch()

Definition at line 108 of file runPythia8.py.

◆ type

runPythia8.type

Definition at line 23 of file runPythia8.py.

◆ X

ROOT runPythia8.X = ROOT.FixedTargetGenerator()

Definition at line 39 of file runPythia8.py.

◆ ybeam

def runPythia8.ybeam = yBeam(pbeam=options.fMom)

Definition at line 112 of file runPythia8.py.

◆ Zstar

dict runPythia8.Zstar = py.event[nmu[ks[0]]]

Definition at line 130 of file runPythia8.py.

Functions

Variables

Function Documentation

◆ debugging()

◆ muflux()

◆ na50()

◆ yBeam()

Variable Documentation

◆ A

◆ action

◆ cosCS

◆ ctime

◆ d0

◆ d1

◆ default

◆ dest

◆ else

◆ f

◆ float

◆ G

◆ generators

◆ h

◆ help

◆ hname

◆ ii

◆ int

◆ k

◆ ks

◆ M

◆ nantilep

◆ nlep

◆ nmu

◆ ntagged

◆ options

◆ P1mi

◆ P1pl

◆ P2mi

◆ P2pl

◆ parser

◆ processes

◆ py

◆ rc

◆ rtime

◆ sigma

◆ signal

◆ str

◆ theSeed

◆ timer

◆ type

◆ X

◆ ybeam

◆ Zstar