void PDF_D_hhpi0_dilution::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case cleo:{
obsValSource = "1410.3964";
setObservable("F_pipipi0_obs", 0.968);
setObservable("F_kkpi0_obs", 0.731);
break;
}
case cleo2015:{
obsValSource = "1504.05878";
setObservable("F_pipipi0_obs", 0.973);
setObservable("F_kkpi0_obs", 0.732);
break;
}
default:{
cout << "PDF_D_hhpi0_dilution::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_GLW_DK_2var::setObservables(config c)
{
switch(c)
{
case truth:
{
setObservablesTruth();
break;
}
case toy:
{
setObservablesToy();
break;
}
case cdf:
obsValSource = "FERMILAB-PUB-09-549-E"; // http://inspirehep.net/search?p=fin%20j%20PR,D81,031105
setObservable("acpp_dk_obs",0.39);
setObservable("rcpp_dk_obs",1.30);
break;
case lhcb:
obsValSource = "LHCb-PAPER-2012-001"; // http://inspirehep.net/search?p=fin+j+PL,B712,203
setObservable("acpp_dk_obs",0.145);
setObservable("rcpp_dk_obs",1.007);
break;
default:
cout << "PDF_GLW_DK_2var::setObservables : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
void PDF_Dmixing::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case lumi1fb:{
obsValSource = "1fb-1 LHCb-PAPER-2012-038";
setObservable("RD_obs", 3.52e-3);
setObservable("yprime_obs", 7.2e-3);
setObservable("xprimesq_obs",-0.09e-3);
break;
}
default:{
cout << "PDF_Dmixing::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_CrossCorTest1::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case manual:{
obsValSource = "manual";
break;
}
case lumi1fbPapercFit:{
obsValSource = "PAPER 25.6.2014 After bug fix, mail by Manuel";
setObservable("c_dsk_obs", 0.526386 );
setObservable("s_dsk_obs", -1.08902 );
setObservable("sb_dsk_obs", 0.356743 );
break;
}
default:{
cout << "PDF_CrossCorTest1::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_ADS_DKstz_Kpi::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case lumi3fb:{
obsValSource = "3fb-1 LHCb-PAPER-2014-028-v2, June 4";
setObservable("afav_dkstz_kpi_obs", -0.032); // -0.00626 0.02810 // A_d^Kpi
setObservable("rp_dkstz_obs", 0.057); // 0.08537 0.06271 // R_d^+
setObservable("rm_dkstz_obs", 0.056); // 0.11856 0.11395 // R_d^-
// truth values corresponding to: ^^ rB=0.3 ^^ rB=0.3
// dB=190 dB=270
break;
}
default:{
cout << "PDF_ADS_DKstz_Kpi::setObservables : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_D_CleoNoK2pi::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case cleo:{
obsValSource = "Cleo 2009 ARXIV:0903.4853";
setObservable("kD_k3pi_obs",0.33);
setObservable("dD_k3pi_obs",DegToRad(114.0));
setObservable("xD_obs", 0.96 /100.);
setObservable("yD_obs", 0.81 /100.);
setObservable("dD_kpi_obs", DegToRad(-151.5+360.));
setObservable("B1_obs", 3.89e-2);
setObservable("B2_obs", 1.47e-4);
setObservable("B3_obs", 7.96e-2);
setObservable("B4_obs", 2.65e-4);
break;
}
default:{
cout << "PDF_D_CleoNoK2pi::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_GGSZ_DstK::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
// for HFAG average
case babar:{
obsValSource = "PRL 105 (2010) 121801";
setObservable("xm_dstk_obs",-0.104);
setObservable("ym_dstk_obs",-0.052);
setObservable("xp_dstk_obs", 0.147);
setObservable("yp_dstk_obs",-0.032);
break;
}
case belle:{
obsValSource = "PRD 81 (2010) 112002";
setObservable("xm_dstk_obs",-0.036);
setObservable("ym_dstk_obs",-0.249);
setObservable("xp_dstk_obs", 0.083);
setObservable("yp_dstk_obs", 0.157);
break;
}
default:{
cout << "PDF_GGSZ_DstK::setObservables() : ERROR : config not found: " << ConfigToTString(c) << endl;
exit(1);
}
}
}
void PDF_GLWADS_DK_hh_newVars::setObservables(config c)
{
switch(c)
{
case truth: {
setObservablesTruth();
break;
}
case toy: {
setObservablesToy();
break;
}
case lumi3fb: {
obsValSource = "3fb-1 ANA v7 - lifted from gammadini "; // https://twiki.cern.ch/twiki/pub/LHCbPhysics/B2D0K/LHCb-ANA-2014-071-v7.pdf (see Vavas email 04/08/15)
setObservable( "aads_dk_kpi_obs", -4.0267e-01 );
setObservable( "acp_dk_kk_obs", 8.6577e-02 );
setObservable( "acp_dk_pipi_obs", 1.2760e-01 );
setObservable( "afav_dk_kpi_obs", -1.9441e-02 );
setObservable( "rads_dk_kpi_obs", 1.8757e-02 );
setObservable( "rcp_kk_obs", 9.6834e-01 );
setObservable( "rcp_pipi_obs", 1.0022e+00 );
break;
}
default: {
cout << "PDF_GLWADS_DK_hh_newVars::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_D_KSKpi_noDmix::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case cleoFullDP:{
obsValSource = "CLEO, 1203.3804, PRD 85 092016, whole Dalitz region";
setObservable("RD_kskpi_obs",0.592);
setObservable("dD_kskpi_obs",DegToRad(8.3));
setObservable("kD_kskpi_obs",0.73);
break;
}
case cleo:{
obsValSource = "CLEO, 1203.3804, PRD 85 092016, K* region of Dalitz Plot";
setObservable("RD_kskpi_obs",0.356);
setObservable("dD_kskpi_obs",DegToRad(26.5));
setObservable("kD_kskpi_obs",1.0);
break;
}
default:{
cout << "PDF_D_KSKpi_noDmix::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_GLW_DK_pipipi0_2var::setObservables(config c)
{
switch(c)
{
case truth:
{
setObservablesTruth();
break;
}
case toy:
{
setObservablesToy();
break;
}
case lhcb:
obsValSource = "arXiv:1504.05442"; // http://arxiv.org/pdf/1504.05442v2.pdf
setObservable("aglw_dk_pipipi0_obs",0.05);
setObservable("rglw_dk_pipipi0_obs",0.98);
break;
default:
cout << "PDF_GLW_DK_pipipi0_2var::setObservables : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
void PDF_GLW_Dpi::setObservables(config c)
{
switch(c)
{
case truth:
{
setObservablesTruth();
break;
}
case toy:
{
setObservablesToy();
break;
}
case lumi1fb:
obsValSource = "1fb-1 paper LHCb-PAPER-2012-001";
setObservable("acpp_dpi_kk_obs",-0.020);
setObservable("acpp_dpi_pipi_obs",-0.001);
break;
default:
cout << "PDF_GLW_DzKstz::setObservables : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
void PDF_ADS_DK_Kpi_2var::setObservables(config c)
{
switch(c)
{
case truth:
{
setObservablesTruth();
break;
}
case toy:
{
setObservablesToy();
break;
}
case babar:
obsValSource = "PRD 82 (2010) 072006";
setObservable("aads_dk_obs", -0.86);
setObservable("rads_dk_obs", 0.011);
break;
case belle:
obsValSource = "PRL 106 (2011) 231803";
setObservable("aads_dk_obs", -0.390 );
setObservable("rads_dk_obs", 0.0163 );
break;
case cdf:
obsValSource = "PRD 84 (2011) 091504";
setObservable("aads_dk_obs", -0.82 );
setObservable("rads_dk_obs", 0.0220 );
break;
case lhcb:
obsValSource = "PLB 712 (2012) 203";
setObservable("aads_dk_obs", -0.52 );
setObservable("rads_dk_obs", 0.0152 );
break;
default:
cout << "PDF_ADS_DK_Kpi_2var::setObservables : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
void PDF_D_HFAGfit_4D::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case hfagFP2014:{
// http://www.slac.stanford.edu/xorg/hfag/charm/FPCP14/results_mix+cpv.html
obsValSource = "HFAG FP2014, CPV allowed";
setObservable("dD_kpi_obs", DegToRad(7.3+180.0));
setObservable("rD_kpi_obs", 0.349e-2);
setObservable("AcpDpipi_obs", 0.14e-2);
setObservable("AcpDKK_obs", -0.11e-2);
break;
}
case hfagCHARM2015:{
// http://www.slac.stanford.edu/xorg/hfag/charm/CHARM15/results_mix_cpv.html
obsValSource = "HFAG CHARM2015, CPV allowed";
setObservable("dD_kpi_obs", DegToRad(11.8+180.0));
setObservable("rD_kpi_obs", 0.349e-2);
setObservable("AcpDpipi_obs", 0.10e-2);
setObservable("AcpDKK_obs", -0.15e-2);
break;
}
default:{
cout << "PDF_D_HFAGfit_4D::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_gamma::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case year2014:{
obsValSource = "year2014";
setObservable("g_obs",70.0/57.3);
break;
}
default:{
cout << "PDF_gamma::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_D_KSKpi_RWS::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case lhcb:{
obsValSource = "arXiv:1509.06628";
setObservable("RD_kskpi_obs",0.370);
//setObservable("kD_kskpi_obs",0.571);
break;
}
default:{
cout << "PDF_D_KSKpi_RWS::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_Dmixing_xy::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case hfagLP2011:{
obsValSource = "HFAG LP2011 CPV allowed http://www.slac.stanford.edu/xorg/hfag/charm/LP11/results_mix+cpv.html";
setObservable("xD_obs", 0.0063);
setObservable("yD_obs", 0.0075);
break;
}
case hfag:{
obsValSource = "HFAG 2014 CPV allowed http://www.slac.stanford.edu/xorg/hfag/charm/CHARM13/results_mix+cpv.html";
setObservable("xD_obs", 0.0039);
setObservable("yD_obs", 0.0067);
break;
}
case lhcb2012:{
//added by sam harnew - the above numbers seem to be the Nov2013 results at http://www.slac.stanford.edu/xorg/hfag/charm/CHARM13/results_mix_cpv.html
obsValSource = "HFAG 2014 CPV allowed"; // (From website, above looks to be 2013???!) http://www.slac.stanford.edu/xorg/hfag/charm/FPCP14/results_mix_cpv.html";
setObservable("xD_obs", 0.00414);
setObservable("yD_obs", 0.00627);
break;
}
default:{
cout << "PDF_Dmixing_xy::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_ADS_DK_Kpi::setObservables(config c)
{
switch(c)
{
case truth:
{
setObservablesTruth();
break;
}
case toy:
{
setObservablesToy();
break;
}
case lumi1fb:
{
obsValSource = "1fb-1 LHCb-PAPER-2012-001";
if ( cTrad == useTradObs )
{
setObservable("rads_dk_kpi_obs", 0.0152);
setObservable("aads_dk_kpi_obs", -0.52);
setObservable("afav_dk_kpi_obs", 0.0044);
}
else
{
setObservable("rp_dk_kpi_obs",0.0232);
setObservable("rm_dk_kpi_obs",0.0073);
setObservable("afav_dk_kpi_obs",0.0044);
}
break;
}
default:
cout << "PDF_ADS_DK_Kpi::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
void PDF_Dmixing_CLEO2D::setObservables(config c)
{
switch(c)
{
case truth: {
setObservablesTruth();
break;
}
case toy: {
setObservablesToy();
break;
}
case cleo2009: {
obsValSource = "arXiv:0903.4853 - CLEO 2009";
setObservable("dD_k3pi_obs",DegToRad(114.0));
setObservable("kD_k3pi_obs",0.33);
break;
}
case cleo2014: {
obsValSource = "arXiv:1401.1904 - CLEO 2013";
setObservable("dD_k3pi_obs",DegToRad(255.0)); // CLEO 2014
setObservable("kD_k3pi_obs",0.32); // CLEO 2014
break;
}
case cleo2015: {
obsValSource = "Guy";
setObservable("dD_k3pi_obs",DegToRad(170.0)); // CLEO 2015
setObservable("kD_k3pi_obs",0.32); // CLEO 2015
break;
}
default: {
cout << "PDF_Dmixing_CLEO2D::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_ADS_DKDpi_K3pi_Dmix_newVars::setObservables(config c)
{
switch(c)
{
case truth:
{
setObservablesTruth();
break;
}
case toy:
{
setObservablesToy();
break;
}
case lumi3fb:
{
obsValSource = "3fb-1 ANA v7 - lifted from gammadini "; // https://twiki.cern.ch/twiki/pub/LHCbPhysics/B2D0K/LHCb-ANA-2014-071-v7.pdf (see Vavas email 04/08/15)
setObservable( "aads_dk_k3pi_obs", -3.1291e-01 );
setObservable( "aads_dpi_k3pi_obs", 2.3228e-02 );
setObservable( "afav_dk_k3pi_obs", -4.1540e-04 );
setObservable( "rads_dk_k3pi_obs", 1.3992e-02 );
setObservable( "rads_dpi_k3pi_obs", 3.7718e-03 );
setObservable( "rkp_k3pi_obs", 7.9306e-02 );
break;
}
case check:
{
obsValSource = "check";
setObservable( "aads_dk_k3pi_obs", -3.1291e-01 );
setObservable( "aads_dpi_k3pi_obs", 2.3228e-02 );
setObservable( "afav_dk_k3pi_obs", -4.1540e-04 );
setObservable( "rads_dk_k3pi_obs", 1.3992e-02 );
//setObservable( "rads_dpi_k3pi_obs", 3.7718e-03 ); // THE REAL VALUE
setObservable( "rads_dpi_k3pi_obs", 3.39e-03 ); // check value
setObservable( "rkp_k3pi_obs", 7.9306e-02 );
break;
}
default:
{
cout << "PDF_ADS_DKDpi_K3pi_Dmix_newVars::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
Observer::Observer(Observable *o)
: _observable(0)
{
setObservable( o );
}
void PDF_GLWADS_DKDpi_KSKpi::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case lumi3fbFullDP:{
obsValSource = "3fb-1, Dan's presentation at LHCb Week Dec 2013, whole Dalitz plot";
setObservable("rfavsup_dpi_kskpi_obs", 1.528);
setObservable("rfav_dkdpi_kskpi_obs", 0.092);
setObservable("rsup_dkdpi_kskpi_obs", 0.066);
setObservable("afav_dk_kskpi_obs", 0.040);
setObservable("asup_dk_kskpi_obs", 0.233);
setObservable("afav_dpi_kskpi_obs", -0.025);
setObservable("asup_dpi_kskpi_obs", -0.052);
break;
}
case lumi3fb:{
obsValSource = "3fb-1, paper draft as of LHCb Week Dec 2013, K* region of Dalitz plot";
setObservable("rfavsup_dpi_kskpi_obs", 2.57);
setObservable("rfav_dkdpi_kskpi_obs", 0.084);
setObservable("rsup_dkdpi_kskpi_obs", 0.056);
setObservable("afav_dk_kskpi_obs", 0.026); // SS
setObservable("asup_dk_kskpi_obs", 0.336); // OS
setObservable("afav_dpi_kskpi_obs", -0.012);
setObservable("asup_dpi_kskpi_obs", -0.054);
break;
}
default:
{
cout << "PDF_GLWADS_DKDpi_KSKpi::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_GGSZ_DKstz::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case lumi3fb:{
obsValSource = "Alexis' estimate"; // pers. comm. 10/6/2015
setObservable("xm_dkstz_obs", -0.03);
setObservable("ym_dkstz_obs", -0.24);
setObservable("xp_dkstz_obs",0.16);
setObservable("yp_dkstz_obs",0.18);
break;
}
case lhcb:{
obsValSource = "MI GGSZ from Sneha's email"; // Sneha's mail 22/01/16
setObservable("xm_dkstz_obs", -0.32);
setObservable("ym_dkstz_obs", 0.30);
setObservable("xp_dkstz_obs", 0.04);
setObservable("yp_dkstz_obs", -0.80);
break;
}
// for HFAG
case belle:{
obsValSource = "arXiv:1509.01098";
setObservable("xm_dkstz_obs", 0.4);
setObservable("ym_dkstz_obs",-0.6);
setObservable("xp_dkstz_obs", 0.1);
setObservable("yp_dkstz_obs", 0.3);
break;
}
default:{
cout << "PDF_GGSZ_DKstz::setObservables() : ERROR : config not found: " << ConfigToTString(c) << endl;
exit(1);
}
}
}
void PDF_GLWADS_DKDpi_hh_2012::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case lumi3fb:{
obsValSource = "3fb-1 ANA v2, ADS values still blind (using Asimov/HFAG/1fb values)";
setObservable("rkp_kpi_obs", 0.0785 );
setObservable("rcp_pipi_obs", 1.045 );
setObservable("rcp_kk_obs", 0.981 );
setObservable("afav_dpi_kpi_obs", -0.0033); // now Aprod
setObservable("afav_dk_kpi_obs", -0.0158);
setObservable("acp_dpi_kk_obs", -0.0162);
setObservable("acp_dpi_pipi_obs", 0.00315);
setObservable("acp_dk_kk_obs", 0.0844 );
setObservable("acp_dk_pipi_obs", 0.117 );
setObservable("rads_dk_kpi_obs", 0.015 ); // HFAG 2014
setObservable("aads_dk_kpi_obs", -0.54 ); // HFAG 2014
setObservable("rads_dpi_kpi_obs", 0.0035 ); // Asimov r_dpi=0.005
setObservable("aads_dpi_kpi_obs", 0.12 ); // Asimov r_dpi=0.005
break;
}
default:{
cout << "PDF_GLWADS_DKDpi_hh_2012::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}
void PDF_GLWADS_DKDpi_hh_Dmix_newVars_uv::setObservables(config c)
{
switch(c)
{
case truth:{
setObservablesTruth();
break;
}
case toy:{
setObservablesToy();
break;
}
case lumi3fb:{
obsValSource = "3fb-1 ANA v7 - lifted from gammadini "; // https://twiki.cern.ch/twiki/pub/LHCbPhysics/B2D0K/LHCb-ANA-2014-071-v7.pdf (see Vavas email 04/08/15)
setObservable( "aads_dk_kpi_obs", -4.0267e-01 );
setObservable( "aads_dpi_kpi_obs", 1.0031e-01 );
setObservable( "acp_dk_kk_obs", 8.6577e-02 );
setObservable( "acp_dk_pipi_obs", 1.2760e-01 );
setObservable( "acp_dpi_kk_obs", -1.4524e-02 );
setObservable( "acp_dpi_pipi_obs", 4.2984e-03 );
setObservable( "afav_dk_kpi_obs", -1.9441e-02 );
setObservable( "rads_dk_kpi_obs", 1.8757e-02 );
setObservable( "rads_dpi_kpi_obs", 3.5977e-03 );
setObservable( "rcp_kk_obs", 9.6834e-01 );
setObservable( "rcp_pipi_obs", 1.0022e+00 );
setObservable( "rkp_kpi_obs", 7.7937e-02 );
break;
}
case check:{
obsValSource = "check";
setObservable( "aads_dk_kpi_obs", -4.0267e-01 );
setObservable( "aads_dpi_kpi_obs", 1.0031e-01 );
//setObservable( "acp_dk_kk_obs", 8.6577e-02 ); // THE REAL VALUE
setObservable( "acp_dk_kk_obs", 0.12497 ); // check
setObservable( "acp_dk_pipi_obs", 1.2760e-01 );
setObservable( "acp_dpi_kk_obs", -1.4524e-02 ); // THE REAL VALUE
setObservable( "acp_dpi_pipi_obs", 4.2984e-03 ); // THE REAL VALUE
//setObservable( "acp_dpi_kk_obs", -0.0087 ); // check
//setObservable( "acp_dpi_pipi_obs",-0.0062 ); // check
setObservable( "afav_dk_kpi_obs", -1.9441e-02 );
setObservable( "rads_dk_kpi_obs", 1.8757e-02 );
setObservable( "rads_dpi_kpi_obs", 3.5977e-03 ); // THE REAL VALUE
//setObservable( "rads_dpi_kpi_obs", 3.75e-03 ); // check
setObservable( "rcp_kk_obs", 9.6834e-01 );
setObservable( "rcp_pipi_obs", 1.0022e+00 );
setObservable( "rkp_kpi_obs", 7.7937e-02 );
break;
}
default:{
cout << "PDF_GLWADS_DKDpi_hh_Dmix_newVars_uv::setObservables() : ERROR : config "+ConfigToTString(c)+" not found." << endl;
exit(1);
}
}
}