diff --git a/test/jsonFiles/BC_VHAD=VLL__Bc_Jpsipipi0_mumu.json b/test/jsonFiles/BC_VHAD=VLL__Bc_Jpsipipi0_mumu.json
index 25ca80a..79b6f10 100644
--- a/test/jsonFiles/BC_VHAD=VLL__Bc_Jpsipipi0_mumu.json
+++ b/test/jsonFiles/BC_VHAD=VLL__Bc_Jpsipipi0_mumu.json
@@ -1,38 +1,38 @@
 {
     "parent" : "B_c+",
     "daughters" : ["J/psi", "pi+", "pi0"],
     "grand_daughters" : [["mu+", "mu-"], [], []],
     "models" : ["BC_VHAD", "VLL", "", ""],
     "parameters" : [["1"], [], [], []],
     "extras" : ["noPhotos"],
     "events" : 10000,
     "histograms" : [
         {"variable" : "prob", "title" : "Prob(BC_VHAD)",  "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 1.0},
         {"variable" : "prob_daug1", "title" : "Prob(VLL)",  "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 1.0},
         {"variable" : "mass", "title" : "m(Jpsi)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 3.0, "xmax" : 3.2},
         {"variable" : "mass", "title" : "m(Jpsi pi+)", "d1" : 1, "d2" : 2, "nbins" : 100, "xmin" : 3.0, "xmax" : 6.5},
         {"variable" : "mass", "title" : "m(Jpsi pi0)", "d1" : 1, "d2" : 3, "nbins" : 100, "xmin" : 3.0, "xmax" : 6.5},
         {"variable" : "mass", "title" : "m(pi+ pi0)", "d1" : 2, "d2" : 3, "nbins" : 100, "xmin" : 0.0, "xmax" : 3.0},
         {"variable" : "pSumSq", "title" : "qSq(pi+ pi0)", "d1" : 2, "d2" : 3, "nbins" : 100, "xmin" : 0.0, "xmax" : 4.0},
         {"variable" : "p", "title" : "p(Jpsi)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.5, "xmax" : 2.5},
         {"variable" : "pz", "title" : "pz(Jpsi)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -3.0, "xmax" : 3.0},
         {"variable" : "cosHel", "title" : "cosHel(Jpsi,pi+)", "d1" : 1, "d2" : 2, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
         {"variable" : "cosHel", "title" : "cosHel(pi+,pi0)", "d1" : 2, "d2" : 3, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
         {"variable" : "cosHel", "title" : "cosHel(Jpsi,pi0)", "d1" : 1, "d2" : 3, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
         {"variable" : "cosHel_daug1", "title" : "cosHel(mu+,Jpsi)", "d1" : 1, "d2" : 2, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
         {"variable" : "cosHel_daug1", "title" : "cosHel(mu-,Jpsi)", "d1" : 2, "d2" : 1, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
         {"variable" : "cosAcoplanarityAngle", "title" : "cosAcopl", "d1" : 1, "d2" : 2, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
         {"variable" : "acoplanarityAngle", "title" : "acopl", "d1" : 1, "d2" : 2, "nbins" : 100, "xmin" : 0, "xmax" : 180.0},
         {"variable" : "cosTheta", "title" : "cosTheta(Jpsi)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
 	{"variable" : "E", "title" : "E(Jpsi)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 3.5, "xmax" : 4.2},
 	{"variable" : "phi", "title" : "phi(Jpsi)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
 	{"variable" : "pLab_daug1", "title" : "pLab(mu+)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 0.5, "xmax" : 3.5},
 	{"variable" : "pLab_daug1", "title" : "pLab(mu-)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 0.5, "xmax" : 3.5},
 	{"variable" : "cosTheta_daug1", "title" : "cosTheta(mu+)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
 	{"variable" : "phi_daug1", "title" : "phi(mu+)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
 	{"variable" : "cosTheta_daug1", "title" : "cosTheta(mu-)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
 	{"variable" : "phi_daug1", "title" : "phi(mu-)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0}
     ],
     "outfile" : "BC_VHAD=VLL__Bc_Jpsipipi0_mumu.root",
-    "reference" : "RefBC_VHAD=VLL__Bc_JpsiKpipi0_mumu.root"
+    "reference" : "RefBC_VHAD=VLL__Bc_Jpsipipi0_mumu.root"
 }