MarlinKinfit
0.4.0
|
Implements a soft constraint with chi^2=((efact*sum(E_i) + pxfact*sum(p_x,i)+pyfact*sum(p_y,i)+pzfact*sum(p_z,i)-value)/sigma)^2. More...
#include <SoftGaussMomentumConstraint.h>
Public Member Functions | |
SoftGaussMomentumConstraint (double sigma_, double efact_, double pxfact_, double pyfact_, double pzfact_, double value_=0) | |
Constructor. More... | |
virtual | ~SoftGaussMomentumConstraint () |
Virtual destructor. | |
virtual double | getValue () const |
Returns the value of the constraint function. | |
virtual void | getDerivatives (int idim, double der[]) const |
Get first order derivatives. More... | |
![]() | |
SoftGaussParticleConstraint (double sigma_) | |
Creates an empty SoftGaussParticleConstraint object. More... | |
virtual | ~SoftGaussParticleConstraint () |
Virtual destructor. | |
virtual void | setFOList (std::vector< ParticleFitObject * > *fitobjects_) |
Adds several ParticleFitObject objects to the list. More... | |
virtual void | addToFOList (ParticleFitObject &fitobject, int flag=1) |
Adds one ParticleFitObject objects to the list. | |
virtual void | resetFOList () |
Resests ParticleFitObject list. | |
virtual double | getChi2 () const |
Returns the chi2. | |
virtual double | getError () const |
Returns the error on the value of the constraint. | |
virtual double | getSigma () const |
Returns the sigma. | |
virtual double | setSigma (double sigma_) |
Sets the sigma. More... | |
virtual void | add2ndDerivativesToMatrix (double *M, int idim) const |
Adds second order derivatives to global covariance matrix M. More... | |
virtual void | addToGlobalChi2DerVector (double *y, int idim) const |
Add derivatives of chi squared to global derivative matrix. More... | |
virtual void | invalidateCache () const |
Invalidates any cached values for the next event. | |
void | test1stDerivatives () |
void | test2ndDerivatives () |
double | num1stDerivative (int ifo, int ilocal, double eps) |
Evaluates numerically the 1st derivative w.r.t. a parameter. More... | |
double | num2ndDerivative (int ifo1, int ilocal1, double eps1, int ifo2, int ilocal2, double eps2) |
Evaluates numerically the 2nd derivative w.r.t. 2 parameters. More... | |
int | getVarBasis () const |
![]() | |
virtual | ~BaseSoftConstraint () |
Virtual destructor. | |
![]() | |
BaseConstraint () | |
Creates an empty BaseConstraint object. | |
BaseConstraint (const BaseConstraint &rhs) | |
Copy constructor. More... | |
BaseConstraint & | operator= (const BaseConstraint &rhs) |
Assignment. More... | |
virtual | ~BaseConstraint () |
Virtual destructor. | |
virtual const char * | getName () const |
Returns the name of the constraint. | |
virtual void | setName (const char *name_) |
Set object's name. | |
virtual std::ostream & | print (std::ostream &os) const |
print object to ostream More... | |
Protected Member Functions | |
virtual bool | secondDerivatives (int i, int j, double *derivatives) const |
Second derivatives with respect to the 4-vectors of Fit objects i and j; result false if all derivatives are zero. More... | |
virtual bool | firstDerivatives (int i, double *derivatives) const |
First derivatives with respect to the 4-vector of Fit objects i; result false if all derivatives are zero. More... | |
Protected Attributes | |
double | efact |
double | pxfact |
double | pyfact |
double | pzfact |
double | value |
![]() | |
FitObjectContainer | fitobjects |
The FitObjectContainer. | |
std::vector< double > | derivatives |
The derivatives. | |
std::vector< int > | flags |
The flags can be used to divide the FitObjectContainer into several subsets used for example to implement an equal mass constraint (see MassConstraint). More... | |
double | sigma |
The sigma of the Gaussian. | |
![]() | |
char * | name |
![]() | |
char * | name |
Additional Inherited Members | |
![]() | |
enum | { VAR_BASIS =BaseDefs::VARBASIS_EPXYZ } |
typedef std::vector < ParticleFitObject * > | FitObjectContainer |
Vector of pointers to ParticleFitObjects. | |
typedef FitObjectContainer::iterator | FitObjectIterator |
Iterator through vector of pointers to ParticleFitObjects. | |
typedef FitObjectContainer::const_iterator | ConstFitObjectIterator |
Constant iterator through vector of pointers to ParticleFitObjects. | |
![]() | |
std::ostream & | operator<< (std::ostream &os, const BaseConstraint &bc) |
Prints out a BaseConstraint, using its print method. More... | |
Implements a soft constraint with chi^2=((efact*sum(E_i) + pxfact*sum(p_x,i)+pyfact*sum(p_y,i)+pzfact*sum(p_z,i)-value)/sigma)^2.
This class implements soft constraints on the sum of energies, px, py, and pz. It can be used to for constraints of the form
Author: Jenny List, Benno List Last update:
by:
SoftGaussMomentumConstraint::SoftGaussMomentumConstraint | ( | double | sigma_, |
double | efact_, | ||
double | pxfact_, | ||
double | pyfact_, | ||
double | pzfact_, | ||
double | value_ = 0 |
||
) |
Constructor.
sigma_ | The sigma value |
efact_ | Factor for energy sum |
pxfact_ | Factor for px sum |
pyfact_ | Factor for py sum |
pzfact_ | Factor for pz sum |
value_ | Target value of sum |
|
protectedvirtual |
First derivatives with respect to the 4-vector of Fit objects i; result false if all derivatives are zero.
i | number of 1st FitObject |
derivatives | The result 4-vector |
Implements SoftGaussParticleConstraint.
|
virtual |
Get first order derivatives.
Call this with a predefined array "der" with the necessary number of entries!
idim | First dimension of the array |
der | Array of derivatives, at least idim x idim |
Implements SoftGaussParticleConstraint.
References SoftGaussParticleConstraint::fitobjects.
|
protectedvirtual |
Second derivatives with respect to the 4-vectors of Fit objects i and j; result false if all derivatives are zero.
i | number of 1st FitObject |
j | number of 2nd FitObject |
derivatives | The result 4x4 matrix |
Implements SoftGaussParticleConstraint.