EvtD0ToKSKppimpi0.cc

Go to the documentation of this file.

//--------------------------------------------------------------------------
// Environment:
//      This software is part of models developed at BES collaboration
//      based on the EvtGen framework.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/BesCopyright
//      Copyright (A) 2006      Ping Rong-Gang @IHEP
//
// Module: EvtD0ToKSKppimpi0.cc
//         the necessary file: EvtD0ToKSKppimpi0.hh
//
// Description: D0 -> KS K+ pi- pi0
//             
// Modification history:
//
//    Zhenxuan Li    Jun. 15, 2023       Module created
//------------------------------------------------------------------------
#include "EvtGenBase/EvtPatches.hh"
#include "EvtGenBase/EvtParticle.hh"
#include "EvtGenBase/EvtGenKine.hh"
#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtReport.hh"
#include "EvtGenModels/EvtD0ToKSKppimpi0.hh"
#include "EvtGenBase/EvtComplex.hh"
#include "EvtGenBase/EvtDecayTable.hh"
#include <stdlib.h>
#include <iostream>
#include <cmath>
using namespace std;
 
EvtD0ToKSKppimpi0::~EvtD0ToKSKppimpi0() {}
 
void EvtD0ToKSKppimpi0::getName(std::string& model_name){
    model_name="D0ToKSKppimpi0";
}
 
EvtDecayBase* EvtD0ToKSKppimpi0::clone(){
    return new EvtD0ToKSKppimpi0;
}
 
