diff --git a/EvtGenModels/EvtBHadronic.hh b/EvtGenModels/EvtBHadronic.hh
deleted file mode 100644
--- a/EvtGenModels/EvtBHadronic.hh
+++ /dev/null
@@ -1,37 +0,0 @@
-
-/***********************************************************************
-* Copyright 1998-2020 CERN for the benefit of the EvtGen authors *
-* *
-* This file is part of EvtGen. *
-* *
-* EvtGen is free software: you can redistribute it and/or modify *
-* it under the terms of the GNU General Public License as published by *
-* the Free Software Foundation, either version 3 of the License, or *
-* (at your option) any later version. *
-* *
-* EvtGen is distributed in the hope that it will be useful, *
-* but WITHOUT ANY WARRANTY; without even the implied warranty of *
-* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
-* GNU General Public License for more details. *
-* *
-* You should have received a copy of the GNU General Public License *
-* along with EvtGen. If not, see . *
-***********************************************************************/
-
-#ifndef EVTBHADRONIC_HH
-#define EVTBHADRONIC_HH
-
-#include "EvtGenBase/EvtDecayAmp.hh"
-
-class EvtParticle;
-
-class EvtBHadronic : public EvtDecayAmp {
- public:
- std::string getName() override;
- EvtBHadronic* clone() override;
-
- void init() override;
- void decay( EvtParticle* p ) override;
-};
-
-#endif
diff --git a/doc/evt_models.tex b/doc/evt_models.tex
--- a/doc/evt_models.tex
+++ b/doc/evt_models.tex
@@ -12,19 +12,6 @@
what the arguments mean and one or more examples are given. For
further examples of use please see the decay table, {\tt DECAY.DEC}
-\Model{BHADRONIC}
-
-\label{bhadronic}
-
-\Auth{Ryd}
-
-\Usage{P1 P2 ... PN}{JH JW;}
-
-\Expl
-This is an experimental model for hadronic $B$ decays. Until
-further developed this is not recommended to be used. For
-questions ask Anders Ryd.
-
\Model{BTO3PI\_CP}
\label{bto3picp}
diff --git a/src/EvtGenModels/EvtBHadronic.cpp b/src/EvtGenModels/EvtBHadronic.cpp
deleted file mode 100644
--- a/src/EvtGenModels/EvtBHadronic.cpp
+++ /dev/null
@@ -1,299 +0,0 @@
-
-/***********************************************************************
-* Copyright 1998-2020 CERN for the benefit of the EvtGen authors *
-* *
-* This file is part of EvtGen. *
-* *
-* EvtGen is free software: you can redistribute it and/or modify *
-* it under the terms of the GNU General Public License as published by *
-* the Free Software Foundation, either version 3 of the License, or *
-* (at your option) any later version. *
-* *
-* EvtGen is distributed in the hope that it will be useful, *
-* but WITHOUT ANY WARRANTY; without even the implied warranty of *
-* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
-* GNU General Public License for more details. *
-* *
-* You should have received a copy of the GNU General Public License *
-* along with EvtGen. If not, see . *
-***********************************************************************/
-
-#include "EvtGenModels/EvtBHadronic.hh"
-
-#include "EvtGenBase/EvtGenKine.hh"
-#include "EvtGenBase/EvtPDL.hh"
-#include "EvtGenBase/EvtParticle.hh"
-#include "EvtGenBase/EvtPatches.hh"
-#include "EvtGenBase/EvtReport.hh"
-#include "EvtGenBase/EvtTensor4C.hh"
-#include "EvtGenBase/EvtVector4C.hh"
-#include "EvtGenBase/EvtVector4R.hh"
-
-#include "EvtGenModels/EvtISGW2FF.hh"
-
-#include
-#include
-using std::endl;
-
-std::string EvtBHadronic::getName()
-{
- return "BHADRONIC";
-}
-
-EvtBHadronic* EvtBHadronic::clone()
-{
- return new EvtBHadronic;
-}
-
-void EvtBHadronic::init()
-{
- // check that there are 2 argument
- checkNArg( 2 );
-}
-
-void EvtBHadronic::decay( EvtParticle* p )
-{
- //added by Lange Jan4,2000
- static EvtId B0 = EvtPDL::getId( "B0" );
- static EvtId D0 = EvtPDL::getId( "D0" );
- static EvtId DST0 = EvtPDL::getId( "D*0" );
- static EvtId D3P10 = EvtPDL::getId( "D'_10" );
- static EvtId D3P20 = EvtPDL::getId( "D_2*0" );
- static EvtId D3P00 = EvtPDL::getId( "D_0*0" );
- static EvtId D1P10 = EvtPDL::getId( "D_10" );
-
- static EvtISGW2FF ffmodel;
-
- p->initializePhaseSpace( getNDaug(), getDaugs() );
-
- EvtVector4R p4[MAX_DAUG];
- p->mass();
-
- int i, j;
-
- for ( i = 0; i < getNDaug(); i++ ) {
- p4[i] = p->getDaug( i )->getP4();
- }
-
- int bcurrent, wcurrent;
- int nbcurrent = 0;
- int nwcurrent = 0;
-
- bcurrent = (int)getArg( 0 );
- wcurrent = (int)getArg( 1 );
-
- EvtVector4C jb[5], jw[5];
- EvtTensor4C g, tds;
-
- EvtVector4R p4b;
- p4b.set( p->mass(), 0.0, 0.0, 0.0 );
-
- EvtVector4R q;
- double q2;
-
- EvtComplex ep_meson_bb[5];
- double f, gf, ap, am;
- double fp, fm;
- double hf, kf, bp, bm;
- EvtVector4R pp, pm;
- EvtVector4C ep_meson_b[5];
-
- switch ( bcurrent ) {
- // D0
- case 1:
- q = p4b - p4[0];
- q2 = q * q;
- nbcurrent = 1;
- ffmodel.getscalarff( B0, D0, EvtPDL::getMeanMass( D0 ),
- EvtPDL::getMeanMass( getDaugs()[1] ), &fp, &fm );
- jb[0] = EvtVector4C( fp * ( p4b + p4[0] ) - fm * q );
- break;
- // D*
- case 2:
- q = p4b - p4[0];
- q2 = q * q;
- nbcurrent = 3;
- ffmodel.getvectorff( B0, DST0, EvtPDL::getMeanMass( DST0 ), q2, &f,
- &gf, &ap, &am );
-
- g.setdiag( 1.0, -1.0, -1.0, -1.0 );
- tds = -f * g -
- ap * ( EvtGenFunctions::directProd( p4b, p4b ) +
- EvtGenFunctions::directProd( p4b, p4[0] ) ) -
- gf * EvtComplex( 0.0, 1.0 ) *
- dual( EvtGenFunctions::directProd( p4[0] + p4b,
- p4b - p4[0] ) ) -
- am * ( ( EvtGenFunctions::directProd( p4b, p4b ) -
- EvtGenFunctions::directProd( p4b, p4[0] ) ) );
- jb[0] = tds.cont1( p->getDaug( 0 )->epsParent( 0 ).conj() );
- jb[1] = tds.cont1( p->getDaug( 0 )->epsParent( 1 ).conj() );
- jb[2] = tds.cont1( p->getDaug( 0 )->epsParent( 2 ).conj() );
- break;
- // D1
- case 3:
- q = p4b - p4[0];
- q2 = q * q;
- nbcurrent = 3;
- ffmodel.getvectorff( B0, D3P10, EvtPDL::getMeanMass( D3P10 ), q2,
- &f, &gf, &ap, &am );
-
- g.setdiag( 1.0, -1.0, -1.0, -1.0 );
- tds = -f * g -
- ap * ( EvtGenFunctions::directProd( p4b, p4b ) +
- EvtGenFunctions::directProd( p4b, p4[0] ) ) -
- gf * EvtComplex( 0.0, 1.0 ) *
- dual( EvtGenFunctions::directProd( p4[0] + p4b,
- p4b - p4[0] ) ) -
- am * ( ( EvtGenFunctions::directProd( p4b, p4b ) -
- EvtGenFunctions::directProd( p4b, p4[0] ) ) );
- jb[0] = tds.cont1( p->getDaug( 0 )->epsParent( 0 ).conj() );
- jb[1] = tds.cont1( p->getDaug( 0 )->epsParent( 1 ).conj() );
- jb[2] = tds.cont1( p->getDaug( 0 )->epsParent( 2 ).conj() );
-
- break;
- // D2*
- case 4:
- q = p4b - p4[0];
- q2 = q * q;
- nbcurrent = 5;
- ffmodel.gettensorff( B0, D3P20, EvtPDL::getMeanMass( D3P20 ), q2,
- &hf, &kf, &bp, &bm );
- g.setdiag( 1.0, -1.0, -1.0, -1.0 );
-
- ep_meson_b[0] =
- ( ( p->getDaug( 0 )->epsTensorParent( 0 ) ).cont2( p4b ) ).conj();
- ep_meson_b[1] =
- ( ( p->getDaug( 0 )->epsTensorParent( 1 ) ).cont2( p4b ) ).conj();
- ep_meson_b[2] =
- ( ( p->getDaug( 0 )->epsTensorParent( 2 ) ).cont2( p4b ) ).conj();
- ep_meson_b[3] =
- ( ( p->getDaug( 0 )->epsTensorParent( 3 ) ).cont2( p4b ) ).conj();
- ep_meson_b[4] =
- ( ( p->getDaug( 0 )->epsTensorParent( 4 ) ).cont2( p4b ) ).conj();
-
- pp = p4b + p4[0];
- pm = p4b - p4[0];
-
- ep_meson_bb[0] = ep_meson_b[0] * ( p4b );
- ep_meson_bb[1] = ep_meson_b[1] * ( p4b );
- ep_meson_bb[2] = ep_meson_b[2] * ( p4b );
- ep_meson_bb[3] = ep_meson_b[3] * ( p4b );
- ep_meson_bb[4] = ep_meson_b[4] * ( p4b );
-
- jb[0] = EvtComplex( 0.0, ( 1.0 ) * hf ) *
- dual( EvtGenFunctions::directProd( pp, pm ) )
- .cont2( ep_meson_b[0] ) -
- kf * ep_meson_b[0] - bp * ep_meson_bb[0] * pp -
- bm * ep_meson_bb[0] * pm;
-
- jb[1] = EvtComplex( 0.0, ( 1.0 ) * hf ) *
- dual( EvtGenFunctions::directProd( pp, pm ) )
- .cont2( ep_meson_b[1] ) -
- kf * ep_meson_b[1] - bp * ep_meson_bb[1] * pp -
- bm * ep_meson_bb[1] * pm;
-
- jb[2] = EvtComplex( 0.0, ( 1.0 ) * hf ) *
- dual( EvtGenFunctions::directProd( pp, pm ) )
- .cont2( ep_meson_b[2] ) -
- kf * ep_meson_b[2] - bp * ep_meson_bb[2] * pp -
- bm * ep_meson_bb[2] * pm;
-
- jb[3] = EvtComplex( 0.0, ( 1.0 ) * hf ) *
- dual( EvtGenFunctions::directProd( pp, pm ) )
- .cont2( ep_meson_b[3] ) -
- kf * ep_meson_b[3] - bp * ep_meson_bb[3] * pp -
- bm * ep_meson_bb[3] * pm;
-
- jb[4] = EvtComplex( 0.0, ( 1.0 ) * hf ) *
- dual( EvtGenFunctions::directProd( pp, pm ) )
- .cont2( ep_meson_b[4] ) -
- kf * ep_meson_b[4] - bp * ep_meson_bb[4] * pp -
- bm * ep_meson_bb[4] * pm;
- break;
- // D_0*
- case 5:
- q = p4b - p4[0];
- q2 = q * q;
- double f, gf, ap, am;
- nbcurrent = 3;
- ffmodel.getvectorff( B0, D1P10, EvtPDL::getMeanMass( D1P10 ), q2,
- &f, &gf, &ap, &am );
- g.setdiag( 1.0, -1.0, -1.0, -1.0 );
- tds = -f * g -
- ap * ( EvtGenFunctions::directProd( p4b, p4b ) +
- EvtGenFunctions::directProd( p4b, p4[0] ) ) +
- gf * EvtComplex( 0.0, 1.0 ) *
- dual( EvtGenFunctions::directProd( p4[0] + p4b,
- p4b - p4[0] ) ) -
- am * ( ( EvtGenFunctions::directProd( p4b, p4b ) -
- EvtGenFunctions::directProd( p4b, p4[0] ) ) );
- jb[0] = tds.cont1( p->getDaug( 0 )->epsParent( 0 ).conj() );
- jb[1] = tds.cont1( p->getDaug( 0 )->epsParent( 1 ).conj() );
- jb[2] = tds.cont1( p->getDaug( 0 )->epsParent( 2 ).conj() );
- break;
- // D_1
- case 6:
- q = p4b - p4[0];
- q2 = q * q;
- nbcurrent = 1;
- ffmodel.getscalarff( B0, D3P00, EvtPDL::getMeanMass( D3P00 ), q2,
- &fp, &fm );
- jb[0] = fp * ( p4b + p4[0] ) + fm * q;
- break;
- default:
- EvtGenReport( EVTGEN_ERROR, "EvtGen" )
- << "In EvtBHadronic, unknown hadronic current." << endl;
- }
-
- double norm;
-
- switch ( wcurrent ) {
- case 1:
- case 3:
- case 4:
- nwcurrent = 1;
- jw[0] = p4[getNDaug() - 1];
- break;
-
- case 2:
- case 5:
- case 6:
- nwcurrent = 3;
- norm = 1.0 /
- sqrt( p4[1].get( 0 ) * p4[1].get( 0 ) / p4[1].mass2() - 1.0 );
- jw[0] = norm * p->getDaug( getNDaug() - 1 )->epsParent( 0 );
- jw[1] = norm * p->getDaug( getNDaug() - 1 )->epsParent( 1 );
- jw[2] = norm * p->getDaug( getNDaug() - 1 )->epsParent( 2 );
- break;
-
- default:
- EvtGenReport( EVTGEN_ERROR, "EvtGen" )
- << "In EvtBHadronic, unknown W current." << endl;
- }
-
- if ( nbcurrent == 1 && nwcurrent == 1 ) {
- vertex( jb[0] * jw[0] );
- return;
- }
-
- if ( nbcurrent == 1 ) {
- for ( j = 0; j < nwcurrent; j++ ) {
- vertex( j, jb[0] * jw[j] );
- }
- return;
- }
-
- if ( nwcurrent == 1 ) {
- for ( j = 0; j < nbcurrent; j++ ) {
- vertex( j, jb[j] * jw[0] );
- }
- return;
- }
-
- for ( j = 0; j < nbcurrent; j++ ) {
- for ( i = 0; i < nwcurrent; i++ ) {
- vertex( j, i, jb[j] * jw[i] );
- }
- }
- return;
-}
diff --git a/src/EvtGenModels/EvtModelReg.cpp b/src/EvtGenModels/EvtModelReg.cpp
--- a/src/EvtGenModels/EvtModelReg.cpp
+++ b/src/EvtGenModels/EvtModelReg.cpp
@@ -26,7 +26,6 @@
#include "EvtGenBase/EvtPatches.hh"
#include "EvtGenModels/EvtBBScalar.hh"
-#include "EvtGenModels/EvtBHadronic.hh"
#include "EvtGenModels/EvtBLLNuL.hh"
#include "EvtGenModels/EvtBTo3piCP.hh"
#include "EvtGenModels/EvtBTo4piCP.hh"
@@ -201,7 +200,6 @@
modelist.registerModel( new EvtHQET );
modelist.registerModel( new EvtHQET2 );
modelist.registerModel( new EvtISGW );
- modelist.registerModel( new EvtBHadronic );
modelist.registerModel( new EvtVSS );
modelist.registerModel( new EvtVSSMix );
modelist.registerModel( new EvtVSSBMixCPT );
diff --git a/test/do_tests b/test/do_tests
--- a/test/do_tests
+++ b/test/do_tests
@@ -71,7 +71,6 @@
time ./evtgenlhc_test1 twobody 1000 exampleFiles/TWOBODY2.DEC twobody2.root
time ./evtgenlhc_test1 twobody 1000 exampleFiles/TWOBODY3.DEC twobody3.root
time ./evtgenlhc_test1 pipipi 1000
-time ./evtgenlhc_test1 bhadronic 1000
time ./evtgenlhc_test1 singleb 1000
time ./evtgenlhc_test1 pipi 1000
time ./evtgenlhc_test1 pipipipi 1000
diff --git a/test/evtgenlhc_test1.cc b/test/evtgenlhc_test1.cc
--- a/test/evtgenlhc_test1.cc
+++ b/test/evtgenlhc_test1.cc
@@ -101,7 +101,6 @@
void runD1( int nevent, EvtGen& myGenerator );
void runGenericCont( int nevent, EvtGen& myGenerator );
void runPiPiPi( int nevent, EvtGen& myGenerator );
-void runBHadronic( int nevent, EvtGen& myGenerator );
void runSingleB( int nevent, EvtGen& myGenerator );
void runA2Pi( int nevent, EvtGen& myGenerator );
void runAlias();
@@ -419,11 +418,6 @@
runPiPiPi( nevent, myGenerator );
}
- if ( !strcmp( argv[1], "bhadronic" ) ) {
- int nevent = atoi( argv[2] );
- runBHadronic( nevent, myGenerator );
- }
-
if ( !strcmp( argv[1], "singleb" ) ) {
int nevent = atoi( argv[2] );
runSingleB( nevent, myGenerator );
@@ -4951,48 +4945,6 @@
EvtGenReport( EVTGEN_INFO, "EvtGen" ) << "SUCCESS\n";
}
-void runBHadronic( int nevent, EvtGen& myGenerator )
-{
- std::ofstream outmix;
- outmix.open( "bhadronic.dat" );
-
- int count;
-
- myGenerator.readUDecay( "exampleFiles/BHADRONIC.DEC" );
-
- static EvtId B0 = EvtPDL::getId( std::string( "B0" ) );
-
- count = 1;
-
- do {
- EvtVector4R p_init( EvtPDL::getMass( B0 ), 0.0, 0.0, 0.0 );
-
- EvtParticle* root_part = EvtParticleFactory::particleFactory( B0, p_init );
- root_part->setDiagonalSpinDensity();
-
- myGenerator.generateDecay( root_part );
-
- EvtParticle* p;
-
- // root_part->printTree();
-
- p = root_part;
-
- do {
- outmix << p->getId().getId() << " " << p->getP4Lab().d3mag()
- << std::endl;
- p = p->nextIter();
-
- } while ( p != nullptr );
-
- root_part->deleteTree();
-
- } while ( count++ < nevent );
-
- outmix.close();
- EvtGenReport( EVTGEN_INFO, "EvtGen" ) << "SUCCESS\n";
-}
-
void runSingleB( int nevent, EvtGen& myGenerator )
{
int count;
diff --git a/test/exampleFiles/BHADRONIC.DEC b/test/exampleFiles/BHADRONIC.DEC
deleted file mode 100644
--- a/test/exampleFiles/BHADRONIC.DEC
+++ /dev/null
@@ -1,50 +0,0 @@
-Decay B0
-1.000 anti-D0 pi0 BHADRONIC 1 1;
-1.000 anti-D0 rho0 BHADRONIC 1 2;
-1.000 anti-D0 eta BHADRONIC 1 3;
-1.000 anti-D0 eta' BHADRONIC 1 4;
-1.000 anti-D0 omega BHADRONIC 1 5;
-1.000 anti-D0 a_10 BHADRONIC 1 6;
-1.000 anti-D*0 pi0 BHADRONIC 2 1;
-1.000 anti-D*0 rho0 BHADRONIC 2 2;
-1.000 anti-D*0 eta BHADRONIC 2 3;
-1.000 anti-D*0 eta' BHADRONIC 2 4;
-1.000 anti-D*0 omega BHADRONIC 2 5;
-1.000 anti-D*0 a_10 BHADRONIC 2 6;
-1.000 anti-D_10 pi0 BHADRONIC 3 1;
-1.000 anti-D_10 rho0 BHADRONIC 3 2;
-1.000 anti-D_10 eta BHADRONIC 3 3;
-1.000 anti-D_10 eta' BHADRONIC 3 4;
-1.000 anti-D_10 omega BHADRONIC 3 5;
-1.000 anti-D_10 a_10 BHADRONIC 3 6;
-1.000 anti-D_2*0 pi0 BHADRONIC 4 1;
-1.000 anti-D_2*0 rho0 BHADRONIC 4 2;
-1.000 anti-D_2*0 eta BHADRONIC 4 3;
-1.000 anti-D_2*0 eta' BHADRONIC 4 4;
-1.000 anti-D_2*0 omega BHADRONIC 4 5;
-1.000 anti-D_2*0 a_10 BHADRONIC 4 6;
-1.000 anti-D'_10 pi0 BHADRONIC 5 1;
-1.000 anti-D'_10 rho0 BHADRONIC 5 2;
-1.000 anti-D'_10 eta BHADRONIC 5 3;
-1.000 anti-D'_10 eta' BHADRONIC 5 4;
-1.000 anti-D'_10 omega BHADRONIC 5 5;
-1.000 anti-D'_10 a_10 BHADRONIC 5 6;
-1.000 anti-D_0*0 pi0 BHADRONIC 6 1;
-1.000 anti-D_0*0 rho0 BHADRONIC 6 2;
-1.000 anti-D_0*0 eta BHADRONIC 6 3;
-1.000 anti-D_0*0 eta' BHADRONIC 6 4;
-1.000 anti-D_0*0 omega BHADRONIC 6 5;
-1.000 anti-D_0*0 a_10 BHADRONIC 6 6;
-Enddecay
-End
-
-
-
-
-
-
-
-
-
-
-
diff --git a/test/exampleFiles/do_tests b/test/exampleFiles/do_tests
--- a/test/exampleFiles/do_tests
+++ b/test/exampleFiles/do_tests
@@ -79,7 +79,6 @@
time ../../testEvtGen twobody 1000 TWOBODY2.DEC
time ../../testEvtGen twobody 1000 TWOBODY3.DEC
time ../../testEvtGen pipipi 1000
-time ../../testEvtGen bhadronic 1000
time ../../testEvtGen singleb 1000
time ../../testEvtGen pipi 1000
time ../../testEvtGen pipipipi 1000
diff --git a/test/jsonFiles/BHADRONIC__Bd_D0pi0.json b/test/jsonFiles/BHADRONIC__Bd_D0pi0.json
deleted file mode 100644
--- a/test/jsonFiles/BHADRONIC__Bd_D0pi0.json
+++ /dev/null
@@ -1,23 +0,0 @@
-{
- "parent" : "B0",
- "daughters" : ["anti-D0", "pi0"],
- "models" : ["BHADRONIC", "", ""],
- "parameters" : [["1", "1"], [], []],
- "do_conjugate_decay" : [ true, false, false ],
- "extras" : ["noPhotos"],
- "events" : 10000,
- "histograms" : [
- {"variable" : "prob", "title" : "Prob", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 1.0},
- {"variable" : "parMass", "title" : "m(B0)", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 5.25, "xmax" : 5.3},
- {"variable" : "mass", "title" : "m(D0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.8, "xmax" : 1.9},
- {"variable" : "mass", "title" : "m(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 0.12, "xmax" : 0.15},
- {"variable" : "p", "title" : "p(D0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 2.0, "xmax" : 2.5},
- {"variable" : "phi", "title" : "phi(D0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(D0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
- {"variable" : "p", "title" : "p(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 2.0, "xmax" : 2.5},
- {"variable" : "phi", "title" : "phi(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0}
- ],
- "outfile" : "BHADRONIC__Bd_D0pi0.root",
- "reference" : "RefBHADRONIC__Bd_D0pi0.root"
-}
diff --git a/test/jsonFiles/BHADRONIC__Bd_D0rho0.json b/test/jsonFiles/BHADRONIC__Bd_D0rho0.json
deleted file mode 100644
--- a/test/jsonFiles/BHADRONIC__Bd_D0rho0.json
+++ /dev/null
@@ -1,23 +0,0 @@
-{
- "parent" : "B0",
- "daughters" : ["anti-D0", "rho0"],
- "models" : ["BHADRONIC", "", ""],
- "parameters" : [["1", "2"], [], []],
- "do_conjugate_decay" : [ true, false, false ],
- "extras" : ["noPhotos"],
- "events" : 10000,
- "histograms" : [
- {"variable" : "prob", "title" : "Prob", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 1.0},
- {"variable" : "parMass", "title" : "m(B0)", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 5.25, "xmax" : 5.3},
- {"variable" : "mass", "title" : "m(D0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.8, "xmax" : 1.9},
- {"variable" : "mass", "title" : "m(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 2.0},
- {"variable" : "p", "title" : "p(D0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.5, "xmax" : 2.5},
- {"variable" : "phi", "title" : "phi(D0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(D0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
- {"variable" : "p", "title" : "p(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 1.5, "xmax" : 2.5},
- {"variable" : "phi", "title" : "phi(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0}
- ],
- "outfile" : "BHADRONIC__Bd_D0rho0.root",
- "reference" : "RefBHADRONIC__Bd_D0rho0.root"
-}
diff --git a/test/jsonFiles/BHADRONIC__Bd_D2star0pi0.json b/test/jsonFiles/BHADRONIC__Bd_D2star0pi0.json
deleted file mode 100644
--- a/test/jsonFiles/BHADRONIC__Bd_D2star0pi0.json
+++ /dev/null
@@ -1,23 +0,0 @@
-{
- "parent" : "B0",
- "daughters" : ["anti-D_2*0", "pi0"],
- "models" : ["BHADRONIC", "", ""],
- "parameters" : [["4", "1"], [], []],
- "do_conjugate_decay" : [ true, false, false ],
- "extras" : ["noPhotos"],
- "events" : 10000,
- "histograms" : [
- {"variable" : "prob", "title" : "Prob", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 1.0},
- {"variable" : "parMass", "title" : "m(B0)", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 5.25, "xmax" : 5.3},
- {"variable" : "mass", "title" : "m(D_2*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 2.2, "xmax" : 3.0},
- {"variable" : "mass", "title" : "m(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 0.12, "xmax" : 0.15},
- {"variable" : "p", "title" : "p(D_2*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.6, "xmax" : 2.2},
- {"variable" : "phi", "title" : "phi(D_2*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(D_2*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
- {"variable" : "p", "title" : "p(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 1.6, "xmax" : 2.2},
- {"variable" : "phi", "title" : "phi(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0}
- ],
- "outfile" : "BHADRONIC__Bd_D2star0pi0.root",
- "reference" : "RefBHADRONIC__Bd_D2star0pi0.root"
-}
diff --git a/test/jsonFiles/BHADRONIC__Bd_D2star0rho0.json b/test/jsonFiles/BHADRONIC__Bd_D2star0rho0.json
deleted file mode 100644
--- a/test/jsonFiles/BHADRONIC__Bd_D2star0rho0.json
+++ /dev/null
@@ -1,23 +0,0 @@
-{
- "parent" : "B0",
- "daughters" : ["anti-D_2*0", "rho0"],
- "models" : ["BHADRONIC", "", ""],
- "parameters" : [["4", "2"], [], []],
- "do_conjugate_decay" : [ true, false, false ],
- "extras" : ["noPhotos"],
- "events" : 10000,
- "histograms" : [
- {"variable" : "prob", "title" : "Prob", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 1.0},
- {"variable" : "parMass", "title" : "m(B0)", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 5.25, "xmax" : 5.3},
- {"variable" : "mass", "title" : "m(D_2*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 2.2, "xmax" : 3.0},
- {"variable" : "mass", "title" : "m(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 2.0},
- {"variable" : "p", "title" : "p(D_2*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.2, "xmax" : 2.2},
- {"variable" : "phi", "title" : "phi(D_2*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(D_2*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
- {"variable" : "p", "title" : "p(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 1.2, "xmax" : 2.2},
- {"variable" : "phi", "title" : "phi(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0}
- ],
- "outfile" : "BHADRONIC__Bd_D2star0rho0.root",
- "reference" : "RefBHADRONIC__Bd_D2star0rho0.root"
-}
diff --git a/test/jsonFiles/BHADRONIC__Bd_Dstar0pi0.json b/test/jsonFiles/BHADRONIC__Bd_Dstar0pi0.json
deleted file mode 100644
--- a/test/jsonFiles/BHADRONIC__Bd_Dstar0pi0.json
+++ /dev/null
@@ -1,23 +0,0 @@
-{
- "parent" : "B0",
- "daughters" : ["anti-D*0", "pi0"],
- "models" : ["BHADRONIC", "", ""],
- "parameters" : [["2", "1"], [], []],
- "do_conjugate_decay" : [ true, false, false ],
- "extras" : ["noPhotos"],
- "events" : 10000,
- "histograms" : [
- {"variable" : "prob", "title" : "Prob", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 1.0},
- {"variable" : "parMass", "title" : "m(B0)", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 5.25, "xmax" : 5.3},
- {"variable" : "mass", "title" : "m(D*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.9, "xmax" : 2.1},
- {"variable" : "mass", "title" : "m(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 0.12, "xmax" : 0.15},
- {"variable" : "p", "title" : "p(D*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 2.0, "xmax" : 2.5},
- {"variable" : "phi", "title" : "phi(D*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(D*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
- {"variable" : "p", "title" : "p(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 2.0, "xmax" : 2.5},
- {"variable" : "phi", "title" : "phi(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(pi0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0}
- ],
- "outfile" : "BHADRONIC__Bd_Dstar0pi0.root",
- "reference" : "RefBHADRONIC__Bd_Dstar0pi0.root"
-}
diff --git a/test/jsonFiles/BHADRONIC__Bd_Dstar0rho0.json b/test/jsonFiles/BHADRONIC__Bd_Dstar0rho0.json
deleted file mode 100644
--- a/test/jsonFiles/BHADRONIC__Bd_Dstar0rho0.json
+++ /dev/null
@@ -1,23 +0,0 @@
-{
- "parent" : "B0",
- "daughters" : ["anti-D*0", "rho0"],
- "models" : ["BHADRONIC", "", ""],
- "parameters" : [["2", "2"], [], []],
- "do_conjugate_decay" : [ true, false, false ],
- "extras" : ["noPhotos"],
- "events" : 10000,
- "histograms" : [
- {"variable" : "prob", "title" : "Prob", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 1.0},
- {"variable" : "parMass", "title" : "m(B0)", "d1" : 0, "d2" : 0, "nbins" : 100, "xmin" : 5.25, "xmax" : 5.3},
- {"variable" : "mass", "title" : "m(D*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.9, "xmax" : 2.1},
- {"variable" : "mass", "title" : "m(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 0.0, "xmax" : 2.0},
- {"variable" : "p", "title" : "p(D*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : 1.5, "xmax" : 2.5},
- {"variable" : "phi", "title" : "phi(D*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(D*0bar)", "d1" : 1, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0},
- {"variable" : "p", "title" : "p(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : 1.5, "xmax" : 2.5},
- {"variable" : "phi", "title" : "phi(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -180.0, "xmax" : 180.0},
- {"variable" : "cosTheta", "title" : "cosTheta(rho0)", "d1" : 2, "d2" : 0, "nbins" : 100, "xmin" : -1.0, "xmax" : 1.0}
- ],
- "outfile" : "BHADRONIC__Bd_Dstar0rho0.root",
- "reference" : "RefBHADRONIC__Bd_Dstar0rho0.root"
-}