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)))