ParticleFitObject.h

Go to the documentation of this file.

#ifndef __PARTICLEFITOBJECT_H
#define __PARTICLEFITOBJECT_H

#include "BaseFitObject.h"
#include "FourVector.h"


// Class ParticleFitObject

class ParticleFitObject: public BaseFitObject {
  public:
    ParticleFitObject();
    
        
    ParticleFitObject (const ParticleFitObject& rhs              
                   );
                   
    ParticleFitObject& operator= (const ParticleFitObject& rhs   
                             );

    virtual ~ParticleFitObject();
    
    virtual ParticleFitObject& assign (const BaseFitObject& source   
                                  );

    virtual bool setMass (double mass_);
    virtual double getMass () const;

    virtual std::ostream& print4Vector (std::ostream& os   
                                       ) const;
        
    virtual FourVector getFourMomentum() const;

    virtual double getE() const;
    virtual double getPx() const;
    virtual double getPy() const;
    virtual double getPz() const;
    
    virtual double getP() const;
    virtual double getP2() const;
    virtual double getPt() const;
    virtual double getPt2() const;
    
    virtual double getDPx (int ilocal              
                           ) const = 0;
    virtual double getDPy (int ilocal              
                           ) const = 0;
    virtual double getDPz (int ilocal              
                           ) const = 0;
    virtual double getDE (int ilocal               
                          ) const = 0;
    
    virtual void getDerivatives (double der[], int idim) const;

    virtual void addToGlobalChi2DerMatrixNum (double *M,   
                                              int idim,    
                                              double eps   
                                              );

    virtual void addToGlobalChi2DerVectorNum (double *y,    
                                              int idim,     
                                              double eps    
                                             );
    
    virtual std::ostream& print (std::ostream& os    
                                ) const;
                                
    void test1stDerivatives ();
    void test2ndDerivatives ();
                                
    double num1stDerivative (int ilocal,  
                             double eps   
                            );

    double num2ndDerivative (int ilocal1, 
                             double eps1, 
                             int ilocal2, 
                             double eps2  
                            );

    virtual double getChi2 () const;

  protected:
    double mass;

    mutable FourVector fourMomentum;

    // this is to flag phi angle, for example
    double paramCycl[BaseDefs::MAXPAR];

};

#endif // __PARTICLEFITOBJECT_H