No OneTemporary
Actions

Size

46 KB

Subscribers

None

View Options

	diff --git a/src/Lau1DCubicSpline.cc b/src/Lau1DCubicSpline.cc
	index 00b55eb..92cf599 100644
	--- a/src/Lau1DCubicSpline.cc
	+++ b/src/Lau1DCubicSpline.cc
	@@ -1,320 +1,326 @@

	// Copyright University of Warwick 2004 - 2015.
	// Distributed under the Boost Software License, Version 1.0.
	// (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

	// Authors:
	// Thomas Latham
	// John Back
	// Paul Harrison

	/*! \file Lau1DCubicSpline.cc
	\brief File containing implementation of Lau1DCubicSpline class.
	*/

	#include <cmath>
	#include <cstdlib>
	#include <iostream>

	#include <TH2.h>
	#include <TSystem.h>

	#include "Lau1DCubicSpline.hh"

	Lau1DCubicSpline::Lau1DCubicSpline(const std::vector<Double_t>& xs, const std::vector<Double_t>& ys,
	LauSplineType type,
	LauSplineBoundaryType leftBound,
	LauSplineBoundaryType rightBound,
	Double_t dydx0, Double_t dydxn) :
	nKnots_(xs.size()),
	x_(xs),
	y_(ys),
	type_(type),
	leftBound_(leftBound),
	rightBound_(rightBound),
	dydx0_(dydx0),
	dydxn_(dydxn)
	{
	init();
	}

	Lau1DCubicSpline::~Lau1DCubicSpline()
	{
	}

	Double_t Lau1DCubicSpline::evaluate(Double_t x) const
	{
	// do not attempt to extrapolate the spline
	if( x<x_[0] \|\| x>x_[nKnots_-1] ) {
	std::cout << "WARNING in Lau1DCubicSpline::evaluate : function is only defined between " << x_[0] << " and " << x_[nKnots_-1] << std::endl;
	std::cout << " value at " << x << " returned as 0" << std::endl;
	return 0.;
	}

	// first determine which 'cell' of the spline x is in
	// cell i runs from knot i to knot i+1
	Int_t cell(0);

	while( x > x_[cell+1] ) {
	++cell;
	}

	// obtain x- and y-values of the neighbouring knots
	Double_t xLow = x_[cell];
	Double_t xHigh = x_[cell+1];
	Double_t yLow = y_[cell];
	Double_t yHigh = y_[cell+1];

	if(type_ == Lau1DCubicSpline::LinearInterpolation) {
	return yHigh(x-xLow)/(xHigh-xLow) + yLow(xHigh-x)/(xHigh-xLow);
	}

	// obtain t, the normalised x-coordinate within the cell,
	// and the coefficients a and b, which are defined in cell i as:
	//
	// a_i = k_i *(x_i+1 - x_i) - (y_i+1 - y_i),
	// b_i = -k_i+1*(x_i+1 - x_i) + (y_i+1 - y_i)
	//
	// where k_i is (by construction) the first derivative at knot i
	Double_t t = (x - xLow) / (xHigh - xLow);
	Double_t a = dydx_[cell] * (xHigh - xLow) - (yHigh - yLow);
	Double_t b = -1.dydx_[cell+1] (xHigh - xLow) + (yHigh - yLow);

	Double_t retVal = (1 - t) * yLow + t * yHigh + t * (1 - t) * ( a * (1 - t) + b * t );

	return retVal;
	}

	void Lau1DCubicSpline::updateYValues(const std::vector<Double_t>& ys)
	{
	y_ = ys;
	- calcDerivatives();
	+ this->calcDerivatives();
	}

	void Lau1DCubicSpline::updateType(LauSplineType type)
	{
	if(type_ != type) {
	type_ = type;
	- calcDerivatives();
	+ this->calcDerivatives();
	}
	}

	void Lau1DCubicSpline::updateBoundaryConditions(LauSplineBoundaryType leftBound,
	LauSplineBoundaryType rightBound,
	Double_t dydx0, Double_t dydxn)
	{
	Bool_t updateDerivatives(kFALSE);

	if(leftBound_ != leftBound \|\| rightBound_ != rightBound) {
	leftBound_ = leftBound;
	rightBound_ = rightBound;
	updateDerivatives = kTRUE;
	}

	if(dydx0_ != dydx0) {
	dydx0_ = dydx0;
	if(leftBound_ == Lau1DCubicSpline::Clamped) updateDerivatives = kTRUE;
	}

	if(dydxn_ != dydxn) {
	dydxn_ = dydxn;
	if(rightBound_ == Lau1DCubicSpline::Clamped) updateDerivatives = kTRUE;
	}

	- if(updateDerivatives) calcDerivatives();
	+ if(updateDerivatives) {
	+ this->calcDerivatives();
	+ }
	}

	void Lau1DCubicSpline::init()
	{
	if( y_.size() != x_.size()) {
	std::cout << "ERROR in Lau1DCubicSpline::init : The number of y-values given does not match the number of x-values" << std::endl;
	std::cout << " Found " << y_.size() << ", expected " << x_.size() << std::endl;
	gSystem->Exit(EXIT_FAILURE);
	}

	dydx_.insert(dydx_.begin(),nKnots_,0.);

	a_.insert(a_.begin(),nKnots_,0.);
	b_.insert(b_.begin(),nKnots_,0.);
	c_.insert(c_.begin(),nKnots_,0.);
	d_.insert(d_.begin(),nKnots_,0.);

	this->calcDerivatives();
	}

	void Lau1DCubicSpline::calcDerivatives()
	{
	- if(type_ == Lau1DCubicSpline::LinearInterpolation) {
	- return; //derivatives not needed for linear interpolation
	- } else if(type_ == Lau1DCubicSpline::AkimaSpline) {
	- this->calcDerivativesAkima();
	- } else {
	- this->calcDerivativesStandard();
	+ switch ( type_ ) {
	+ case Lau1DCubicSpline::StandardSpline :
	+ this->calcDerivativesStandard();
	+ break;
	+ case Lau1DCubicSpline::AkimaSpline :
	+ this->calcDerivativesAkima();
	+ break;
	+ case Lau1DCubicSpline::LinearInterpolation :
	+ //derivatives not needed for linear interpolation
	+ break;
	}
	}

	void Lau1DCubicSpline::calcDerivativesStandard()
	{
	// derivatives are determined such that the second derivative is continuous at internal knots

	// derivatives, k_i, are the solutions to a set of linear equations of the form:
	// a_i * k_i-1 + b_i * k+i + c_i * k_i+1 = d_i with a_0 = 0, c_n-1 = 0

	// this is solved using the tridiagonal matrix algorithm as on en.wikipedia.org/wiki/Tridiagonal_matrix_algorithm

	// first and last equations give boundary conditions
	// - for natural boundary, require f''(x) = 0 at end knot
	// - for 'not a knot' boundary, require f'''(x) continuous at second knot
	// - for clamped boundary, require predefined value of f'(x) at end knot

	// non-zero values of a_0 and c_n-1 would give cyclic boundary conditions
	a_[0] = 0.;
	c_[nKnots_-1] = 0.;

	// set left boundary condition
	if(leftBound_ == Lau1DCubicSpline::Natural) {
	b_[0] = 2./(x_[1]-x_[0]);
	c_[0] = 1./(x_[1]-x_[0]);
	d_[0] = 3.(y_[1]-y_[0])/((x_[1]-x_[0])(x_[1]-x_[0]));

	} else if(leftBound_ == Lau1DCubicSpline::NotAKnot) {
	// define the width, h, and the 'slope', delta, of the first cell
	Double_t h1(x_[1]-x_[0]), h2(x_[2]-x_[0]);
	Double_t delta1((y_[1]-y_[0])/h1), delta2((y_[2]-y_[1])/h2);

	// these coefficients can be determined by requiring f'''_0(x_1) = f'''_1(x_1)
	// the requirement f''_0(x_1) = f''_1(x_1) has been used to remove the dependence on k_2
	b_[0] = h2;
	c_[0] = h1+h2;
	d_[0] = delta1(2.h2h2 + 3.h1h2)/(h1+h2) + delta25.h1h1/(h1+h2);

	} else { //Clamped
	b_[0] = 1.;
	c_[0] = 0.;
	d_[0] = dydx0_;
	}

	// set right boundary condition
	if(rightBound_ == Lau1DCubicSpline::Natural) {
	a_[nKnots_-1] = 1./(x_[nKnots_-1]-x_[nKnots_-2]);
	b_[nKnots_-1] = 2./(x_[nKnots_-1]-x_[nKnots_-2]);
	d_[nKnots_-1] = 3.(y_[nKnots_-1]-y_[nKnots_-2])/((x_[nKnots_-1]-x_[nKnots_-2])(x_[nKnots_-1]-x_[nKnots_-2]));

	} else if(rightBound_ == Lau1DCubicSpline::NotAKnot) {
	// define the width, h, and the 'slope', delta, of the last cell
	Double_t hnm1(x_[nKnots_-1]-x_[nKnots_-2]), hnm2(x_[nKnots_-2]-x_[nKnots_-3]);
	Double_t deltanm1((y_[nKnots_-1]-y_[nKnots_-2])/hnm1), deltanm2((y_[nKnots_-2]-y_[nKnots_-3])/hnm2);

	// these coefficients can be determined by requiring f'''_n-3(x_n-2) = f'''_n-2(x_n-2)
	// the requirement f''_n-3(x_n-2) = f''_n-2(x_n-2) has been used to remove
	// the dependence on k_n-3
	a_[nKnots_-1] = hnm2 + hnm1;
	b_[nKnots_-1] = hnm1;
	d_[nKnots_-1] = deltanm2hnm1hnm1/(hnm2+hnm1) + deltanm1(2.hnm2hnm2 + 3.hnm2*hnm1)/(hnm2+hnm1);

	} else { //Clamped
	a_[nKnots_-1] = 0.;
	b_[nKnots_-1] = 1.;
	d_[nKnots_-1] = dydxn_;
	}

	// the remaining equations ensure that f_i-1''(x_i) = f''_i(x_i) for all internal knots
	for(UInt_t i=1; i<nKnots_-1; ++i) {
	a_[i] = 1./(x_[i]-x_[i-1]);
	b_[i] = 2./(x_[i]-x_[i-1]) + 2./(x_[i+1]-x_[i]);
	c_[i] = 1./(x_[i+1]-x_[i]);
	d_[i] = 3.(y_[i]-y_[i-1])/((x_[i]-x_[i-1])(x_[i]-x_[i-1])) + 3.(y_[i+1]-y_[i])/((x_[i+1]-x_[i])(x_[i+1]-x_[i]));
	}

	// forward sweep - replace c_i and d_i with
	//
	// c'_i = c_i / b_i for i = 0
	// c_i / (b_i - a_i * c_i-1) otherwise
	//
	// and
	//
	// d'_i = d_i / b_i for i = 0
	// (d_i - a_i * d_i-1) / (b_i - a_i * c_i-1) otherwise

	c_[0] /= b_[0];
	d_[0] /= b_[0];

	for(UInt_t i=1; i<nKnots_-1; ++i) {
	c_[i] = c_[i] / (b_[i] - a_[i]*c_[i-1]);
	d_[i] = (d_[i] - a_[i]d_[i-1]) / (b_[i] - a_[i]c_[i-1]);
	}

	d_[nKnots_-1] = (d_[nKnots_-1] - a_[nKnots_-1]d_[nKnots_-2]) / (b_[nKnots_-1] - a_[nKnots_-1]c_[nKnots_-2]);

	// back substitution - calculate k_i = dy/dx at each knot
	//
	// k_i = d'_i for i = n-1
	// d'_i - c'_i * k_i+1 otherwise

	dydx_[nKnots_-1] = d_[nKnots_-1];

	for(Int_t i=nKnots_-2; i>=0; --i) {
	dydx_[i] = d_[i] - c_[i]*dydx_[i+1];
	}
	}

	void Lau1DCubicSpline::calcDerivativesAkima() {
	//derivatives are calculated according to the Akima method
	// J.ACM vol. 17 no. 4 pp 589-602

	Double_t am1(0.), an(0.), anp1(0.);

	// a[i] is the slope of the segment from point i-1 to point i
	//
	// n.b. segment 0 is before point 0 and segment n is after point n-1
	// internal segments are numbered 1 - n-1
	for(UInt_t i=1; i<nKnots_; ++i) {
	a_[i] = (y_[i]-y_[i-1])/(x_[i]-x_[i-1]);
	}

	// calculate slopes between additional points on each end according to the Akima method
	// method assumes that the additional points follow a quadratic defined by the last three points
	// this leads to the relations a[2] - a[1] = a[1] - a[0] = a[0] - a[-1] and a[n-1] - a[n-2] = a[n] - a[n-1] = a[n+1] - a[n]
	a_[0] = 2*a_[1] - a_[2];
	am1 = 2*a_[0] - a_[1];

	an = 2*a_[nKnots_-1] - a_[nKnots_-2];
	anp1 = 2*an - a_[nKnots_-1];

	// b[i] is the weight of a[i] towards dydx[i]
	// c[i] is the weight of a[i+1] towards dydx[i]
	// See Appendix A of J.ACM vol. 17 no. 4 pp 589-602 for a justification of these weights
	b_[0] = TMath::Abs(a_[1] - a_[0]);
	b_[1] = TMath::Abs(a_[2] - a_[1]);

	c_[0] = TMath::Abs(a_[0] - am1 );
	c_[1] = TMath::Abs(a_[1] - a_[0]);

	for(UInt_t i=2; i<nKnots_-2; ++i) {
	b_[i] = TMath::Abs(a_[i+2] - a_[i+1]);
	c_[i] = TMath::Abs(a_[i] - a_[i-1]);
	}

	b_[nKnots_-2] = TMath::Abs(an - a_[nKnots_-1]);
	b_[nKnots_-1] = TMath::Abs(anp1 - an );

	c_[nKnots_-2] = TMath::Abs(a_[nKnots_-2] - a_[nKnots_-3]);
	c_[nKnots_-1] = TMath::Abs(a_[nKnots_-1] - a_[nKnots_-2]);

	// dy/dx calculated as the weighted average of a[i] and a[i+1]:
	// dy/dx_i = ( \| a_i+2 - a_i+1 \| a_i + \| a_i - a_i-1 \| a_i+1 ) / ( \| a_i+2 - a_i+1 \| + \| a_i - a_i-1 \| )
	// in the special case a_i-1 == a_i != a_i+1 == a_i+2 this function is undefined so dy/dx is then defined as (a_i + a_i+1) / 2
	for(UInt_t i=0; i<nKnots_-2; ++i) {
	if(b_[i]==0 && c_[i]==0) {
	dydx_[i] = ( a_[i] + a_[i+1] ) / 2.;
	} else {
	dydx_[i] = ( b_[i] * a_[i] + c_[i] * a_[i+1] ) / ( b_[i] + c_[i] );
	}
	}

	if(b_[nKnots_-1]==0 && c_[nKnots_-1]==0) {
	dydx_[nKnots_-1] = ( a_[nKnots_-1] + an ) / 2.;
	} else {
	dydx_[nKnots_-1] = ( b_[nKnots_-1] * a_[nKnots_-1] + c_[nKnots_-1] * an ) / ( b_[nKnots_-1] + c_[nKnots_-1] );
	}
	}
	diff --git a/src/LauResonanceMaker.cc b/src/LauResonanceMaker.cc
	index 0cf89c9..0262d5d 100644
	--- a/src/LauResonanceMaker.cc
	+++ b/src/LauResonanceMaker.cc
	@@ -1,723 +1,723 @@

	// Copyright University of Warwick 2004 - 2014.
	// Distributed under the Boost Software License, Version 1.0.
	// (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

	// Authors:
	// Thomas Latham
	// John Back
	// Paul Harrison

	/*! \file LauResonanceMaker.cc
	\brief File containing implementation of LauResonanceMaker class.
	*/

	#include <iostream>

	#include "LauAbsResonance.hh"
	#include "LauBelleNR.hh"
	#include "LauBelleSymNR.hh"
	#include "LauBreitWignerRes.hh"
	#include "LauDabbaRes.hh"
	#include "LauDaughters.hh"
	#include "LauEFKLLMRes.hh"
	#include "LauFlatteRes.hh"
	#include "LauFlatNR.hh"
	#include "LauGaussIncohRes.hh"
	#include "LauGounarisSakuraiRes.hh"
	#include "LauKappaRes.hh"
	#include "LauLASSRes.hh"
	#include "LauLASSBWRes.hh"
	#include "LauLASSNRRes.hh"
	#include "LauModIndPartWaveMagPhase.hh"
	#include "LauModIndPartWaveRealImag.hh"
	#include "LauModIndPartWaveSymMagPhase.hh"
	#include "LauModIndPartWaveSymRealImag.hh"
	#include "LauNRAmplitude.hh"
	#include "LauPolNR.hh"
	#include "LauRelBreitWignerRes.hh"
	#include "LauResonanceInfo.hh"
	#include "LauResonanceMaker.hh"
	#include "LauSigmaRes.hh"

	ClassImp(LauResonanceMaker);


	LauResonanceMaker* LauResonanceMaker::resonanceMaker_ = 0;


	LauResonanceMaker::LauResonanceMaker() :
	nResDefMax_(0)
	{
	this->createResonanceVector();
	this->setDefaultBWRadius( LauBlattWeisskopfFactor::Parent, 4.0 );
	}

	LauResonanceMaker::~LauResonanceMaker()
	{
	for ( std::vector<LauBlattWeisskopfFactor*>::iterator iter = bwIndepFactors_.begin(); iter != bwIndepFactors_.end(); ++iter ) {
	delete *iter;
	}
	bwIndepFactors_.clear();
	for ( BWFactorCategoryMap::iterator iter = bwFactors_.begin(); iter != bwFactors_.end(); ++iter ) {
	delete iter->second;
	}
	bwFactors_.clear();
	}

	LauResonanceMaker& LauResonanceMaker::get()
	{
	if ( resonanceMaker_ == 0 ) {
	resonanceMaker_ = new LauResonanceMaker();
	}

	return *resonanceMaker_;
	}

	void LauResonanceMaker::createResonanceVector()
	{
	// Function to create all possible resonances that this class supports.
	// Also add in the sigma and kappa - but a special paramterisation is used
	// instead of the PDG "pole mass and width" values.

	std::cout << "INFO in LauResonanceMaker::createResonanceVector : Setting up possible resonance states..." << std::endl;

	LauResonanceInfo* neutral(0);
	LauResonanceInfo* positve(0);
	LauResonanceInfo* negatve(0);

	// Define the resonance names and store them in the array
	resInfo_.clear();
	resInfo_.reserve(100);
	// rho resonances name, mass, width, spin, charge, default BW category, BW radius parameter (defaults to 4.0)
	// rho(770)
	neutral = new LauResonanceInfo("rho0(770)", 0.77526, 0.1478, 1, 0, LauBlattWeisskopfFactor::Light, 5.3);
	positve = new LauResonanceInfo("rho+(770)", 0.77511, 0.1491, 1, 1, LauBlattWeisskopfFactor::Light, 5.3);
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// rho(1450)
	neutral = new LauResonanceInfo("rho0(1450)", 1.465, 0.400, 1, 0, LauBlattWeisskopfFactor::Light );
	positve = new LauResonanceInfo("rho+(1450)", 1.465, 0.400, 1, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// rho(1700)
	neutral = new LauResonanceInfo("rho0(1700)", 1.720, 0.250, 1, 0, LauBlattWeisskopfFactor::Light );
	positve = new LauResonanceInfo("rho+(1700)", 1.720, 0.250, 1, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );

	// K* resonances name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// K*(892)
	neutral = new LauResonanceInfo("K*0(892)", 0.89581, 0.0474, 1, 0, LauBlattWeisskopfFactor::Kstar, 3.0);
	positve = new LauResonanceInfo("K*+(892)", 0.89166, 0.0508, 1, 1, LauBlattWeisskopfFactor::Kstar, 3.0);
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// K*(1410)
	neutral = new LauResonanceInfo("K*0(1410)", 1.414, 0.232, 1, 0, LauBlattWeisskopfFactor::Kstar );
	positve = new LauResonanceInfo("K*+(1410)", 1.414, 0.232, 1, 1, LauBlattWeisskopfFactor::Kstar );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// K*_0(1430)
	neutral = new LauResonanceInfo("K*0_0(1430)", 1.425, 0.270, 0, 0, LauBlattWeisskopfFactor::Kstar );
	positve = new LauResonanceInfo("K*+_0(1430)", 1.425, 0.270, 0, 1, LauBlattWeisskopfFactor::Kstar );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// LASS nonresonant model
	neutral = neutral->createSharedParameterRecord("LASSNR0");
	positve = positve->createSharedParameterRecord("LASSNR+");
	negatve = negatve->createSharedParameterRecord("LASSNR-");
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// K*_2(1430)
	neutral = new LauResonanceInfo("K*0_2(1430)", 1.4324, 0.109, 2, 0, LauBlattWeisskopfFactor::Kstar );
	positve = new LauResonanceInfo("K*+_2(1430)", 1.4256, 0.0985, 2, 1, LauBlattWeisskopfFactor::Kstar );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// K*(1680)
	neutral = new LauResonanceInfo("K*0(1680)", 1.717, 0.322, 1, 0, LauBlattWeisskopfFactor::Kstar );
	positve = new LauResonanceInfo("K*+(1680)", 1.717, 0.322, 1, 1, LauBlattWeisskopfFactor::Kstar );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// K*(1950)
	neutral = new LauResonanceInfo("K*0_0(1950)", 1.945, 0.201, 0, 0, LauBlattWeisskopfFactor::Kstar );
	positve = new LauResonanceInfo("K*+_0(1950)", 1.945, 0.201, 0, 1, LauBlattWeisskopfFactor::Kstar );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );

	// phi resonances name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// phi(1020)
	neutral = new LauResonanceInfo("phi(1020)", 1.019461, 0.004266, 1, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// phi(1680)
	neutral = new LauResonanceInfo("phi(1680)", 1.680, 0.150, 1, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );

	// f resonances name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// f_0(980)
	neutral = new LauResonanceInfo("f_0(980)", 0.990, 0.070, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// f_2(1270)
	neutral = new LauResonanceInfo("f_2(1270)", 1.2751, 0.1851, 2, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// f_0(1370)
	neutral = new LauResonanceInfo("f_0(1370)", 1.370, 0.350, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// f'_0(1300), from Belle's Kspipi paper
	neutral = new LauResonanceInfo("f'_0(1300)", 1.449, 0.126, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// f_0(1500)
	neutral = new LauResonanceInfo("f_0(1500)", 1.505, 0.109, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// f'_2(1525)
	neutral = new LauResonanceInfo("f'_2(1525)", 1.525, 0.073, 2, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// f_0(1710)
	neutral = new LauResonanceInfo("f_0(1710)", 1.722, 0.135, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// f_2(2010)
	neutral = new LauResonanceInfo("f_2(2010)", 2.011, 0.202, 2, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );

	// omega resonances name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// omega(782)
	neutral = new LauResonanceInfo("omega(782)", 0.78265, 0.00849, 1, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );

	// a resonances name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// a_0(980)
	neutral = new LauResonanceInfo("a0_0(980)", 0.980, 0.092, 0, 0, LauBlattWeisskopfFactor::Light );
	positve = new LauResonanceInfo("a+_0(980)", 0.980, 0.092, 0, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// a_0(1450)
	neutral = new LauResonanceInfo("a0_0(1450)", 1.474, 0.265, 0, 0, LauBlattWeisskopfFactor::Light );
	positve = new LauResonanceInfo("a+_0(1450)", 1.474, 0.265, 0, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// a_2(1320)
	neutral = new LauResonanceInfo("a0_2(1320)", 1.3190, 0.1050, 2, 0, LauBlattWeisskopfFactor::Light );
	positve = new LauResonanceInfo("a+_2(1320)", 1.3190, 0.1050, 2, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );

	// charmonium resonances name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// chi_c0
	neutral = new LauResonanceInfo("chi_c0", 3.41475, 0.0105, 0, 0, LauBlattWeisskopfFactor::Charmonium );
	resInfo_.push_back( neutral );
	// chi_c1
	neutral = new LauResonanceInfo("chi_c1", 3.51066, 0.00084, 0, 0, LauBlattWeisskopfFactor::Charmonium );
	resInfo_.push_back( neutral );
	// chi_c2
	neutral = new LauResonanceInfo("chi_c2", 3.55620, 0.00193, 2, 0, LauBlattWeisskopfFactor::Charmonium );
	resInfo_.push_back( neutral );
	// X(3872)
	neutral = new LauResonanceInfo("X(3872)", 3.87169, 0.0012, 1, 0, LauBlattWeisskopfFactor::Charmonium );
	resInfo_.push_back( neutral );

	// unknown scalars name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// sigma
	neutral = new LauResonanceInfo("sigma0", 0.475, 0.550, 0, 0, LauBlattWeisskopfFactor::Light );
	positve = new LauResonanceInfo("sigma+", 0.475, 0.550, 0, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// kappa
	neutral = new LauResonanceInfo("kappa0", 0.682, 0.547, 0, 0, LauBlattWeisskopfFactor::Kstar );
	positve = new LauResonanceInfo("kappa+", 0.682, 0.547, 0, 1, LauBlattWeisskopfFactor::Kstar );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// dabba
	neutral = new LauResonanceInfo("dabba0", 2.098, 0.520, 0, 0, LauBlattWeisskopfFactor::Charm );
	positve = new LauResonanceInfo("dabba+", 2.098, 0.520, 0, 1, LauBlattWeisskopfFactor::Charm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );

	// excited charm states name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// D*
	neutral = new LauResonanceInfo("D*0", 2.00696, 0.0021, 1, 0, LauBlattWeisskopfFactor::Charm );
	positve = new LauResonanceInfo("D*+", 2.01026, 83.4e-6, 1, 1, LauBlattWeisskopfFactor::Charm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// D*_0
	neutral = new LauResonanceInfo("D*0_0", 2.318, 0.267, 0, 0, LauBlattWeisskopfFactor::Charm );
	positve = new LauResonanceInfo("D*+_0", 2.403, 0.283, 0, 1, LauBlattWeisskopfFactor::Charm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// D*_2
	//AVERAGE--neutral = new LauResonanceInfo("D*0_2", 2.4618, 0.049, 2, 0 );
	neutral = new LauResonanceInfo("D*0_2", 2.4626, 0.049, 2, 0, LauBlattWeisskopfFactor::Charm );
	positve = new LauResonanceInfo("D*+_2", 2.4643, 0.037, 2, 1, LauBlattWeisskopfFactor::Charm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// D1(2420)
	neutral = new LauResonanceInfo("D0_1(2420)", 2.4214, 0.0274, 1, 0, LauBlattWeisskopfFactor::Charm );
	positve = new LauResonanceInfo("D+_1(2420)", 2.4232, 0.025, 1, 1, LauBlattWeisskopfFactor::Charm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// D(2600)
	//OLD--neutral = new LauResonanceInfo("D0(2600)", 2.6087, 0.093, 0, 0 );
	//OLD--positve = new LauResonanceInfo("D+(2600)", 2.6213, 0.093, 0, 1 );
	neutral = new LauResonanceInfo("D0(2600)", 2.612, 0.093, 0, 0, LauBlattWeisskopfFactor::Charm );
	positve = new LauResonanceInfo("D+(2600)", 2.612, 0.093, 0, 1, LauBlattWeisskopfFactor::Charm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// D(2760)
	//OLD-- neutral = new LauResonanceInfo("D0(2760)", 2.7633, 0.061, 1, 0 );
	//OLD-- positve = new LauResonanceInfo("D+(2760)", 2.7697, 0.061, 1, 1 );
	neutral = new LauResonanceInfo("D0(2760)", 2.761, 0.063, 1, 0, LauBlattWeisskopfFactor::Charm );
	positve = new LauResonanceInfo("D+(2760)", 2.761, 0.063, 1, 1, LauBlattWeisskopfFactor::Charm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// D(2900)
	neutral = new LauResonanceInfo("D0(3000)", 3.0, 0.15, 0, 0, LauBlattWeisskopfFactor::Charm );
	resInfo_.push_back( neutral );
	// D(3400)
	neutral = new LauResonanceInfo("D0(3400)", 3.4, 0.15, 0, 0, LauBlattWeisskopfFactor::Charm );
	resInfo_.push_back( neutral );

	// excited strange charm name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// Ds*
	positve = new LauResonanceInfo("Ds*+", 2.1121, 0.0019, 1, 1, LauBlattWeisskopfFactor::StrangeCharm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// Ds0*(2317)
	positve = new LauResonanceInfo("Ds*+_0(2317)", 2.3177, 0.0038, 0, 1, LauBlattWeisskopfFactor::StrangeCharm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// Ds2*(2573)
	positve = new LauResonanceInfo("Ds*+_2(2573)", 2.5719, 0.017, 2, 1, LauBlattWeisskopfFactor::StrangeCharm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// Ds1*(2700)
	positve = new LauResonanceInfo("Ds*+_1(2700)", 2.709, 0.117, 1, 1, LauBlattWeisskopfFactor::StrangeCharm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// Ds1*(2860)
	positve = new LauResonanceInfo("Ds*+_1(2860)", 2.862, 0.180, 1, 1, LauBlattWeisskopfFactor::StrangeCharm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// Ds3*(2860)
	positve = new LauResonanceInfo("Ds*+_3(2860)", 2.862, 0.058, 3, 1, LauBlattWeisskopfFactor::StrangeCharm );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );

	// excited bottom states name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// B*
	neutral = new LauResonanceInfo("B*0", 5.3252, 0.00, 1, 0, LauBlattWeisskopfFactor::Beauty, 6.0);
	positve = new LauResonanceInfo("B*+", 5.3252, 0.00, 1, 1, LauBlattWeisskopfFactor::Beauty, 6.0);
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( neutral );
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );

	// excited strange bottom name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// Bs*
	neutral = new LauResonanceInfo("Bs*0", 5.4154, 0.00, 1, 0, LauBlattWeisskopfFactor::StrangeBeauty, 6.0);
	resInfo_.push_back( neutral );

	// nonresonant models name, mass, width, spin, charge, BW category, BW radius parameter (defaults to 4.0)
	// Phase-space nonresonant model
	neutral = new LauResonanceInfo("NonReson", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// Theory-based nonresonant model
	neutral = new LauResonanceInfo("NRModel", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// Belle nonresonant models
	neutral = new LauResonanceInfo("BelleSymNR", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	neutral = new LauResonanceInfo("BelleNR", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	positve = new LauResonanceInfo("BelleNR+", 0.0, 0.0, 0, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	neutral = new LauResonanceInfo("BelleNR_Swave", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	positve = new LauResonanceInfo("BelleNR_Swave+",0.0, 0.0, 0, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	neutral = new LauResonanceInfo("BelleNR_Pwave", 0.0, 0.0, 1, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	positve = new LauResonanceInfo("BelleNR_Pwave+",0.0, 0.0, 1, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	neutral = new LauResonanceInfo("BelleNR_Dwave", 0.0, 0.0, 2, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	positve = new LauResonanceInfo("BelleNR_Dwave+",0.0, 0.0, 2, 1, LauBlattWeisskopfFactor::Light );
	negatve = positve->createChargeConjugate();
	resInfo_.push_back( positve );
	resInfo_.push_back( negatve );
	// Taylor expansion nonresonant model
	neutral = new LauResonanceInfo("NRTaylor", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	// Polynomial nonresonant models
	neutral = new LauResonanceInfo("PolNR_S0", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	neutral = new LauResonanceInfo("PolNR_S1", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	neutral = new LauResonanceInfo("PolNR_S2", 0.0, 0.0, 0, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	neutral = new LauResonanceInfo("PolNR_P0", 0.0, 0.0, 1, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	neutral = new LauResonanceInfo("PolNR_P1", 0.0, 0.0, 1, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );
	neutral = new LauResonanceInfo("PolNR_P2", 0.0, 0.0, 1, 0, LauBlattWeisskopfFactor::Light );
	resInfo_.push_back( neutral );

	nResDefMax_ = resInfo_.size();
	}

	void LauResonanceMaker::setDefaultBWRadius(const LauBlattWeisskopfFactor::BlattWeisskopfCategory bwCategory, const Double_t bwRadius)
	{
	if ( bwCategory == LauBlattWeisskopfFactor::Default \|\| bwCategory == LauBlattWeisskopfFactor::Indep ) {
	std::cerr << "WARNING in LauResonanceMaker::setDefaultBWRadius : cannot set radius values for Default or Indep categories" << std::endl;
	return;
	}

	// Check if a LauBlattWeisskopfFactor object has been created for this category
	// If it has then we can set it directly
	// If not then we just store the value to be used when it is created later
	BWFactorCategoryMap::iterator factor_iter = bwFactors_.find( bwCategory );
	if ( factor_iter != bwFactors_.end() ) {
	LauBlattWeisskopfFactor* bwFactor = factor_iter->second;
	LauParameter* radius = bwFactor->getRadiusParameter();
	radius->value(bwRadius);
	radius->initValue(bwRadius);
	radius->genValue(bwRadius);
	}

	bwDefaultRadii_[bwCategory] = bwRadius;
	}

	void LauResonanceMaker::fixBWRadius(const LauBlattWeisskopfFactor::BlattWeisskopfCategory bwCategory, const Bool_t fixRadius)
	{
	if ( bwCategory == LauBlattWeisskopfFactor::Default \|\| bwCategory == LauBlattWeisskopfFactor::Indep ) {
	std::cerr << "WARNING in LauResonanceMaker::fixBWRadius : cannot fix/float radius values for Default or Indep categories" << std::endl;
	return;
	}

	// Check if a LauBlattWeisskopfFactor object has been created for this category
	// If it has then we can set it directly
	// If not then we just store the value to be used when it is created later
	BWFactorCategoryMap::iterator factor_iter = bwFactors_.find( bwCategory );
	if ( factor_iter != bwFactors_.end() ) {
	LauBlattWeisskopfFactor* bwFactor = factor_iter->second;
	LauParameter* radius = bwFactor->getRadiusParameter();
	radius->fixed(fixRadius);
	}

	bwFixRadii_[bwCategory] = fixRadius;
	}

	LauBlattWeisskopfFactor* LauResonanceMaker::getBWFactor( const LauBlattWeisskopfFactor::BlattWeisskopfCategory bwCategory, const LauResonanceInfo* resInfo, const LauBlattWeisskopfFactor::BarrierType bwType )
	{
	LauBlattWeisskopfFactor* bwFactor(0);

	BWFactorCategoryMap::iterator factor_iter = bwFactors_.find( bwCategory );
	BWRadiusCategoryMap::const_iterator radius_iter = bwDefaultRadii_.find( bwCategory );

	if ( bwCategory == LauBlattWeisskopfFactor::Indep ) {
	bwFactor = new LauBlattWeisskopfFactor( *resInfo, bwType, bwCategory );
	bwIndepFactors_.push_back(bwFactor);
	} else if ( factor_iter == bwFactors_.end() ) {
	if ( radius_iter == bwDefaultRadii_.end() ) {
	bwFactor = new LauBlattWeisskopfFactor( *resInfo, bwType, bwCategory );
	} else {
	bwFactor = new LauBlattWeisskopfFactor( *resInfo, radius_iter->second, bwType, bwCategory );
	}
	bwFactors_[bwCategory] = bwFactor;
	} else {
	const UInt_t resSpin = resInfo->getSpin();
	bwFactor = factor_iter->second->createClone( resSpin );
	if ( bwFactor->getBarrierType() != bwType ) {
	std::cerr << "WARNING in LauResonanceMaker::getBWFactor : A barrier factor already exists for the specified category but does not have the specified barrier type.\n";
	std::cerr << " : If you want to use that type you will need to use a different category for this resonance." << std::endl;
	}
	}

	BWRadiusFixedCategoryMap::const_iterator radius_fixed_iter = bwFixRadii_.find( bwCategory );

	if ( radius_fixed_iter != bwFixRadii_.end() ) {
	const Bool_t fixed = radius_fixed_iter->second;

	LauParameter* radius = bwFactor->getRadiusParameter();
	radius->fixed( fixed );
	}

	return bwFactor;
	}

	LauAbsResonance* LauResonanceMaker::getResonance(const LauDaughters* daughters, const TString& resName, const Int_t resPairAmpInt, const LauAbsResonance::LauResonanceModel resType, const LauBlattWeisskopfFactor::BlattWeisskopfCategory bwCategory, const LauBlattWeisskopfFactor::BarrierType bwType)
	{
	// Routine to return the appropriate LauAbsResonance object given the resonance
	// name (resName), which daughter is the bachelor track (resPairAmpInt = 1,2 or 3),
	// and the resonance type ("BW" = Breit-Wigner, "Flatte" = Flatte distribution).

	// Loop over all possible resonance states we have defined in
	// createResonanceVector() until we get a match with the name of the resonance

	LauResonanceInfo* resInfo(0);
	for (std::vector<LauResonanceInfo*>::const_iterator iter=resInfo_.begin(); iter!=resInfo_.end(); ++iter) {

	if (resName == (*iter)->getName()) {
	// We have recognised the resonance name.
	std::cout<<"INFO in LauResonanceMaker::getResonance : Creating resonance: "<<resName<<std::endl;

	resInfo = (*iter);

	// stop looping
	break;
	}
	}

	// if we couldn't find the right resonance then we should return a null pointer
	if ( resInfo == 0 ) {
	- std::cout<<"WARNING in LauResonanceMaker::getResonance : Creating resonance: "<<resName<<std::endl;
	+ std::cout<<"ERROR in LauResonanceMaker::getResonance : Unable to locate resonance info for: "<<resName<<std::endl;
	return 0;
	}

	LauAbsResonance* theResonance(0);

	// Now construct the resonnace using the right type.
	// If we don't recognise the resonance model name, just use a simple Breit-Wigner.

	switch ( resType ) {

	case LauAbsResonance::BW :
	// Simple non-relativistic Breit-Wigner
	std::cout<<" : Using simple Breit-Wigner lineshape. "<<std::endl;
	theResonance = new LauBreitWignerRes(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::RelBW :
	{
	// Relativistic Breit-Wigner with Blatt-Weisskopf factors.
	std::cout<<" : Using relativistic Breit-Wigner lineshape. "<<std::endl;
	theResonance = new LauRelBreitWignerRes(resInfo, resPairAmpInt, daughters);
	LauBlattWeisskopfFactor::BlattWeisskopfCategory parCategory = LauBlattWeisskopfFactor::Parent;
	LauBlattWeisskopfFactor::BlattWeisskopfCategory resCategory = bwCategory;
	if ( bwCategory == LauBlattWeisskopfFactor::Default ) {
	resCategory = resInfo->getBWCategory();
	}
	LauBlattWeisskopfFactor* resBWFactor = this->getBWFactor( resCategory, resInfo, bwType );
	LauBlattWeisskopfFactor* parBWFactor = this->getBWFactor( parCategory, resInfo, bwType );
	theResonance->setBarrierRadii( resBWFactor, parBWFactor );
	break;
	}

	case LauAbsResonance::GS :
	{
	// Gounaris-Sakurai function to try and model the rho(770) better
	std::cout<<" : Using Gounaris-Sakurai lineshape. "<<std::endl;
	theResonance = new LauGounarisSakuraiRes(resInfo, resPairAmpInt, daughters);
	LauBlattWeisskopfFactor::BlattWeisskopfCategory parCategory = LauBlattWeisskopfFactor::Parent;
	LauBlattWeisskopfFactor::BlattWeisskopfCategory resCategory = bwCategory;
	if ( bwCategory == LauBlattWeisskopfFactor::Default ) {
	resCategory = resInfo->getBWCategory();
	}
	LauBlattWeisskopfFactor* resBWFactor = this->getBWFactor( resCategory, resInfo, bwType );
	LauBlattWeisskopfFactor* parBWFactor = this->getBWFactor( parCategory, resInfo, bwType );
	theResonance->setBarrierRadii( resBWFactor, parBWFactor );
	break;
	}

	case LauAbsResonance::Flatte :
	// Flatte distribution - coupled channel Breit-Wigner
	std::cout<<" : Using Flatte lineshape. "<<std::endl;
	theResonance = new LauFlatteRes(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::Sigma :
	// Sigma model - should only be used for the pi-pi system
	std::cout<<" : Using Sigma lineshape. "<<std::endl;
	theResonance = new LauSigmaRes(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::Kappa :
	// Kappa model - should only be used for the K-pi system
	std::cout<<" : Using Kappa lineshape. "<<std::endl;
	theResonance = new LauKappaRes(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::Dabba :
	// Dabba model - should only be used for the D-pi system
	std::cout<<" : Using Dabba lineshape. "<<std::endl;
	theResonance = new LauDabbaRes(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::LASS :
	// LASS function to try and model the K-pi S-wave better
	std::cout<<" : Using LASS lineshape. "<<std::endl;
	theResonance = new LauLASSRes(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::LASS_BW :
	// LASS function to try and model the K-pi S-wave better
	std::cout<<" : Using LASS lineshape (resonant part only). "<<std::endl;
	theResonance = new LauLASSBWRes(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::LASS_NR :
	// LASS function to try and model the K-pi S-wave better
	std::cout<<" : Using LASS lineshape (nonresonant part only). "<<std::endl;
	theResonance = new LauLASSNRRes(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::EFKLLM :
	// EFKLLM form-factor description of the K-pi S-wave
	std::cout<<" : Using EFKLLM lineshape. "<<std::endl;
	theResonance = new LauEFKLLMRes(resInfo, resPairAmpInt, daughters);
	break;

	+ case LauAbsResonance::KMatrix :
	+ // K-matrix description of S-wave
	+ std::cerr<<"ERROR in LauResonanceMaker::getResonance : K-matrix type specified, which should be separately handled."<<std::endl;
	+ break;
	+
	case LauAbsResonance::FlatNR :
	// uniform NR amplitude - arguments are there to preserve the interface
	std::cout<<" : Using uniform NR lineshape. "<<std::endl;
	// we override resPairAmpInt here
	theResonance = new LauFlatNR(resInfo, 0, daughters);
	break;

	case LauAbsResonance::NRModel :
	// NR amplitude model - arguments are there to preserve the interface
	std::cout<<" : Using NR-model lineshape. "<<std::endl;
	// we override resPairAmpInt here
	theResonance = new LauNRAmplitude(resInfo, 0, daughters);
	break;

	case LauAbsResonance::BelleNR :
	case LauAbsResonance::PowerLawNR :
	// Belle NR amplitude model - arguments are there to preserve the interface
	std::cout<<" : Using Belle NR lineshape. "<<std::endl;
	theResonance = new LauBelleNR(resInfo, resType, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::BelleSymNR :
	case LauAbsResonance::BelleSymNRNoInter :
	case LauAbsResonance::TaylorNR :
	// Belle NR amplitude model - arguments are there to preserve the interface
	std::cout<<" : Using Belle symmetric NR lineshape. "<<std::endl;
	theResonance = new LauBelleSymNR(resInfo, resType, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::PolNR :
	// Polynomial NR amplitude model - arguments are there to preserve the interface
	std::cout<<" : Using polynomial NR lineshape. "<<std::endl;
	theResonance = new LauPolNR(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::MIPW_MagPhase :
	// Model independent partial wave
	std::cout<<" : Using model independent partial wave lineshape (magnitude and phase). "<<std::endl;
	theResonance = new LauModIndPartWaveMagPhase(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::MIPW_RealImag :
	// Model independent partial wave
	std::cout<<" : Using model independent partial wave lineshape (real and imaginary part). "<<std::endl;
	theResonance = new LauModIndPartWaveRealImag(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::MIPW_Sym_MagPhase :
	// Model independent partial wave
	std::cout<<" : Using model independent partial wave lineshape for symmetric DPs (magnitude and phase). "<<std::endl;
	theResonance = new LauModIndPartWaveSymMagPhase(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::MIPW_Sym_RealImag :
	// Model independent partial wave
	std::cout<<" : Using model independent partial wave lineshape for symmetric DPs (real and imaginary part). "<<std::endl;
	theResonance = new LauModIndPartWaveSymRealImag(resInfo, resPairAmpInt, daughters);
	break;

	case LauAbsResonance::GaussIncoh :
	// Incoherent Gaussian
	std::cout<<" : Using incoherent Gaussian lineshape. "<<std::endl;
	theResonance = new LauGaussIncohRes(resInfo, resPairAmpInt, daughters);
	break;
	-
	- default :
	- std::cerr << "ERROR in LauResonanceMaker::getResonance : Could not match resonace type \"" << resType << "\"." << std::endl;
	- break;
	}

	return theResonance;
	-
	}

	Int_t LauResonanceMaker::resTypeInt(const TString& name) const
	{
	// Internal function that returns the resonance integer, specified by the
	// order of the resonance vector defined in createResonanceVector(),
	// for a given resonance name.
	Int_t resTypInt(-99);
	Int_t i(0);

	for (std::vector<LauResonanceInfo*>::const_iterator iter=resInfo_.begin(); iter!=resInfo_.end(); ++iter) {

	if (name.BeginsWith((*iter)->getName(), TString::kExact) == kTRUE) {
	// We have recognised the resonance from those that are available
	resTypInt = i;
	return resTypInt;
	}
	++i;
	}

	return resTypInt;
	}

	void LauResonanceMaker::printAll( std::ostream& stream ) const
	{
	for ( std::vector<LauResonanceInfo*>::const_iterator iter = resInfo_.begin(); iter != resInfo_.end(); ++iter ) {
	stream << (**iter) << std::endl;
	}
	}

File Metadata

Mime Type: text/x-diff
Expires: Tue, Nov 19, 6:17 PM (1 d, 19 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 3796684
Default Alt Text: (46 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions