Changeset View
Changeset View
Standalone View
Standalone View
examples/GenFit3pi.cc
Show First 20 Lines • Show All 148 Lines • ▼ Show 20 Lines | int main( int argc, char** argv ) | ||||
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 | // Reset the maximum signal DP ASq value | ||||
// This will be automatically adjusted to avoid bias or extreme | // This will be automatically adjusted to avoid bias or extreme | ||||
// inefficiency if you get the value wrong but best to set this by | // inefficiency if you get the value wrong but best to set this by | ||||
// hand once you've found the right value through some trial and error. | // hand once you've found the right value through some trial and error. | ||||
sigModel->setASqMaxValue(0.35); | //sigModel->setASqMaxValue(0.35); | ||||
sigModel->autoSetASqMaxValue(); | |||||
// 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) | ||||
std::vector<LauAbsCoeffSet*> coeffset; | std::vector<LauAbsCoeffSet*> coeffset; | ||||
▲ Show 20 Lines • Show All 94 Lines • Show Last 20 Lines |