laura is hosted by Hepforge, IPPP Durham
Laura++  v1r1p1
A maximum likelihood fitting package for performing Dalitz-plot analysis.
LauIsobarDynamics.hh
Go to the documentation of this file.
1 
2 // Copyright University of Warwick 2005 - 2013.
3 // Distributed under the Boost Software License, Version 1.0.
4 // (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
5 
6 // Authors:
7 // Thomas Latham
8 // John Back
9 // Paul Harrison
10 
19 #ifndef LAU_ISOBAR_DYNAMICS
20 #define LAU_ISOBAR_DYNAMICS
21 
22 #include <vector>
23 
24 #include "TString.h"
25 
26 #include "LauAbsDPDynamics.hh"
27 #include "LauAbsResonance.hh"
28 #include "LauComplex.hh"
29 
30 class LauDaughters;
31 class LauEffModel;
32 class LauFitDataTree;
34 
36 
37  public:
39 
44  LauIsobarDynamics(LauDaughters* daughters, LauEffModel* effModel, LauEffModel* scfFractionModel = 0);
45 
47 
52  LauIsobarDynamics(LauDaughters* daughters, LauEffModel* effModel, LauTagCatScfFractionModelMap scfFractionModel);
53 
55  virtual ~LauIsobarDynamics();
56 
58 
61  virtual void initialise(const std::vector<LauComplex>& coeffs);
62 
64 
67  inline void setIntFileName(const TString& fileName) {intFileName_ = fileName;}
68 
69  // Integration
71 
75  virtual void setIntegralBinWidths(Double_t m13BinWidth, Double_t m23BinWidth);
76 
78 
86  virtual void addResonance(const TString& resName, Int_t resPairAmpInt, const TString& resType = "RelBW",
87  Double_t newMass = -1.0, Double_t newWidth = -1.0, Int_t newSpin = -1);
88 
90 
98  virtual void defineKMatrixPropagator(const TString& propName, const TString& paramFileName,
99  Int_t resPairAmpInt, Int_t nChannels, Int_t nPoles, Int_t rowIndex = 1);
100 
102 
107  virtual void addKMatrixProdPole(const TString& poleName, const TString& propName, Int_t poleIndex);
108 
110 
115  virtual void addKMatrixProdSVP(const TString& SVPName, const TString& propName, Int_t channelIndex);
116 
118 
125  virtual void changeResonance(const TString& resName, Double_t newMass = -1.0, Double_t newWidth = -1.0, Int_t newSpin = -1);
126 
128 
131  virtual Bool_t generate();
132 
134 
139  virtual ToyMCStatus checkToyMC(Bool_t printErrorMessages = kTRUE, Bool_t printInfoMessages = kFALSE);
140 
142 
145  inline Int_t maxGenIterations() const {return iterationsMax_;}
146 
148 
151  virtual void calcLikelihoodInfo(UInt_t iEvt);
152 
154 
158  virtual void calcLikelihoodInfo(Double_t m13Sq, Double_t m23Sq);
159 
161 
167  virtual void calcLikelihoodInfo(Double_t m13Sq, Double_t m23Sq, Int_t tagCat);
168 
170 
173  virtual const LauComplex& getEvtDPAmp() const {return totAmp_;}
174 
176 
179  virtual Double_t getEvtm13Sq() const {return m13Sq_;}
180 
182 
185  virtual Double_t getEvtm23Sq() const {return m23Sq_;}
186 
188 
191  virtual Double_t getEvtmPrime() const {return mPrime_;}
192 
194 
197  virtual Double_t getEvtthPrime() const {return thPrime_;}
198 
200 
203  virtual Double_t getEvtEff() const {return eff_;}
204 
206 
209  virtual Double_t getEvtScfFraction() const {return scfFraction_;}
210 
212 
215  virtual Double_t getEvtJacobian() const {return jacobian_;}
216 
218 
223  virtual Double_t getEvtLikelihood() const {return evtLike_;}
224 
226 
229  virtual void calcExtraInfo(Bool_t init = kFALSE);
230 
232 
236  virtual Bool_t gotReweightedEvent();
237 
239 
242  virtual Double_t getEventWeight();
243 
245 
248  inline void setASqMaxValue(Double_t value) {aSqMaxSet_ = value;}
249 
251 
254  inline Double_t getASqMaxSetValue() const { return aSqMaxSet_; }
255 
257 
260  inline Double_t getASqMaxVarValue() const { return aSqMaxVar_; }
261 
263 
267  virtual LauComplex getDynamicAmp(Int_t resID) const {return ff_[resID].scale(fNorm_[resID]);}
268 
270 
273  inline const std::vector< std::vector<LauComplex> >& getFiFjSum() const {return fifjSum_;}
274 
276 
279  inline const std::vector< std::vector<LauComplex> >& getFiFjEffSum() const {return fifjEffSum_;}
280 
282 
285  inline const std::vector<Double_t>& getFNorm() const {return fNorm_;}
286 
288 
291  virtual void fillDataTree(const LauFitDataTree& fitDataTree);
292 
294 
297  virtual void setDataEventNo(UInt_t iEvt);
298 
300 
303  inline void setBelleNRAlpha(Double_t alpha) {BelleNRAlpha_ = alpha;}
304 
306 
309  inline Double_t getBelleNRAlpha() const {return BelleNRAlpha_;}
310 
312 
321  inline void setLASSParameters(Double_t a, Double_t r, Double_t R = 1.0, Double_t phiR = 0.0, Double_t B = 1.0, Double_t phiB = 0.0, Double_t cutOff = 1.8)
322  {
323  this->setLASSScatteringLength(a);
324  this->setLASSEffectiveRange(r);
325  this->setLASSResonanceMag(R);
326  this->setLASSResonancePhase(phiR);
327  this->setLASSBackgroundMag(B);
328  this->setLASSBackgroundPhase(phiB);
329  this->setLASSCutOff(cutOff);
330  }
331 
333 
337 
339 
343 
345 
348  inline void setLASSResonanceMag(Double_t R) {LASSResonanceMag_ = R; changeLASSResonanceMag_ = kTRUE;}
349 
351 
354  inline void setLASSResonancePhase(Double_t phiR) {LASSResonancePhase_ = phiR; changeLASSResonancePhase_ = kTRUE;}
355 
357 
360  inline void setLASSBackgroundMag(Double_t B) {LASSBackgroundMag_ = B; changeLASSBackgroundMag_ = kTRUE;}
361 
363 
366  inline void setLASSBackgroundPhase(Double_t phiB) {LASSBackgroundPhase_ = phiB; changeLASSBackgroundPhase_ = kTRUE;}
367 
369 
372  inline void setLASSCutOff(Double_t cutOff) {LASSCutOff_ = cutOff; changeLASSCutOff_ = kTRUE;}
373 
375 
378  inline Double_t getLASSScatteringLength() const {return LASSScatteringLength_;}
379 
381 
384  inline Double_t getLASSEffectiveRange() const {return LASSEffectiveRange_;}
385 
387 
390  inline Double_t getLASSResonanceMag() const {return LASSResonanceMag_;}
391 
393 
396  inline Double_t getLASSResonancePhase() const {return LASSResonancePhase_;}
397 
399 
402  inline Double_t getLASSBackgroundMag() const {return LASSBackgroundMag_;}
403 
405 
408  inline Double_t getLASSBackgroundPhase() const {return LASSBackgroundPhase_;}
409 
411 
414  inline Double_t getLASSCutOff() const {return LASSCutOff_;}
415 
417 
421  inline void setFlatteParameters(Double_t g1, Double_t g2) {this->setFlatteg1(g1); this->setFlatteg2(g2);}
422 
424 
427  inline void setFlatteg1(Double_t g1) {FlatteParameterg1_ = g1; changeFlatteParameterg1_ = kTRUE;}
428 
430 
433  inline void setFlatteg2(Double_t g2) {FlatteParameterg2_ = g2; changeFlatteParameterg2_ = kTRUE;}
434 
436 
439  inline Double_t getFlatteg1() const {return FlatteParameterg1_;}
440 
442 
445  inline Double_t getFlatteg2() const {return FlatteParameterg2_;}
446 
448 
453  inline void setBarrierRadii( Double_t resRadius, Double_t parRadius = 4.0,
455  {
456  this->setResBarrierRadius(resRadius);
457  this->setParBarrierRadius(parRadius);
458  this->setBarrierType(type);
459  }
460 
462 
465  inline void setResBarrierRadius(Double_t radius) {resBarrierRadius_ = radius;}
466 
468 
471  inline void setParBarrierRadius(Double_t radius) {parBarrierRadius_ = radius;}
472 
474 
478 
480 
483  inline Double_t getResBarrierRadius() const {return resBarrierRadius_;}
484 
486 
489  inline Double_t getParBarrierRadius() const {return parBarrierRadius_;}
490 
492 
495  inline void flipHelicityForCPEigenstates(Bool_t boolean) {flipHelicity_ = boolean;}
496 
497  protected:
499  virtual void initSummary();
500 
502  virtual void initialiseVectors();
503 
505  virtual void calcDPNormalisation();
506 
508 
516  virtual void calcDPPartialIntegral(Double_t minm13, Double_t maxm13, Double_t minm23,
517  Double_t maxm23, Double_t m13BinWidth, Double_t m23BinWidth);
518 
520  virtual void writeIntegralsFile();
521 
523 
528  virtual void setFFTerm(UInt_t index, Double_t realPart, Double_t imagPart);
529 
531 
536  virtual void dynamics(Bool_t cacheResData = kTRUE, Double_t weight = 1.0, Bool_t useEff = kTRUE);
537 
539 
542  inline void setASqMaxVarValue(Double_t value) {aSqMaxVar_ = value;}
543 
545 
548  virtual Double_t calcSigDPNorm();
549 
551 
554  virtual LauComplex resAmp(Int_t index);
555 
557 
561  virtual LauAbsResonance* findResonance(const TString& name);
562 
564 
568  virtual const LauAbsResonance* findResonance(const TString& name) const;
569 
571 
574  virtual void removeCharge(TString& string) const;
575 
577 
582  Bool_t gotKMatrixMatch(UInt_t resAmpInt, const TString& propName) const;
583 
584  private:
586  typedef std::map<TString, LauKMatrixPropagator*> KMPropMap;
587 
589  typedef std::map<TString, TString> KMStringMap;
590 
592  std::vector<LauAbsResonance*> sigResonances_;
593 
596 
599 
601  std::vector<TString> resTypAmp_;
602 
604  std::vector<Int_t> resIntAmp_;
605 
607  std::vector<Int_t> resPairAmp_;
608 
610  std::vector<Int_t> typDaug_;
611 
614 
617 
619  TString intFileName_;
620 
622  Double_t m13BinWidth_;
623 
625  Double_t m23BinWidth_;
626 
628  Double_t m13Sq_;
629 
631  Double_t m23Sq_;
632 
634  Double_t mPrime_;
635 
637  Double_t thPrime_;
638 
640  Double_t eff_;
641 
643  Double_t scfFraction_;
644 
646  Double_t jacobian_;
647 
649  Double_t ASq_;
650 
652  Double_t evtLike_;
653 
656 
658 
661  std::vector< std::vector<LauComplex> > fifjEffSum_;
662 
664 
667  std::vector< std::vector<LauComplex> > fifjSum_;
668 
670  std::vector<LauComplex> ff_;
671 
673  std::vector<Double_t> fSqSum_;
674 
676  std::vector<Double_t> fNorm_;
677 
680 
683 
685  Double_t aSqMaxSet_;
686 
688  Double_t aSqMaxVar_;
689 
691  Double_t BelleNRAlpha_;
692 
695 
698 
701 
704 
707 
710 
712  Double_t LASSCutOff_;
713 
716 
719 
722 
725 
728 
731 
734 
737 
740 
743 
746 
749 
752 
755 
758 
759  ClassDef(LauIsobarDynamics,0)
760 
761 };
762 
763 #endif
void setParBarrierRadius(Double_t radius)
Set the radius of the barrier factor due to the parent to use for new amplitude components.
KMStringMap kMatrixPropSet_
The names of the M-matrix components in the model mapped to their propagators.
Bool_t changeLASSEffectiveRange_
Whether the default value of the LASS effective range has been changed.
Double_t getFlatteg1() const
Retrieve the g1 parameter for new Flatte components.
virtual Double_t getEvtJacobian() const
Retrieve the Jacobian, for the transformation into square DP coordinates, for the current event...
void setLASSCutOff(Double_t cutOff)
Set the cutoff value for new LASS components.
Int_t nSigGenLoop_
The number of unsucessful attempts to generate an event so far.
Class for defining the abstract interface for signal Dalitz plot dynamics.
Double_t getASqMaxVarValue() const
Retrieve the maximum of A squared that has been found while generating.
void setFlatteParameters(Double_t g1, Double_t g2)
Set the parameters for new Flatte components.
Bool_t changeLASSResonancePhase_
Whether the default value of the LASS resonance phase has been changed.
Double_t resBarrierRadius_
The radius of the resonance barrier factor for new amplitude components.
virtual LauComplex resAmp(Int_t index)
Calculate the dynamic part of the amplitude for a given component at the current point in the Dalitz ...
virtual ~LauIsobarDynamics()
Destructor.
const std::vector< Double_t > & getFNorm() const
Retrieve the normalisation factors for the dynamic parts of the amplitudes for all of the amplitude c...
Double_t FlatteParameterg1_
The constant parameter g1 for new Flatte components.
Double_t m13BinWidth_
The bin width to use when integrating over m13.
Bool_t flipHelicity_
The helicity flip flag for new amplitude components.
Double_t LASSBackgroundMag_
The background magnitude for new LASS components.
virtual void calcDPPartialIntegral(Double_t minm13, Double_t maxm13, Double_t minm23, Double_t maxm23, Double_t m13BinWidth, Double_t m23BinWidth)
Calculate the Dalitz plot normalisation integrals over a given range.
virtual Double_t getEvtLikelihood() const
Retrieve the likelihood for the current event.
Double_t getParBarrierRadius() const
Retrieve the radius of the barrier factor due to the parent to use for new amplitude components...
void setResBarrierRadius(Double_t radius)
Set the radius of the barrier factor due to the resonance to use for new amplitude components...
std::vector< Int_t > resPairAmp_
The index of the daughter not produced by the resonance for each amplitude component.
Double_t getLASSBackgroundPhase() const
Retrieve the background phase for new LASS components.
virtual void addResonance(const TString &resName, Int_t resPairAmpInt, const TString &resType="RelBW", Double_t newMass=-1.0, Double_t newWidth=-1.0, Int_t newSpin=-1)
Add a resonance to the Dalitz plot.
Class that defines the particular 3-body decay under study.
Definition: LauDaughters.hh:33
virtual void defineKMatrixPropagator(const TString &propName, const TString &paramFileName, Int_t resPairAmpInt, Int_t nChannels, Int_t nPoles, Int_t rowIndex=1)
Define a new K-matrix Propagator.
Double_t LASSScatteringLength_
The scattering length for new LASS components.
Bool_t changeLASSResonanceMag_
Whether the default value of the LASS resonance magnitude has been changed.
virtual void dynamics(Bool_t cacheResData=kTRUE, Double_t weight=1.0, Bool_t useEff=kTRUE)
Calculate the total Dalitz plot amplitude at the current point in the Dalitz plot.
void setASqMaxValue(Double_t value)
Set the maximum value of A squared to be used in the accept/reject.
Double_t aSqMaxVar_
The maximum value of A squared that has been seen so far while generating.
Double_t getLASSResonancePhase() const
Retrieve the resonance phase for new LASS components.
Double_t scfFraction_
The fraction of events that are poorly reconstructed (the self cross feed fraction) at the current po...
Double_t m13Sq_
The invariant mass squared of the first and third daughters.
Bool_t changeLASSBackgroundMag_
Whether the default value of the LASS background magnitude has been changed.
Int_t iterationsMax_
The maximum allowed number of attempts when generating an event.
LauComplex totAmp_
The total amplitude for the current event.
Double_t getASqMaxSetValue() const
Retrieve the maximum value of A squared to be used in the accept/reject.
LauAbsResonance::BarrierType barrierType_
The type of the barrier factor for new amplitude components.
Double_t aSqMaxSet_
The maximum allowed value of A squared.
Bool_t symmetricalDP_
Whether the Dalitz plot is symmetrical.
void setLASSResonancePhase(Double_t phiR)
Set the resonance phase for new LASS components.
Double_t parBarrierRadius_
The radius of the parent barrier factor for new amplitude components.
virtual void calcDPNormalisation()
Calculate the Dalitz plot normalisation integrals across the whole Dalitz plot.
virtual void setFFTerm(UInt_t index, Double_t realPart, Double_t imagPart)
Set the dynamic part of the amplitude for a given amplitude component at the current point in the Dal...
void setASqMaxVarValue(Double_t value)
Set the maximum of A squared that has been found.
KMPropMap kMatrixPropagators_
The K-matrix propagators.
std::vector< Double_t > fSqSum_
The event-by-event running total of the dynamical amplitude squared for each amplitude component...
virtual Double_t getEvtEff() const
Retrieve the efficiency for the current event.
Double_t m23Sq_
The invariant mass squared of the second and third daughters.
Double_t thPrime_
The square Dalitz plot coordinate theta&#39;.
Double_t evtLike_
The normalised likelihood for the current event.
std::vector< std::vector< LauComplex > > fifjEffSum_
The event-by-event running total of efficiency corrected amplitude cross terms for each pair of ampli...
virtual void initSummary()
Print a summary of the model to be used.
virtual void fillDataTree(const LauFitDataTree &fitDataTree)
Fill the internal data structure that caches the resonance dynamics.
void setLASSBackgroundPhase(Double_t phiB)
Set the background phase for new LASS components.
virtual LauAbsResonance * findResonance(const TString &name)
Retrieve the named resonance.
Double_t ASq_
The value of A squared for the current event.
Double_t getLASSBackgroundMag() const
Retrieve the background magnitude for new LASS components.
virtual const LauComplex & getEvtDPAmp() const
Retrieve the total amplitude for the current event.
virtual LauComplex getDynamicAmp(Int_t resID) const
Retrieve the normalised dynamic part of the amplitude of the given amplitude component at the current...
virtual void changeResonance(const TString &resName, Double_t newMass=-1.0, Double_t newWidth=-1.0, Int_t newSpin=-1)
Change the properties of a resonance particle within this model.
const std::vector< std::vector< LauComplex > > & getFiFjSum() const
Retrieve the event-by-event running totals of amplitude cross terms for all pairs of amplitude compon...
Double_t m23BinWidth_
The bin width to use when integrating over m23.
TString intFileName_
The name of the file to save integrals to.
void setLASSScatteringLength(Double_t a)
Set the scattering length for new LASS components.
virtual Double_t getEvtm13Sq() const
Retrieve the invariant mass squared of the first and third daughters in the current event...
void setFlatteg1(Double_t g1)
Set the g1 parameter for new Flatte components.
virtual void addKMatrixProdPole(const TString &poleName, const TString &propName, Int_t poleIndex)
Add a K-matrix production pole term to the model.
Bool_t changeFlatteParameterg2_
Whether the default value of the Flatte parameter g2 has been changed.
std::vector< TString > resTypAmp_
The resonance types of all of the amplitude components.
virtual void calcLikelihoodInfo(UInt_t iEvt)
Calculate the likelihood (and all associated information) for the given event number.
virtual Double_t getEvtm23Sq() const
Retrieve the invariant mass squared of the second and third daughters in the current event...
std::vector< Int_t > resIntAmp_
The index within the resonance maker for each amplitude component.
virtual void initialiseVectors()
Initialise the internal storage for this model.
Double_t jacobian_
The Jacobian, for the transformation into square DP coordinates at the current point in the Dalitz pl...
Double_t LASSEffectiveRange_
The effective range for new LASS components.
Double_t getFlatteg2() const
Retrieve the g2 parameter for new Flatte components.
File containing declaration of LauComplex class.
Double_t getLASSCutOff() const
Retrieve the cutoff value for new LASS components.
void setLASSResonanceMag(Double_t R)
Set the resonance magnitude for new LASS components.
std::vector< std::vector< LauComplex > > fifjSum_
The event-by-event running total of the amplitude cross terms for each pair of amplitude components...
virtual Double_t getEvtthPrime() const
Retrieve the square Dalitz plot coordinate, theta&#39;, for the current event.
std::map< Int_t, LauEffModel * > LauTagCatScfFractionModelMap
The type used for containing multiple self cross feed fraction models for different categories (e...
Bool_t changeLASSCutOff_
Whether the default value of the LASS cutoff has been changed.
virtual Double_t calcSigDPNorm()
Calculate the normalisation factor for the log-likelihood function.
virtual Double_t getEvtScfFraction() const
Retrieve the fraction of events that are poorly reconstructed (the self cross feed fraction) for the ...
void setLASSBackgroundMag(Double_t B)
Set the background magnitude for new LASS components.
void setBarrierRadii(Double_t resRadius, Double_t parRadius=4.0, LauAbsResonance::BarrierType type=LauAbsResonance::BWPrimeBarrier)
Set the parameters for the barrier factors for new resonances.
Double_t FlatteParameterg2_
The constant parameter g2 for new Flatte components.
void setBelleNRAlpha(Double_t alpha)
Set the alpha parameter for new Belle non-resonant components.
Bool_t changeFlatteParameterg1_
Whether the default value of the Flatte parameter g1 has been changed.
virtual Bool_t generate()
Generate a toy MC signal event.
File containing declaration of LauAbsDPDynamics class.
virtual void removeCharge(TString &string) const
Remove the charge from the given particle name.
std::vector< Int_t > typDaug_
The PDG codes of the daughters.
void setLASSParameters(Double_t a, Double_t r, Double_t R=1.0, Double_t phiR=0.0, Double_t B=1.0, Double_t phiB=0.0, Double_t cutOff=1.8)
Set the parameters for new LASS components.
void setFlatteg2(Double_t g2)
Set the g2 parameter for new Flatte components.
void flipHelicityForCPEigenstates(Bool_t boolean)
Set the helicity flip flag for new amplitude components.
Double_t eff_
The efficiency at the current point in the Dalitz plot.
Double_t LASSBackgroundPhase_
The background phase for new LASS components.
Double_t getResBarrierRadius() const
Retrieve the radius of the barrier factor due to the resonance to use for new amplitude components...
Class that implements the efficiency description across the signal Dalitz plot.
Definition: LauEffModel.hh:37
Int_t maxGenIterations() const
Retrieve the maximum number of iterations allowed when generating an event.
void setBarrierType(LauAbsResonance::BarrierType type)
Set the type of barrier factor to use for new amplitude components.
virtual Bool_t gotReweightedEvent()
Calculates whether an event with the current kinematics should be accepted in order to produce a dist...
Abstract class for defining type for resonance amplitude models (Breit-Wigner, Flatte etc...
Double_t LASSCutOff_
The cutoff value for new LASS components.
virtual void addKMatrixProdSVP(const TString &SVPName, const TString &propName, Int_t channelIndex)
Add a K-matrix slowly-varying part (SVP) term to the model.
Bool_t changeLASSScatteringLength_
Whether the default value of the LASS scattering length has been changed.
Class for defining signal dynamics using the isobar model.
virtual Double_t getEventWeight()
Calculate the acceptance rate, for events with the current kinematics, when generating events accordi...
void setIntFileName(const TString &fileName)
Set the name of the file to which to save the results of the integrals.
Double_t getBelleNRAlpha() const
Retreive the alpha parameter for new Belle non-resonant components.
ToyMCStatus
The possible statuses for toy MC generation.
virtual void writeIntegralsFile()
Write the results of the integrals (and related information) to a file.
virtual ToyMCStatus checkToyMC(Bool_t printErrorMessages=kTRUE, Bool_t printInfoMessages=kFALSE)
Check the status of the toy MC generation.
File containing declaration of LauAbsResonance class.
virtual void setIntegralBinWidths(Double_t m13BinWidth, Double_t m23BinWidth)
Set the widths of the bins to use when integrating across the Dalitz plot.
Class for defining a complex number.
Definition: LauComplex.hh:47
LauIsobarDynamics(LauDaughters *daughters, LauEffModel *effModel, LauEffModel *scfFractionModel=0)
Constructor.
Bool_t changeLASSBackgroundPhase_
Whether the default value of the LASS background phase has been changed.
Bool_t gotKMatrixMatch(UInt_t resAmpInt, const TString &propName) const
Check whether a resonance is a K-matrix component of a given propagator.
BarrierType
Define the allowed types of barrier factors.
std::map< TString, TString > KMStringMap
The type used for mapping K-matrix components to their propagators.
virtual void setDataEventNo(UInt_t iEvt)
Load the data for a given event.
std::map< TString, LauKMatrixPropagator * > KMPropMap
The type used for containing the K-matrix propagators.
virtual Double_t getEvtmPrime() const
Retrieve the square Dalitz plot coordinate, m&#39;, for the current event.
Double_t getLASSEffectiveRange() const
Retrieve the effective range for new LASS components.
Bool_t integralsDone_
Whether the integrals have been performed.
std::vector< LauComplex > ff_
The dynamic part of the amplitude for each amplitude component at the current point in the Dalitz plo...
Double_t getLASSScatteringLength() const
Retrieve the scattering length for new LASS components.
Double_t LASSResonancePhase_
The resonance phase for new LASS components.
void setLASSEffectiveRange(Double_t r)
Set the effective range for new LASS components.
Double_t getLASSResonanceMag() const
Retrieve the resonance magnitude for new LASS components.
Double_t LASSResonanceMag_
The resonance magnitude for new LASS components.
virtual void initialise(const std::vector< LauComplex > &coeffs)
Initialise the Dalitz plot dynamics.
const std::vector< std::vector< LauComplex > > & getFiFjEffSum() const
Retrieve the event-by-event running totals of efficiency corrected amplitude cross terms for all pair...
Double_t mPrime_
The square Dalitz plot coordinate, m&#39;.
std::vector< Double_t > fNorm_
The normalisation factors for the dynamic parts of the amplitude for each amplitude component...
Class to store the input fit variables.
Double_t BelleNRAlpha_
The alpha parameter for new Belle non-resonant components.
std::vector< LauAbsResonance * > sigResonances_
The resonances in the model.
virtual void calcExtraInfo(Bool_t init=kFALSE)
Calculate the fit fractions, mean efficiency and total DP rate.
Class for defining a K-matrix propagator.