MINT2
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
Rho0Omega Class Reference

#include <Rho0Omega.h>

Inheritance diagram for Rho0Omega:
BW_BW ILineshape ILineshape MINT::FitParDependent MINT::IFitParDependent MINT::IFitParDependent MINT::IFitParRegister MINT::PolymorphVector< FitParRef > MINT::IFitParDependent

Public Member Functions

 Rho0Omega (const AssociatedDecayTree &tree, const std::string &namePrefix="")
 
virtual std::complex< double > getVal (IDalitzEvent &evt)
 
virtual void print (std::ostream &out=std::cout) const
 
virtual void print (IDalitzEvent &evt, std::ostream &out=std::cout)
 
virtual std::string name () const
 
virtual ~Rho0Omega ()
 
- Public Member Functions inherited from BW_BW
virtual double prSq () const
 
virtual double prSqForGofM () const
 
virtual double pABSq ()
 
virtual double prSqMax () const
 
 BW_BW (const AssociatedDecayTree &decay, const std::string &lineshapePrefix="", MINT::MinuitParameterSet *mps=0)
 
 BW_BW (const BW_BW &other)
 
virtual ~BW_BW ()
 
virtual DalitzCoordinate getDalitzCoordinate (double nSigma=3) const
 
virtual MINT::counted_ptr< IGenFct > generatingFunction () const
 
virtual void setGenerationLimits (double mi, double ma)
 
virtual int twoLPlusOne () const
 
virtual std::complex< double > EvtGenValue (IDalitzEvent &evt)
 
- Public Member Functions inherited from ILineshape
virtual ~ILineshape ()
 
- Public Member Functions inherited from MINT::FitParDependent
virtual unsigned int size () const
 
virtual const FitParRef & operator[] (unsigned int i) const
 
virtual FitParRef & operator[] (unsigned int i)
 
virtual bool changedSinceLastCall () const
 
virtual void rememberFitParValues ()
 
virtual bool registerFitParDependence (const IFitParDependent &fpd)
 
bool registerFitParDependence (const FitParRef &fpr)
 
void removeAllFitParDependencies ()
 
 FitParDependent (IFitParRegister *daddy=0)
 
 FitParDependent (const FitParDependent &other, IFitParRegister *newDaddy=0)
 
void listFitParDependencies (std::ostream &os=std::cout) const
 
- Public Member Functions inherited from MINT::PolymorphVector< FitParRef >
 PolymorphVector ()
 
 PolymorphVector (unsigned int N)
 
 PolymorphVector (unsigned int N, const FitParRef &c)
 
 PolymorphVector (const PolymorphVector &other)
 
 PolymorphVector (const typename std::vector< FitParRef > &other)
 
virtual ~PolymorphVector ()
 
std::vector< FitParRef > & theVector ()
 
const std::vector< FitParRef > & theVector () const
 
FitParRef & operator[] (unsigned int i)
 
const FitParRef & operator[] (unsigned int i) const
 
FitParRef & at (unsigned int i)
 
const FitParRef & at (unsigned int i) const
 
std::vector< FitParRef >::iterator begin ()
 
std::vector< FitParRef >::const_iterator begin () const
 
std::vector< FitParRef >::iterator end ()
 
std::vector< FitParRef >::const_iterator end () const
 
std::vector< FitParRef >::iterator find (const FitParRef &c)
 
std::vector< FitParRef >::const_iterator find (const FitParRef &c) const
 
FitParRef & front ()
 
const FitParRef & front () const
 
FitParRef & back ()
 
const FitParRef & back () const
 
unsigned int size () const
 
bool empty () const
 
void push_back (const FitParRef &c)
 
void pop_back ()
 
void erase (typename std::vector< FitParRef >::iterator pos)
 
void erase (typename std::vector< FitParRef >::iterator first, typename std::vector< FitParRef >::iterator last)
 
PolymorphVector< FitParRef > & operator= (const PolymorphVector< FitParRef > &other)
 
void clear ()
 
void resize (unsigned int N)
 
