MomentumConstraint.h

Go to the documentation of this file.

// Class MomentumConstraint
//
// Author: Benno List, Jenny List
// Last update: $Date: 2008/02/18 09:59:34 $
//          by: $Author: blist $
// 
// Description: constraint 
// efact*sum(E_i) + pxfact*sum(p_x,i)+pyfact*sum(p_y,i)+pzfact*sum(p_z,i)-value = 0
//               

#ifndef __MOMENTUMCONSTRAINT_H
#define __MOMENTUMCONSTRAINT_H

#include "ParticleConstraint.h"

class ParticleFitObject;

//  Class PConstraint:

class MomentumConstraint : public ParticleConstraint {
  public:
    MomentumConstraint (double efact_=0,      
                        double pxfact_=0,     
                        double pyfact_=0,     
                        double pzfact_=0,     
                        double value_ = 0     
                       );
    virtual ~MomentumConstraint();
    virtual double getValue() const;
    virtual void getDerivatives(int idim,      
                                double der[]   
                               ) const;
                                   
    // not used at all
    //    virtual void addToGlobalDerMatrix (double lambda, int idim, double *M) const;
    
    virtual void invalidateCache() const;

    virtual int getVarBasis() const;
  
  protected:
    void updateCache() const;
  
  
    double efact;
    double pxfact;
    double pyfact;
    double pzfact;
    double value;
    
    mutable bool cachevalid;
    mutable int  nparams;
  
    virtual bool secondDerivatives (int i,                        
                                    int j,                        
                                    double *derivatives           
                                   ) const;
    virtual bool firstDerivatives (int i,                        
                                   double *derivatives           
                                  ) const;

    enum { VAR_BASIS=0 }; // this means that the constraint knows about E,px,py,pz

};

#endif // __MOMENTUMCONSTRAINT_H