Changeset View
Changeset View
Standalone View
Standalone View
src/LauKMatrixPropagator.cc
| Show First 20 Lines • Show All 909 Lines • ▼ Show 20 Lines | void LauKMatrixPropagator::calcGammaMatrix(const Double_t s) | ||||
| // Initialise all entries to zero | // Initialise all entries to zero | ||||
| GammaMatrix_.Zero(); | GammaMatrix_.Zero(); | ||||
| Double_t gamma(0.0); | Double_t gamma(0.0); | ||||
| for (Int_t iChannel (0); iChannel < nChannels_; ++iChannel) { | for (Int_t iChannel (0); iChannel < nChannels_; ++iChannel) { | ||||
| LauKMatrixPropagator::KMatrixChannels phaseSpaceIndex = phaseSpaceTypes_[iChannel]; | |||||
| if ( L_[iChannel] != 0 ) { | if ( L_[iChannel] != 0 ) { | ||||
| gamma = this->calcGamma(iChannel,s,phaseSpaceIndex); | gamma = this->calcGamma(iChannel,s); | ||||
| } else { | } else { | ||||
| gamma = 1.0; // S-wave | gamma = 1.0; // S-wave | ||||
| } | } | ||||
| if (verbose_) { | if (verbose_) { | ||||
| cout<<"GammaMatrix("<<iChannel<<", "<<iChannel<<") = "<<gamma<<endl; | cout<<"GammaMatrix("<<iChannel<<", "<<iChannel<<") = "<<gamma<<endl; | ||||
| } | } | ||||
| GammaMatrix_(iChannel, iChannel) = gamma; | GammaMatrix_(iChannel, iChannel) = gamma; | ||||
| } | } | ||||
| } | } | ||||
| Double_t LauKMatrixPropagator::calcGamma(const Int_t iCh, const Double_t s, const LauKMatrixPropagator::KMatrixChannels phaseSpaceIndex) const | Double_t LauKMatrixPropagator::calcGamma(const Int_t iCh, const Double_t s) const | ||||
| { | { | ||||
| // Calculate the barrier factor | // Calculate the barrier factor | ||||
| Double_t gamma(0.0); | Double_t gamma(0.0); | ||||
| LauKMatrixPropagator::KMatrixChannels phaseSpaceIndex = phaseSpaceTypes_[iCh]; | |||||
| LauComplex rho = getRho(s,phaseSpaceIndex); | LauComplex rho = getRho(s,phaseSpaceIndex); | ||||
| Double_t q = 0.5 * sqrt(s) * rho.abs(); | Double_t q = 0.5 * sqrt(s) * rho.abs(); | ||||
| gamma = pow(q,L_[iCh]); | gamma = pow(q,L_[iCh]); | ||||
| if (includeBWBarrierFactor_) | if (includeBWBarrierFactor_) | ||||
| { | { | ||||
| gamma /= pow( q*q + gamAInvRadSq_[iCh] , L_[iCh]/2. ); | gamma /= pow( q*q + gamAInvRadSq_[iCh] , L_[iCh]/2. ); | ||||
| } | } | ||||
| ▲ Show 20 Lines • Show All 477 Lines • Show Last 20 Lines | |||||