dpProductionRates Namespace Reference

Functions
	isDP (pdg)
	pbremProdRateVDM (mass, epsilon, doprint=True)
	pbremProdRateDipole (mass, epsilon, doprint=False)
	getAverageMesonRate (mumPdg)
	mesonBRtoPhoton (mumPdg, doprint=False)
	brMesonToGammaDP (mass, epsilon, mumPdg, doprint=False)
	brMesonToMesonDP (mass, epsilon, mumPdg, dauPdg, doprint=False)
	brMesonToDP (mass, epsilon, mumPdg, doprint=False)
	mesonProdRate (mass, epsilon, mumPdg, doprint=False)
	qcdprodRate (mass, epsilon, doprint=False)
	getDPprodRate (mass, epsilon, prodMode, mumPdg, doprint=False)

Variables
	PDG = ROOT.TDatabasePDG.Instance()
int	protonFlux = 2e20

Function Documentation

brMesonToDP()

dpProductionRates.brMesonToDP	(		mass,
			epsilon,
			mumPdg,
			doprint = False
	)

Definition at line 74 of file dpProductionRates.py.

def brMesonToDP(mass,epsilon,mumPdg,doprint=False):
    if mumPdg==223: return brMesonToMesonDP(mass,epsilon,mumPdg,111,doprint)
    elif (mumPdg==111 or mumPdg==221): return brMesonToGammaDP(mass,epsilon,mumPdg,doprint)
    elif mumPdg==331: return brMesonToMesonDP(mass,epsilon,mumPdg,113,doprint),brMesonToGammaDP(mass,epsilon,mumPdg,doprint)
    else: 
        print("Warning! Unknown mother pdgId %d, not implemented. Setting br to 0."%mumPdg)
        return 1
 

brMesonToGammaDP()

dpProductionRates.brMesonToGammaDP	(		mass,
			epsilon,
			mumPdg,
			doprint = False
	)

Definition at line 52 of file dpProductionRates.py.

def brMesonToGammaDP(mass,epsilon,mumPdg,doprint=False):
    mMeson = PDG.GetParticle(mumPdg).Mass()
    if (doprint==True): print("Mass of mother %d meson is %3.3f"%(mumPdg,mMeson))
    if (mass<mMeson): br = 2*epsilon**2*pow((1-mass**2/mMeson**2),3)*mesonBRtoPhoton(mumPdg,doprint)
    else: br = 0
    if (doprint==True): print("Branching ratio of %d meson to DP is %.8g"%(mumPdg,br))
    return br
 

brMesonToMesonDP()

dpProductionRates.brMesonToMesonDP	(		mass,
			epsilon,
			mumPdg,
			dauPdg,
			doprint = False
	)

Definition at line 60 of file dpProductionRates.py.

def brMesonToMesonDP(mass,epsilon,mumPdg,dauPdg,doprint=False):
    mMeson = PDG.GetParticle(mumPdg).Mass()
    mDaughterMeson = PDG.GetParticle(dauPdg).Mass()
    if (doprint==True): print("Mass of mother %d meson is %3.3f"%(mumPdg,mMeson))
    if (doprint==True): print("Mass of daughter %d meson is %3.3f"%(dauPdg,mDaughterMeson))
    if (mass<(mMeson-mDaughterMeson)):
        fac1 = pow(mMeson**2.-mDaughterMeson**2.,-3.)
        fac2 = pow((mass**2.-(mMeson+mDaughterMeson)**2.)*(mass**2.-(mMeson-mDaughterMeson)**2.),1.5)
        massfactor = fac1*fac2
        br = (epsilon**2.)*mesonBRtoPhoton(mumPdg,doprint)*massfactor
    else: br = 0
    if (doprint==True): print("Branching ratio of %d meson to DP is %.8g"%(mumPdg,br))
    return br
 

getAverageMesonRate()

dpProductionRates.getAverageMesonRate

(

mumPdg

)

Definition at line 34 of file dpProductionRates.py.

def getAverageMesonRate(mumPdg):
    if (mumPdg==111): return 6.166
    if (mumPdg==221): return 0.7012
    if (mumPdg==223): return 0.8295
    if (mumPdg==331): return 0.07825
    print(" -- ERROR, unknown mother pdgId %d"%mumPdg)
    return 0
 
#from the PDG, decay to photon channels available for mixing with DP

getDPprodRate()

dpProductionRates.getDPprodRate	(		mass,
			epsilon,
			prodMode,
			mumPdg,
			doprint = False
	)

Definition at line 103 of file dpProductionRates.py.

