Changeset View
Changeset View
Standalone View
Standalone View
examples/GenFitBelleCPKpipi.cc
Show First 20 Lines • Show All 124 Lines • ▼ Show 20 Lines | int main( int argc, char** argv ) | ||||
res = posSigModel->addResonance("NonReson", 1, LauAbsResonance::BelleNR); | res = posSigModel->addResonance("NonReson", 1, LauAbsResonance::BelleNR); | ||||
res->setResonanceParameter("alpha", 0.50); | res->setResonanceParameter("alpha", 0.50); | ||||
res->floatResonanceParameter("alpha"); | res->floatResonanceParameter("alpha"); | ||||
// Set the file names for the integrals information (can be useful for debugging) | // Set the file names for the integrals information (can be useful for debugging) | ||||
negSigModel->setIntFileName("integ_neg.dat"); | negSigModel->setIntFileName("integ_neg.dat"); | ||||
posSigModel->setIntFileName("integ_pos.dat"); | posSigModel->setIntFileName("integ_pos.dat"); | ||||
// 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, | ||||
Double_t aSqMaxValue = 1.62; | // 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. | |||||
const Double_t aSqMaxValue{1.62}; | |||||
negSigModel->setASqMaxValue(aSqMaxValue); | negSigModel->setASqMaxValue(aSqMaxValue); | ||||
posSigModel->setASqMaxValue(aSqMaxValue); | posSigModel->setASqMaxValue(aSqMaxValue); | ||||
// Create the fit model, giving it both isobar models | // Create the fit model, giving it both isobar models | ||||
LauCPFitModel* fitModel = new LauCPFitModel( negSigModel, posSigModel ); | LauCPFitModel* fitModel = new LauCPFitModel( negSigModel, posSigModel ); | ||||
// Create the complex coefficients for the isobar model | // Create the complex coefficients for the isobar model | ||||
// Here we're using the "Belle" form: | // Here we're using the "Belle" form: | ||||
▲ Show 20 Lines • Show All 103 Lines • Show Last 20 Lines |