def pbremProdRateNoFF(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
def pbremProdRate(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
0., 10.)
h['xsMeson'].Sumw2()
ut.bookHist(h, 'xsQCD',
';m_{#gamma^{D}} (GeV); #sigma_{qcd}/#varepsilon^2 (mb)', 1000, 0.,
10.)
h['xsMeson'].Sumw2()
print proton_bremsstrahlung.rhoFormFactor(
0), proton_bremsstrahlung.rhoFormFactor(
0.77), proton_bremsstrahlung.rhoFormFactor(2.)
for step in range(0, 1000):
mass = step * 0.0100000
xswg = proton_bremsstrahlung.prodRate(400, mass, epsilon) * 10.7
penalty = proton_bremsstrahlung.penaltyFactor(mass)
rhoff = proton_bremsstrahlung.rhoFormFactor(mass)
#if (rhoff<1): rhoff = 1
h['xsPbrem'].SetBinContent(step + 1, xswg * penalty / (epsilon**2))
h['xsPbremNew'].SetBinContent(step + 1,
xswg * penalty * rhoff**2 / (epsilon**2))
h['xsPbremFF'].SetBinContent(step + 1, xswg * rhoff**2 / (epsilon**2))
for step in range(0, 900):
mass = math.sqrt(step * 0.0100000)
rhoff = proton_bremsstrahlung.rhoFormFactor(mass)
h['ff'].SetBinContent(step + 1, rhoff)
for step in range(1, 1000):