void resize (unsigned int N, const FitParRef &c)
 
 operator const typename std::vector< FitParRef > & () const
 
 operator typename std::vector< FitParRef > & ()
 
bool operator== (const MINT::PolymorphVector< FitParRef > &v2) const
 
bool operator!= (const MINT::PolymorphVector< FitParRef > &v2) const
 
bool operator< (const MINT::PolymorphVector< FitParRef > &v2) const
 
bool operator> (const MINT::PolymorphVector< FitParRef > &v2) const
 

Protected Member Functions

double GetAOmegatoRho0 () const
 
double GetEpsCoh () const
 
std::complex< double > GetExpiBeta () const
 
double GetDelta () const
 
double GetRho0Mass () const
 
double GetOmegaMass () const
 
double GetRho0Width () const
 
double GetOmegaWidth () const
 
double GetPipMass () const
 
double Getp (const double &mpipi) const
 
double Getq (const double &mpipi) const
 
double GetD1 (const double &q) const
 
double GetB1 (const double &q, const double &q0) const
 
double Getrhos (const double &mpipi) const
 
double GetGammas (const double &mpipi, const double &m0, const double &Gamma0) const
 
std::complex< double > BW (const double &mpipi, const double &m0, const double &Gamma0, const double &q, const double &q0) const
 
- Protected Member Functions inherited from BW_BW
const MINT::MinuitParameterSet * getMinuitParameterSet () const
 
MINT::MinuitParameterSet * getMinuitParameterSet ()
 
const std::string & prefix () const
 
bool substitutePDGForReco () const
 
virtual double twoBody_dgtPsq_in_MumsFrame (double mumsMass, double mA, double mB) const
 
virtual double twoBody_dgtPsq_in_MumsPDGFrame () const
 
virtual double twoBody_recodgtPsq_in_MumsPDGFrame () const
 
virtual double twoBody_dgtPsq_in_MumsRecoFrame ()
 
bool setEventPtr (IDalitzEvent &evt) const
 
IDalitzEvent * getEvent () const
 
bool setAllFitParameters ()
 
ResonancePropertiesList * resonancePropertiesList () const
 
const ResonanceProperties * resonanceProperties () const
 
virtual const ParticleProperties * mumsProperties () const
 
virtual ResonancePropertiesFitRef & mumsFittableProperties () const
 
virtual int mumsPID () const
 
virtual double mumsMass () const
 
virtual double mumsWidth () const
 
virtual double mumsRadius () const
 
virtual double globalRadius () const
 
virtual double Radius () const
 
double GetAlpha () const
 
virtual int lowestPossibleTwoLPlusOne () const
 
virtual int maxDaughterSpinSum () const
 
virtual int minDaughterSpinSum () const
 
virtual int minDaughterSpinSum2 () const
 
virtual int minDaughterSpinSum3 () const
 
virtual int minDaughterSpinSum4 () const
 
virtual int maxDaughterPairSpinSum (int i, int j) const
 
virtual int minDaughterPairSpinSum (int i, int j) const
 
virtual std::string mumsSpin () const
 
virtual int mumsSpinValue () const
 
virtual int mumsParity () const
 
virtual bool startOfDecayChain () const
 
virtual bool parityConservingL (int L) const
 
virtual double mumsRecoMass2 () const
 
virtual double mumsRecoMass () const
 
virtual MultiQuarkContent mumsQuarkContent () const
 
bool nonResonant () const
 
virtual int numDaughters () const
 
virtual TLorentzVector daughterP4 (int i) const
 
virtual double daughterPDGMass (const int &i) const
 
virtual double daughterWidth (int i) const
 
virtual double daughterRecoMass2 (int i) const
 
virtual double daughterRecoMass (int i) const
 
virtual std::string daughterSpin (int i) const
 
virtual int daughterSpinValue (int i) const
 
virtual MultiQuarkContent daughterQuarkContent (int i) const
 
virtual int daughterP (int i) const
 
virtual int dgtrsInternalParity () const
 
virtual MultiQuarkContent dgtrsQuarkContent () const
 
virtual bool isWeakDecay () const
 
virtual double Fr_PDG_BL ()
 
