Changeset View
Changeset View
Standalone View
Standalone View
examples/SimFitTask.cc
Show First 20 Lines • Show All 131 Lines • ▼ Show 20 Lines | int main( int argc, char** argv ) | ||||
LauIsobarDynamics* sigModel = new LauIsobarDynamics(daughters, effModel); | LauIsobarDynamics* sigModel = new LauIsobarDynamics(daughters, effModel); | ||||
LauAbsResonance* res(0); | LauAbsResonance* res(0); | ||||
res = sigModel->addResonance("f_0(980)", 3, LauAbsResonance::Flatte); // resPairAmpInt = 3 => resonance mass is m12. | res = sigModel->addResonance("f_0(980)", 3, LauAbsResonance::Flatte); // resPairAmpInt = 3 => resonance mass is m12. | ||||
res = sigModel->addResonance("f_2(1270)", 3, LauAbsResonance::RelBW); | res = sigModel->addResonance("f_2(1270)", 3, LauAbsResonance::RelBW); | ||||
res = sigModel->addResonance("K*0(892)", 1, LauAbsResonance::RelBW); | res = sigModel->addResonance("K*0(892)", 1, LauAbsResonance::RelBW); | ||||
res->fixMass(kFALSE); | res->fixMass(kFALSE); | ||||
res = sigModel->addResonance("K*0_0(1430)", 1, LauAbsResonance::LASS); | res = sigModel->addResonance("K*0_0(1430)", 1, LauAbsResonance::LASS); | ||||
// 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 | // be optimal. | ||||
// error. | // 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(1.25); | sigModel->setASqMaxValue(1.25); | ||||
TString integralsFileName("integ-"); | TString integralsFileName("integ-"); | ||||
integralsFileName += category; | integralsFileName += category; | ||||
integralsFileName += ".dat"; | integralsFileName += ".dat"; | ||||
sigModel->setIntFileName( integralsFileName ); | sigModel->setIntFileName( integralsFileName ); | ||||
// Create the fit model | // Create the fit model | ||||
▲ Show 20 Lines • Show All 100 Lines • Show Last 20 Lines |