Changeset View
Changeset View
Standalone View
Standalone View
examples/GenFitEFKLLM.cc
| Show First 20 Lines • Show All 112 Lines • ▼ Show 20 Lines | int main( int argc, char** argv ) | ||||
| LauAbsResonance* res(0); | LauAbsResonance* res(0); | ||||
| //addResonance arguments: resName, resPairAmpInt, resType | //addResonance arguments: resName, resPairAmpInt, resType | ||||
| res = sigModel->addResonance("kappa0", 1, LauAbsResonance::EFKLLM); | res = sigModel->addResonance("kappa0", 1, LauAbsResonance::EFKLLM); | ||||
| res->setResonanceParameter("massFactor", 0.0); | res->setResonanceParameter("massFactor", 0.0); | ||||
| res = sigModel->addResonance("K*0_0(1430)", 1, LauAbsResonance::EFKLLM); | res = sigModel->addResonance("K*0_0(1430)", 1, LauAbsResonance::EFKLLM); | ||||
| res->setResonanceParameter("massFactor", -2.0); | res->setResonanceParameter("massFactor", -2.0); | ||||
| // 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(14.5); | sigModel->setASqMaxValue(14.5); | ||||
| // 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 69 Lines • Show Last 20 Lines | |||||