virtual double Fr_BELLE (double prSquared)
 
virtual double Fr_BELLE_Max ()
 
virtual double Fr ()
 
virtual double FrForGofM ()
 
virtual double FrMax ()
 
virtual double GofM ()
 
virtual std::complex< double > BreitWigner ()
 
virtual void resetInternals ()
 
virtual void resetPDG ()
 

Protected Attributes

const MINT::NamedParameter< double > _pr
 
const MINT::NamedParameter< double > _aRatioSq
 
const MINT::NamedParameter< double > _eps_coh
 
const MINT::NamedParameter< double > _beta
 
const MINT::NamedParameter< double > _delta
 
double _mRho0
 
double _mOmega
 
double _GRho0
 
double _GOmega
 
const double _mPip
 
- Protected Attributes inherited from BW_BW
MINT::MinuitParameterSet * _mps
 
std::string _prefix
 
MINT::NamedParameter< int > _normBF
 
bool _useGlobalRadius
 
const AssociatedDecayTree & _theDecay
 
double _gen_s_mi
 
double _gen_s_ma
 
ResonancePropertiesList * _RPL
 
ResonancePropertiesFitRef * _fittableResonancePropertiesPtr
 
MINT::FitParRef * _fittableGlobalRadiusPtr
 
- Protected Attributes inherited from MINT::PolymorphVector< FitParRef >
std::vector< FitParRef > _vec
 

Detailed Description

Definition at line 13 of file Rho0Omega.h.

Constructor & Destructor Documentation

◆ Rho0Omega()

Rho0Omega::Rho0Omega ( const AssociatedDecayTree &  tree,
const std::string &  namePrefix = "" 
)

Constructor Lineshape taking rho0-omega interference into account Default parameters taken from the Crystal Barrel Collaboration Phys. Lett. B 411, 354 (1997)

Definition at line 14 of file Rho0Omega.cpp.

15  : BW_BW(tree, namePrefix)
16  , _pr("Rho0Omega::pr", 0.1973)
17  , _aRatioSq("Rho0Omega::aRatioSq", 3.9)
18  , _eps_coh("Rho0Omega::eps_coh", 1.14)
19  , _beta("Rho0Omega::beta", -4.3*pi/180.0)
20  , _delta("Rho0Omega::delta", 2.5/GeV)
21  , _mRho0(-9999), _mOmega(-9999), _GRho0(-9999), _GOmega(-9999)
22  , _mPip(ParticlePropertiesList::mass(211)/GeV)
23 {
24 }
Rho0Omega::_pr
const MINT::NamedParameter< double > _pr
Definition: Rho0Omega.h:155
Rho0Omega::_delta
const MINT::NamedParameter< double > _delta
Definition: Rho0Omega.h:155
Rho0Omega::_mPip
const double _mPip
Definition: Rho0Omega.h:158
pi
static const double pi
Definition: CLHEPPhysicalConstants.h:44
BW_BW::BW_BW
BW_BW(const AssociatedDecayTree &decay, const std::string &lineshapePrefix="", MINT::MinuitParameterSet *mps=0)
Definition: BW_BW.cpp:26
Rho0Omega::_beta
const MINT::NamedParameter< double > _beta
Definition: Rho0Omega.h:155
Rho0Omega::_mOmega
double _mOmega
Definition: Rho0Omega.h:157
Rho0Omega::_GRho0
double _GRho0
Definition: Rho0Omega.h:157
Rho0Omega::_eps_coh
const MINT::NamedParameter< double > _eps_coh
Definition: Rho0Omega.h:155
GeV
static const double GeV
Definition: CLHEPSystemOfUnits.h:152
Rho0Omega::_mRho0
double _mRho0
Definition: Rho0Omega.h:157
ParticlePropertiesList::mass
static double mass(const std::string &name)
Definition: ParticlePropertiesList.cpp:50
Rho0Omega::_aRatioSq
const MINT::NamedParameter< double > _aRatioSq
Definition: Rho0Omega.h:155
Rho0Omega::_GOmega
double _GOmega
Definition: Rho0Omega.h:157

◆ ~Rho0Omega()

