Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenModels/EvtBsMuMuKK.cpp
Show First 20 Lines • Show All 547 Lines • ▼ Show 20 Lines | |||||
// Integral | // Integral | ||||
double EvtBsMuMuKK::Integral( const double Gamma0, const double m0, const int JR, | double EvtBsMuMuKK::Integral( const double Gamma0, const double m0, const int JR, | ||||
const int JB, const double q0, const double M_KK_ll, | const int JB, const double q0, const double M_KK_ll, | ||||
const double M_KK_ul, const int fcntype ) const | const double M_KK_ul, const int fcntype ) const | ||||
{ | { | ||||
int bins = 1000; | int bins = 1000; | ||||
double bin_width = ( M_KK_ul - M_KK_ll ) / static_cast<double>( bins ); | double bin_width = ( M_KK_ul - M_KK_ll ) / static_cast<double>( bins ); | ||||
EvtComplex integral( 0.0, 0.0 ); | EvtComplex integral( 0.0, 0.0 ); | ||||
double MKpiKm2 = pow( MKp + MKm, 2 ); | double sumMKpKm2 = pow( MKp + MKm, 2 ); | ||||
double diffMKpKm2 = pow( MKp - MKm, 2 ); | |||||
double MBs2 = pow( MBs, 2 ); | double MBs2 = pow( MBs, 2 ); | ||||
for ( int i = 0; i < bins; i++ ) { | for ( int i = 0; i < bins; i++ ) { | ||||
double M_KK_i = M_KK_ll + static_cast<double>( i ) * bin_width; | double M_KK_i = M_KK_ll + static_cast<double>( i ) * bin_width; | ||||
double M_KK_f = M_KK_ll + static_cast<double>( i + 1 ) * bin_width; | double M_KK_f = M_KK_ll + static_cast<double>( i + 1 ) * bin_width; | ||||
double M_KK_i_sq = M_KK_i * M_KK_i; | |||||
double M_KK_f_sq = M_KK_f * M_KK_f; | |||||
double p3Kp_KK_CMS_i = ( pow( M_KK_i, 2 ) - MKpiKm2 ) / ( 2.0 * M_KK_i ); | double p3Kp_KK_CMS_i = sqrt( ( M_KK_i_sq - sumMKpKm2 ) * | ||||
double p3Kp_KK_CMS_f = ( pow( M_KK_f, 2 ) - MKpiKm2 ) / ( 2.0 * M_KK_f ); | ( M_KK_i_sq - diffMKpKm2 ) ) / | ||||
( 2.0 * M_KK_i ); | |||||
double p3Kp_KK_CMS_f = sqrt( ( M_KK_f_sq - sumMKpKm2 ) * | |||||
( M_KK_f_sq - diffMKpKm2 ) ) / | |||||
( 2.0 * M_KK_f ); | |||||
double p3Jpsi_Bs_CMS_i = sqrt( ( MBs2 - pow( M_KK_i + MJpsi, 2 ) ) * | double p3Jpsi_Bs_CMS_i = sqrt( ( MBs2 - pow( M_KK_i + MJpsi, 2 ) ) * | ||||
( MBs2 - pow( M_KK_i - MJpsi, 2 ) ) ) / | ( MBs2 - pow( M_KK_i - MJpsi, 2 ) ) ) / | ||||
( 2.0 * MBs ); | ( 2.0 * MBs ); | ||||
double p3Jpsi_Bs_CMS_f = sqrt( ( MBs2 - pow( M_KK_f + MJpsi, 2 ) ) * | double p3Jpsi_Bs_CMS_f = sqrt( ( MBs2 - pow( M_KK_f + MJpsi, 2 ) ) * | ||||
( MBs2 - pow( M_KK_f - MJpsi, 2 ) ) ) / | ( MBs2 - pow( M_KK_f - MJpsi, 2 ) ) ) / | ||||
( 2.0 * MBs ); | ( 2.0 * MBs ); | ||||
▲ Show 20 Lines • Show All 103 Lines • Show Last 20 Lines |