calcuex Class Reference

class to do the calculus of a constraint based on a GiNaC symbolic expression More...

#include <model.h>

Inheritance diagram for calcuex:

Collaboration diagram for calcuex:

Public Member Functions
	calcuex (observable *ob, const ex &e0)
	~calcuex ()
double	operator() (const parameters &p) const
double	error (const parameters &p) const

Public Attributes
shared_ptr< observable >	o
ex	e

Detailed Description

class to do the calculus of a constraint based on a GiNaC symbolic expression

Definition at line 282 of file model.h.

Constructor & Destructor Documentation

calcuex::calcuex ( observable *  ob, const ex &  e0  )

inline

Definition at line 284 of file model.h.

: calcu(), o(ob), e(e0){}

calcuex::~calcuex ( )

inline

Definition at line 285 of file model.h.

{}

Member Function Documentation

double calcuex::error ( const parameters &  p ) const

inline

Definition at line 310 of file model.h.

References parameters::p.

                                            {
                 double ret=1000;
         int pass=1;
         try{
                 cout<<e<<endl;
                 cout<<e.subs(p.p)<<endl;
                 
                 ret=ex_to<numeric>(e.subs(p.p,subs_options::no_pattern).evalf()).to_double();
         }
         catch(GiNaC::pole_error er){
          pass=0;
          cout<<"Pole error"<<endl;
         }
         catch(exception er){
          pass=0;
          
          cout<<er.what()<<endl;
          cout<<e.subs(p.p,subs_options::no_pattern).evalf()<<endl;
          }
         catch(...){
          cout<<"Other exception"<<endl;
          exit(1);
         }
         if(pass) ret=o->error(ret);
         else ret=1000;
         
         return ret;
    }

double calcuex::operator() ( const parameters &  p ) const

inlinevirtual

Parameters

hipothesis

the theoretical hypothesis

Returns

the logarithm of the probability of measuring what was measured, assuming that the hypothesis is true

Implements calcu.

Definition at line 287 of file model.h.

References parameters::p.

                                                     {
                 double ret=1000;
         int pass=1;
         try{
                 ret=ex_to<numeric>(e.subs(p.p,subs_options::no_pattern).evalf()).to_double();
         }
         catch(GiNaC::pole_error e){
          pass=0;
          cout<<"Pole error"<<endl;
         }
         catch(exception e){
          cout<<e.what()<<endl;
          }
         catch(...){
          cout<<"Other exception"<<endl;
          exit(1);
         }
         if(pass) ret=o->loglikelihood(ret);
         else ret=1000;
         
         return ret;
    }

Member Data Documentation

ex calcuex::e

Definition at line 340 of file model.h.

shared_ptr<observable> calcuex::o

Definition at line 339 of file model.h.

---

The documentation for this class was generated from the following file:

• model.h