virtual Rho0Omega::~Rho0Omega ( )
inlinevirtual

Destructor

Definition at line 48 of file Rho0Omega.h.

48 {}

Member Function Documentation

◆ BW()

std::complex< double > Rho0Omega::BW ( const double &  mpipi,
const double &  m0,
const double &  Gamma0,
const double &  q,
const double &  q0 
) const
protected

Return Crystal Barrel measurement definition of the Breit-Wigner Note that the production barrier factor of Eq. 2, D1(p), has been removed as it is equivalent to the decay barrier factor of the rho0 mother which is already included higher up in the isobar chain

Definition at line 205 of file Rho0Omega.cpp.

207  {
208  const double s = mpipi*mpipi;
209  const double s0 = m0*m0;
210 
211  const double rhos = Getrhos(mpipi);
212  const double rhos0 = Getrhos(m0);
213  const double B1 = GetB1(q,q0);
214 
215  const double Gamma = Gamma0 * rhos / rhos0 * B1 * B1;
216 
217  const std::complex<double> bw =
218  m0 * Gamma0 / rhos0 * B1 /
219  std::complex<double>(s0 - s,-m0*Gamma);
220 
221  return bw;
222 }
s
static const double s
Definition: CLHEPSystemOfUnits.h:126
Rho0Omega::Getrhos
double Getrhos(const double &mpipi) const
Definition: Rho0Omega.cpp:184
Rho0Omega::GetB1
double GetB1(const double &q, const double &q0) const
Definition: Rho0Omega.cpp:179

◆ GetAOmegatoRho0()

double Rho0Omega::GetAOmegatoRho0 ( ) const
protected

Production amplitude ratio of omega to rho0 Amplitude ratio-squared determined by the Crystal Barrel measurement Default value is 3.9 +/- 0.4

Definition at line 94 of file Rho0Omega.cpp.

94  {
95  return sqrt(_aRatioSq);
96 }
Rho0Omega::_aRatioSq
const MINT::NamedParameter< double > _aRatioSq
Definition: Rho0Omega.h:155

◆ GetB1()

double Rho0Omega::GetB1 ( const double &  q,
const double &  q0 
) const
protected

Return normalised barrier factors

Definition at line 179 of file Rho0Omega.cpp.

179  {
180 
181  return GetD1(q)/GetD1(q0);
182 }
Rho0Omega::GetD1
double GetD1(const double &q) const
Definition: Rho0Omega.cpp:171

◆ GetD1()

double Rho0Omega::GetD1 ( const double &  q) const
protected

Return barrier factor pi+pi- radius set by the Crystal Barrel measurement Default value is 197.3 MeV/c

Definition at line 171 of file Rho0Omega.cpp.

171  {
172  const double qoverpr = q/_pr;
173 
174  const double D1 = qoverpr * sqrt(2.0/((qoverpr*qoverpr)+1.0));
175 
176  return D1;
177 }
Rho0Omega::_pr
const MINT::NamedParameter< double > _pr
Definition: Rho0Omega.h:155

◆ GetDelta()

double Rho0Omega::GetDelta ( ) const
protected

Magnitude of the off-diagonal element in the rho0-omega mixing matrix Derived from rho-omega interference Fixed parameter in the Crystal Barrel measurement Default value is 2.5 +/- 0.2 MeV

Definition at line 108 of file Rho0Omega.cpp.

108  {
109  return _delta;
110 }
Rho0Omega::_delta
const MINT::NamedParameter< double > _delta
Definition: Rho0Omega.h:155

◆ GetEpsCoh()

double Rho0Omega::GetEpsCoh ( ) const
protected

Magnitude of the off-diagonal element in the rho0-omega mixing matrix Free parameter in the Crystal Barrel measurement Default value is 1.14 +/- 0.15 +/- 0.13

Definition at line 98 of file Rho0Omega.cpp.

98  {
99  return _eps_coh;
100 }
Rho0Omega::_eps_coh
const MINT::NamedParameter< double > _eps_coh
Definition: Rho0Omega.h:155

◆ GetExpiBeta()

