Changeset View
Changeset View
Standalone View
Standalone View
examples/GenFit3pi.cc
Show First 20 Lines • Show All 144 Lines • ▼ Show 20 Lines | int main( int argc, char** argv ) | ||||
reson = sigModel->addResonance("f_0(980)", 1, LauAbsResonance::Flatte); | reson = sigModel->addResonance("f_0(980)", 1, LauAbsResonance::Flatte); | ||||
reson->setResonanceParameter("g1",0.2); | reson->setResonanceParameter("g1",0.2); | ||||
reson->setResonanceParameter("g2",1.0); | reson->setResonanceParameter("g2",1.0); | ||||
reson = sigModel->addResonance("f_2(1270)", 1, LauAbsResonance::RelBW); | reson = sigModel->addResonance("f_2(1270)", 1, LauAbsResonance::RelBW); | ||||
const TString nrName = "BelleNR_Swave"; | const TString nrName = "BelleNR_Swave"; | ||||
reson = sigModel->addResonance(nrName, 0, LauAbsResonance::BelleSymNRNoInter); | reson = sigModel->addResonance(nrName, 0, LauAbsResonance::BelleSymNRNoInter); | ||||
reson->setResonanceParameter("alpha", 0.2); | reson->setResonanceParameter("alpha", 0.2); | ||||
// Reset the maximum signal DP ASq value | // Set the maximum signal DP ASq value | ||||
// This will be automatically adjusted to avoid bias or extreme | // If you do not provide a value, one will be determined automatically, | ||||
// inefficiency if you get the value wrong but best to set this by | // which should be close to the true maximum but is not guaranteed to | ||||
// hand once you've found the right value through some trial and error. | // be optimal. | ||||
// Any value, whether manually provided or automatically determined, | |||||
// will be automatically adjusted to avoid bias or extreme inefficiency | |||||
// but it is best to set this by hand once you've found the right value | |||||
// through some trial and error. | |||||
sigModel->setASqMaxValue(0.35); | sigModel->setASqMaxValue(0.35); | ||||
// Create the fit model | // Create the fit model | ||||
LauSimpleFitModel* fitModel = new LauSimpleFitModel(sigModel); | LauSimpleFitModel* fitModel = new LauSimpleFitModel(sigModel); | ||||
// Create the complex coefficients for the isobar model | // Create the complex coefficients for the isobar model | ||||
// Here we're using the magnitude and phase form: | // Here we're using the magnitude and phase form: | ||||
// c_j = a_j exp(i*delta_j) | // c_j = a_j exp(i*delta_j) | ||||
▲ Show 20 Lines • Show All 95 Lines • Show Last 20 Lines |