Abstract class for a model. More...

#include <model.h>

Inheritance diagram for Model:

Collaboration diagram for Model:

Public Member Functions
	Model ()

virtual	~Model ()

virtual parameters	getlist (const parameters &p) const =0

virtual parameters	generateparameters (int max=0) const =0

virtual int	veto (const parameters &p, int max=0) const

double	likelihood (const parameters &p, bool check=1, int max=0) const
	calculates the probability of getting all the experimental measures if the model describes the reality More...

double	loglike (const parameters &p, bool check=1, int max=0) const

Public Attributes
TRandom3 *	r

Detailed Description

Abstract class for a model.

Definition at line 361 of file model.h.

Constructor & Destructor Documentation

Model::Model ( )

inline

Definition at line 364 of file model.h.

364 : r(new TRandom3(0)){}

Model::r

TRandom3 * r

Definition: model.h:405

virtual Model::~Model ( )

inlinevirtual

Definition at line 365 of file model.h.

365 {delete r;};

Model::r

TRandom3 * r

Definition: model.h:405

Member Function Documentation

virtual parameters Model::generateparameters ( int max = 0 ) const

pure virtual

Implemented in BGLmodels::BGL, and BGL2.

virtual parameters Model::getlist ( const parameters & p ) const

pure virtual

Implemented in BGLmodels::BGL, BGL2, and Juca.

double Model::likelihood	(	const parameters &	p,
		bool	check = `1`,
		int	max = `0`
	)		const

inline

calculates the probability of getting all the experimental measures if the model describes the reality

Parameters

p	vector with the values of the free parameters

Definition at line 374 of file model.h.

                                                                       {
         if(veto(p,max) && check) return 0;
         double total=loglike(p,0);
         if(total<-1000) return 0;
         return exp(total);
         }

double Model::loglike	(	const parameters &	p,
		bool	check = `1`,
		int	max = `0`
	)		const

inline

Definition at line 381 of file model.h.

Referenced by main().

                                                                    {
         if(veto(p,max) && check) return -1000;
         parameters pp(getlist(p));
         
         double total=0; 
         int n=0;
         for(const_iterator i=begin();i!=end();i++) {
         try{
                 n++;
                 total+=i->loglikelihood(pp);
                                 }
     catch(exception e){
           cout<<n<<e.what()<<endl;
           exit(1);
          }
         //catch(...){
                 //cout<<"DD "<<n<<endl;
                 //exit(1);
                 //}
 }
 
         return -total;
         }