std::complex< double > Rho0Omega::GetExpiBeta ( ) const
protected

Phase of the off-diagonal element in the rho0-omega mixing matrix Free parameter in the Crystal Barrel measurement Default value is -4.3 +/- 1.8 +/- 3.9 degrees

Definition at line 102 of file Rho0Omega.cpp.

102  {
103  const double beta = _beta;
104 
105  return std::polar(1.0, beta);
106 }
Rho0Omega::_beta
const MINT::NamedParameter< double > _beta
Definition: Rho0Omega.h:155

◆ GetGammas()

double Rho0Omega::GetGammas ( const double &  mpipi,
const double &  m0,
const double &  Gamma0 
) const
protected

Return mass-dependent width Not given in the Crystal Barrel reference Given in Phys. Lett. B 354, 14 (1995)

Definition at line 193 of file Rho0Omega.cpp.

194  {
195  const double s = mpipi*mpipi;
196  const double s0 = m0*m0;
197  const double m2pi = GetPipMass()*GetPipMass();
198 
199  const double Gamma =
200  Gamma0*pow(((s-(4.0*m2pi))/(s0-(4.0*m2pi))), 1.5)*m0/mpipi;
201 
202  return Gamma;
203 }
Rho0Omega::GetPipMass
double GetPipMass() const
Definition: Rho0Omega.cpp:132
s
static const double s
Definition: CLHEPSystemOfUnits.h:126

◆ GetOmegaMass()

double Rho0Omega::GetOmegaMass ( ) const
protected

Return omega mass Not possible to vary this in the fit in the current implementation

Definition at line 117 of file Rho0Omega.cpp.

117  {
118  if(_mOmega < 0) _mOmega = ParticlePropertiesList::mass(223)/GeV;
119  return _mOmega;
120 }
Rho0Omega::_mOmega
double _mOmega
Definition: Rho0Omega.h:157
GeV
static const double GeV
Definition: CLHEPSystemOfUnits.h:152
ParticlePropertiesList::mass
static double mass(const std::string &name)
Definition: ParticlePropertiesList.cpp:50

◆ GetOmegaWidth()

double Rho0Omega::GetOmegaWidth ( ) const
protected

Return omega width Not possible to vary this in the fit in the current implementation

Definition at line 127 of file Rho0Omega.cpp.

127  {
128  if(_GOmega < 0) _GOmega = ParticlePropertiesList::width(223)/GeV;
129  return _GOmega;
130 }
ParticlePropertiesList::width
static double width(const std::string &name)
Definition: ParticlePropertiesList.cpp:66
GeV
static const double GeV
Definition: CLHEPSystemOfUnits.h:152
Rho0Omega::_GOmega
double _GOmega
Definition: Rho0Omega.h:157

◆ Getp()

double Rho0Omega::Getp ( const double &  mpipi) const
protected

Return pi+pi- breakup momentum in the rest frame of the mother of the pi+pi-

Definition at line 136 of file Rho0Omega.cpp.

136  {
137  const TLorentzVector p4pipi =
138  (BW_BW::daughterP4(0)+BW_BW::daughterP4(1))*(1.0/GeV);
139  const TLorentzVector p4init =
140  _theDecay.getParent()->getVal().p(getEvent())*(1.0/GeV);
141  const TLorentzVector p4rec = p4init-p4pipi;
142 
143  const double minit = p4init.M();
144  const double mrec = p4rec.M();
145 
146  const double m2init = minit*minit;
147  const double mpipi_p2_mrec = (mpipi + mrec)*(mpipi + mrec);
148  const double mpipi_m2_mrec = (mpipi - mrec)*(mpipi - mrec);
149 
150  const double p =
151  sqrt((m2init-mpipi_p2_mrec)*(m2init-mpipi_m2_mrec))/
152  (2.0*minit);
153 
154  if( (m2init-mpipi_p2_mrec)*(m2init-mpipi_m2_mrec) < 0.0 ){
155  std::cout << "p is negative. This should never happen." << std::endl;
156  exit(1);
157  }
158 
159  return p;
160 }
BW_BW::_theDecay
const AssociatedDecayTree & _theDecay
Definition: BW_BW.h:69
DDTree::getVal
const ValueType & getVal() const
Definition: DDTree.h:102
BW_BW::getEvent
IDalitzEvent * getEvent() const
Definition: BW_BW.cpp:147
BW_BW::daughterP4
virtual TLorentzVector daughterP4(int i) const
Definition: BW_BW.cpp:566
GeV
static const double GeV
Definition: CLHEPSystemOfUnits.h:152
DDTree::getParent
const DDTree< ValueType > * getParent() const
Definition: DDTree.h:261
AssociatedDecayTreeItem::p
TLorentzVector p(const IDalitzEvent *evt) const
Definition: AssociatedDecayTreeItem.cpp:48

