Struct rustitude_core::manager::ExtendedLogLikelihood

pub struct ExtendedLogLikelihood {
    pub data_manager: Manager,
    pub mc_manager: Manager,
}

The ExtendedLogLikelihood stores two Managers, one for data and one for a Monte-Carlo dataset used for acceptance correction. These should probably have the same Manager in practice, but this is left to the user.

Fields

data_manager: Manager

Manager for data

mc_manager: Manager

Manager for Monte-Carlo

Implementations

impl ExtendedLogLikelihood

pub const fn new(data_manager: Manager, mc_manager: Manager) -> Self

Create a new ExtendedLogLikelihood from a data and Monte-Carlo Managers.

pub fn evaluate(&self, parameters: &[f64]) -> Result<f64, RustitudeError>

Evaluate the ExtendedLogLikelihood over the Dataset with the given free parameters.

Errors

This method will return a RustitudeError if the amplitude calculation fails. See Model::compute for more information.

pub fn par_evaluate( &self, parameters: &[f64], num_threads: usize, ) -> Result<f64, RustitudeError>

Evaluate the ExtendedLogLikelihood over the Dataset with the given free parameters This method also allows the user to input a maximum number of threads to use in the calculation, as it uses a parallel loop over events.

Errors

This method will return a RustitudeError if the amplitude calculation fails. See Model::compute for more information.

pub fn intensity( &self, parameters: &[f64], dataset: &Dataset, ) -> Result<Vec<f64>, RustitudeError>

Evaluate the unnormalized intensity function over the given Dataset with the given free parameters.

Errors

This method will return a RustitudeError if the amplitude calculation fails. See Model::compute for more information.

pub fn par_intensity( &self, parameters: &[f64], dataset: &Dataset, num_threads: usize, ) -> Result<Vec<f64>, RustitudeError>

Evaluate the unnormalized intensity function over the given Dataset with the given free parameters. This method also allows the user to input a maximum number of threads to use in the calculation. This version uses a parallel loop over events.

Errors

This method will return a RustitudeError if the amplitude calculation fails. See Model::compute for more information.

pub fn norm_int( &self, parameters: &[f64], weighted: bool, ) -> Result<f64, RustitudeError>

Find the normalization integral for the Model over the Dataset with the given free parameters.

Errors

This method will return a RustitudeError if the amplitude calculation fails. See Model::norm_int for more information.

pub fn par_norm_int( &self, parameters: &[f64], num_threads: usize, weighted: bool, ) -> Result<f64, RustitudeError>

Find the normalization integral for the Model over the Dataset with the given free parameters.

Errors

This method will return a RustitudeError if the amplitude calculation fails. See Model::norm_int for more information.

pub fn get_amplitude( &self, amplitude_name: &str, ) -> Result<Amplitude, RustitudeError>

Get a copy of an Amplitude in the Model by name.

Errors

This method will return a RustitudeError if there is no amplitude found with the given name in the parent Model. See Model::get_amplitude for more information.

pub fn get_parameter( &self, amplitude_name: &str, parameter_name: &str, ) -> Result<Parameter, RustitudeError>

Get a copy of a Parameter in a Model by name and the name of the parent Amplitude.

Errors

This method will return a RustitudeError if there is no parameter found with the given name in the parent Model. It will also first check if the given amplitude exists, and this method can also fail in the same way (see Model::get_amplitude and Model::get_parameter).

pub fn print_parameters(&self)

Print the free parameters in the Model. See Model::print_parameters for more information.

pub fn constrain( &mut self, amplitude_1: &str, parameter_1: &str, amplitude_2: &str, parameter_2: &str, ) -> Result<(), RustitudeError>

Constrain two parameters by name, reducing the number of free parameters by one.

Errors

This method will fail if any of the given amplitude or parameter names don’t correspond to a valid amplitude-parameter pair. See Model::constrain for more information.

pub fn fix( &mut self, amplitude: &str, parameter: &str, value: f64, ) -> Result<(), RustitudeError>

Fix a parameter by name to the given value.

Errors

This method will fail if the given amplitude-parameter pair does not exist. See Model::fix for more information.

pub fn free( &mut self, amplitude: &str, parameter: &str, ) -> Result<(), RustitudeError>

Free a fixed parameter by name.

Errors

This method will fail if the given amplitude-parameter pair does not exist. See Model::free for more information.

pub fn set_bounds( &mut self, amplitude: &str, parameter: &str, bounds: (f64, f64), ) -> Result<(), RustitudeError>

Set the bounds of a parameter by name.

Errors

This method will fail if the given amplitude-parameter pair does not exist. See Model::set_bounds for more information.

pub fn set_initial( &mut self, amplitude: &str, parameter: &str, initial: f64, ) -> Result<(), RustitudeError>

Set the initial value of a parameter by name.

Errors

This method will fail if the given amplitude-parameter pair does not exist. See Model::set_initial for more information.

pub fn get_bounds(&self) -> Vec<(f64, f64)>

Get a list of bounds for all free parameters in the Model. See Model::get_bounds for more information.

pub fn get_initial(&self) -> Vec<f64>

Get a list of initial values for all free parameters in the Model. See Model::get_initial for more information.

pub fn get_n_free(&self) -> usize

Get the number of free parameters in the Model See Model::get_n_free for more information.

pub fn activate(&mut self, amplitude: &str)

Activate an Amplitude by name. See Model::activate for more information.

pub fn deactivate(&mut self, amplitude: &str)

Deactivate an Amplitude by name. See Model::deactivate for more information.