def magV(mass, resonance, wave, waves, normint):
"""
Returns the magnitude of the production amplitued V, for a specified wave,
resonance, wave in the set of waves, and normalization integral.
"""
return numpy.sqrt(
(1. / normint[wave.epsilon, wave.epsilon,
waves.index(wave),
waves.index(wave)]) * resonance.wR[waves.index(wave)] *
resonance.cR * breitWigner(mass, resonance.w0, resonance.r0) *
numpy.conjugate(breitWigner(mass, resonance.w0, resonance.r0)))
def V(self,mass,beta,k,eps):
BW=numpy.complex(0.,0.)
g=self.getStrength(beta)
if(k<1):
if eps <=1:
return {
1:breitWigner(mass,self.MASS1,self.GAMMA1)*g,
2:breitWigner(mass,self.MASS2,self.GAMMA2)*g,
3:breitWigner(mass,self.MASS3,self.GAMMA3)*g,
4:breitWigner(mass,self.MASS4,self.GAMMA4)*g,
5:breitWigner(mass,self.MASS5,self.GAMMA5)*g,
6:breitWigner(mass,self.MASS_A1,self.WIDTH_A1)*g,
7:breitWigner(mass,self.MASS_PI2,self.WIDTH_PI2)*g,
8:breitWigner(mass,self.MASS_PI1,self.WIDTH_PI1)*g
}[beta]
def V(self, mass, beta, k, eps):
BW = numpy.complex(0., 0.)
g = self.getStrength(beta)
if (k < 1):
if eps <= 1:
return {
1: breitWigner(mass, self.MASS1, self.GAMMA1) * g,
2: breitWigner(mass, self.MASS2, self.GAMMA2) * g,
3: breitWigner(mass, self.MASS3, self.GAMMA3) * g,
4: breitWigner(mass, self.MASS4, self.GAMMA4) * g,
5: breitWigner(mass, self.MASS5, self.GAMMA5) * g,
6: breitWigner(mass, self.MASS_A1, self.WIDTH_A1) * g,
7: breitWigner(mass, self.MASS_PI2, self.WIDTH_PI2) * g,
8: breitWigner(mass, self.MASS_PI1, self.WIDTH_PI1) * g
}[beta]
def magV(mass,resonance,wave,waves,normint):
"""
Returns the magnitude of the production amplitued V, for a specified wave,
resonance, wave in the set of waves, and normalization integral.
"""
return numpy.sqrt((1./normint[wave.epsilon,wave.epsilon,waves.index(wave),waves.index(wave)])*resonance.wR[waves.index(wave)]*resonance.cR*breitWigner(mass,resonance.w0,resonance.r0)*numpy.conjugate(breitWigner(mass,resonance.w0,resonance.r0)))