◆ GetPipMass()

double Rho0Omega::GetPipMass ( ) const
protected

Return pi+ mass

Definition at line 132 of file Rho0Omega.cpp.

132  {
133  return _mPip;
134 }
Rho0Omega::_mPip
const double _mPip
Definition: Rho0Omega.h:158

◆ Getq()

double Rho0Omega::Getq ( const double &  mpipi) const
protected

Return pi+ breakup momentum in the rest frame of the pi+pi- system

Definition at line 162 of file Rho0Omega.cpp.

162  {
163  const double s = mpipi*mpipi;
164  const double m2pi = GetPipMass()*GetPipMass();
165 
166  const double q = sqrt(s-(4.0*m2pi))/2.0;
167 
168  return q;
169 }
Rho0Omega::GetPipMass
double GetPipMass() const
Definition: Rho0Omega.cpp:132
s
static const double s
Definition: CLHEPSystemOfUnits.h:126

◆ GetRho0Mass()

double Rho0Omega::GetRho0Mass ( ) const
protected

Return rho0 mass Not possible to vary this in the fit in the current implementation

Definition at line 112 of file Rho0Omega.cpp.

112  {
113  if(_mRho0 < 0) _mRho0 = ParticlePropertiesList::mass(113)/GeV;
114  return _mRho0;
115 }
GeV
static const double GeV
Definition: CLHEPSystemOfUnits.h:152
Rho0Omega::_mRho0
double _mRho0
Definition: Rho0Omega.h:157
ParticlePropertiesList::mass
static double mass(const std::string &name)
Definition: ParticlePropertiesList.cpp:50

◆ GetRho0Width()

double Rho0Omega::GetRho0Width ( ) const
protected

Return rho0 width Not possible to vary this in the fit in the current implementation

Definition at line 122 of file Rho0Omega.cpp.

122  {
123  if(_GRho0 < 0) _GRho0 = ParticlePropertiesList::width(113)/GeV;
124  return _GRho0;
125 }
ParticlePropertiesList::width
static double width(const std::string &name)
Definition: ParticlePropertiesList.cpp:66
Rho0Omega::_GRho0
double _GRho0
Definition: Rho0Omega.h:157
GeV
static const double GeV
Definition: CLHEPSystemOfUnits.h:152

◆ Getrhos()

double Rho0Omega::Getrhos ( const double &  mpipi) const
protected

Return pi+pi- phase space factor

Definition at line 184 of file Rho0Omega.cpp.

184  {
185  const double s = mpipi*mpipi;
186  const double m2pi = GetPipMass()*GetPipMass();
187 
188  const double rhos = 1.0/mpipi*sqrt(s - (4.0*m2pi));
189 
190  return rhos;
191 }
Rho0Omega::GetPipMass
double GetPipMass() const
Definition: Rho0Omega.cpp:132
s
static const double s
Definition: CLHEPSystemOfUnits.h:126

◆ getVal()

std::complex< double > Rho0Omega::getVal ( IDalitzEvent &  evt)
virtual

Evaluate rho0-omega lineshape from the Crystal Barrel measurement

Reimplemented from BW_BW.

Definition at line 26 of file Rho0Omega.cpp.