void EvtD0ToKSKppimpi0::init(){
    checkNArg(0);
    checkNDaug(4);
    checkSpinParent(EvtSpinType::SCALAR);
    //cout << "Initializing EvtD0ToKSKppimpi0" << endl;
 
    mD = 1.86484; mPi  = 0.13957; mKa  = 0.493677;
    math_pi = 3.1415926;
    mKstr0 = 0.89555; mrhop = 0.77511;
    GKstr0 = 0.0473;  Grhop = 0.1491;
    mKstrp = 0.89167; mrho0 = 0.77526;
    GKstrp = 0.0514;  Grho0 = 0.1474;
    ma980  = 0.980;
    Ga980  = 0.0972; 
    ma1_1260 = 1.230; mK1_1270 = 1.272;
    Ga1_1260 = 0.420; GK1_1270 = 0.087;
    meta1475 = 1.475; mK1_1400 = 1.403;
    Geta1475 = 0.09;  GK1_1400 = 0.174;
 
    rho[1] = 1.0;   //2S
    phi[1] = 0.0;   
    rho[0] = 1.1688955687;  //1S  
    phi[0] = 5.4715577597;  
    rho[2] = 1.2403434033;  //2P
    phi[2] =-4.2080261661;  
    rho[3] = 4.9828746114;  //16S
    phi[3] = 5.6470624082;  
    rho[4] = 2.8452580959;  //201
    phi[4] = 0.1889705326;  
    rho[5] = 2.8452580959;  //201
    phi[5] = 0.1889705326;  
    rho[6] = 8.7553053301;  //101
    phi[6] =-5.9775594415;  
 
    modetype[0]=1001;
    modetype[1]=1002;
    modetype[2]=1002;
    modetype[3]=2008;
    modetype[4]=4001;
    modetype[5]=4002;
    modetype[6]=3001;
 
    int GG[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
    int EE[4][4][4][4] =
    { { {{0,0,0,0},  {0,0,0,0},  {0,0,0,0},  {0,0,0,0} },
          {{0,0,0,0},  {0,0,0,0},  {0,0,0,1},  {0,0,-1,0}},
          {{0,0,0,0},  {0,0,0,-1}, {0,0,0,0},  {0,1,0,0} },
          {{0,0,0,0},  {0,0,1,0},  {0,-1,0,0}, {0,0,0,0} } },
    { {{0,0,0,0},  {0,0,0,0},  {0,0,0,-1}, {0,0,1,0} },
        {{0,0,0,0},  {0,0,0,0},  {0,0,0,0},  {0,0,0,0} },
        {{0,0,0,1},  {0,0,0,0},  {0,0,0,0},  {-1,0,0,0}},
        {{0,0,-1,0}, {0,0,0,0},  {1,0,0,0},  {0,0,0,0} } },
    { {{0,0,0,0},  {0,0,0,1},  {0,0,0,0},  {0,-1,0,0}},
        {{0,0,0,-1}, {0,0,0,0},  {0,0,0,0},  {1,0,0,0} },
        {{0,0,0,0},  {0,0,0,0},  {0,0,0,0},  {0,0,0,0} },
        {{0,1,0,0},  {-1,0,0,0}, {0,0,0,0},  {0,0,0,0} } },
    { {{0,0,0,0},  {0,0,-1,0}, {0,1,0,0},  {0,0,0,0} },
        {{0,0,1,0},  {0,0,0,0},  {-1,0,0,0}, {0,0,0,0} },
        {{0,-1,0,0}, {1,0,0,0},  {0,0,0,0},  {0,0,0,0} },
        {{0,0,0,0},  {0,0,0,0},  {0,0,0,0},  {0,0,0,0} } } };
    for (int i=0; i<4; i++) {
        for (int j=0; j<4; j++) {
            G[i][j] = GG[i][j];
            for (int k=0; k<4; k++) {
                for (int l=0; l<4; l++) {
                    E[i][j][k][l] = EE[i][j][k][l];
                }
            }
        }
    }
}
 
void EvtD0ToKSKppimpi0::initProbMax() {
    setProbMax(172500);
}
 
void EvtD0ToKSKppimpi0::decay( EvtParticle *p){
    //-----------for max value------------------
/*  double maxprob = 0.0;
    for(int ir=0;ir<=60000000;ir++){
        p->initializePhaseSpace(getNDaug(),getDaugs());
        EvtVector4R ks  = p->getDaug(0)->getP4();
        EvtVector4R kp  = p->getDaug(1)->getP4();
        EvtVector4R pim = p->getDaug(2)->getP4();
        EvtVector4R pi0 = p->getDaug(3)->getP4();
        mother_c=EvtPDL::getStdHep(p->getId());
        int cc;
        if(mother_c==421){ 
            cc=1;
        } else if(mother_c==-421){
            cc=-1;
        }
        double value;
        double KS[4],Kp[4],Pim[4],Pi0[4];
        KS[0] = ks.get(0);    Kp[0] = kp.get(0);    Pim[0] = pim.get(0);    Pi0[0] = pi0.get(0);
        KS[1] = cc*ks.get(1); Kp[1] = cc*kp.get(1); Pim[1] = cc*pim.get(1); Pi0[1] = cc*pi0.get(1);
        KS[2] = cc*ks.get(2); Kp[2] = cc*kp.get(2); Pim[2] = cc*pim.get(2); Pi0[2] = cc*pi0.get(2);
        KS[3] = cc*ks.get(3); Kp[3] = cc*kp.get(3); Pim[3] = cc*pim.get(3); Pi0[3] = cc*pi0.get(3);
        calPDF(KS, Kp, Pim, Pi0, value);
        if(value>maxprob) {
            maxprob=value;
            std::cout << "Max PDF = " << ir << " prob= " << value << std::endl;
        }
    } 
    std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
    return;*/
    //-----------------------------------------------
        p->initializePhaseSpace(getNDaug(),getDaugs());
        EvtVector4R ks  = p->getDaug(0)->getP4();
        EvtVector4R kp  = p->getDaug(1)->getP4();
        EvtVector4R pim = p->getDaug(2)->getP4();
        EvtVector4R pi0 = p->getDaug(3)->getP4();
 
        mother_c=EvtPDL::getStdHep(p->getId());
        int cc;
        if(mother_c==421){ 
        cc=1;
        } else if(mother_c==-421){
        cc=-1;
        }
 
        double KS[4],Kp[4],Pim[4],Pi0[4];
        KS[0] = ks.get(0);    Kp[0] = kp.get(0);    Pim[0] = pim.get(0);    Pi0[0] = pi0.get(0);
        KS[1] = cc*ks.get(1); Kp[1] = cc*kp.get(1); Pim[1] = cc*pim.get(1); Pi0[1] = cc*pi0.get(1);
        KS[2] = cc*ks.get(2); Kp[2] = cc*kp.get(2); Pim[2] = cc*pim.get(2); Pi0[2] = cc*pi0.get(2);
        KS[3] = cc*ks.get(3); Kp[3] = cc*kp.get(3); Pim[3] = cc*pim.get(3); Pi0[3] = cc*pi0.get(3);
    //KS[0] = 0.581505; Kp[0] = 0.654953; Pim[0] = 0.344256; Pi0[0] = 0.297861;
    //KS[1] = 0.302649; Kp[1] =-0.369050; Pim[1] =-0.164686; Pi0[1] = 0.027881;
    //KS[2] = 0.038261; Kp[2] =-0.221134; Pim[2] =-0.247290; Pi0[2] = 0.203629;
    //KS[3] =-0.032410; Kp[3] =-0.012204; Pim[3] =-0.103723; Pi0[3] = 0.168109;
 
    double value;
    calPDF(KS, Kp, Pim, Pi0, value);
    setProb(value);
    return;
 
}
 
double EvtD0ToKSKppimpi0::calPDF(double KS[], double Kp[], double Pim[], double Pi0[], double & Result) {
    double cof[2], amp_tmp[2], amp_PDF[2], PDF[2];
    double flag[3], mass_R[2], width_R[2];  
    double sa[3], sb[3], sc[3], B[3];
    double t1D[4], t1V[4];//, t1A[4];
    double pS[4], pV[4], pD[4], pA[4], pP[4];
    double pro1[2], pro2[2], proKPi_S[2], pro[2];
    double temp_PDF = 0;
 
    double rD2 = 25.0;
    double rRes2 = 9.0;
    double mass1[7] ={mKstr0, mKstrp, mKstrp,    mrhop,   mKstrp,   mKstrp, ma980};
    double mass2[7] ={mKstr0, mKstrp, mKstrp, mK1_1270, meta1475, meta1475, mrhop};
    double width1[7]={GKstr0, GKstrp, GKstrp,    Grhop,   GKstrp,   GKstrp, Ga980};
    double width2[7]={GKstr0, GKstrp, GKstrp, GK1_1270, Geta1475, Geta1475, Grhop};
    double g0[7]    ={     1,      1,      1,        1,        1,        1,     1};
    double g1[7]    ={     1,      1,      1,        1,        1,        1,     1};
    double g2[7]    ={     0,      0,      1,        0,        0,        0,     1};
 
    PDF[0]=0; PDF[1]=0;
 
    for(int i=0; i<7; i++) {
        flag[0] = g0[i]; flag[1] = g1[i];flag[2] = g2[i];
        mass_R[0] = mass1[i]; mass_R[1] = mass2[i];
        width_R[0] = width1[i]; width_R[1] = width2[i];
 
        amp_tmp[0] = 0;
        amp_tmp[1] = 0;
 
        cof[0] = rho[i]*cos(phi[i]);
        cof[1] = rho[i]*sin(phi[i]);
 
        if( modetype[i] == 1001 ){//D0 -> K*0 K*0
            double t1V1[4], t1V2[4], t2D[4][4];
            temp_PDF = 0;
            double pV1[4], pV2[4];
            for(int ii=0; ii!=4; ii++){
                pV1[ii] = KS[ii] + Pi0[ii];
                pV2[ii] = Kp[ii] + Pim[ii];
                pD[ii] = pV1[ii] + pV2[ii];
            }
            sa[0] = SCADot(pV1,pV1);
            sb[0] = SCADot(KS,KS);
            sc[0] = SCADot(Pi0,Pi0);
            sa[1] = SCADot(pV2,pV2);
            sb[1] = SCADot(Kp,Kp);
            sc[1] = SCADot(Pim,Pim);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[0];
            sc[2] = sa[1];
            if(flag[0] == 1){
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0],width_R[0],sa[0],sb[0],sc[0],rRes2,1, pro1);
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
            if(flag[1] == 1){
                propagatorRBW(mass_R[1]*mass_R[1],mass_R[1],width_R[1],sa[1],sb[1],sc[1],rRes2,1, pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
            B[0] = Barrier(mass_R[0]*mass_R[0],1,sa[0],sb[0],sc[0],rRes2);
            B[1] = Barrier(mass_R[1]*mass_R[1],1,sa[1],sb[1],sc[1],rRes2);
            calt1(KS,Pi0,t1V1);
            calt1(Kp,Pim,t1V2);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    temp_PDF += (G[ii][ii])*t1V1[ii]*t1V2[ii];
                }
                B[2] = 1;
            }
            if(flag[2] == 1){
                calt1(pV1,pV2,t1D);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        for(int k=0; k!=4; k++){
                            for(int l=0; l!=4; l++){
                                temp_PDF += E[ii][j][k][l]*pD[ii]*t1D[j]*t1V1[k]*t1V2[l]*
                                    (G[ii][ii])*(G[j][j])*(G[l][l])*(G[k][k]);
                            }
                        }
                    }
                }
                B[2] = Barrier(mD*mD,1,sa[2],sb[2],sc[2],rD2);
            }
            if(flag[2] == 2){
                calt2(pV1,pV2,t2D);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t2D[ii][j]*t1V1[ii]*t1V2[j]*(G[ii][ii])*(G[j][j]);
                    }
                }
                B[2] = Barrier(mD*mD,2,sa[2],sb[2],sc[2],rD2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if( modetype[i] == 1002 ){//D0 -> K*- K*+
            double t1V1[4], t1V2[4], t2D[4][4];
            temp_PDF = 0;
            double pV1[4], pV2[4];
            for(int ii=0; ii!=4; ii++){
                pV1[ii] = KS[ii] + Pim[ii];
                pV2[ii] = Kp[ii] + Pi0[ii];
                pD[ii] = pV1[ii] + pV2[ii];
            }
            sa[0] = SCADot(pV1,pV1);
            sb[0] = SCADot(KS,KS);
            sc[0] = SCADot(Pim,Pim);
            sa[1] = SCADot(pV2,pV2);
            sb[1] = SCADot(Kp,Kp);
            sc[1] = SCADot(Pi0,Pi0);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[0];
            sc[2] = sa[1];
            if(flag[0] == 1){
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0],width_R[0],sa[0],sb[0],sc[0],rRes2,1, pro1);
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
            if(flag[1] == 1){
                propagatorRBW(mass_R[1]*mass_R[1],mass_R[1],width_R[1],sa[1],sb[1],sc[1],rRes2,1, pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
            B[0] = Barrier(mass_R[0]*mass_R[0],1,sa[0],sb[0],sc[0],rRes2);
            B[1] = Barrier(mass_R[1]*mass_R[1],1,sa[1],sb[1],sc[1],rRes2);
            calt1(KS,Pim,t1V1);
            calt1(Kp,Pi0,t1V2);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    temp_PDF += (G[ii][ii])*t1V1[ii]*t1V2[ii];
                }
                B[2] = 1;
            }
            if(flag[2] == 1){
                calt1(pV1,pV2,t1D);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        for(int k=0; k!=4; k++){
                            for(int l=0; l!=4; l++){
                                temp_PDF += E[ii][j][k][l]*pD[ii]*t1D[j]*t1V1[k]*t1V2[l]*
                                    (G[ii][ii])*(G[j][j])*(G[l][l])*(G[k][k]);
                            }
                        }
                    }
                }
                B[2] = Barrier(mD*mD,1,sa[2],sb[2],sc[2],rD2);
            }
            if(flag[2] == 2){
                calt2(pV1,pV2,t2D);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t2D[ii][j]*t1V1[ii]*t1V2[j]*(G[ii][ii])*(G[j][j]);
                    }
                }
                B[2] = Barrier(mD*mD,2,sa[2],sb[2],sc[2],rD2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2001){//D0 -> a1()+ pi-, a1()+ -> K*+ KS
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = Kp[ii] + Pi0[ii];
                pA[ii] = pV[ii] + KS[ii];
                pD[ii] = pA[ii] + Pim[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(Kp,Kp);
            sc[0] = SCADot(Pi0,Pi0);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(KS,KS);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Pim,Pim);
 
            if(flag[0] == 1){
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(Kp,Pi0,t1V);
            calt1(pA,Pim,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,KS,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2002){//D0 -> a1()+ pi-, a1()+ -> K*0 K+
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = KS[ii] + Pi0[ii];
                pA[ii] = pV[ii] + Kp[ii];
                pD[ii] = pA[ii] + Pim[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(KS,KS);
            sc[0] = SCADot(Pi0,Pi0);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(Kp,Kp);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Pim,Pim);
 
            if(flag[0] == 1){
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(KS,Pi0,t1V);
            calt1(pA,Pim,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,Kp,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2003){//D0 -> a1()0 pi0, a1()0 -> K*0 KS
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = Kp[ii] + Pim[ii];
                pA[ii] = pV[ii] + KS[ii];
                pD[ii] = pA[ii] + Pi0[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(Kp,Kp);
            sc[0] = SCADot(Pim,Pim);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(KS,KS);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Pi0,Pi0);
 
            if(flag[0] == 1){
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(Kp,Pim,t1V);
            calt1(pA,Pi0,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,KS,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2004){//D0 -> a1()0 pi0, a1()0 -> K*- K+
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = KS[ii] + Pim[ii];
                pA[ii] = pV[ii] + Kp[ii];
                pD[ii] = pA[ii] + Pi0[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(KS,KS);
            sc[0] = SCADot(Pim,Pim);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(Kp,Kp);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Pi0,Pi0);
 
            if(flag[0] == 1){
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(KS,Pim,t1V);
            calt1(pA,Pi0,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,Kp,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2005){//D0 -> K1(1270)- K+, K1()- -> KS rho-
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = Pim[ii] + Pi0[ii];
                pA[ii] = pV[ii] + KS[ii];
                pD[ii] = pA[ii] + Kp[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(Pim,Pim);
            sc[0] = SCADot(Pi0,Pi0);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(KS,KS);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Kp,Kp);
 
            if(flag[0] == 1){
                propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, pro1);//rho
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW_K1(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(Pim,Pi0,t1V);
            calt1(pA,Kp,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,KS,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2006){//D0 -> K1(1270/1400)- K+, K1()- -> K*0 pi-
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = KS[ii] + Pi0[ii];
                pA[ii] = pV[ii] + Pim[ii];
                pD[ii] = pA[ii] + Kp[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(KS,KS);
            sc[0] = SCADot(Pi0,Pi0);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(Pim,Pim);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Kp,Kp);
 
            if(flag[0] == 1){
                //propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, pro1);//rho
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW_K1(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(KS,Pi0,t1V);
            calt1(pA,Kp,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,Pim,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2007){//D0 -> K1(1270/1400)- K+, K1()- -> K*- pi0
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = KS[ii] + Pim[ii];
                pA[ii] = pV[ii] + Pi0[ii];
                pD[ii] = pA[ii] + Kp[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(KS,KS);
            sc[0] = SCADot(Pim,Pim);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(Pi0,Pi0);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Kp,Kp);
 
            if(flag[0] == 1){
                //propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, pro1);//rho
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW_K1(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(KS,Pim,t1V);
            calt1(pA,Kp,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,Pi0,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = -1.0*temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = -1.0*temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2008){//D0 -> K1(1270)0 KS, K1()0 -> K+ rho-
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = Pim[ii] + Pi0[ii];
                pA[ii] = pV[ii] + Kp[ii];
                pD[ii] = pA[ii] + KS[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(Pim,Pim);
            sc[0] = SCADot(Pi0,Pi0);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(Kp,Kp);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(KS,KS);
 
            if(flag[0] == 1){
                propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, pro1);//rho
                //propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW_K1(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(Pim,Pi0,t1V);
            calt1(pA,KS,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,Kp,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2009){//D0 -> K1(1270/1400)0 KS, K1()0 -> K*+ pi-
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = Kp[ii] + Pi0[ii];
                pA[ii] = pV[ii] + Pim[ii];
                pD[ii] = pA[ii] + KS[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(Kp,Kp);
            sc[0] = SCADot(Pi0,Pi0);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(Pim,Pim);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(KS,KS);
 
            if(flag[0] == 1){
                //propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, pro1);//rho
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW_K1(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(Kp,Pi0,t1V);
            calt1(pA,KS,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,Pim,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if(modetype[i] == 2010){//D0 -> K1(1270/1400)0 KS, K1()0 -> K*0 pi0
            temp_PDF = 0;
            double t2A[4][4];
            for(int ii=0; ii!=4; ii++){
                pV[ii] = Kp[ii] + Pim[ii];
                pA[ii] = pV[ii] + Pi0[ii];
                pD[ii] = pA[ii] + KS[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(Kp,Kp);
            sc[0] = SCADot(Pim,Pim);
            sa[1] = SCADot(pA,pA);
            sb[1] = sa[0];
            sc[1] = SCADot(Pi0,Pi0);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(KS,KS);
 
            if(flag[0] == 1){
                //propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, pro1);//rho
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW_K1(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, flag[2], pro2);
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1,sa[0],sb[0],sc[0],rRes2);
            B[2] = Barrier(mD*mD, 1,sa[2],sb[2],sc[2],rD2);
            calt1(Kp,Pim,t1V);
            calt1(pA,KS,t1D);
            if(flag[2] == 0){
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*(G[ii][j] - pA[ii]*pA[j]/sa[1])*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = 1;
            }
            else if(flag[2] == 2){
                calt2(pV,Pi0,t2A);
                for(int ii=0; ii!=4; ii++){
                    for(int j=0; j!=4; j++){
                        temp_PDF += t1D[ii]*t2A[ii][j]*t1V[j]*G[ii][ii]*G[j][j];
                    }
                }
                B[1] = Barrier(mass_R[1]*mass_R[1], 2,sa[1],sb[1],sc[1],rRes2);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = -1.0*temp_PDF*B[0]*B[1]*B[2]*pro[0];
            amp_tmp[1] = -1.0*temp_PDF*B[0]*B[1]*B[2]*pro[1];
 
        }else if( modetype[i] == 3001){//D0 to a0(980)+ rho-
            temp_PDF = 0;
            for(int ii=0; ii!=4; ii++){
                pS[ii] = KS[ii] + Kp[ii];
                pV[ii] = Pim[ii] + Pi0[ii];
                pD[ii] = pS[ii] + pV[ii];
            }
            sa[0] = SCADot(pS,pS);
            sb[0] = SCADot(KS,KS);
            sc[0] = SCADot(Kp,Kp);
            sa[1] = SCADot(pV,pV);
            sb[1] = SCADot(Pim,Pim);
            sc[1] = SCADot(Pi0,Pi0);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[0];
            sc[2] = sa[1];
            if(flag[0] == 1){
                propagatorFlatte(mass_R[0],width_R[0],sa[0],pro1);
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorGS(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, pro2);//rho
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[1] = Barrier(mass_R[1]*mass_R[1], 1, sa[1], sb[1], sc[1], rRes2);
            B[2] = Barrier(mD*mD, 1, sa[2], sb[2], sc[2], rD2);
 
            calt1(Pim, Pi0, t1V);
            calt1(pS, pV, t1D);
            for(int ii=0; ii!=4; ii++){
                temp_PDF += G[ii][ii]*t1D[ii]*t1V[ii];
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[1]*B[2]*pro[0];
            amp_tmp[1] = temp_PDF*B[1]*B[2]*pro[1];
 
        }else if( modetype[i] == 4001){//D0 to eta() pi0, eta() -> K*- K+
            temp_PDF = 0;
            for(int ii=0; ii!=4; ii++){
                pV[ii] = KS[ii] + Pim[ii];
                pP[ii] = pV[ii] + Kp[ii];
                pD[ii] = pP[ii] + Pi0[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(KS,KS);
            sc[0] = SCADot(Pim,Pim);
            sa[1] = SCADot(pP,pP);
            sb[1] = sa[0];
            sc[1] = SCADot(Kp,Kp);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Pi0,Pi0);;
            if(flag[0] == 1){
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, 1, pro2);//eta
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1, sa[0], sb[0], sc[0], rRes2);
            B[1] = Barrier(mass_R[1]*mass_R[1], 1, sa[1], sb[1], sc[1], rRes2);
 
            calt1(KS, Pim, t1V);
            for(int ii=0; ii!=4; ii++){
                temp_PDF += Kp[ii]*t1V[ii]*(G[ii][ii]);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*pro[1];
 
        }else if( modetype[i] == 4002){//D0 to eta() pi0, eta() -> K*0 KS
            temp_PDF = 0;
            for(int ii=0; ii!=4; ii++){
                pV[ii] = Kp[ii] + Pim[ii];
                pP[ii] = pV[ii] + KS[ii];
                pD[ii] = pP[ii] + Pi0[ii];
            }
            sa[0] = SCADot(pV,pV);
            sb[0] = SCADot(Kp,Kp);
            sc[0] = SCADot(Pim,Pim);
            sa[1] = SCADot(pP,pP);
            sb[1] = sa[0];
            sc[1] = SCADot(KS,KS);
            sa[2] = SCADot(pD,pD);
            sb[2] = sa[1];
            sc[2] = SCADot(Pi0,Pi0);;
            if(flag[0] == 1){
                propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2, 1, pro1);//K*
            }
            else if(flag[0] == 0){pro1[0] = 1; pro1[1] = 0;}
 
            if(flag[1] == 1){
                propagatorRBW(mass_R[1]*mass_R[1],mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2, 1, pro2);//eta
            }
            else if(flag[1] == 0){pro2[0] = 1; pro2[1] = 0;}
 
            B[0] = Barrier(mass_R[0]*mass_R[0], 1, sa[0], sb[0], sc[0], rRes2);
            B[1] = Barrier(mass_R[1]*mass_R[1], 1, sa[1], sb[1], sc[1], rRes2);
 
            calt1(Kp, Pim, t1V);
            for(int ii=0; ii!=4; ii++){
                temp_PDF += KS[ii]*t1V[ii]*(G[ii][ii]);
            }
            Com_Multi(pro1,pro2,pro);
            amp_tmp[0] = temp_PDF*B[0]*B[1]*pro[0];
            amp_tmp[1] = temp_PDF*B[0]*B[1]*pro[1];
 
        }
 
        Com_Multi(amp_tmp,cof,amp_PDF);
        //printf("pdf: %d  %.10f  %.10f\n",i,amp_PDF[0],amp_PDF[1]);
        PDF[0] += amp_PDF[0];
        PDF[1] += amp_PDF[1];
    }
    double value = PDF[0]*PDF[0] + PDF[1]*PDF[1];
    Result = value;
    //printf("PDF[0] :  %.10f \n",PDF[0]);
    //printf("PDF[1] :  %.10f \n",PDF[1]);
    //printf("Value :  %.10f \n",value);
}
 
void EvtD0ToKSKppimpi0::Com_Multi(double a1[2], double a2[2], double res[2])
{
    res[0] = a1[0]*a2[0]-a1[1]*a2[1];
    res[1] = a1[1]*a2[0]+a1[0]*a2[1];
}
void EvtD0ToKSKppimpi0::Com_Divide(double a1[2], double a2[2], double res[2])
{
    double tmp = a2[0]*a2[0]+a2[1]*a2[1];
    res[0] = (a1[0]*a2[0]+a1[1]*a2[1])/tmp;
    res[1] = (a1[1]*a2[0]-a1[0]*a2[1])/tmp;
}
double EvtD0ToKSKppimpi0::SCADot(double a1[4], double a2[4])
{
    double _cal = a1[0]*a2[0]-a1[1]*a2[1]-a1[2]*a2[2]-a1[3]*a2[3];
    return _cal;
}
double EvtD0ToKSKppimpi0::Barrier(double mass2, int l, double sa, double sb, double sc, double r2)
{
    double F;
    double tmp = sa+sb-sc;
    double q = fabs(0.25*tmp*tmp/sa-sb);
    //if (q < 0) q = 1e-16;
    double tmp2 = mass2+sb-sc;
    double q0 = fabs(0.25*tmp2*tmp2/mass2-sb);
    //if (q0 < 0) q0 = 1e-16;
    double z = q*r2;
    double z0 = q0*r2;
    if (l==1) {
        F = sqrt((1.0+z0)/(1.0+z));
    }
    else if (l==2) {
        double z2 = z*z; double z02 = z0*z0;
        F = sqrt((9.0+3.0*z0+z02)/(9.0+3.0*z+z2));
    } else {
        F = 1.0;
    }
    return F;
}
void EvtD0ToKSKppimpi0::calt1(double daug1[4], double daug2[4], double t1[4])
{
    double p, pq, tmp;
    double pa[4], qa[4];
    for(int i=0; i<4; i++) {
        pa[i] = daug1[i] + daug2[i];
        qa[i] = daug1[i] - daug2[i];
    }
    p = SCADot(pa,pa);
    pq = SCADot(pa,qa);
    tmp = pq/p;
    for(int i=0; i<4; i++) {
        t1[i] = qa[i] - tmp*pa[i];
    }
}
void EvtD0ToKSKppimpi0::calt2(double daug1[4], double daug2[4], double t2[4][4])
{
    double p, r;
    double pa[4], t1[4];
    calt1(daug1,daug2,t1);
    r = SCADot(t1,t1)/3.0;
    for(int i=0; i<4; i++) {
        pa[i] = daug1[i] + daug2[i];
    }
    p = SCADot(pa,pa);
    for(int i=0; i<4; i++) {
        for(int j=0; j<4; j++) {
            t2[i][j] = t1[i]*t1[j] - r*(G[i][j]-pa[i]*pa[j]/p);
        }
    }
}
void EvtD0ToKSKppimpi0::propagator(double mass2, double mass, double width, double sx, double prop[2])
{
    double a[2], b[2];
    a[0] = 1;
    a[1] = 0;
    b[0] = mass2-sx;
    b[1] = -mass*width;
    Com_Divide(a,b,prop);
}
double EvtD0ToKSKppimpi0::wid(double mass2, double mass, double sa, double sb, double sc, double r2, int l)
{
    double widm = 0.;
    double m = sqrt(sa);
    double tmp  = sb-sc;
    double tmp1 = sa+tmp;
    double q = fabs(0.25*tmp1*tmp1/sa-sb);
    //if(q<0) q = 1e-16;
    double tmp2 = mass2+tmp;
    double q0 = fabs(0.25*tmp2*tmp2/mass2-sb);
    //if(q0<0) q0 = 1e-16;
    double z = q*r2;
    double z0 = q0*r2;
    double t = q/q0;
    if(l == 0)      {widm = sqrt(t)*mass/m;}
    else if(l == 1) {widm = t*sqrt(t)*mass/m*(1+z0)/(1+z);}
    else if(l == 2) {widm = t*t*sqrt(t)*mass/m*(9+3*z0+z0*z0)/(9+3*z+z*z);}
    return widm;
}
double EvtD0ToKSKppimpi0::widl1(double mass2, double mass, double sa, double sb, double sc, double r2)
{
    double widm = 0.;
    double m = sqrt(sa);
    double tmp  = sb-sc;
    double tmp1 = sa+tmp;
    double q = fabs(0.25*tmp1*tmp1/sa-sb);
    //if(q<0) q = 1e-16;
    double tmp2 = mass2+tmp;
    double q0 = fabs(0.25*tmp2*tmp2/mass2-sb);
    //if(q0<0) q0 = 1e-16;
    double z = q*r2;
    double z0 = q0*r2;
    double F = (1+z0)/(1+z);
    double t = q/q0;
    widm = t*sqrt(t)*mass/m*F;
    return widm;
}
void EvtD0ToKSKppimpi0::propagatorRBW(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2])
{
    double a[2], b[2];
    a[0] = 1;
    a[1] = 0;
    b[0] = mass2-sa;
    b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r2,l);
    Com_Divide(a,b,prop);
}
void EvtD0ToKSKppimpi0::propagatorRBWl1(double mass2, double mass, double width, double sa, double sb, double sc, double r2, double prop[2])
{
    double a[2], b[2];
    a[0] = 1;
    a[1] = 0;
    b[0] = mass2-sa;
    b[1] = -mass*width*widl1(mass2,mass,sa,sb,sc,r2);
    Com_Divide(a,b,prop);
}
void EvtD0ToKSKppimpi0::propagatorRBW_a1(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2])
{
    double a[2], b[2];
    int iii=int(sa*1000/1.55)-115;
    if(iii<0){ iii=0; }
    double a1_width[2000]={//1885 bin, s = (0.17825, 3.1)
        0.0000000172, 0.0000000487, 0.0000001049, 0.0000001927, 0.0000003185, 
        0.0000004888, 0.0000007099, 0.0000009878, 0.0000013285, 0.0000017378, 
        0.0000022212, 0.0000027845, 0.0000034328, 0.0000041717, 0.0000050061, 
        0.0000059414, 0.0000069824, 0.0000081342, 0.0000094014, 0.0000107890, 
        0.0000123016, 0.0000139437, 0.0000157200, 0.0000176350, 0.0000196931, 
        0.0000218986, 0.0000242559, 0.0000267692, 0.0000294428, 0.0000322808, 
        0.0000352875, 0.0000384668, 0.0000418228, 0.0000453596, 0.0000490811, 
        0.0000529913, 0.0000570941, 0.0000613934, 0.0000658931, 0.0000705969, 
        0.0000755087, 0.0000806324, 0.0000859716, 0.0000915300, 0.0000973115, 
        0.0001033197, 0.0001095583, 0.0001160310, 0.0001227415, 0.0001296933, 
        0.0001368901, 0.0001443355, 0.0001520332, 0.0001599867, 0.0001681995, 
        0.0001766753, 0.0001854177, 0.0001944301, 0.0002037162, 0.0002132795, 
        0.0002231235, 0.0002332517, 0.0002436677, 0.0002543750, 0.0002653771, 
        0.0002766775, 0.0002882799, 0.0003001876, 0.0003124041, 0.0003249331, 
        0.0003377780, 0.0003509423, 0.0003644296, 0.0003782433, 0.0003923870, 
        0.0004068643, 0.0004216785, 0.0004368333, 0.0004523323, 0.0004681788, 
        0.0004843765, 0.0005009290, 0.0005178397, 0.0005351122, 0.0005527502, 
        0.0005707571, 0.0005891365, 0.0006078922, 0.0006270275, 0.0006465462, 
        0.0006664518, 0.0006867480, 0.0007074384, 0.0007285267, 0.0007500165, 
        0.0007719114, 0.0007942151, 0.0008169314, 0.0008400639, 0.0008636163, 
        0.0008875923, 0.0009119957, 0.0009368303, 0.0009620997, 0.0009878078, 
        0.0010139580, 0.0010405550, 0.0010676020, 0.0010951030, 0.0011230610, 
        0.0011514820, 0.0011803670, 0.0012097230, 0.0012395510, 0.0012698570, 
        0.0013006440, 0.0013319160, 0.0013636770, 0.0013959320, 0.0014286840, 
        0.0014619370, 0.0014956950, 0.0015299630, 0.0015647450, 0.0016000440, 
        0.0016358650, 0.0016722130, 0.0017090900, 0.0017465020, 0.0017844530, 
        0.0018229470, 0.0018619890, 0.0019015820, 0.0019417310, 0.0019824400, 
        0.0020237150, 0.0020655590, 0.0021079760, 0.0021509730, 0.0021945510, 
        0.0022387180, 0.0022834760, 0.0023288310, 0.0023747880, 0.0024213500, 
        0.0024685230, 0.0025163120, 0.0025647210, 0.0026137560, 0.0026634200, 
        0.0027137200, 0.0027646590, 0.0028162430, 0.0028684780, 0.0029213670, 
        0.0029749170, 0.0030291320, 0.0030840170, 0.0031395780, 0.0031958200, 
        0.0032527490, 0.0033103700, 0.0033686880, 0.0034277080, 0.0034874370, 
        0.0035478800, 0.0036090420, 0.0036709290, 0.0037335480, 0.0037969030, 
        0.0038610000, 0.0039258460, 0.0039914470, 0.0040578080, 0.0041249360, 
        0.0041928370, 0.0042615160, 0.0043309810, 0.0044012380, 0.0044722930, 
        0.0045441520, 0.0046168220, 0.0046903100, 0.0047646230, 0.0048397670, 
        0.0049157480, 0.0049925750, 0.0050702540, 0.0051487920, 0.0052281970, 
        0.0053084750, 0.0053896340, 0.0054716810, 0.0055546240, 0.0056384710, 
        0.0057232280, 0.0058089050, 0.0058955090, 0.0059830470, 0.0060715290, 
        0.0061609620, 0.0062513540, 0.0063427140, 0.0064350500, 0.0065283720, 
        0.0066226870, 0.0067180050, 0.0068143350, 0.0069116850, 0.0070100650, 
        0.0071094850, 0.0072099530, 0.0073114800, 0.0074140740, 0.0075177470, 
        0.0076225070, 0.0077283660, 0.0078353330, 0.0079434180, 0.0080526330, 
        0.0081629880, 0.0082744940, 0.0083871620, 0.0085010030, 0.0086160290, 
        0.0087322510, 0.0088496800, 0.0089683290, 0.0090882100, 0.0092093350, 
        0.0093317160, 0.0094553660, 0.0095802980, 0.0097065250, 0.0098340590, 
        0.0099629140, 0.0100931000, 0.0102246400, 0.0103575500, 0.0104918200, 
        0.0106274900, 0.0107645700, 0.0109030700, 0.0110430000, 0.0111843800, 
        0.0113272300, 0.0114715700, 0.0116174000, 0.0117647500, 0.0119136300, 
        0.0120640600, 0.0122160600, 0.0123696300, 0.0125248100, 0.0126816100, 
        0.0128400500, 0.0130001400, 0.0131619100, 0.0133253800, 0.0134905500, 
        0.0136574600, 0.0138261300, 0.0139965600, 0.0141688000, 0.0143428500, 
        0.0145187400, 0.0146964900, 0.0148761200, 0.0150576500, 0.0152411200, 
        0.0154265300, 0.0156139200, 0.0158033100, 0.0159947200, 0.0161881800, 
        0.0163837200, 0.0165813500, 0.0167811100, 0.0169830300, 0.0171871200, 
        0.0173934200, 0.0176019600, 0.0178127600, 0.0180258500, 0.0182412700, 
        0.0184590400, 0.0186791900, 0.0189017500, 0.0191267600, 0.0193542500, 
        0.0195842500, 0.0198167900, 0.0200519000, 0.0202896300, 0.0205300000, 
        0.0207730500, 0.0210188200, 0.0212673300, 0.0215186400, 0.0217727700, 
        0.0220297700, 0.0222896700, 0.0225525200, 0.0228183400, 0.0230871900, 
        0.0233591000, 0.0236341200, 0.0239122800, 0.0241936400, 0.0244782300, 
        0.0247661000, 0.0250572900, 0.0253518600, 0.0256498400, 0.0259512800, 
        0.0262562300, 0.0265647400, 0.0268768600, 0.0271926300, 0.0275121100, 
        0.0278353500, 0.0281624000, 0.0284933100, 0.0288281300, 0.0291669200, 
        0.0295097300, 0.0298566200, 0.0302076400, 0.0305628500, 0.0309223000, 
        0.0312860500, 0.0316541600, 0.0320266900, 0.0324036900, 0.0327852300, 
        0.0331713600, 0.0335621500, 0.0339576500, 0.0343579300, 0.0347630500, 
        0.0351730600, 0.0355880400, 0.0360080500, 0.0364331400, 0.0368633900, 
        0.0372988500, 0.0377395900, 0.0381856600, 0.0386371500, 0.0390941000, 
        0.0395565900, 0.0400246700, 0.0404984100, 0.0409778800, 0.0414631300, 
        0.0419542300, 0.0424512400, 0.0429542200, 0.0434632400, 0.0439783500, 
        0.0444996200, 0.0450271000, 0.0455608500, 0.0461009400, 0.0466474000, 
        0.0472003200, 0.0477597200, 0.0483256700, 0.0488982300, 0.0494774300, 
        0.0500633300, 0.0506559700, 0.0512554000, 0.0518616600, 0.0524748000, 
        0.0530948400, 0.0537218300, 0.0543558000, 0.0549967800, 0.0556448000, 
        0.0562998900, 0.0569620700, 0.0576313600, 0.0583077800, 0.0589913500, 
        0.0596820600, 0.0603799400, 0.0610849900, 0.0617972100, 0.0625166000, 
        0.0632431500, 0.0639768500, 0.0647176900, 0.0654656600, 0.0662207300, 
        0.0669828800, 0.0677520900, 0.0685283100, 0.0693115200, 0.0701016800, 
        0.0708987400, 0.0717026500, 0.0725133700, 0.0733308400, 0.0741550000, 
        0.0749858000, 0.0758231500, 0.0766670000, 0.0775172700, 0.0783738800, 
        0.0792367600, 0.0801058200, 0.0809809800, 0.0818621400, 0.0827492100, 
        0.0836421000, 0.0845407100, 0.0854449400, 0.0863546800, 0.0872698400, 
        0.0881902900, 0.0891159400, 0.0900466700, 0.0909823700, 0.0919229100, 
        0.0928682000, 0.0938181000, 0.0947725000, 0.0957312800, 0.0966943100, 
        0.0976614900, 0.0986326700, 0.0996077500, 0.1005866000, 0.1015691000, 
        0.1025551000, 0.1035445000, 0.1045372000, 0.1055331000, 0.1065320000, 
        0.1075338000, 0.1085384000, 0.1095456000, 0.1105555000, 0.1115677000, 
        0.1125823000, 0.1135991000, 0.1146180000, 0.1156389000, 0.1166616000, 
        0.1176862000, 0.1187123000, 0.1197401000, 0.1207692000, 0.1217997000, 
        0.1228315000, 0.1238644000, 0.1248983000, 0.1259332000, 0.1269690000, 
        0.1280055000, 0.1290427000, 0.1300805000, 0.1311188000, 0.1321575000, 
        0.1331966000, 0.1342359000, 0.1352754000, 0.1363150000, 0.1373546000, 
        0.1383942000, 0.1394337000, 0.1404730000, 0.1415120000, 0.1425508000, 
        0.1435891000, 0.1446271000, 0.1456645000, 0.1467014000, 0.1477377000, 
        0.1487733000, 0.1498082000, 0.1508423000, 0.1518756000, 0.1529081000, 
        0.1539396000, 0.1549702000, 0.1559998000, 0.1570284000, 0.1580558000, 
        0.1590822000, 0.1601074000, 0.1611314000, 0.1621541000, 0.1631757000, 
        0.1641959000, 0.1652148000, 0.1662323000, 0.1672484000, 0.1682632000, 
        0.1692765000, 0.1702883000, 0.1712986000, 0.1723075000, 0.1733147000, 
        0.1743205000, 0.1753246000, 0.1763271000, 0.1773281000, 0.1783273000, 
        0.1793250000, 0.1803209000, 0.1813152000, 0.1823078000, 0.1832986000, 
        0.1842877000, 0.1852751000, 0.1862607000, 0.1872445000, 0.1882266000, 
        0.1892068000, 0.1901853000, 0.1911619000, 0.1921368000, 0.1931098000, 
        0.1940809000, 0.1950503000, 0.1960177000, 0.1969834000, 0.1979471000, 
        0.1989090000, 0.1998690000, 0.2008271000, 0.2017834000, 0.2027378000, 
        0.2036903000, 0.2046408000, 0.2055895000, 0.2065363000, 0.2074812000, 
        0.2084242000, 0.2093653000, 0.2103045000, 0.2112418000, 0.2121772000, 
        0.2131106000, 0.2140422000, 0.2149719000, 0.2158996000, 0.2168255000, 
        0.2177494000, 0.2186715000, 0.2195916000, 0.2205098000, 0.2214262000, 
        0.2223406000, 0.2232532000, 0.2241638000, 0.2250726000, 0.2259795000, 
        0.2268844000, 0.2277875000, 0.2286888000, 0.2295881000, 0.2304856000, 
        0.2313812000, 0.2322749000, 0.2331668000, 0.2340568000, 0.2349449000, 
        0.2358312000, 0.2367157000, 0.2375983000, 0.2384790000, 0.2393580000, 
        0.2402351000, 0.2411103000, 0.2419837000, 0.2428554000, 0.2437252000, 
        0.2445931000, 0.2454593000, 0.2463237000, 0.2471863000, 0.2480471000, 
        0.2489061000, 0.2497633000, 0.2506187000, 0.2514724000, 0.2523243000, 
        0.2531744000, 0.2540228000, 0.2548694000, 0.2557143000, 0.2565574000, 
        0.2573988000, 0.2582385000, 0.2590764000, 0.2599126000, 0.2607471000, 
        0.2615799000, 0.2624110000, 0.2632404000, 0.2640681000, 0.2648941000, 
        0.2657184000, 0.2665410000, 0.2673620000, 0.2681813000, 0.2689989000, 
        0.2698149000, 0.2706292000, 0.2714419000, 0.2722529000, 0.2730623000, 
        0.2738701000, 0.2746762000, 0.2754807000, 0.2762836000, 0.2770849000, 
        0.2778846000, 0.2786827000, 0.2794792000, 0.2802741000, 0.2810675000, 
        0.2818592000, 0.2826494000, 0.2834380000, 0.2842251000, 0.2850106000, 
        0.2857946000, 0.2865770000, 0.2873578000, 0.2881372000, 0.2889150000, 
        0.2896913000, 0.2904661000, 0.2912393000, 0.2920111000, 0.2927813000, 
        0.2935501000, 0.2943173000, 0.2950831000, 0.2958474000, 0.2966102000, 
        0.2973715000, 0.2981314000, 0.2988898000, 0.2996468000, 0.3004023000, 
        0.3011564000, 0.3019090000, 0.3026602000, 0.3034099000, 0.3041583000, 
        0.3049052000, 0.3056507000, 0.3063948000, 0.3071374000, 0.3078787000, 
        0.3086186000, 0.3093571000, 0.3100942000, 0.3108299000, 0.3115643000, 
        0.3122973000, 0.3130289000, 0.3137591000, 0.3144880000, 0.3152155000, 
        0.3159417000, 0.3166666000, 0.3173901000, 0.3181122000, 0.3188331000, 
        0.3195526000, 0.3202708000, 0.3209877000, 0.3217033000, 0.3224175000, 
        0.3231305000, 0.3238421000, 0.3245525000, 0.3252616000, 0.3259694000, 
        0.3266759000, 0.3273811000, 0.3280851000, 0.3287878000, 0.3294893000, 
        0.3301894000, 0.3308884000, 0.3315861000, 0.3322825000, 0.3329777000, 
        0.3336716000, 0.3343644000, 0.3350559000, 0.3357461000, 0.3364352000, 
        0.3371230000, 0.3378096000, 0.3384951000, 0.3391793000, 0.3398623000, 
        0.3405441000, 0.3412247000, 0.3419042000, 0.3425824000, 0.3432595000, 
        0.3439354000, 0.3446101000, 0.3452837000, 0.3459561000, 0.3466273000, 
        0.3472974000, 0.3479663000, 0.3486341000, 0.3493008000, 0.3499663000, 
        0.3506306000, 0.3512939000, 0.3519560000, 0.3526169000, 0.3532768000, 
        0.3539355000, 0.3545931000, 0.3552497000, 0.3559051000, 0.3565594000, 
        0.3572126000, 0.3578647000, 0.3585157000, 0.3591656000, 0.3598144000, 
        0.3604622000, 0.3611089000, 0.3617545000, 0.3623990000, 0.3630425000, 
        0.3636849000, 0.3643263000, 0.3649665000, 0.3656058000, 0.3662440000, 
        0.3668811000, 0.3675172000, 0.3681523000, 0.3687863000, 0.3694193000, 
        0.3700512000, 0.3706822000, 0.3713121000, 0.3719409000, 0.3725688000, 
        0.3731957000, 0.3738215000, 0.3744464000, 0.3750702000, 0.3756930000, 
        0.3763149000, 0.3769357000, 0.3775556000, 0.3781745000, 0.3787924000, 
        0.3794093000, 0.3800252000, 0.3806401000, 0.3812541000, 0.3818671000, 
        0.3824792000, 0.3830902000, 0.3837004000, 0.3843095000, 0.3849178000, 
        0.3855250000, 0.3861313000, 0.3867367000, 0.3873411000, 0.3879446000, 
        0.3885472000, 0.3891488000, 0.3897495000, 0.3903492000, 0.3909481000, 
        0.3915460000, 0.3921430000, 0.3927391000, 0.3933343000, 0.3939285000, 
        0.3945219000, 0.3951143000, 0.3957059000, 0.3962966000, 0.3968863000, 
        0.3974752000, 0.3980631000, 0.3986502000, 0.3992364000, 0.3998218000, 
        0.4004062000, 0.4009898000, 0.4015725000, 0.4021543000, 0.4027352000, 
        0.4033153000, 0.4038946000, 0.4044729000, 0.4050504000, 0.4056271000, 
        0.4062029000, 0.4067778000, 0.4073519000, 0.4079252000, 0.4084976000, 
        0.4090692000, 0.4096399000, 0.4102098000, 0.4107789000, 0.4113472000, 
        0.4119146000, 0.4124812000, 0.4130470000, 0.4136119000, 0.4141760000, 
        0.4147394000, 0.4153019000, 0.4158636000, 0.4164245000, 0.4169846000, 
        0.4175439000, 0.4181024000, 0.4186601000, 0.4192170000, 0.4197731000, 
        0.4203284000, 0.4208830000, 0.4214367000, 0.4219897000, 0.4225419000, 
        0.4230933000, 0.4236439000, 0.4241938000, 0.4247429000, 0.4252912000, 
        0.4258388000, 0.4263856000, 0.4269317000, 0.4274769000, 0.4280215000, 
        0.4285652000, 0.4291083000, 0.4296505000, 0.4301921000, 0.4307329000, 
        0.4312729000, 0.4318122000, 0.4323508000, 0.4328886000, 0.4334257000, 
        0.4339621000, 0.4344977000, 0.4350326000, 0.4355668000, 0.4361003000, 
        0.4366331000, 0.4371651000, 0.4376964000, 0.4382270000, 0.4387569000, 
        0.4392861000, 0.4398146000, 0.4403424000, 0.4408695000, 0.4413959000, 
        0.4419216000, 0.4424466000, 0.4429709000, 0.4434945000, 0.4440174000, 
        0.4445396000, 0.4450612000, 0.4455821000, 0.4461023000, 0.4466218000, 
        0.4471407000, 0.4476589000, 0.4481764000, 0.4486932000, 0.4492094000, 
        0.4497249000, 0.4502398000, 0.4507540000, 0.4512675000, 0.4517804000, 
        0.4522926000, 0.4528042000, 0.4533152000, 0.4538255000, 0.4543351000, 
        0.4548442000, 0.4553525000, 0.4558603000, 0.4563674000, 0.4568739000, 
        0.4573797000, 0.4578849000, 0.4583895000, 0.4588935000, 0.4593969000, 
        0.4598996000, 0.4604018000, 0.4609033000, 0.4614042000, 0.4619045000, 
        0.4624042000, 0.4629032000, 0.4634017000, 0.4638996000, 0.4643969000, 
        0.4648936000, 0.4653897000, 0.4658852000, 0.4663802000, 0.4668745000, 
        0.4673683000, 0.4678614000, 0.4683541000, 0.4688461000, 0.4693375000, 
        0.4698284000, 0.4703188000, 0.4708085000, 0.4712977000, 0.4717863000, 
        0.4722744000, 0.4727619000, 0.4732489000, 0.4737353000, 0.4742212000, 
        0.4747065000, 0.4751913000, 0.4756756000, 0.4761593000, 0.4766424000, 
        0.4771251000, 0.4776072000, 0.4780888000, 0.4785699000, 0.4790504000, 
        0.4795304000, 0.4800099000, 0.4804889000, 0.4809674000, 0.4814454000, 
        0.4819229000, 0.4823999000, 0.4828764000, 0.4833524000, 0.4838279000, 
        0.4843029000, 0.4847774000, 0.4852515000, 0.4857250000, 0.4861981000, 
        0.4866708000, 0.4871429000, 0.4876146000, 0.4880858000, 0.4885566000, 
        0.4890269000, 0.4894968000, 0.4899662000, 0.4904351000, 0.4909037000, 
        0.4913718000, 0.4918394000, 0.4923067000, 0.4927735000, 0.4932398000, 
        0.4937058000, 0.4941713000, 0.4946365000, 0.4951012000, 0.4955655000, 
        0.4960295000, 0.4964930000, 0.4969562000, 0.4974189000, 0.4978813000, 
        0.4983434000, 0.4988050000, 0.4992663000, 0.4997272000, 0.5001878000, 
        0.5006480000, 0.5011079000, 0.5015674000, 0.5020266000, 0.5024855000, 
        0.5029440000, 0.5034023000, 0.5038602000, 0.5043178000, 0.5047751000, 
        0.5052321000, 0.5056889000, 0.5061453000, 0.5066015000, 0.5070574000, 
        0.5075131000, 0.5079685000, 0.5084237000, 0.5088786000, 0.5093333000, 
        0.5097878000, 0.5102420000, 0.5106961000, 0.5111499000, 0.5116036000, 
        0.5120571000, 0.5125104000, 0.5129636000, 0.5134166000, 0.5138694000, 
        0.5143221000, 0.5147747000, 0.5152272000, 0.5156796000, 0.5161319000, 
        0.5165842000, 0.5170363000, 0.5174884000, 0.5179405000, 0.5183925000, 
        0.5188445000, 0.5192965000, 0.5197485000, 0.5202005000, 0.5206526000, 
        0.5211048000, 0.5215570000, 0.5220092000, 0.5224616000, 0.5229141000, 
        0.5233668000, 0.5238195000, 0.5242725000, 0.5247257000, 0.5251790000, 
        0.5256326000, 0.5260864000, 0.5265405000, 0.5269949000, 0.5274496000, 
        0.5279047000, 0.5283601000, 0.5288159000, 0.5292721000, 0.5297287000, 
        0.5301858000, 0.5306434000, 0.5311015000, 0.5315601000, 0.5320194000, 
        0.5324792000, 0.5329396000, 0.5334008000, 0.5338626000, 0.5343252000, 
        0.5347885000, 0.5352526000, 0.5357176000, 0.5361834000, 0.5366502000, 
        0.5371179000, 0.5375866000, 0.5380563000, 0.5385271000, 0.5389989000, 
        0.5394720000, 0.5399462000, 0.5404217000, 0.5408984000, 0.5413764000, 
        0.5418558000, 0.5423367000, 0.5428189000, 0.5433026000, 0.5437879000, 
        0.5442747000, 0.5447631000, 0.5452532000, 0.5457450000, 0.5462385000, 
        0.5467337000, 0.5472307000, 0.5477296000, 0.5482303000, 0.5487328000, 
        0.5492373000, 0.5497437000, 0.5502520000, 0.5507623000, 0.5512745000, 
        0.5517887000, 0.5523048000, 0.5528230000, 0.5533430000, 0.5538651000, 
        0.5543890000, 0.5549149000, 0.5554426000, 0.5559722000, 0.5565036000, 
        0.5570367000, 0.5575716000, 0.5581082000, 0.5586464000, 0.5591862000, 
        0.5597274000, 0.5602702000, 0.5608143000, 0.5613597000, 0.5619063000, 
        0.5624541000, 0.5630030000, 0.5635529000, 0.5641037000, 0.5646553000, 
        0.5652078000, 0.5657609000, 0.5663146000, 0.5668689000, 0.5674236000, 
        0.5679787000, 0.5685341000, 0.5690898000, 0.5696456000, 0.5702015000, 
        0.5707575000, 0.5713134000, 0.5718692000, 0.5724249000, 0.5729803000, 
        0.5735355000, 0.5740903000, 0.5746448000, 0.5751989000, 0.5757525000, 
        0.5763056000, 0.5768581000, 0.5774101000, 0.5779615000, 0.5785121000, 
        0.5790622000, 0.5796114000, 0.5801600000, 0.5807078000, 0.5812548000, 
        0.5818009000, 0.5823462000, 0.5828907000, 0.5834343000, 0.5839770000, 
        0.5845188000, 0.5850597000, 0.5855996000, 0.5861386000, 0.5866767000, 
        0.5872137000, 0.5877498000, 0.5882850000, 0.5888191000, 0.5893522000, 
        0.5898844000, 0.5904155000, 0.5909457000, 0.5914748000, 0.5920029000, 
        0.5925300000, 0.5930560000, 0.5935811000, 0.5941051000, 0.5946282000, 
        0.5951502000, 0.5956711000, 0.5961911000, 0.5967100000, 0.5972280000, 
        0.5977449000, 0.5982608000, 0.5987757000, 0.5992896000, 0.5998025000, 
        0.6003144000, 0.6008253000, 0.6013352000, 0.6018441000, 0.6023520000, 
        0.6028589000, 0.6033649000, 0.6038699000, 0.6043739000, 0.6048769000, 
        0.6053790000, 0.6058801000, 0.6063803000, 0.6068795000, 0.6073777000, 
        0.6078750000, 0.6083714000, 0.6088668000, 0.6093613000, 0.6098549000, 
        0.6103476000, 0.6108393000, 0.6113301000, 0.6118201000, 0.6123091000, 
        0.6127972000, 0.6132844000, 0.6137707000, 0.6142562000, 0.6147407000, 
        0.6152244000, 0.6157072000, 0.6161892000, 0.6166703000, 0.6171505000, 
        0.6176298000, 0.6181084000, 0.6185860000, 0.6190629000, 0.6195388000, 
        0.6200140000, 0.6204883000, 0.6209618000, 0.6214345000, 0.6219064000, 
        0.6223774000, 0.6228476000, 0.6233171000, 0.6237857000, 0.6242536000, 
        0.6247206000, 0.6251869000, 0.6256523000, 0.6261170000, 0.6265809000, 
        0.6270441000, 0.6275064000, 0.6279681000, 0.6284289000, 0.6288890000, 
        0.6293483000, 0.6298069000, 0.6302647000, 0.6307218000, 0.6311782000, 
        0.6316338000, 0.6320887000, 0.6325428000, 0.6329963000, 0.6334490000, 
        0.6339009000, 0.6343522000, 0.6348028000, 0.6352526000, 0.6357018000, 
        0.6361502000, 0.6365979000, 0.6370450000, 0.6374913000, 0.6379370000, 
        0.6383820000, 0.6388262000, 0.6392699000, 0.6397128000, 0.6401550000, 
        0.6405966000, 0.6410375000, 0.6414778000, 0.6419174000, 0.6423563000, 
        0.6427946000, 0.6432322000, 0.6436692000, 0.6441055000, 0.6445412000, 
        0.6449762000, 0.6454106000, 0.6458443000, 0.6462774000, 0.6467099000, 
        0.6471418000, 0.6475730000, 0.6480036000, 0.6484336000, 0.6488629000, 
        0.6492917000, 0.6497198000, 0.6501473000, 0.6505742000, 0.6510005000, 
        0.6514262000, 0.6518513000, 0.6522758000, 0.6526997000, 0.6531230000, 
        0.6535457000, 0.6539678000, 0.6543893000, 0.6548103000, 0.6552306000, 
        0.6556504000, 0.6560696000, 0.6564882000, 0.6569062000, 0.6573237000, 
        0.6577406000, 0.6581569000, 0.6585727000, 0.6589879000, 0.6594025000, 
        0.6598166000, 0.6602301000, 0.6606431000, 0.6610555000, 0.6614673000, 
        0.6618787000, 0.6622894000, 0.6626996000, 0.6631093000, 0.6635184000, 
        0.6639270000, 0.6643350000, 0.6647425000, 0.6651495000, 0.6655560000, 
        0.6659619000, 0.6663672000, 0.6667721000, 0.6671764000, 0.6675802000, 
        0.6679835000, 0.6683863000, 0.6687885000, 0.6691902000, 0.6695914000, 
        0.6699921000, 0.6703923000, 0.6707920000, 0.6711912000, 0.6715898000, 
        0.6719880000, 0.6723856000, 0.6727828000, 0.6731794000, 0.6735756000, 
        0.6739712000, 0.6743664000, 0.6747610000, 0.6751552000, 0.6755489000, 
        0.6759421000, 0.6763348000, 0.6767270000, 0.6771187000, 0.6775100000, 
        0.6779008000, 0.6782910000, 0.6786809000, 0.6790702000, 0.6794590000, 
        0.6798474000, 0.6802353000, 0.6806228000, 0.6810097000, 0.6813962000, 
        0.6817823000, 0.6821678000, 0.6825530000, 0.6829376000, 0.6833218000, 
        0.6837055000, 0.6840888000, 0.6844716000, 0.6848539000, 0.6852358000, 
        0.6856172000, 0.6859982000, 0.6863788000, 0.6867588000, 0.6871385000, 
        0.6875177000, 0.6878964000, 0.6882747000, 0.6886526000, 0.6890300000, 
        0.6894070000, 0.6897835000, 0.6901596000, 0.6905353000, 0.6909105000, 
        0.6912853000, 0.6916596000, 0.6920336000, 0.6924071000, 0.6927801000, 
        0.6931528000, 0.6935250000, 0.6938967000, 0.6942681000, 0.6946390000, 
        0.6950095000, 0.6953796000, 0.6957493000, 0.6961185000, 0.6964874000, 
        0.6968558000, 0.6972238000, 0.6975914000, 0.6979585000, 0.6983253000, 
        0.6986916000, 0.6990576000, 0.6994231000, 0.6997882000, 0.7001529000, 
        0.7005172000, 0.7008811000, 0.7012446000, 0.7016077000, 0.7019704000, 
        0.7023326000, 0.7026945000, 0.7030560000, 0.7034171000, 0.7037778000, 
        0.7041381000, 0.7044980000, 0.7048575000, 0.7052166000, 0.7055753000, 
        0.7059336000, 0.7062916000, 0.7066491000, 0.7070063000, 0.7073630000, 
        0.7077194000, 0.7080754000, 0.7084310000, 0.7087863000, 0.7091411000, 
        0.7094956000, 0.7098497000, 0.7102034000, 0.7105567000, 0.7109096000, 
        0.7112622000, 0.7116144000, 0.7119662000, 0.7123177000, 0.7126687000, 
        0.7130194000, 0.7133697000, 0.7137197000, 0.7140693000, 0.7144185000, 
        0.7147673000, 0.7151158000, 0.7154639000, 0.7158117000, 0.7161591000, 
        0.7165061000, 0.7168527000, 0.7171990000, 0.7175450000, 0.7178905000, 
        0.7182357000, 0.7185806000, 0.7189251000, 0.7192692000, 0.7196130000, 
        0.7199564000, 0.7202995000, 0.7206422000, 0.7209845000, 0.7213266000, 
        0.7216682000, 0.7220095000, 0.7223505000, 0.7226911000, 0.7230313000, 
        0.7233712000, 0.7237108000, 0.7240500000, 0.7243889000, 0.7247274000, 
        0.7250656000, 0.7254034000, 0.7257409000, 0.7260780000, 0.7264148000, 
        0.7267513000, 0.7270874000, 0.7274232000, 0.7277587000, 0.7280938000, 
        0.7284286000, 0.7287630000, 0.7290971000, 0.7294309000, 0.7297643000, 
        0.7300974000, 0.7304302000, 0.7307626000, 0.7310947000, 0.7314265000, 
        0.7317579000, 0.7320891000, 0.7324198000, 0.7327503000, 0.7330804000, 
        0.7334102000, 0.7337397000, 0.7340689000, 0.7343977000, 0.7347262000, 
        0.7350544000, 0.7353823000, 0.7357098000, 0.7360370000, 0.7363639000, 
        0.7366905000, 0.7370168000, 0.7373427000, 0.7376683000, 0.7379936000, 
        0.7383186000, 0.7386433000, 0.7389677000, 0.7392917000, 0.7396155000, 
        0.7399389000, 0.7402620000, 0.7405848000, 0.7409073000, 0.7412294000, 
        0.7415513000, 0.7418728000, 0.7421941000, 0.7425150000, 0.7428357000, 
        0.7431560000, 0.7434760000, 0.7437957000, 0.7441151000, 0.7444342000, 
        0.7447530000, 0.7450715000, 0.7453897000, 0.7457076000, 0.7460251000, 
        0.7463424000, 0.7466594000, 0.7469761000, 0.7472925000, 0.7476085000, 
        0.7479243000, 0.7482398000, 0.7485550000, 0.7488699000, 0.7491845000, 
        0.7494988000, 0.7498128000, 0.7501265000, 0.7504399000, 0.7507531000, 
        0.7510659000, 0.7513784000, 0.7516907000, 0.7520027000, 0.7523143000, 
        0.7526257000, 0.7529368000, 0.7532476000, 0.7535581000, 0.7538683000, 
        0.7541783000, 0.7544879000, 0.7547973000, 0.7551064000, 0.7554152000, 
        0.7557237000, 0.7560319000, 0.7563398000, 0.7566475000, 0.7569549000, 
        0.7572620000, 0.7575688000, 0.7578753000, 0.7581816000, 0.7584875000, 
        0.7587932000, 0.7590986000, 0.7594038000, 0.7597086000, 0.7600132000, 
        0.7603175000, 0.7606215000, 0.7609253000, 0.7612288000, 0.7615320000, 
        0.7618349000, 0.7621376000, 0.7624399000, 0.7627420000, 0.7630439000, 
        0.7633454000, 0.7636467000, 0.7639477000, 0.7642485000, 0.7645490000, 
        0.7648492000, 0.7651491000, 0.7654488000, 0.7657482000, 0.7660473000, 
        0.7663462000, 0.7666448000, 0.7669431000, 0.7672412000, 0.7675390000, 
        0.7678366000, 0.7681338000, 0.7684309000, 0.7687276000, 0.7690241000, 
        0.7693203000, 0.7696163000, 0.7699120000, 0.7702074000, 0.7705026000, 
        0.7707975000, 0.7710922000, 0.7713866000, 0.7716807000, 0.7719746000, 
        0.7722682000, 0.7725616000, 0.7728547000, 0.7731475000, 0.7734401000, 
        0.7737325000, 0.7740245000, 0.7743164000, 0.7746079000, 0.7748993000, 
        0.7751903000, 0.7754811000, 0.7757717000, 0.7760620000, 0.7763521000, 
        0.7766419000, 0.7769314000, 0.7772207000, 0.7775098000, 0.7777986000, 
        0.7780871000, 0.7783754000, 0.7786635000, 0.7789513000, 0.7792389000, 
        0.7795262000, 0.7798132000, 0.7801000000, 0.7803866000, 0.7806729000, 
        0.7809590000, 0.7812449000, 0.7815304000, 0.7818158000, 0.7821009000, 
        0.7823858000, 0.7826704000, 0.7829547000, 0.7832389000, 0.7835228000, 
        0.7838064000, 0.7840898000, 0.7843730000, 0.7846559000, 0.7849386000, 
        0.7852210000, 0.7855032000, 0.7857852000, 0.7860669000, 0.7863484000, 
        0.7866297000, 0.7869107000, 0.7871915000, 0.7874720000, 0.7877523000, 
        0.7880324000, 0.7883122000, 0.7885919000, 0.7888712000, 0.7891504000, 
        0.7894293000, 0.7897079000, 0.7899864000, 0.7902646000, 0.7905425000, 
        0.7908203000, 0.7910978000, 0.7913750000, 0.7916521000, 0.7919289000, 
        0.7922055000, 0.7924818000, 0.7927579000, 0.7930338000, 0.7933095000, 
        0.7935849000, 0.7938601000, 0.7941351000, 0.7944099000, 0.7946844000, 
        0.7949587000, 0.7952328000, 0.7955066000, 0.7957803000, 0.7960537000, 
        0.7963268000, 0.7965998000, 0.7968725000, 0.7971450000, 0.7974173000, 
        0.7976893000, 0.7979612000, 0.7982328000, 0.7985041000, 0.7987753000, 
        0.7990463000, 0.7993170000, 0.7995875000, 0.7998578000, 0.8001278000, 
        0.8003977000, 0.8006673000, 0.8009367000, 0.8012059000, 0.8014748000, 
        0.8017436000, 0.8020121000, 0.8022804000, 0.8025485000, 0.8028164000, 
        0.8030841000, 0.8033515000, 0.8036187000, 0.8038857000, 0.8041525000, 
        0.8044191000, 0.8046855000, 0.8049517000, 0.8052176000, 0.8054833000, 
        0.8057488000, 0.8060141000, 0.8062792000, 0.8065441000, 0.8068088000, 
        0.8070732000, 0.8073375000, 0.8076015000, 0.8078653000, 0.8081289000, 
        0.8083923000, 0.8086555000, 0.8089185000, 0.8091813000, 0.8094438000, 
        0.8097062000, 0.8099683000, 0.8102303000, 0.8104920000, 0.8107535000, 
        0.8110148000, 0.8112760000, 0.8115369000, 0.8117976000, 0.8120580000, 
        0.8123183000, 0.8125784000, 0.8128383000, 0.8130980000, 0.8133574000, 
        0.8136167000, 0.8138757000, 0.8141346000, 0.8143933000, 0.8146517000, 
        0.8149099000, 0.8151680000, 0.8154258000, 0.8156835000, 0.8159409000, 
        0.8161982000, 0.8164552000, 0.8167120000, 0.8169687000, 0.8172251000, 
        0.8174813000, 0.8177374000, 0.8179932000, 0.8182488000, 0.8185043000, 
        0.8187595000, 0.8190146000, 0.8192694000, 0.8195241000, 0.8197785000, 
        0.8200328000, 0.8202868000, 0.8205407000, 0.8207944000, 0.8210478000, 
        0.8213011000, 0.8215542000, 0.8218071000, 0.8220598000, 0.8223122000, 
        0.8225645000, 0.8228167000, 0.8230686000, 0.8233203000, 0.8235718000, 
        0.8238232000, 0.8240743000, 0.8243252000, 0.8245760000, 0.8248266000, 
        0.8250769000, 0.8253271000, 0.8255771000, 0.8258269000, 0.8260765000, 
        0.8263259000, 0.8265752000, 0.8268242000, 0.8270731000, 0.8273217000, 
        0.8275702000, 0.8278185000, 0.8280666000, 0.8283145000, 0.8285622000, 
        0.8288097000, 0.8290571000, 0.8293043000, 0.8295512000, 0.8297980000, 
        0.8300446000, 0.8302910000, 0.8305373000, 0.8307833000, 0.8310292000, 
        0.8312748000, 0.8315203000, 0.8317656000, 0.8320107000, 0.8322557000, 
        0.8325004000, 0.8327450000, 0.8329894000, 0.8332336000, 0.8334776000, 
        0.8337214000, 0.8339651000, 0.8342085000, 0.8344518000, 0.8346949000 
    };
    double width_a1=a1_width[iii];
    a[0] = 1;
    a[1] = 0; 
    b[0] = mass2-sa;
    //b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r2,l); 
    b[1] = -mass*width_a1; 
    Com_Divide(a,b,prop);
}
void EvtD0ToKSKppimpi0::propagatorRBW_K1(double mass2, double mass, double width, double sa, double sb, double sc, double r2, int l, double prop[2])
{
    double a[2], b[2];
    int iii=int(sa*1000/2)-300;//(4000-0)/2000=2
    if(iii<0){ iii=0; }
    double K1270_width[2000]={//1700 bin, s = (0.6, 4)
        0.0000000014, 0.0000000052, 0.0000000127, 0.0000000252, 0.0000000442, 
        0.0000000709, 0.0000001067, 0.0000001528, 0.0000002108, 0.0000002819, 
        0.0000003675, 0.0000004689, 0.0000005876, 0.0000007249, 0.0000008822, 
        0.0000010609, 0.0000012624, 0.0000014882, 0.0000017396, 0.0000020181, 
        0.0000023251, 0.0000026620, 0.0000030303, 0.0000034314, 0.0000038668, 
        0.0000043379, 0.0000048463, 0.0000053933, 0.0000059805, 0.0000066094, 
        0.0000072815, 0.0000079983, 0.0000087612, 0.0000095719, 0.0000104318, 
        0.0000113424, 0.0000123054, 0.0000133223, 0.0000143946, 0.0000155239, 
        0.0000167117, 0.0000179597, 0.0000192695, 0.0000206426, 0.0000220807, 
        0.0000235853, 0.0000251581, 0.0000268008, 0.0000285149, 0.0000303022, 
        0.0000321643, 0.0000341028, 0.0000361195, 0.0000382160, 0.0000403941, 
        0.0000426554, 0.0000450017, 0.0000474347, 0.0000499561, 0.0000525678, 
        0.0000552715, 0.0000580689, 0.0000609619, 0.0000639522, 0.0000670417, 
        0.0000702323, 0.0000735256, 0.0000769237, 0.0000804283, 0.0000840414, 
        0.0000877648, 0.0000916004, 0.0000955503, 0.0000996163, 0.0001038003, 
        0.0001081045, 0.0001125306, 0.0001170808, 0.0001217571, 0.0001265615, 
        0.0001314961, 0.0001365630, 0.0001417641, 0.0001471018, 0.0001525781, 
        0.0001581951, 0.0001639551, 0.0001698602, 0.0001759127, 0.0001821148, 
        0.0001884687, 0.0001949769, 0.0002016415, 0.0002084650, 0.0002154497, 
        0.0002225981, 0.0002299124, 0.0002373952, 0.0002450491, 0.0002528764, 
        0.0002608798, 0.0002690618, 0.0002774250, 0.0002859721, 0.0002947058, 
        0.0003036288, 0.0003127438, 0.0003220537, 0.0003315613, 0.0003412695, 
        0.0003511812, 0.0003612993, 0.0003716269, 0.0003821671, 0.0003929230, 
        0.0004038976, 0.0004150944, 0.0004265164, 0.0004381671, 0.0004500499, 
        0.0004621681, 0.0004745254, 0.0004871253, 0.0004999714, 0.0005130674, 
        0.0005264172, 0.0005400246, 0.0005538936, 0.0005680282, 0.0005824325, 
        0.0005971108, 0.0006120672, 0.0006273063, 0.0006428325, 0.0006586504, 
        0.0006747647, 0.0006911803, 0.0007079021, 0.0007249351, 0.0007422845, 
        0.0007599558, 0.0007779542, 0.0007962855, 0.0008149555, 0.0008339699, 
        0.0008533351, 0.0008730571, 0.0008931425, 0.0009135979, 0.0009344302, 
        0.0009556463, 0.0009772537, 0.0009992598, 0.0010216720, 0.0010444990, 
        0.0010677490, 0.0010914300, 0.0011155510, 0.0011401210, 0.0011651510, 
        0.0011906480, 0.0012166250, 0.0012430910, 0.0012700570, 0.0012975350, 
        0.0013255370, 0.0013540740, 0.0013831600, 0.0014128090, 0.0014430330, 
        0.0014738470, 0.0015052670, 0.0015373090, 0.0015699870, 0.0016033210, 
        0.0016373270, 0.0016720240, 0.0017074330, 0.0017435730, 0.0017804660, 
        0.0018181340, 0.0018566010, 0.0018958920, 0.0019360320, 0.0019770490, 
        0.0020189710, 0.0020618280, 0.0021056510, 0.0021504740, 0.0021963300, 
        0.0022432570, 0.0022912920, 0.0023404760, 0.0023908500, 0.0024424590, 
        0.0024953480, 0.0025495680, 0.0026051670, 0.0026622000, 0.0027207220, 
        0.0027807900, 0.0028424670, 0.0029058130, 0.0029708960, 0.0030377810, 
        0.0031065390, 0.0031772420, 0.0032499620, 0.0033247750, 0.0034017560, 
        0.0034809810, 0.0035625250, 0.0036464650, 0.0037328720, 0.0038218190, 
        0.0039133710, 0.0040075910, 0.0041045360, 0.0042042540, 0.0043067850, 
        0.0044121600, 0.0045203990, 0.0046315090, 0.0047454840, 0.0048623040, 
        0.0049819360, 0.0051043310, 0.0052294280, 0.0053571500, 0.0054874080, 
        0.0056201030, 0.0057551230, 0.0058923500, 0.0060316590, 0.0061729200, 
        0.0063159990, 0.0064607640, 0.0066070800, 0.0067548180, 0.0069038480, 
        0.0070540500, 0.0072053060, 0.0073575050, 0.0075105450, 0.0076643270, 
        0.0078187630, 0.0079737720, 0.0081292780, 0.0082852130, 0.0084415170, 
        0.0085981350, 0.0087550200, 0.0089121270, 0.0090694200, 0.0092268660, 
        0.0093844380, 0.0095421120, 0.0096998670, 0.0098576870, 0.0100155600, 
        0.0101734700, 0.0103314200, 0.0104894000, 0.0106474000, 0.0108054200, 
        0.0109634700, 0.0111215500, 0.0112796500, 0.0114377900, 0.0115959700, 
        0.0117541900, 0.0119124600, 0.0120707900, 0.0122292000, 0.0123876800, 
        0.0125462500, 0.0127049200, 0.0128637100, 0.0130226100, 0.0131816400, 
        0.0133408200, 0.0135001600, 0.0136596700, 0.0138193500, 0.0139792400, 
        0.0141393300, 0.0142996400, 0.0144601800, 0.0146209700, 0.0147820200, 
        0.0149433400, 0.0151049500, 0.0152668600, 0.0154290800, 0.0155916300, 
        0.0157545200, 0.0159177700, 0.0160813900, 0.0162453900, 0.0164097800, 
        0.0165745900, 0.0167398200, 0.0169054800, 0.0170716000, 0.0172381900, 
        0.0174052600, 0.0175728200, 0.0177408900, 0.0179094800, 0.0180786100, 
        0.0182483000, 0.0184185500, 0.0185893800, 0.0187608000, 0.0189328400, 
        0.0191054900, 0.0192787900, 0.0194527400, 0.0196273600, 0.0198026700, 
        0.0199786700, 0.0201553900, 0.0203328300, 0.0205110200, 0.0206899700, 
        0.0208696900, 0.0210502100, 0.0212315300, 0.0214136700, 0.0215966500, 
        0.0217804800, 0.0219651900, 0.0221507800, 0.0223372700, 0.0225246900, 
        0.0227130400, 0.0229023400, 0.0230926100, 0.0232838700, 0.0234761300, 
        0.0236694200, 0.0238637400, 0.0240591200, 0.0242555800, 0.0244531300, 
        0.0246517900, 0.0248515800, 0.0250525200, 0.0252546300, 0.0254579200, 
        0.0256624200, 0.0258681400, 0.0260751100, 0.0262833500, 0.0264928700, 
        0.0267036900, 0.0269158400, 0.0271293400, 0.0273442100, 0.0275604600, 
        0.0277781300, 0.0279972300, 0.0282177800, 0.0284398100, 0.0286633400, 
        0.0288884000, 0.0291149900, 0.0293431600, 0.0295729200, 0.0298042900, 
        0.0300373100, 0.0302719900, 0.0305083600, 0.0307464500, 0.0309862700, 
        0.0312278600, 0.0314712500, 0.0317164500, 0.0319635000, 0.0322124200, 
        0.0324632400, 0.0327159900, 0.0329706900, 0.0332273800, 0.0334860800, 
        0.0337468200, 0.0340096300, 0.0342745400, 0.0345415800, 0.0348107800, 
        0.0350821700, 0.0353557800, 0.0356316400, 0.0359097900, 0.0361902600, 
        0.0364730800, 0.0367582700, 0.0370458900, 0.0373359500, 0.0376284900, 
        0.0379235500, 0.0382211500, 0.0385213500, 0.0388241600, 0.0391296300, 
        0.0394377800, 0.0397486700, 0.0400623100, 0.0403787600, 0.0406980400, 
        0.0410202000, 0.0413452600, 0.0416732700, 0.0420042700, 0.0423382900, 
        0.0426753700, 0.0430155500, 0.0433588700, 0.0437053700, 0.0440550800, 
        0.0444080500, 0.0447643100, 0.0451239000, 0.0454868700, 0.0458532400, 
        0.0462230700, 0.0465963900, 0.0469732300, 0.0473536500, 0.0477376700, 
        0.0481253400, 0.0485166900, 0.0489117600, 0.0493106000, 0.0497132400, 
        0.0501197100, 0.0505300600, 0.0509443200, 0.0513625300, 0.0517847200, 
        0.0522109400, 0.0526412100, 0.0530755700, 0.0535140600, 0.0539567000, 
        0.0544035400, 0.0548545900, 0.0553099000, 0.0557694900, 0.0562334000, 
        0.0567016400, 0.0571742500, 0.0576512600, 0.0581326800, 0.0586185500, 
        0.0591088800, 0.0596036900, 0.0601030200, 0.0606068600, 0.0611152500, 
        0.0616281900, 0.0621457100, 0.0626678100, 0.0631945000, 0.0637258000, 
        0.0642617000, 0.0648022200, 0.0653473700, 0.0658971300, 0.0664515200, 
        0.0670105300, 0.0675741500, 0.0681423900, 0.0687152300, 0.0692926600, 
        0.0698746700, 0.0704612600, 0.0710524000, 0.0716480700, 0.0722482600, 
        0.0728529400, 0.0734620900, 0.0740756900, 0.0746937100, 0.0753161100, 
        0.0759428600, 0.0765739400, 0.0772093000, 0.0778489100, 0.0784927300, 
        0.0791407200, 0.0797928300, 0.0804490200, 0.0811092400, 0.0817734600, 
        0.0824416000, 0.0831136300, 0.0837895000, 0.0844691400, 0.0851525100, 
        0.0858395400, 0.0865301900, 0.0872243800, 0.0879220700, 0.0886231900, 
        0.0893276700, 0.0900354700, 0.0907465100, 0.0914607200, 0.0921780600, 
        0.0928984400, 0.0936218200, 0.0943481100, 0.0950772500, 0.0958091900, 
        0.1023867000, 0.1055505000, 0.1081475000, 0.1104523000, 0.1125711000, 
        0.1145581000, 0.1164458000, 0.1182550000, 0.1200004000, 0.1216926000, 
        0.1233395000, 0.1249474000, 0.1265212000, 0.1280649000, 0.1295817000, 
        0.1310743000, 0.1325452000, 0.1339961000, 0.1354288000, 0.1368448000, 
        0.1382454000, 0.1396316000, 0.1410044000, 0.1423647000, 0.1437133000, 
        0.1450510000, 0.1463782000, 0.1476956000, 0.1490038000, 0.1503031000, 
        0.1515941000, 0.1528770000, 0.1541523000, 0.1554203000, 0.1566812000, 
        0.1579354000, 0.1591832000, 0.1604247000, 0.1616602000, 0.1628899000, 
        0.1641140000, 0.1653326000, 0.1665460000, 0.1677543000, 0.1689576000, 
        0.1701561000, 0.1713499000, 0.1725391000, 0.1737239000, 0.1749043000, 
        0.1760805000, 0.1772526000, 0.1784206000, 0.1795847000, 0.1807449000, 
        0.1819014000, 0.1830541000, 0.1842031000, 0.1853487000, 0.1864907000, 
        0.1876292000, 0.1887644000, 0.1898963000, 0.1910250000, 0.1921504000, 
        0.1932727000, 0.1943919000, 0.1955080000, 0.1966211000, 0.1977313000, 
        0.1988386000, 0.1999430000, 0.2010446000, 0.2021434000, 0.2032395000, 
        0.2043329000, 0.2054236000, 0.2065117000, 0.2075972000, 0.2086801000, 
        0.2097605000, 0.2108384000, 0.2119139000, 0.2129870000, 0.2140576000, 
        0.2151259000, 0.2161919000, 0.2172555000, 0.2183169000, 0.2193761000, 
        0.2204330000, 0.2214878000, 0.2225404000, 0.2235909000, 0.2246392000, 
        0.2256855000, 0.2267297000, 0.2277719000, 0.2288121000, 0.2298503000, 
        0.2308866000, 0.2319210000, 0.2329534000, 0.2339840000, 0.2350126000, 
        0.2360395000, 0.2370646000, 0.2380878000, 0.2391093000, 0.2401291000, 
        0.2411471000, 0.2421634000, 0.2431781000, 0.2441911000, 0.2452024000, 
        0.2462122000, 0.2472203000, 0.2482269000, 0.2492319000, 0.2502354000, 
        0.2512374000, 0.2522379000, 0.2532369000, 0.2542345000, 0.2552306000, 
        0.2562253000, 0.2572187000, 0.2582106000, 0.2592012000, 0.2601905000, 
        0.2611784000, 0.2621651000, 0.2631504000, 0.2641345000, 0.2651174000, 
        0.2660991000, 0.2670795000, 0.2680587000, 0.2690368000, 0.2700137000, 
        0.2709895000, 0.2719642000, 0.2729378000, 0.2739103000, 0.2748817000, 
        0.2758521000, 0.2768215000, 0.2777898000, 0.2787572000, 0.2797235000, 
        0.2806890000, 0.2816534000, 0.2826170000, 0.2835796000, 0.2845414000, 
        0.2855022000, 0.2864623000, 0.2874214000, 0.2883798000, 0.2893373000, 
        0.2902940000, 0.2912500000, 0.2922052000, 0.2931597000, 0.2941134000, 
        0.2950664000, 0.2960187000, 0.2969704000, 0.2979213000, 0.2988716000, 
        0.2998213000, 0.3007704000, 0.3017189000, 0.3026667000, 0.3036140000, 
        0.3045608000, 0.3055070000, 0.3064527000, 0.3073978000, 0.3083425000, 
        0.3092867000, 0.3102304000, 0.3111737000, 0.3121166000, 0.3130590000, 
        0.3140010000, 0.3149426000, 0.3158839000, 0.3168248000, 0.3177653000, 
        0.3187056000, 0.3196455000, 0.3205851000, 0.3215244000, 0.3224634000, 
        0.3234022000, 0.3243408000, 0.3252791000, 0.3262172000, 0.3271551000, 
        0.3280928000, 0.3290303000, 0.3299677000, 0.3309050000, 0.3318421000, 
        0.3327791000, 0.3337160000, 0.3346528000, 0.3355896000, 0.3365263000, 
        0.3374629000, 0.3383995000, 0.3393362000, 0.3402728000, 0.3412094000, 
        0.3421460000, 0.3430827000, 0.3440194000, 0.3449562000, 0.3458931000, 
        0.3468301000, 0.3477672000, 0.3487044000, 0.3496418000, 0.3505793000, 
        0.3515170000, 0.3524548000, 0.3533929000, 0.3543311000, 0.3552696000, 
        0.3562083000, 0.3571472000, 0.3580864000, 0.3590259000, 0.3599657000, 
        0.3609057000, 0.3618461000, 0.3627868000, 0.3637279000, 0.3646693000, 
        0.3656111000, 0.3665532000, 0.3674958000, 0.3684387000, 0.3693821000, 
        0.3703259000, 0.3712701000, 0.3722149000, 0.3731600000, 0.3741057000, 
        0.3750519000, 0.3759986000, 0.3769458000, 0.3778936000, 0.3788419000, 
        0.3797908000, 0.3807402000, 0.3816903000, 0.3826409000, 0.3835922000, 
        0.3845441000, 0.3854966000, 0.3864498000, 0.3874037000, 0.3883583000, 
        0.3893135000, 0.3902695000, 0.3912262000, 0.3921836000, 0.3931418000, 
        0.3941007000, 0.3950604000, 0.3960209000, 0.3969822000, 0.3979443000, 
        0.3989072000, 0.3998710000, 0.4008356000, 0.4018011000, 0.4027675000, 
        0.4037347000, 0.4047029000, 0.4056720000, 0.4066420000, 0.4076129000, 
        0.4085848000, 0.4095576000, 0.4105315000, 0.4115063000, 0.4124821000, 
        0.4134590000, 0.4144368000, 0.4154157000, 0.4163957000, 0.4173767000, 
        0.4183588000, 0.4193420000, 0.4203264000, 0.4213118000, 0.4222983000, 
        0.4232860000, 0.4242749000, 0.4252649000, 0.4262561000, 0.4272484000, 
        0.4282420000, 0.4292368000, 0.4302328000, 0.4312301000, 0.4322286000, 
        0.4332284000, 0.4342295000, 0.4352318000, 0.4362354000, 0.4372404000, 
        0.4382467000, 0.4392543000, 0.4402633000, 0.4412736000, 0.4422853000, 
        0.4432984000, 0.4443129000, 0.4453288000, 0.4463461000, 0.4473648000, 
        0.4483850000, 0.4494067000, 0.4504298000, 0.4514544000, 0.4524804000, 
        0.4535080000, 0.4545371000, 0.4555678000, 0.4565999000, 0.4576337000, 
        0.4586689000, 0.4597058000, 0.4607442000, 0.4617843000, 0.4628259000, 
        0.4638692000, 0.4649141000, 0.4659606000, 0.4670088000, 0.4680586000, 
        0.4691101000, 0.4701633000, 0.4712183000, 0.4722749000, 0.4733332000, 
        0.4743933000, 0.4754551000, 0.4765186000, 0.4775840000, 0.4786511000, 
        0.4797199000, 0.4807906000, 0.4818631000, 0.4829374000, 0.4840135000, 
        0.4850915000, 0.4861713000, 0.4872530000, 0.4883366000, 0.4894220000, 
        0.4905093000, 0.4915985000, 0.4926897000, 0.4937827000, 0.4948777000, 
        0.4959747000, 0.4970736000, 0.4981744000, 0.4992772000, 0.5003820000, 
        0.5014888000, 0.5025976000, 0.5037085000, 0.5048213000, 0.5059362000, 
        0.5070531000, 0.5081720000, 0.5092931000, 0.5104162000, 0.5115413000, 
        0.5126686000, 0.5137980000, 0.5149294000, 0.5160630000, 0.5171987000, 
        0.5183366000, 0.5194766000, 0.5206187000, 0.5217630000, 0.5229095000, 
        0.5240581000, 0.5252090000, 0.5263620000, 0.5275173000, 0.5286748000, 
        0.5298344000, 0.5309964000, 0.5321605000, 0.5333269000, 0.5344956000, 
        0.5356665000, 0.5368397000, 0.5380152000, 0.5391930000, 0.5403731000, 
        0.5415554000, 0.5427401000, 0.5439272000, 0.5451165000, 0.5463082000, 
        0.5475022000, 0.5486986000, 0.5498973000, 0.5510984000, 0.5523019000, 
        0.5535077000, 0.5547159000, 0.5559266000, 0.5571396000, 0.5583551000, 
        0.5595729000, 0.5607932000, 0.5620159000, 0.5632411000, 0.5644686000, 
        0.5656987000, 0.5669312000, 0.5681661000, 0.5694036000, 0.5706435000, 
        0.5718859000, 0.5731307000, 0.5743781000, 0.5756280000, 0.5768804000, 
        0.5781352000, 0.5793926000, 0.5806526000, 0.5819150000, 0.5831800000, 
        0.5844476000, 0.5857176000, 0.5869903000, 0.5882655000, 0.5895432000, 
        0.5908235000, 0.5921064000, 0.5933919000, 0.5946799000, 0.5959706000, 
        0.5972638000, 0.5985596000, 0.5998581000, 0.6011591000, 0.6024627000, 
        0.6037690000, 0.6050779000, 0.6063894000, 0.6077036000, 0.6090203000, 
        0.6103397000, 0.6116618000, 0.6129865000, 0.6143139000, 0.6156439000, 
        0.6169765000, 0.6183119000, 0.6196499000, 0.6209905000, 0.6223339000, 
        0.6236799000, 0.6250286000, 0.6263800000, 0.6277340000, 0.6290908000, 
        0.6304503000, 0.6318124000, 0.6331773000, 0.6345448000, 0.6359151000, 
        0.6372881000, 0.6386638000, 0.6400422000, 0.6414233000, 0.6428071000, 
        0.6441937000, 0.6455830000, 0.6469750000, 0.6483698000, 0.6497673000, 
        0.6511675000, 0.6525705000, 0.6539762000, 0.6553846000, 0.6567958000, 
        0.6582098000, 0.6596264000, 0.6610459000, 0.6624681000, 0.6638930000, 
        0.6653208000, 0.6667512000, 0.6681845000, 0.6696205000, 0.6710592000, 
        0.6725007000, 0.6739450000, 0.6753921000, 0.6768419000, 0.6782945000, 
        0.6797499000, 0.6812080000, 0.6826690000, 0.6841327000, 0.6855992000, 
        0.6870684000, 0.6885404000, 0.6900153000, 0.6914929000, 0.6929732000, 
        0.6944564000, 0.6959423000, 0.6974311000, 0.6989226000, 0.7004169000, 
        0.7019140000, 0.7034139000, 0.7049165000, 0.7064220000, 0.7079302000, 
        0.7094412000, 0.7109550000, 0.7124716000, 0.7139910000, 0.7155132000, 
        0.7170382000, 0.7185659000, 0.7200965000, 0.7216298000, 0.7231660000, 
        0.7247049000, 0.7262466000, 0.7277911000, 0.7293383000, 0.7308884000, 
        0.7324413000, 0.7339969000, 0.7355553000, 0.7371166000, 0.7386806000, 
        0.7402474000, 0.7418169000, 0.7433893000, 0.7449644000, 0.7465424000, 
        0.7481231000, 0.7497066000, 0.7512929000, 0.7528819000, 0.7544737000, 
        0.7560684000, 0.7576657000, 0.7592659000, 0.7608689000, 0.7624746000, 
        0.7640831000, 0.7656943000, 0.7673084000, 0.7689252000, 0.7705448000, 
        0.7721671000, 0.7737922000, 0.7754201000, 0.7770507000, 0.7786841000, 
        0.7803203000, 0.7819592000, 0.7836009000, 0.7852454000, 0.7868926000, 
        0.7885425000, 0.7901952000, 0.7918507000, 0.7935089000, 0.7951698000, 
        0.7968335000, 0.7985000000, 0.8001692000, 0.8018411000, 0.8035158000, 
        0.8051932000, 0.8068733000, 0.8085562000, 0.8102418000, 0.8119301000, 
        0.8136212000, 0.8153150000, 0.8170115000, 0.8187108000, 0.8204127000, 
        0.8221174000, 0.8238248000, 0.8255349000, 0.8272477000, 0.8289633000, 
        0.8306815000, 0.8324025000, 0.8341261000, 0.8358525000, 0.8375815000, 
        0.8393133000, 0.8410477000, 0.8427849000, 0.8445247000, 0.8462672000, 
        0.8480124000, 0.8497603000, 0.8515108000, 0.8532641000, 0.8550200000, 
        0.8567786000, 0.8585399000, 0.8603038000, 0.8620704000, 0.8638396000, 
        0.8656116000, 0.8673861000, 0.8691634000, 0.8709433000, 0.8727258000, 
        0.8745110000, 0.8762988000, 0.8780893000, 0.8798824000, 0.8816782000, 
        0.8834766000, 0.8852776000, 0.8870812000, 0.8888875000, 0.8906964000, 
        0.8925079000, 0.8943221000, 0.8961388000, 0.8979582000, 0.8997801000, 
        0.9016047000, 0.9034319000, 0.9052617000, 0.9070941000, 0.9089291000, 
        0.9107666000, 0.9126068000, 0.9144495000, 0.9162949000, 0.9181428000, 
        0.9199933000, 0.9218463000, 0.9237020000, 0.9255602000, 0.9274210000, 
        0.9292843000, 0.9311502000, 0.9330186000, 0.9348896000, 0.9367632000, 
        0.9386393000, 0.9405179000, 0.9423991000, 0.9442828000, 0.9461691000, 
        0.9480579000, 0.9499492000, 0.9518431000, 0.9537394000, 0.9556383000, 
        0.9575397000, 0.9594436000, 0.9613501000, 0.9632590000, 0.9651704000, 
        0.9670844000, 0.9690008000, 0.9709197000, 0.9728411000, 0.9747650000, 
        0.9766914000, 0.9786203000, 0.9805516000, 0.9824854000, 0.9844217000, 
        0.9863605000, 0.9883017000, 0.9902454000, 0.9921915000, 0.9941401000, 
        0.9960911000, 0.9980446000, 1.0000010000, 1.0019590000, 1.0039200000, 
        1.0058830000, 1.0078490000, 1.0098170000, 1.0117870000, 1.0137600000, 
        1.0157360000, 1.0177130000, 1.0196940000, 1.0216760000, 1.0236610000, 
        1.0256490000, 1.0276380000, 1.0296310000, 1.0316250000, 1.0336220000, 
        1.0356210000, 1.0376230000, 1.0396270000, 1.0416340000, 1.0436430000, 
        1.0456540000, 1.0476670000, 1.0496830000, 1.0517020000, 1.0537220000, 
        1.0557450000, 1.0577710000, 1.0597980000, 1.0618280000, 1.0638610000, 
        1.0658950000, 1.0679320000, 1.0699720000, 1.0720130000, 1.0740570000, 
        1.0761040000, 1.0781520000, 1.0802030000, 1.0822560000, 1.0843120000, 
        1.0863690000, 1.0884300000, 1.0904920000, 1.0925570000, 1.0946240000, 
        1.0966930000, 1.0987640000, 1.1008380000, 1.1029140000, 1.1049920000, 
        1.1070730000, 1.1091560000, 1.1112410000, 1.1133280000, 1.1154180000, 
        1.1175090000, 1.1196030000, 1.1217000000, 1.1237980000, 1.1258990000, 
        1.1280020000, 1.1301070000, 1.1322140000, 1.1343240000, 1.1364360000, 
        1.1385500000, 1.1406660000, 1.1427850000, 1.1449050000, 1.1470280000, 
        1.1491530000, 1.1512800000, 1.1534100000, 1.1555410000, 1.1576750000, 
        1.1598110000, 1.1619490000, 1.1640890000, 1.1662310000, 1.1683760000, 
        1.1705230000, 1.1726710000, 1.1748220000, 1.1769750000, 1.1791310000, 
        1.1812880000, 1.1834480000, 1.1856090000, 1.1877730000, 1.1899390000, 
        1.1921070000, 1.1942770000, 1.1964490000, 1.1986240000, 1.2008000000, 
        1.2029780000, 1.2051590000, 1.2073420000, 1.2095270000, 1.2117130000, 
        1.2139020000, 1.2160930000, 1.2182860000, 1.2204820000, 1.2226790000, 
        1.2248780000, 1.2270790000, 1.2292830000, 1.2314880000, 1.2336960000, 
        1.2359050000, 1.2381170000, 1.2403300000, 1.2425460000, 1.2447640000, 
        1.2469830000, 1.2492050000, 1.2514290000, 1.2536540000, 1.2558820000, 
        1.2581120000, 1.2603430000, 1.2625770000, 1.2648130000, 1.2670510000, 
        1.2692900000, 1.2715320000, 1.2737760000, 1.2760210000, 1.2782690000, 
        1.2805180000, 1.2827700000, 1.2850230000, 1.2872790000, 1.2895360000, 
        1.2917950000, 1.2940560000, 1.2963200000, 1.2985850000, 1.3008520000, 
        1.3031210000, 1.3053920000, 1.3076650000, 1.3099390000, 1.3122160000, 
        1.3144940000, 1.3167750000, 1.3190570000, 1.3213420000, 1.3236280000, 
        1.3259160000, 1.3282060000, 1.3304970000, 1.3327910000, 1.3350870000, 
        1.3373840000, 1.3396830000, 1.3419850000, 1.3442880000, 1.3465920000, 
        1.3488990000, 1.3512080000, 1.3535180000, 1.3558300000, 1.3581440000, 
        1.3604600000, 1.3627780000, 1.3650970000, 1.3674190000, 1.3697420000, 
        1.3720670000, 1.3743940000, 1.3767220000, 1.3790530000, 1.3813850000, 
        1.3837190000, 1.3860540000, 1.3883920000, 1.3907310000, 1.3930720000, 
        1.3954150000, 1.3977600000, 1.4001060000, 1.4024540000, 1.4048040000, 
        1.4071560000, 1.4095090000, 1.4118640000, 1.4142210000, 1.4165800000, 
        1.4189400000, 1.4213020000, 1.4236660000, 1.4260310000, 1.4283990000, 
        1.4307670000, 1.4331380000, 1.4355100000, 1.4378840000, 1.4402600000, 
        1.4426370000, 1.4450160000, 1.4473970000, 1.4497800000, 1.4521640000, 
        1.4545490000, 1.4569370000, 1.4593260000, 1.4617170000, 1.4641090000, 
        1.4665030000, 1.4688990000, 1.4712960000, 1.4736950000, 1.4760960000, 
        1.4784980000, 1.4809020000, 1.4833070000, 1.4857140000, 1.4881230000, 
        1.4905330000, 1.4929450000, 1.4953580000, 1.4977730000, 1.5001900000, 
        1.5026080000, 1.5050280000, 1.5074490000, 1.5098720000, 1.5122970000, 
        1.5147230000, 1.5171500000, 1.5195790000, 1.5220100000, 1.5244420000, 
        1.5268760000, 1.5293110000, 1.5317480000, 1.5341870000, 1.5366270000, 
        1.5390680000, 1.5415110000, 1.5439550000, 1.5464010000, 1.5488490000, 
        1.5512980000, 1.5537480000, 1.5562000000, 1.5586530000, 1.5611080000, 
        1.5635650000, 1.5660230000, 1.5684820000, 1.5709430000, 1.5734050000, 
        1.5758680000, 1.5783340000, 1.5808000000, 1.5832680000, 1.5857380000, 
        1.5882080000, 1.5906810000, 1.5931540000, 1.5956290000, 1.5981060000, 
        1.6005840000, 1.6030630000, 1.6055440000, 1.6080260000, 1.6105100000, 
        1.6129950000, 1.6154810000, 1.6179680000, 1.6204570000, 1.6229480000, 
        1.6254400000, 1.6279330000, 1.6304270000, 1.6329230000, 1.6354200000, 
        1.6379190000, 1.6404190000, 1.6429200000, 1.6454220000, 1.6479260000, 
        1.6504310000, 1.6529380000, 1.6554450000, 1.6579540000, 1.6604650000, 
        1.6629760000, 1.6654890000, 1.6680040000, 1.6705190000, 1.6730360000, 
        1.6755540000, 1.6780730000, 1.6805940000, 1.6831160000, 1.6856390000, 
        1.6881630000, 1.6906890000, 1.6932160000, 1.6957440000, 1.6982730000, 
        1.7008030000, 1.7033350000, 1.7058680000, 1.7084020000, 1.7109380000, 
        1.7134740000, 1.7160120000, 1.7185510000, 1.7210910000, 1.7236330000, 
        1.7261750000, 1.7287190000, 1.7312640000, 1.7338100000, 1.7363570000, 
        1.7389050000, 1.7414550000, 1.7440050000, 1.7465570000, 1.7491100000, 
        1.7516640000, 1.7542190000, 1.7567760000, 1.7593330000, 1.7618920000, 
        1.7644520000, 1.7670120000, 1.7695740000, 1.7721370000, 1.7747010000, 
        1.7772660000, 1.7798330000, 1.7824000000, 1.7849680000, 1.7875380000, 
        1.7901080000, 1.7926800000, 1.7952530000, 1.7978260000, 1.8004010000, 
        1.8029770000, 1.8055540000, 1.8081320000, 1.8107100000, 1.8132900000, 
        1.8158710000, 1.8184530000, 1.8210360000, 1.8236200000, 1.8262050000, 
        1.8287910000, 1.8313780000, 1.8339660000, 1.8365550000, 1.8391450000, 
        1.8417360000, 1.8443270000, 1.8469200000, 1.8495140000, 1.8521090000, 
        1.8547040000, 1.8573010000, 1.8598980000, 1.8624970000, 1.8650960000, 
        1.8676970000, 1.8702980000, 1.8729000000, 1.8755030000, 1.8781070000, 
        1.8807120000, 1.8833180000, 1.8859250000, 1.8885320000, 1.8911410000, 
        1.8937500000, 1.8963600000, 1.8989710000, 1.9015830000, 1.9041960000, 
        1.9068100000, 1.9094240000, 1.9120400000, 1.9146560000, 1.9172730000, 
        1.9198910000, 1.9225100000, 1.9251290000, 1.9277500000, 1.9303710000, 
        1.9329930000, 1.9356160000, 1.9382400000, 1.9408640000, 1.9434890000, 
        1.9461150000, 1.9487420000, 1.9513700000, 1.9539980000, 1.9566270000, 
        1.9592570000, 1.9618880000, 1.9645190000, 1.9671520000, 1.9697850000, 
        1.9724180000, 1.9750530000, 1.9776880000, 1.9803240000, 1.9829610000, 
        1.9855980000, 1.9882360000, 1.9908750000, 1.9935140000, 1.9961550000, 
        1.9987960000, 2.0014370000, 2.0040800000, 2.0067230000, 2.0093660000, 
        2.0120110000, 2.0146560000, 2.0173010000, 2.0199480000, 2.0225950000, 
        2.0252430000, 2.0278910000, 2.0305400000, 2.0331900000, 2.0358400000, 
        2.0384910000, 2.0411430000, 2.0437950000, 2.0464480000, 2.0491010000, 
        2.0517550000, 2.0544100000, 2.0570650000, 2.0597210000, 2.0623780000, 
        2.0650350000, 2.0676930000, 2.0703510000, 2.0730100000, 2.0756690000, 
        2.0783290000, 2.0809900000, 2.0836510000, 2.0863130000, 2.0889750000, 
        2.0916380000, 2.0943010000, 2.0969650000, 2.0996290000, 2.1022940000, 
        2.1049600000, 2.1076260000, 2.1102930000, 2.1129600000, 2.1156270000, 
        2.1182950000, 2.1209640000, 2.1236330000, 2.1263030000, 2.1289730000, 
        2.1316440000, 2.1343150000, 2.1369860000, 2.1396580000, 2.1423310000, 
        2.1450040000, 2.1476780000, 2.1503510000, 2.1530260000, 2.1557010000
    };
    double K1400_width[2000]={//1700 bin, s = (0.6, 4)
        0.0000000007, 0.0000000025, 0.0000000061, 0.0000000121, 0.0000000213, 
        0.0000000344, 0.0000000519, 0.0000000747, 0.0000001034, 0.0000001389, 
        0.0000001818, 0.0000002329, 0.0000002930, 0.0000003630, 0.0000004435, 
        0.0000005356, 0.0000006400, 0.0000007576, 0.0000008892, 0.0000010359, 
        0.0000011985, 0.0000013780, 0.0000015753, 0.0000017914, 0.0000020274, 
        0.0000022842, 0.0000025629, 0.0000028646, 0.0000031903, 0.0000035413, 
        0.0000039185, 0.0000043232, 0.0000047566, 0.0000052199, 0.0000057143, 
        0.0000062411, 0.0000068015, 0.0000073970, 0.0000080288, 0.0000086983, 
        0.0000094069, 0.0000101561, 0.0000109474, 0.0000117821, 0.0000126619, 
        0.0000135884, 0.0000145631, 0.0000155876, 0.0000166637, 0.0000177930, 
        0.0000189773, 0.0000202184, 0.0000215182, 0.0000228784, 0.0000243010, 
        0.0000257879, 0.0000273412, 0.0000289629, 0.0000306551, 0.0000324199, 
        0.0000342596, 0.0000361763, 0.0000381724, 0.0000402502, 0.0000424121, 
        0.0000446606, 0.0000469982, 0.0000494275, 0.0000519511, 0.0000545718, 
        0.0000572923, 0.0000601155, 0.0000630443, 0.0000660816, 0.0000692306, 
        0.0000724943, 0.0000758761, 0.0000793791, 0.0000830068, 0.0000867626, 
        0.0000906501, 0.0000946730, 0.0000988349, 0.0001031398, 0.0001075916, 
        0.0001121943, 0.0001169520, 0.0001218691, 0.0001269499, 0.0001321989, 
        0.0001376208, 0.0001432202, 0.0001490021, 0.0001549715, 0.0001611335, 
        0.0001674933, 0.0001740566, 0.0001808287, 0.0001878156, 0.0001950231, 
        0.0002024573, 0.0002101244, 0.0002180310, 0.0002261837, 0.0002345892, 
        0.0002432546, 0.0002521872, 0.0002613945, 0.0002708840, 0.0002806638, 
        0.0002907420, 0.0003011271, 0.0003118276, 0.0003228526, 0.0003342113, 
        0.0003459132, 0.0003579682, 0.0003703863, 0.0003831781, 0.0003963544, 
        0.0004099263, 0.0004239055, 0.0004383037, 0.0004531333, 0.0004684072, 
        0.0004841383, 0.0005003405, 0.0005170277, 0.0005342145, 0.0005519161, 
        0.0005701480, 0.0005889265, 0.0006082683, 0.0006281907, 0.0006487119, 
        0.0006698504, 0.0006916256, 0.0007140575, 0.0007371670, 0.0007609758, 
        0.0007855062, 0.0008107815, 0.0008368260, 0.0008636648, 0.0008913240, 
        0.0009198309, 0.0009492137, 0.0009795018, 0.0010107260, 0.0010429180, 
        0.0010761110, 0.0011103390, 0.0011456390, 0.0011820480, 0.0012196060, 
        0.0012583530, 0.0012983320, 0.0013395870, 0.0013821650, 0.0014261150, 
        0.0014714870, 0.0015183340, 0.0015667120, 0.0016166800, 0.0016682970, 
        0.0017216280, 0.0017767390, 0.0018337000, 0.0018925850, 0.0019534700, 
        0.0020164360, 0.0020815670, 0.0021489530, 0.0022186860, 0.0022908640, 
        0.0023655910, 0.0024429750, 0.0025231290, 0.0026061740, 0.0026922350, 
        0.0027814460, 0.0028739440, 0.0029698790, 0.0030694050, 0.0031726840, 
        0.0032798900, 0.0033912030, 0.0035068140, 0.0036269240, 0.0037517470, 
        0.0038815040, 0.0040164320, 0.0041567790, 0.0043028050, 0.0044547850, 
        0.0046130090, 0.0047777790, 0.0049494150, 0.0051282500, 0.0053146360, 
        0.0055089370, 0.0057115350, 0.0059228300, 0.0061432340, 0.0063731760, 
        0.0066130990, 0.0068634590, 0.0071247200, 0.0073973600, 0.0076818580, 
        0.0079786990, 0.0082883670, 0.0086113370, 0.0089480750, 0.0092990290, 
        0.0096646210, 0.0100452400, 0.0104412500, 0.0108529300, 0.0112805500, 
        0.0117242700, 0.0121842100, 0.0126603900, 0.0131527500, 0.0136611300, 
        0.0141852900, 0.0147248800, 0.0152794700, 0.0158485100, 0.0164314100, 
        0.0170274700, 0.0176359300, 0.0182559700, 0.0188867400, 0.0195273500, 
        0.0201768900, 0.0208344400, 0.0214991000, 0.0221699600, 0.0228461600, 
        0.0235268400, 0.0242112100, 0.0248985100, 0.0255880100, 0.0262790500, 
        0.0269710000, 0.0276632900, 0.0283554000, 0.0290468400, 0.0297371800, 
        0.0304260300, 0.0311130400, 0.0317978800, 0.0324802900, 0.0331600000, 
        0.0338368100, 0.0345105300, 0.0351809800, 0.0358480300, 0.0365115500, 
        0.0371714500, 0.0378276300, 0.0384800200, 0.0391285600, 0.0397732200, 
        0.0404139500, 0.0410507400, 0.0416835500, 0.0423124000, 0.0429372700, 
        0.0435581700, 0.0441751100, 0.0447881100, 0.0453971800, 0.0460023400, 
        0.0466036300, 0.0472010600, 0.0477946700, 0.0483845000, 0.0489705600, 
        0.0495529100, 0.0501315800, 0.0507066000, 0.0512780200, 0.0518458700, 
        0.0524101800, 0.0529710200, 0.0535284000, 0.0540823700, 0.0546329800, 
        0.0551802500, 0.0557242400, 0.0562649800, 0.0568025200, 0.0573368800, 
        0.0578681100, 0.0583962500, 0.0589213300, 0.0594434000, 0.0599624900, 
        0.0604786400, 0.0609918800, 0.0615022500, 0.0620097800, 0.0625145200, 
        0.0630164900, 0.0635157300, 0.0640122700, 0.0645061500, 0.0649973900, 
        0.0654860300, 0.0659721100, 0.0664556400, 0.0669366600, 0.0674152100, 
        0.0678913000, 0.0683649800, 0.0688362600, 0.0693051800, 0.0697717600, 
        0.0702360300, 0.0706980200, 0.0711577500, 0.0716152500, 0.0720705400, 
        0.0725236500, 0.0729746100, 0.0734234300, 0.0738701400, 0.0743147600, 
        0.0747573300, 0.0751978500, 0.0756363500, 0.0760728500, 0.0765073800, 
        0.0769399600, 0.0773706000, 0.0777993300, 0.0782261600, 0.0786511200, 
        0.0790742300, 0.0794955000, 0.0799149600, 0.0803326100, 0.0807484900, 
        0.0811626100, 0.0815749800, 0.0819856300, 0.0823945700, 0.0828018200, 
        0.0832073900, 0.0836113000, 0.0840135700, 0.0844142200, 0.0848132500, 
        0.0852106900, 0.0856065500, 0.0860008400, 0.0863935900, 0.0867848000, 
        0.0871744900, 0.0875626700, 0.0879493600, 0.0883345700, 0.0887183200, 
        0.0891006200, 0.0894814800, 0.0898609100, 0.0902389300, 0.0906155600, 
        0.0909908000, 0.0913646700, 0.0917371700, 0.0921083300, 0.0924781500, 
        0.0928466500, 0.0932138400, 0.0935797200, 0.0939443200, 0.0943076400, 
        0.0946696900, 0.0950304800, 0.0953900400, 0.0957483600, 0.0961054500, 
        0.0964613400, 0.0968160300, 0.0971695200, 0.0975218300, 0.0978729800, 
        0.0982229600, 0.0985718000, 0.0989194900, 0.0992660600, 0.0996115100, 
        0.0999558400, 0.1002991000, 0.1006412000, 0.1009823000, 0.1013223000, 
        0.1016612000, 0.1019991000, 0.1023359000, 0.1026717000, 0.1030064000, 
        0.1033402000, 0.1036729000, 0.1040047000, 0.1043354000, 0.1046652000, 
        0.1049939000, 0.1053218000, 0.1056486000, 0.1059745000, 0.1062995000, 
        0.1066235000, 0.1069466000, 0.1072688000, 0.1075901000, 0.1079104000, 
        0.1082299000, 0.1085485000, 0.1088662000, 0.1091831000, 0.1094990000, 
        0.1098142000, 0.1101284000, 0.1104418000, 0.1107544000, 0.1110662000, 
        0.1113771000, 0.1116873000, 0.1119966000, 0.1123051000, 0.1126128000, 
        0.1129198000, 0.1132259000, 0.1135313000, 0.1138360000, 0.1141398000, 
        0.1144429000, 0.1147453000, 0.1150469000, 0.1153478000, 0.1156480000, 
        0.1159474000, 0.1162461000, 0.1165442000, 0.1168415000, 0.1171381000, 
        0.1174340000, 0.1177292000, 0.1180238000, 0.1183177000, 0.1186109000, 
        0.1189034000, 0.1191953000, 0.1194865000, 0.1197771000, 0.1200670000, 
        0.1203563000, 0.1206449000, 0.1209329000, 0.1212203000, 0.1215071000, 
        0.1217933000, 0.1220788000, 0.1223637000, 0.1226480000, 0.1229318000, 
        0.1232149000, 0.1234974000, 0.1237794000, 0.1240607000, 0.1243415000, 
        0.1246217000, 0.1249013000, 0.1251803000, 0.1254588000, 0.1257367000, 
        0.1260140000, 0.1262908000, 0.1265670000, 0.1268427000, 0.1271178000, 
        0.1273923000, 0.1276663000, 0.1279398000, 0.1282127000, 0.1284850000, 
        0.1287568000, 0.1290281000, 0.1292988000, 0.1295690000, 0.1298387000, 
        0.1301078000, 0.1303764000, 0.1306445000, 0.1309120000, 0.1311790000, 
        0.1314455000, 0.1317114000, 0.1319768000, 0.1322417000, 0.1325060000, 
        0.1327699000, 0.1330332000, 0.1332960000, 0.1335582000, 0.1338199000, 
        0.1340812000, 0.1343418000, 0.1346020000, 0.1348617000, 0.1351208000, 
        0.1355387000, 0.1358629000, 0.1361711000, 0.1364708000, 0.1367648000, 
        0.1370546000, 0.1373411000, 0.1376249000, 0.1379063000, 0.1381858000, 
        0.1384634000, 0.1387394000, 0.1390139000, 0.1392870000, 0.1395588000, 
        0.1398294000, 0.1400988000, 0.1403671000, 0.1406344000, 0.1409007000, 
        0.1411659000, 0.1414302000, 0.1416936000, 0.1419562000, 0.1422178000, 
        0.1424786000, 0.1427386000, 0.1429977000, 0.1432561000, 0.1435136000, 
        0.1437704000, 0.1440265000, 0.1442818000, 0.1445364000, 0.1447902000, 
        0.1450434000, 0.1452958000, 0.1455475000, 0.1457986000, 0.1460490000, 
        0.1462986000, 0.1465477000, 0.1467960000, 0.1470437000, 0.1472908000, 
        0.1475372000, 0.1477830000, 0.1480281000, 0.1482726000, 0.1485165000, 
        0.1487597000, 0.1490024000, 0.1492444000, 0.1494858000, 0.1497266000, 
        0.1499668000, 0.1502065000, 0.1504455000, 0.1506839000, 0.1509218000, 
        0.1511590000, 0.1513957000, 0.1516318000, 0.1518674000, 0.1521023000, 
        0.1523367000, 0.1525706000, 0.1528038000, 0.1530366000, 0.1532687000, 
        0.1535003000, 0.1537314000, 0.1539619000, 0.1541918000, 0.1544213000, 
        0.1546501000, 0.1548785000, 0.1551063000, 0.1553336000, 0.1555603000, 
        0.1557866000, 0.1560122000, 0.1562374000, 0.1564621000, 0.1566862000, 
        0.1569098000, 0.1571329000, 0.1573555000, 0.1575776000, 0.1577992000, 
        0.1580203000, 0.1582409000, 0.1584609000, 0.1586805000, 0.1588996000, 
        0.1591182000, 0.1593363000, 0.1595539000, 0.1597710000, 0.1599876000, 
        0.1602038000, 0.1604194000, 0.1606346000, 0.1608493000, 0.1610636000, 
        0.1612773000, 0.1614906000, 0.1617035000, 0.1619158000, 0.1621277000, 
        0.1623391000, 0.1625501000, 0.1627606000, 0.1629706000, 0.1631802000, 
        0.1633894000, 0.1635981000, 0.1638063000, 0.1640141000, 0.1642214000, 
        0.1644283000, 0.1646347000, 0.1648407000, 0.1650463000, 0.1652514000, 
        0.1654561000, 0.1656604000, 0.1658642000, 0.1660676000, 0.1662706000, 
        0.1664731000, 0.1666752000, 0.1668769000, 0.1670782000, 0.1672790000, 
        0.1674794000, 0.1676794000, 0.1678790000, 0.1680782000, 0.1682769000, 
        0.1684753000, 0.1686732000, 0.1688707000, 0.1690679000, 0.1692646000, 
        0.1694609000, 0.1696568000, 0.1698523000, 0.1700474000, 0.1702421000, 
        0.1704364000, 0.1706304000, 0.1708239000, 0.1710170000, 0.1712098000, 
        0.1714022000, 0.1715941000, 0.1717857000, 0.1719769000, 0.1721677000, 
        0.1723582000, 0.1725482000, 0.1727379000, 0.1729272000, 0.1731162000, 
        0.1733047000, 0.1734929000, 0.1736807000, 0.1738682000, 0.1740552000, 
        0.1742420000, 0.1744283000, 0.1746143000, 0.1747999000, 0.1749851000, 
        0.1751700000, 0.1753546000, 0.1755388000, 0.1757226000, 0.1759060000, 
        0.1760892000, 0.1762719000, 0.1764543000, 0.1766364000, 0.1768181000, 
        0.1769995000, 0.1771805000, 0.1773612000, 0.1775415000, 0.1777215000, 
        0.1779011000, 0.1780804000, 0.1782594000, 0.1784380000, 0.1786163000, 
        0.1787943000, 0.1789719000, 0.1791492000, 0.1793262000, 0.1795028000, 
        0.1796791000, 0.1798551000, 0.1800308000, 0.1802061000, 0.1803811000, 
        0.1805558000, 0.1807301000, 0.1809042000, 0.1810779000, 0.1812513000, 
        0.1814244000, 0.1815971000, 0.1817696000, 0.1819417000, 0.1821136000, 
        0.1822851000, 0.1824563000, 0.1826272000, 0.1827978000, 0.1829681000, 
        0.1831380000, 0.1833077000, 0.1834771000, 0.1836461000, 0.1838149000, 
        0.1839834000, 0.1841515000, 0.1843194000, 0.1844869000, 0.1846542000, 
        0.1848212000, 0.1849879000, 0.1851543000, 0.1853203000, 0.1854862000, 
        0.1856517000, 0.1858169000, 0.1859818000, 0.1861465000, 0.1863108000, 
        0.1864749000, 0.1866387000, 0.1868022000, 0.1869654000, 0.1871284000, 
        0.1872910000, 0.1874534000, 0.1876155000, 0.1877774000, 0.1879389000, 
        0.1881002000, 0.1882612000, 0.1884219000, 0.1885824000, 0.1887426000, 
        0.1889025000, 0.1890621000, 0.1892215000, 0.1893806000, 0.1895395000, 
        0.1896980000, 0.1898563000, 0.1900144000, 0.1901722000, 0.1903297000, 
        0.1904869000, 0.1906439000, 0.1908007000, 0.1909571000, 0.1911133000, 
        0.1912693000, 0.1914250000, 0.1915805000, 0.1917357000, 0.1918906000, 
        0.1920453000, 0.1921997000, 0.1923539000, 0.1925078000, 0.1926615000, 
        0.1928149000, 0.1929681000, 0.1931211000, 0.1932738000, 0.1934262000, 
        0.1935784000, 0.1937304000, 0.1938821000, 0.1940336000, 0.1941848000, 
        0.1943358000, 0.1944865000, 0.1946370000, 0.1947873000, 0.1949374000, 
        0.1950872000, 0.1952367000, 0.1953860000, 0.1955351000, 0.1956840000, 
        0.1958326000, 0.1959810000, 0.1961292000, 0.1962771000, 0.1964248000, 
        0.1965723000, 0.1967196000, 0.1968666000, 0.1970134000, 0.1971600000, 
        0.1973063000, 0.1974524000, 0.1975983000, 0.1977440000, 0.1978894000, 
        0.1980347000, 0.1981797000, 0.1983245000, 0.1984690000, 0.1986134000, 
        0.1987575000, 0.1989015000, 0.1990452000, 0.1991886000, 0.1993319000, 
        0.1994750000, 0.1996178000, 0.1997604000, 0.1999029000, 0.2000451000, 
        0.2001871000, 0.2003289000, 0.2004704000, 0.2006118000, 0.2007530000, 
        0.2008939000, 0.2010347000, 0.2011752000, 0.2013156000, 0.2014557000, 
        0.2015956000, 0.2017353000, 0.2018749000, 0.2020142000, 0.2021533000, 
        0.2022922000, 0.2024310000, 0.2025695000, 0.2027078000, 0.2028459000, 
        0.2029839000, 0.2031216000, 0.2032591000, 0.2033965000, 0.2035336000, 
        0.2036706000, 0.2038073000, 0.2039439000, 0.2040803000, 0.2042165000, 
        0.2043524000, 0.2044883000, 0.2046239000, 0.2047593000, 0.2048945000, 
        0.2050296000, 0.2051644000, 0.2052991000, 0.2054336000, 0.2055679000, 
        0.2057020000, 0.2058360000, 0.2059697000, 0.2061033000, 0.2062367000, 
        0.2063699000, 0.2065029000, 0.2066358000, 0.2067684000, 0.2069009000, 
        0.2070332000, 0.2071653000, 0.2072973000, 0.2074291000, 0.2075607000, 
        0.2076921000, 0.2078233000, 0.2079544000, 0.2080853000, 0.2082161000, 
        0.2083466000, 0.2084770000, 0.2086072000, 0.2087372000, 0.2088671000, 
        0.2089968000, 0.2091263000, 0.2092557000, 0.2093849000, 0.2095139000, 
        0.2096428000, 0.2097715000, 0.2099000000, 0.2100284000, 0.2101566000, 
        0.2102846000, 0.2104125000, 0.2105402000, 0.2106677000, 0.2107951000, 
        0.2109223000, 0.2110494000, 0.2111763000, 0.2113030000, 0.2114296000, 
        0.2115560000, 0.2116823000, 0.2118084000, 0.2119343000, 0.2120601000, 
        0.2121858000, 0.2123112000, 0.2124366000, 0.2125617000, 0.2126868000, 
        0.2128116000, 0.2129363000, 0.2130609000, 0.2131853000, 0.2133095000, 
        0.2134337000, 0.2135576000, 0.2136814000, 0.2138051000, 0.2139286000, 
        0.2140519000, 0.2141751000, 0.2142982000, 0.2144211000, 0.2145439000, 
        0.2146665000, 0.2147890000, 0.2149113000, 0.2150335000, 0.2151555000, 
        0.2152774000, 0.2153992000, 0.2155208000, 0.2156423000, 0.2157636000, 
        0.2158848000, 0.2160059000, 0.2161268000, 0.2162475000, 0.2163682000, 
        0.2164887000, 0.2166090000, 0.2167292000, 0.2168493000, 0.2169693000, 
        0.2170891000, 0.2172088000, 0.2173283000, 0.2174477000, 0.2175670000, 
        0.2176861000, 0.2178051000, 0.2179240000, 0.2180427000, 0.2181613000, 
        0.2182798000, 0.2183982000, 0.2185164000, 0.2186345000, 0.2187524000, 
        0.2188703000, 0.2189880000, 0.2191055000, 0.2192230000, 0.2193403000, 
        0.2194575000, 0.2195745000, 0.2196915000, 0.2198083000, 0.2199250000, 
        0.2200416000, 0.2201580000, 0.2202743000, 0.2203905000, 0.2205066000, 
        0.2206225000, 0.2207384000, 0.2208541000, 0.2209697000, 0.2210851000, 
        0.2212005000, 0.2213157000, 0.2214308000, 0.2215458000, 0.2216607000, 
        0.2217754000, 0.2218901000, 0.2220046000, 0.2221190000, 0.2222333000, 
        0.2223475000, 0.2224615000, 0.2225755000, 0.2226893000, 0.2228030000, 
        0.2229166000, 0.2230301000, 0.2231435000, 0.2232568000, 0.2233699000, 
        0.2234830000, 0.2235959000, 0.2237088000, 0.2238215000, 0.2239341000, 
        0.2240466000, 0.2241590000, 0.2242713000, 0.2243834000, 0.2244955000, 
        0.2246075000, 0.2247193000, 0.2248311000, 0.2249427000, 0.2250543000, 
        0.2251657000, 0.2252770000, 0.2253883000, 0.2254994000, 0.2256104000, 
        0.2257213000, 0.2258321000, 0.2259429000, 0.2260535000, 0.2261640000, 
        0.2262744000, 0.2263847000, 0.2264949000, 0.2266050000, 0.2267151000, 
        0.2268250000, 0.2269348000, 0.2270445000, 0.2271542000, 0.2272637000, 
        0.2273731000, 0.2274825000, 0.2275917000, 0.2277009000, 0.2278099000, 
        0.2279189000, 0.2280277000, 0.2281365000, 0.2282452000, 0.2283538000, 
        0.2284623000, 0.2285707000, 0.2286790000, 0.2287872000, 0.2288954000, 
        0.2290034000, 0.2291114000, 0.2292192000, 0.2293270000, 0.2294347000, 
        0.2295423000, 0.2296498000, 0.2297572000, 0.2298646000, 0.2299718000, 
        0.2300790000, 0.2301861000, 0.2302931000, 0.2304000000, 0.2305068000, 
        0.2306136000, 0.2307202000, 0.2308268000, 0.2309333000, 0.2310397000, 
        0.2311460000, 0.2312523000, 0.2313584000, 0.2314645000, 0.2315705000, 
        0.2316765000, 0.2317823000, 0.2318881000, 0.2319938000, 0.2320994000, 
        0.2322049000, 0.2323104000, 0.2324157000, 0.2325210000, 0.2326263000, 
        0.2327314000, 0.2328365000, 0.2329415000, 0.2330464000, 0.2331513000, 
        0.2332560000, 0.2333607000, 0.2334654000, 0.2335699000, 0.2336744000, 
        0.2337788000, 0.2338831000, 0.2339874000, 0.2340916000, 0.2341957000, 
        0.2342998000, 0.2344038000, 0.2345077000, 0.2346115000, 0.2347153000, 
        0.2348190000, 0.2349227000, 0.2350262000, 0.2351298000, 0.2352332000, 
        0.2353366000, 0.2354399000, 0.2355431000, 0.2356463000, 0.2357494000, 
        0.2358525000, 0.2359555000, 0.2360584000, 0.2361612000, 0.2362640000, 
        0.2363668000, 0.2364694000, 0.2365721000, 0.2366746000, 0.2367771000, 
        0.2368795000, 0.2369819000, 0.2370842000, 0.2371864000, 0.2372886000, 
        0.2373908000, 0.2374928000, 0.2375948000, 0.2376968000, 0.2377987000, 
        0.2379005000, 0.2380023000, 0.2381041000, 0.2382057000, 0.2383074000, 
        0.2384089000, 0.2385104000, 0.2386119000, 0.2387133000, 0.2388146000, 
        0.2389159000, 0.2390172000, 0.2391184000, 0.2392195000, 0.2393206000, 
        0.2394216000, 0.2395226000, 0.2396236000, 0.2397245000, 0.2398253000, 
        0.2399261000, 0.2400268000, 0.2401275000, 0.2402281000, 0.2403287000, 
        0.2404293000, 0.2405298000, 0.2406302000, 0.2407306000, 0.2408310000, 
        0.2409313000, 0.2410316000, 0.2411318000, 0.2412320000, 0.2413321000, 
        0.2414322000, 0.2415323000, 0.2416323000, 0.2417323000, 0.2418322000, 
        0.2419321000, 0.2420319000, 0.2421317000, 0.2422315000, 0.2423312000, 
        0.2424309000, 0.2425305000, 0.2426301000, 0.2427297000, 0.2428292000, 
        0.2429287000, 0.2430281000, 0.2431275000, 0.2432269000, 0.2433262000, 
        0.2434255000, 0.2435248000, 0.2436240000, 0.2437232000, 0.2438224000, 
        0.2439215000, 0.2440206000, 0.2441196000, 0.2442187000, 0.2443176000, 
        0.2444166000, 0.2445155000, 0.2446144000, 0.2447133000, 0.2448121000, 
        0.2449109000, 0.2450097000, 0.2451084000, 0.2452071000, 0.2453058000, 
        0.2454044000, 0.2455031000, 0.2456017000, 0.2457002000, 0.2457988000, 
        0.2458973000, 0.2459957000, 0.2460942000, 0.2461926000, 0.2462910000, 
        0.2463894000, 0.2464878000, 0.2465861000, 0.2466844000, 0.2467827000, 
        0.2468809000, 0.2469792000, 0.2470774000, 0.2471756000, 0.2472737000, 
        0.2473719000, 0.2474700000, 0.2475681000, 0.2476662000, 0.2477642000, 
        0.2478623000, 0.2479603000, 0.2480583000, 0.2481563000, 0.2482542000, 
        0.2483522000, 0.2484501000, 0.2485480000, 0.2486459000, 0.2487437000, 
        0.2488416000, 0.2489394000, 0.2490372000, 0.2491351000, 0.2492328000, 
        0.2493306000, 0.2494284000, 0.2495261000, 0.2496239000, 0.2497216000, 
        0.2498193000, 0.2499170000, 0.2500146000, 0.2501123000, 0.2502100000, 
        0.2503076000, 0.2504052000, 0.2505029000, 0.2506005000, 0.2506981000, 
        0.2507956000, 0.2508932000, 0.2509908000, 0.2510883000, 0.2511859000, 
        0.2512834000, 0.2513810000, 0.2514785000, 0.2515760000, 0.2516735000, 
        0.2517710000, 0.2518685000, 0.2519660000, 0.2520635000, 0.2521610000, 
        0.2522584000, 0.2523559000, 0.2524534000, 0.2525508000, 0.2526483000, 
        0.2527457000, 0.2528432000, 0.2529406000, 0.2530381000, 0.2531355000, 
        0.2532330000, 0.2533304000, 0.2534278000, 0.2535252000, 0.2536227000, 
        0.2537201000, 0.2538175000, 0.2539150000, 0.2540124000, 0.2541098000, 
        0.2542073000, 0.2543047000, 0.2544021000, 0.2544996000, 0.2545970000, 
        0.2546944000, 0.2547919000, 0.2548893000, 0.2549868000, 0.2550842000, 
        0.2551817000, 0.2552791000, 0.2553766000, 0.2554740000, 0.2555715000, 
        0.2556690000, 0.2557665000, 0.2558639000, 0.2559614000, 0.2560589000, 
        0.2561564000, 0.2562539000, 0.2563515000, 0.2564490000, 0.2565465000, 
        0.2566441000, 0.2567416000, 0.2568392000, 0.2569367000, 0.2570343000, 
        0.2571319000, 0.2572295000, 0.2573271000, 0.2574247000, 0.2575223000, 
        0.2576199000, 0.2577176000, 0.2578152000, 0.2579129000, 0.2580106000, 
        0.2581082000, 0.2582059000, 0.2583036000, 0.2584014000, 0.2584991000, 
        0.2585968000, 0.2586946000, 0.2587924000, 0.2588902000, 0.2589879000, 
        0.2590858000, 0.2591836000, 0.2592814000, 0.2593793000, 0.2594771000, 
        0.2595750000, 0.2596729000, 0.2597708000, 0.2598688000, 0.2599667000, 
        0.2600647000, 0.2601627000, 0.2602606000, 0.2603587000, 0.2604567000, 
        0.2605547000, 0.2606528000, 0.2607509000, 0.2608490000, 0.2609471000, 
        0.2610452000, 0.2611434000, 0.2612415000, 0.2613397000, 0.2614379000, 
        0.2615361000, 0.2616344000, 0.2617326000, 0.2618309000, 0.2619292000, 
        0.2620275000, 0.2621259000, 0.2622242000, 0.2623226000, 0.2624210000, 
        0.2625194000, 0.2626179000, 0.2627164000, 0.2628148000, 0.2629133000, 
        0.2630119000, 0.2631104000, 0.2632090000, 0.2633076000, 0.2634062000, 
        0.2635048000, 0.2636035000, 0.2637022000, 0.2638009000, 0.2638996000, 
        0.2639984000, 0.2640971000, 0.2641959000, 0.2642947000, 0.2643936000, 
        0.2644925000, 0.2645913000, 0.2646903000, 0.2647892000, 0.2648882000, 
        0.2649871000, 0.2650861000, 0.2651852000, 0.2652842000, 0.2653833000, 
        0.2654824000, 0.2655816000, 0.2656807000, 0.2657799000, 0.2658791000, 
        0.2659784000, 0.2660776000, 0.2661769000, 0.2662762000, 0.2663755000, 
        0.2664749000, 0.2665743000, 0.2666737000, 0.2667731000, 0.2668726000, 
        0.2669721000, 0.2670716000, 0.2671712000, 0.2672707000, 0.2673703000, 
        0.2674700000, 0.2675696000, 0.2676693000, 0.2677690000, 0.2678687000, 
        0.2679685000, 0.2680683000, 0.2681681000, 0.2682679000, 0.2683678000, 
        0.2684677000, 0.2685676000, 0.2686675000, 0.2687675000, 0.2688675000, 
        0.2689676000, 0.2690676000, 0.2691677000, 0.2692678000, 0.2693680000, 
        0.2694681000, 0.2695683000, 0.2696685000, 0.2697688000, 0.2698691000, 
        0.2699694000, 0.2700697000, 0.2701701000, 0.2702705000, 0.2703709000, 
        0.2704713000, 0.2705718000, 0.2706723000, 0.2707729000, 0.2708734000, 
        0.2709740000, 0.2710746000, 0.2711753000, 0.2712759000, 0.2713766000, 
        0.2714774000, 0.2715781000, 0.2716789000, 0.2717797000, 0.2718806000, 
        0.2719814000, 0.2720823000, 0.2721833000, 0.2722842000, 0.2723852000, 
        0.2724862000, 0.2725873000, 0.2726883000, 0.2727894000, 0.2728905000, 
        0.2729917000, 0.2730929000, 0.2731941000, 0.2732953000, 0.2733966000, 
        0.2734979000, 0.2735992000, 0.2737005000, 0.2738019000, 0.2739033000, 
        0.2740048000, 0.2741062000, 0.2742077000, 0.2743092000, 0.2744108000, 
        0.2745123000, 0.2746139000, 0.2747156000, 0.2748172000, 0.2749189000, 
        0.2750206000, 0.2751223000, 0.2752241000, 0.2753259000, 0.2754277000, 
        0.2755295000, 0.2756314000, 0.2757333000, 0.2758352000, 0.2759372000, 
        0.2760392000, 0.2761412000, 0.2762432000, 0.2763453000, 0.2764474000, 
        0.2765495000, 0.2766516000, 0.2767538000, 0.2768560000, 0.2769582000, 
        0.2770604000, 0.2771627000, 0.2772650000, 0.2773673000, 0.2774697000, 
        0.2775720000, 0.2776745000, 0.2777769000, 0.2778793000, 0.2779818000, 
        0.2780843000, 0.2781868000, 0.2782894000, 0.2783920000, 0.2784946000, 
        0.2785972000, 0.2786999000, 0.2788026000, 0.2789053000, 0.2790080000, 
        0.2791107000, 0.2792135000, 0.2793163000, 0.2794192000, 0.2795220000, 
        0.2796249000, 0.2797278000, 0.2798307000, 0.2799337000, 0.2800366000, 
        0.2801396000, 0.2802427000, 0.2803457000, 0.2804488000, 0.2805519000, 
        0.2806550000, 0.2807581000, 0.2808613000, 0.2809645000, 0.2810677000, 
        0.2811709000, 0.2812742000, 0.2813774000, 0.2814807000, 0.2815841000, 
        0.2816874000, 0.2817908000, 0.2818942000, 0.2819976000, 0.2821010000, 
        0.2822044000, 0.2823079000, 0.2824114000, 0.2825149000, 0.2826185000, 
        0.2827220000, 0.2828256000, 0.2829292000, 0.2830328000, 0.2831365000, 
        0.2832401000, 0.2833438000, 0.2834475000, 0.2835512000, 0.2836550000, 
        0.2837588000, 0.2838625000, 0.2839663000, 0.2840702000, 0.2841740000, 
        0.2842779000, 0.2843817000, 0.2844856000, 0.2845896000, 0.2846935000, 
        0.2847975000, 0.2849014000, 0.2850054000, 0.2851094000, 0.2852135000, 
        0.2853175000, 0.2854216000, 0.2855257000, 0.2856298000, 0.2857339000, 
        0.2858380000, 0.2859422000, 0.2860464000, 0.2861505000, 0.2862548000, 
        0.2863590000, 0.2864632000, 0.2865675000, 0.2866717000, 0.2867760000, 
        0.2868803000, 0.2869847000, 0.2870890000, 0.2871933000, 0.2872977000, 
        0.2874021000, 0.2875065000, 0.2876109000, 0.2877153000, 0.2878198000, 
        0.2879242000, 0.2880287000, 0.2881332000, 0.2882377000, 0.2883422000, 
        0.2884468000, 0.2885513000, 0.2886559000, 0.2887604000, 0.2888650000, 
        0.2889696000, 0.2890742000, 0.2891789000, 0.2892835000, 0.2893882000, 
        0.2894928000, 0.2895975000, 0.2897022000, 0.2898069000, 0.2899116000, 
        0.2900163000, 0.2901211000, 0.2902258000, 0.2903306000, 0.2904354000, 
        0.2905402000, 0.2906449000, 0.2907498000, 0.2908546000, 0.2909594000, 
        0.2910642000, 0.2911691000, 0.2912740000, 0.2913788000, 0.2914837000 
    };
    double width_K1;
    if(mass<1.35){
        width_K1=K1270_width[iii];
    }else{
        width_K1=K1400_width[iii];
    }
    a[0] = 1;
    a[1] = 0;
    b[0] = mass2-sa;
    //b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r2,l); 
    b[1] = -mass*width_K1;
    Com_Divide(a,b,prop);
}
void EvtD0ToKSKppimpi0::propagatorGS(double mass2, double mass, double width, double sa, double sb, double sc, double r2, double prop[2])
{
 
    double GS1 = 0.636619783;
    double GS2 = 0.01860182466;
    double GS3 = 0.1591549458;
    double GS4 = 0.00620060822;
    double a[2], b[2];
    double tmp  = sb-sc;
    double tmp1 = sa+tmp;
    double q2 = fabs(0.25*tmp1*tmp1/sa-sb);
    //if(q2<0) q2 = 1e-16;
 
    double tmp2 = mass2+tmp;
    double q02 = fabs(0.25*tmp2*tmp2/mass2-sb);
    //if(q02<0) q02 = 1e-16;
 
    double q  = sqrt(q2);
    double q0 = sqrt(q02);
    double m  = sqrt(sa);
    double q03 = q0*q02;
    double tmp3 = log(mass+2*q0)+1.2760418309;  // log(mass_2Pion) = 1.2760418309;
 
    double h  = GS1*q/m*(log(m+2*q)+1.2760418309);
    double h0 = GS1*q0/mass*tmp3;
    double dh = h0*(0.125/q02-0.5/mass2)+GS3/mass2;
    double d  = GS2/q02*tmp3+GS3*mass/q0-GS4*mass/q03;
    double f  = mass2/q03*(q2*(h-h0)+(mass2-sa)*q02*dh);
 
    a[0] = 1.0+d*width/mass;
    a[1] = 0.0;
    b[0] = mass2-sa+width*f;
    b[1] = -mass*width*widl1(mass2,mass,sa,sb,sc,r2);
    Com_Divide(a,b,prop);
}
void EvtD0ToKSKppimpi0::rhoab(double sa, double sb, double sc, double res[2]) {
    double tmp = sa+sb-sc;
    double q = 0.25*tmp*tmp/sa-sb;
    if(q>=0) {
        res[0]=2.0*sqrt(q/sa);
        res[1]=0.0;
    } else {
        res[0]=0.0;
        res[1]=2.0*sqrt(-q/sa);
    }
}
void EvtD0ToKSKppimpi0::rho4Pi(double sa, double res[2]) {
    double temp = 1.0-0.3116765584/sa;    // 0.3116765584=0.13957*0.13957*16
    if(temp>=0) {
        res[0]=sqrt(temp)/(1.0+exp(9.8-3.5*sa));
        res[1]=0.0;
    } else {
        res[0]=0.0;
        res[1]=sqrt(-temp)/(1.0+exp(9.8-3.5*sa));
    }
}
void EvtD0ToKSKppimpi0::propagatorsigma500(double sa, double sb, double sc, double prop[2]) {
    double f = 0.5843+1.6663*sa;
    const double M = 0.9264;
    const double mass2 = 0.85821696;  // M*M
    const double mpi2d2 = 0.00973989245;
    double g1 = f*(sa-mpi2d2)/(mass2-mpi2d2)*exp((mass2-sa)/1.082);
    double rho1s[2], rho1M[2], rho2s[2], rho2M[2], rho1[2], rho2[2];
    rhoab(sa,sb,sc,rho1s);
    rhoab(mass2,sb,sc,rho1M);
    rho4Pi(sa,rho2s);
    rho4Pi(mass2,rho2M);
    Com_Divide(rho1s,rho1M,rho1);
    Com_Divide(rho2s,rho2M,rho2);
    double a[2], b[2];
    a[0] = 1.0;
    a[1] = 0.0;
    b[0] = mass2-sa+M*(g1*rho1[1]+0.0024*rho2[1]);
    b[1] =         -M*(g1*rho1[0]+0.0024*rho2[0]);
    Com_Divide(a,b,prop);
}
void  EvtD0ToKSKppimpi0::Flatte_rhoab(double sa, double sb, double sc, double rho[2])
{
    double q = (sa+sb-sc)*(sa+sb-sc)/(4*sa)-sb;
    if(q>0) {
        rho[0]=2* sqrt(q/sa);
        rho[1]=0;
    }
    else if(q<0){
        rho[0]=0;
        rho[1]=2*sqrt(-q/sa);
    }
}
void EvtD0ToKSKppimpi0::propagator980(double mass, double sx, double *sb, double *sc, double prop[2])
{
    double unit[2]={1.0};
    double ci[2]={0,1};
    double rho1[2];
    Flatte_rhoab(sx,sb[0],sc[0],rho1);
    double rho2[2];
    Flatte_rhoab(sx,sb[1],sc[1],rho2);
    double gK_f980=0.69465, gPi_f980=0.165;
    double tmp1[2]={gK_f980,0};
    double tmp11[2];
    double tmp2[2]={gPi_f980,0};
    double tmp22[2];
    Com_Multi(tmp1,rho1,tmp11);
    Com_Multi(tmp2,rho2,tmp22);
    double tmp3[2]={tmp11[0]+tmp22[0],tmp11[1]+tmp22[1]};
    double tmp31[2];
    Com_Multi(tmp3, ci,tmp31);
    double tmp4[2]={mass*mass-sx-tmp31[0], -1.0*tmp31[1]};
    Com_Divide( unit,tmp4, prop);
}
void EvtD0ToKSKppimpi0::propagatorFlatte(double mass, double width, double sa, double prop[2]){
    double q2_Pi, q2_Ka;
    double rhoPi[2], rhoKa[2];
    q2_Pi = 0.25*sa-mPi*mPi;
    q2_Ka = 0.25*sa-mKa*mKa;
    if (q2_Pi > 0) {
        rhoPi[0] = 2.0*sqrt(q2_Pi/sa);
        rhoPi[1] = 0.0;
    }
    if (q2_Pi <= 0) {
        rhoPi[0] = 0.0;
        rhoPi[1] = 2.0*sqrt(-q2_Pi/sa);
    }
    if (q2_Ka > 0) {
        rhoKa[0] = 2.0*sqrt(q2_Ka/sa);
        rhoKa[1] = 0.0;
    }
    if (q2_Ka <= 0) {
        rhoKa[0] = 0.0;
        rhoKa[1] = 2.0*sqrt(-q2_Ka/sa);
    }
    double a[2], b[2];
    a[0] = 1;
    a[1] = 0;
    b[0] = mass*mass - sa + 0.165*rhoPi[1] + 0.69465*rhoKa[1];
    b[1] = - (0.165*rhoPi[0] + 0.69465*rhoKa[0]);
    Com_Divide(a,b,prop);
 
}
void EvtD0ToKSKppimpi0::KPiSLASS(double sa, double sb, double sc, double prop[2]) {
    const double m1430 = 1.441;
    const double sa0   = 2.076481;   // m1430*m1430;
    const double w1430 = 0.193;
    const double Lass1 = 0.25/sa0;
    double tmp = sb-sc;
    double tmp1 = sa0+tmp;
    double q0 = fabs(Lass1*tmp1*tmp1-sb);
    //if(q0<0) q0 = 1e-16;
    double tmp2 = sa+tmp;
    double qs = fabs(0.25*tmp2*tmp2/sa-sb);
    double q = sqrt(qs);
    double width = w1430*q*m1430/sqrt(sa*q0);
    double temp_R = atan(m1430*width/(sa0-sa));
    if(temp_R<0) temp_R += math_pi;
    double deltaR = -109.7*math_pi/180.0 + temp_R;
    double temp_F = atan(0.226*q/(2.0-3.8194*qs));  // 2.0*0.113 = 0.226; -33.8*0.113 = -3.8194
    if(temp_F<0) temp_F += math_pi;
    double deltaF = 0.1*math_pi/180.0 + temp_F;
    double deltaS = deltaR + 2.0*deltaF;
    double t1 = 0.96*sin(deltaF);
    double t2 = sin(deltaR);
    double CF[2], CS[2];
    CF[0] = cos(deltaF);
    CF[1] = sin(deltaF);
    CS[0] = cos(deltaS);
    CS[1] = sin(deltaS);
    prop[0] = t1*CF[0] + t2*CS[0];
    prop[1] = t1*CF[1] + t2*CS[1];
 
}