def getDPprodRate(mass,epsilon,prodMode,mumPdg,doprint=False):
    if ('pbrem' in prodMode):
        print("VDM")
        return pbremProdRateVDM(mass,epsilon,doprint)
    elif ('pbrem1' in prodMode):
        print("Dipole")
        return pbremProdRateDipole(mass,epsilon,doprint)
    elif ('meson' in prodMode):
        return mesonProdRate(mass,epsilon,mumPdg,doprint)
    elif ('qcd' in prodMode):
        return qcdprodRate(mass,epsilon,doprint)
    else:
        print("Unknown production mode! Choose among pbrem, meson or qcd.")
        return 1

isDP()

dpProductionRates.isDP

(

pdg

)

Definition at line 12 of file dpProductionRates.py.

def isDP(pdg):
    if (pdg==9900015 or pdg==4900023): 
        return True
    return False
 

mesonBRtoPhoton()

dpProductionRates.mesonBRtoPhoton	(		mumPdg,
			doprint = False
	)

Definition at line 43 of file dpProductionRates.py.

def mesonBRtoPhoton(mumPdg,doprint=False):
    br = 1
    if (mumPdg==111): br = 0.9879900
    if (mumPdg==221): br = 0.3931181
    if (mumPdg==223): br = 0.0834941
    if (mumPdg==331): br = 0.0219297
    if (doprint==True): print("BR of %d meson to photons: %.8g"%(mumPdg,br))
    return br 
 

mesonProdRate()

dpProductionRates.mesonProdRate	(		mass,
			epsilon,
			mumPdg,
			doprint = False
	)

Definition at line 82 of file dpProductionRates.py.

def mesonProdRate(mass,epsilon,mumPdg,doprint=False):
    brM2DP=brMesonToDP(mass,epsilon,mumPdg,doprint)
    if mumPdg==331:
        avgMeson = getAverageMesonRate(mumPdg)*brM2DP[0] 
        avgMeson1 = getAverageMesonRate(mumPdg)*brM2DP[1]
        return avgMeson*0.6, avgMeson1*0.6
        #return avgMeson + avgMeson1
    if not mumPdg==331:
        avgMeson = getAverageMesonRate(mumPdg)*brM2DP
        return avgMeson*0.6
 
#from interpolation of Pythia XS, normalised to epsilon^2

pbremProdRateDipole()

dpProductionRates.pbremProdRateDipole	(		mass,
			epsilon,
			doprint = False
	)

Definition at line 25 of file dpProductionRates.py.

def pbremProdRateDipole(mass,epsilon,doprint=False):
    xswg = proton_bremsstrahlung.prodRate(mass, epsilon)
    if doprint: print("A' production rate per p.o.t: \t %.8g"%(xswg))
    penalty = proton_bremsstrahlung.penaltyFactor(mass)
    if doprint: print("A' penalty factor: \t %.8g"%penalty)
    if doprint: print("A' rescaled production rate per p.o.t:\t %.8g"%(xswg*penalty))
    return xswg*penalty
 
#obtained with Pythia8: average number of meson expected per p.o.t from inclusive pp to X production, 100k events produced

pbremProdRateVDM()

dpProductionRates.pbremProdRateVDM	(		mass,
			epsilon,
			doprint = True
	)

Definition at line 17 of file dpProductionRates.py.

def pbremProdRateVDM(mass,epsilon,doprint=True):
    xswg = proton_bremsstrahlung.prodRate(mass, epsilon)
    if doprint: print("A' production rate per p.o.t: \t %.8g"%(xswg))
    rhoff = proton_bremsstrahlung.rhoFormFactor(mass)**2
    if doprint: print("A' rho form factor: \t %.8g"%rhoff)
    if doprint: print("A' rescaled production rate per p.o.t:\t %.8g"%(xswg*rhoff))
    return xswg*rhoff
 

qcdprodRate()

dpProductionRates.qcdprodRate	(		mass,
			epsilon,
			doprint = False
	)

Definition at line 94 of file dpProductionRates.py.

def qcdprodRate(mass,epsilon,doprint=False):
    if (mass > 3.):
        xs = math.exp(-5.928-0.8669*mass)
    elif (mass > 1.4):
        xs = math.exp(-4.1477-1.4745*mass)
    else:
        xs = 0.
    return xs*epsilon*epsilon
 

Variable Documentation

PDG

dpProductionRates.PDG = ROOT.TDatabasePDG.Instance()

Definition at line 8 of file dpProductionRates.py.

protonFlux

int dpProductionRates.protonFlux = 2e20

Definition at line 9 of file dpProductionRates.py.