26  {
27  setEventPtr(evt);
28  resetInternals();
29 
30  const double mpipi = (BW_BW::daughterP4(0)+BW_BW::daughterP4(1)).M()/GeV;
31 
32  const double rho0_mass = GetRho0Mass();
33  const double rho0_width = GetRho0Width();
34  const double omega_mass = GetOmegaMass();
35  const double omega_width = GetOmegaWidth();
36 
37  const double p = Getp(mpipi);
38  const double q = Getq(mpipi);
39 
40  // Rho0 BW
41  const double q0_rho0 = Getq(rho0_mass);
42  const std::complex<double> bw_rho0 =
43  BW( mpipi, rho0_mass, rho0_width, q, q0_rho0);
44 
45  // Omega BW
46  const double q0_omega = Getq(omega_mass);
47  const std::complex<double> bw_omega =
48  BW( mpipi, omega_mass, omega_width, q, q0_omega);
49 
50  //Off-diagonal term
51  const double tanphi =
52  ((omega_mass*omega_width)-
53  (rho0_mass*GetGammas(omega_mass,rho0_mass,rho0_width)))/
54  ((omega_mass*omega_mass)-(rho0_mass*rho0_mass));
55  const double phi = atan(tanphi);
56 
57  const double Gammas_rho0 = GetGammas(mpipi, rho0_mass, rho0_width);
58 
59  const std::complex<double> off_diag =
60  std::polar(1.0, phi) * GetDelta() * (rho0_mass+omega_mass) /
61  std::complex<double>((omega_mass*omega_mass)-(rho0_mass*rho0_mass),
62  -((omega_mass*omega_width) - (rho0_mass*Gammas_rho0)));
63 
64  //Rho0-Omega function
65  const std::complex<double> rho0_omega =
66  Fr_PDG_BL()*sqrt(p*q) *
67  (bw_rho0 +
68  (GetAOmegatoRho0()*bw_omega*GetEpsCoh()*GetExpiBeta()*off_diag));
69 
70  return rho0_omega;
71 }
BW_BW::resetInternals
virtual void resetInternals()
Definition: BW_BW.cpp:1114
Rho0Omega::GetRho0Width
double GetRho0Width() const
Definition: Rho0Omega.cpp:122
BW_BW::setEventPtr
bool setEventPtr(IDalitzEvent &evt) const
Definition: BW_BW.cpp:141
Rho0Omega::GetDelta
double GetDelta() const
Definition: Rho0Omega.cpp:108
Rho0Omega::GetGammas
double GetGammas(const double &mpipi, const double &m0, const double &Gamma0) const
Definition: Rho0Omega.cpp:193
Rho0Omega::GetRho0Mass
double GetRho0Mass() const
Definition: Rho0Omega.cpp:112
Rho0Omega::GetOmegaWidth
double GetOmegaWidth() const
Definition: Rho0Omega.cpp:127
Rho0Omega::GetEpsCoh
double GetEpsCoh() const
Definition: Rho0Omega.cpp:98
BW_BW::Fr_PDG_BL
virtual double Fr_PDG_BL()
Definition: BW_BW.cpp:830
Rho0Omega::BW
std::complex< double > BW(const double &mpipi, const double &m0, const double &Gamma0, const double &q, const double &q0) const
Definition: Rho0Omega.cpp:205
BW_BW::daughterP4
virtual TLorentzVector daughterP4(int i) const
Definition: BW_BW.cpp:566
Rho0Omega::Getq
double Getq(const double &mpipi) const
Definition: Rho0Omega.cpp:162
GeV
static const double GeV
Definition: CLHEPSystemOfUnits.h:152
Rho0Omega::GetExpiBeta
std::complex< double > GetExpiBeta() const
Definition: Rho0Omega.cpp:102
Rho0Omega::GetOmegaMass
double GetOmegaMass() const
Definition: Rho0Omega.cpp:117
Rho0Omega::Getp
double Getp(const double &mpipi) const
Definition: Rho0Omega.cpp:136
Rho0Omega::GetAOmegatoRho0
double GetAOmegatoRho0() const
Definition: Rho0Omega.cpp:94

◆ name()

virtual std::string Rho0Omega::name ( ) const
inlinevirtual

Decay Name

Reimplemented from BW_BW.

Definition at line 41 of file Rho0Omega.h.

