예제 #1
0
def tanbForXsecLow(xSecAtLimit,mHp,tanbRef,mu):

    tanb = tanbRef
        
    accuracy = 0.001
    
    counter = 0
    xSec = crosssection.whTauNuCrossSectionMSSM(mHp,tanb,mu)
    while(abs(xSecAtLimit - xSec)/xSecAtLimit > accuracy and xSec > 0 and tanb > 1.1 ):
        tanb = tanb - 0.1*(xSecAtLimit - xSec)/xSecAtLimit*tanb
        xSec = crosssection.whTauNuCrossSectionMSSM(mHp,tanb,mu)
        counter = counter+1
        if counter > 10: # to prevent infinite loops
            counter = 0
            accuracy = accuracy*2
    if tanb < 1.1 or xSec <= 0:
        return -1
    return tanb
def tanbForXsecLow(xSecAtLimit, mHp, tanbRef, mu):

    tanb = tanbRef

    accuracy = 0.001

    counter = 0
    xSec = crosssection.whTauNuCrossSectionMSSM(mHp, tanb, mu)
    while (abs(xSecAtLimit - xSec) / xSecAtLimit > accuracy and xSec > 0
           and tanb > 1.1):
        tanb = tanb - 0.1 * (xSecAtLimit - xSec) / xSecAtLimit * tanb
        xSec = crosssection.whTauNuCrossSectionMSSM(mHp, tanb, mu)
        counter = counter + 1
        if counter > 10:  # to prevent infinite loops
            counter = 0
            accuracy = accuracy * 2
    if tanb < 1.1 or xSec <= 0:
        return -1
    return tanb
예제 #3
0
def tanbForXsec(xSecAtLimit,mHp,tanbRef,mu):

    tanb = tanbRef

    accuracy = 0.001    

    counter = 0
    xSec = crosssection.whTauNuCrossSectionMSSM(mHp,tanb,mu)
    while(abs(xSecAtLimit - xSec)/xSecAtLimit > accuracy and xSec > 0 and tanb < 99 ):
        tanb = tanb + 0.1*(xSecAtLimit - xSec)/xSecAtLimit*tanb
	xSec = crosssection.whTauNuCrossSectionMSSM(mHp,tanb,mu)
	counter = counter+1
	if counter > 10: # to prevent infinite loops
	    counter = 0
	    accuracy = accuracy*2
#        print "       tanbForXsec loop, tanb, xsec ",tanb,xSec
#    print "check tanbForXsec ",tanb
    if tanb > 99 or xSec <= 0:
	return -1
    return tanb
def tanbForXsec(xSecAtLimit, mHp, tanbRef, mu):

    tanb = tanbRef

    accuracy = 0.001

    counter = 0
    xSec = crosssection.whTauNuCrossSectionMSSM(mHp, tanb, mu)
    while (abs(xSecAtLimit - xSec) / xSecAtLimit > accuracy and xSec > 0
           and tanb < 99):
        tanb = tanb + 0.1 * (xSecAtLimit - xSec) / xSecAtLimit * tanb
        xSec = crosssection.whTauNuCrossSectionMSSM(mHp, tanb, mu)
        counter = counter + 1
        if counter > 10:  # to prevent infinite loops
            counter = 0
            accuracy = accuracy * 2


#        print "       tanbForXsec loop, tanb, xsec ",tanb,xSec
#    print "check tanbForXsec ",tanb
    if tanb > 99 or xSec <= 0:
        return -1
    return tanb
예제 #5
0
def signalXsecAtNsigma(nSignal,tanbRef,nSignalAtLimit,mHp,mu):
    xSec_atRef = crosssection.whTauNuCrossSectionMSSM(mHp, tanbRef, mu)
    xSec_atLimit = xSec_atRef*nSignalAtLimit/nSignal
#    print "check signalXsecAtNsigma xSec_atRef,xSec_atLimit,nSignalAtLimit,nSignal ",xSec_atRef,xSec_atLimit,nSignalAtLimit,nSignal
    return xSec_atLimit
예제 #6
0
def signalXsecAtNsigma(nSignal, tanbRef, nSignalAtLimit, mHp, mu):
    xSec_atRef = crosssection.whTauNuCrossSectionMSSM(mHp, tanbRef, mu)
    xSec_atLimit = xSec_atRef * nSignalAtLimit / nSignal
    #    print "check signalXsecAtNsigma xSec_atRef,xSec_atLimit,nSignalAtLimit,nSignal ",xSec_atRef,xSec_atLimit,nSignalAtLimit,nSignal
    return xSec_atLimit