//==========================================================================
// Class member functions for calculating the matrix elements for
%(process_lines)s

//--------------------------------------------------------------------------
// Initialize process. 
  
void %(process_class_name)s::initProc() {
// Instantiate the model class and set parameters that stay fixed during run
    pars = Parameters_%(model_name)s::getInstance();
    pars->setIndependentParameters(particleDataPtr, couplingsPtr, slhaPtr);
    pars->setIndependentCouplings();
    %(initProc_lines)s
} 

//--------------------------------------------------------------------------
// Evaluate |M|^2, part independent of incoming flavour. 

void %(process_class_name)s::sigmaKin() { 
    // Set the parameters which change event by event
    pars->setDependentParameters(particleDataPtr, couplingsPtr, slhaPtr, alpS);
    pars->setDependentCouplings();
    // Reset color flows
    %(reset_jamp_lines)s
    %(sigmaKin_lines)s
}

//--------------------------------------------------------------------------
// Evaluate |M|^2, including incoming flavour dependence. 

double %(process_class_name)s::sigmaHat() {  
    %(sigmaHat_lines)s
}

//--------------------------------------------------------------------------
// Select identity, colour and anticolour.

void %(process_class_name)s::setIdColAcol() {
    %(setIdColAcol_lines)s
}

//--------------------------------------------------------------------------
// Evaluate weight for angles of decay products in process 

double %(process_class_name)s::weightDecay(Event& process, int iResBeg, int iResEnd) {
    %(weightDecay_lines)s
}

//==========================================================================
// Private class member functions

//--------------------------------------------------------------------------
// Evaluate |M|^2 for each subprocess

%(all_sigmaKin)s