41  {
42  return "Rho0Omega("+_theDecay.oneLiner() +")";
43  }
BW_BW::_theDecay
const AssociatedDecayTree & _theDecay
Definition: BW_BW.h:69
DDTree::oneLiner
void oneLiner(std::stringstream &seam, int generation=0) const
Definition: DDTree.h:375

◆ print() [1/2]

void Rho0Omega::print ( std::ostream &  out = std::cout) const
virtual

Print decay

Reimplemented from BW_BW.

Definition at line 73 of file Rho0Omega.cpp.

73  {
74  out << name();
75 }
Rho0Omega::name
virtual std::string name() const
Definition: Rho0Omega.h:41

◆ print() [2/2]

void Rho0Omega::print ( IDalitzEvent &  evt,
std::ostream &  out = std::cout 
)
virtual

Print decay

Reimplemented from BW_BW.

Definition at line 77 of file Rho0Omega.cpp.

77  {
78  setEventPtr(evt);
79  resetInternals();
80  out << name()
81  << "\n\t> co-ordinate: " << getDalitzCoordinate()
82  << "\n\t> This is the decay I'm looking at:"
83  << "\n" << _theDecay
84  << "\n\t> These are a few values: "
85  << " startOfDecayChain? " << startOfDecayChain()
86  << ", BreitWigner = "
87  << BreitWigner()
88  << ", Blatt-Weisskopf penetration factor: "
89  << Fr()
90  << ", total BW_BW: "
91  << getVal(evt);
92 }
BW_BW::resetInternals
virtual void resetInternals()
Definition: BW_BW.cpp:1114
BW_BW::Fr
virtual double Fr()
Definition: BW_BW.cpp:710
BW_BW::_theDecay
const AssociatedDecayTree & _theDecay
Definition: BW_BW.h:69
BW_BW::setEventPtr
bool setEventPtr(IDalitzEvent &evt) const
Definition: BW_BW.cpp:141
BW_BW::startOfDecayChain
virtual bool startOfDecayChain() const
Definition: BW_BW.cpp:532
Rho0Omega::name
virtual std::string name() const
Definition: Rho0Omega.h:41
Rho0Omega::getVal
virtual std::complex< double > getVal(IDalitzEvent &evt)
Definition: Rho0Omega.cpp:26
BW_BW::BreitWigner
virtual std::complex< double > BreitWigner()
Definition: BW_BW.cpp:1093
BW_BW::getDalitzCoordinate
virtual DalitzCoordinate getDalitzCoordinate(double nSigma=3) const
Definition: BW_BW.cpp:537

Member Data Documentation

◆ _aRatioSq

const MINT::NamedParameter<double> Rho0Omega::_aRatioSq
protected

Definition at line 155 of file Rho0Omega.h.

◆ _beta

const MINT::NamedParameter<double> Rho0Omega::_beta
protected

Definition at line 155 of file Rho0Omega.h.

◆ _delta

const MINT::NamedParameter<double> Rho0Omega::_delta
protected

Definition at line 155 of file Rho0Omega.h.

◆ _eps_coh

const MINT::NamedParameter<double> Rho0Omega::_eps_coh
protected

Definition at line 155 of file Rho0Omega.h.

◆ _GOmega

double Rho0Omega::_GOmega
mutableprotected

Definition at line 157 of file Rho0Omega.h.

◆ _GRho0

double Rho0Omega::_GRho0
mutableprotected

Definition at line 157 of file Rho0Omega.h.

◆ _mOmega

double Rho0Omega::_mOmega
mutableprotected

Definition at line 157 of file Rho0Omega.h.

◆ _mPip

const double Rho0Omega::_mPip
protected

Definition at line 158 of file Rho0Omega.h.

◆ _mRho0

double Rho0Omega::_mRho0
mutableprotected

Definition at line 157 of file Rho0Omega.h.

◆ _pr

const MINT::NamedParameter<double> Rho0Omega::_pr
protected

Definition at line 155 of file Rho0Omega.h.

The documentation for this class was generated from the following files:
Generated by   doxygen 1.8.15