New method resetCoefficient() defined to set Wilson coefficients to 0. setMu(double mu) still sets the Wilson coefficients to 0 but now uses resetCoefficient().
Some minor bugs in Y1(x, mu) fixed. WC for B to mu mu now includes sin(\theta_W)^2.
IMPORTANT: UPDATE YOUR CONF FILES. QCD.cpp now reads the lifetime of B_s and B_d as input from conf files and not the width. Conversion to width done after reading conf file.
B_s to mu mu and B_d to mu mu implemented in full glory. Waiting for completion of Delta B = 2. Both time integrated and t = 0 rate implemented in same observable file BR_B(s/d)mumu and separately defined in ThFactory.
Minor modification due to modifications in Y1(x,mu) in StandardModelMatching
One last modification to make B_q to mu mu work