Millepede.h

Go to the documentation of this file.

#ifndef MILLEPEDE_H
#define MILLEPEDE_H 1
 
/** @class Millepede Millepede.h AlignmentTools/Millepede.h
 *
 *
 *  @author Sebastien Viret
 *  @date   2005-07-29
 */
 
#include "MdcAlignAlg/Alignment.h"
#include <vector>
 
using namespace Alignment;
 
class Millepede {
public:
  /// Standard constructor
  Millepede();
  ~Millepede(); ///< Destructor
 
  /// Initialization
  bool initialize();
 
  bool InitMille( bool DOF[], double Sigm[], int nglo, int nloc, double startfact, int nstd,
                  double res_cut, double res_cut_init );
  bool MakeGlobalFit( double par[], double error[], double pull[] );
  bool ParGlo( int index, double param );
  bool ParSig( int index, double sigma );
  bool ConstF( double dercs[], double rhs );
  bool EquLoc( double dergb[], double derlc[], double dernl[], double rmeas, double sigma );
  bool ZerLoc( double dergb[], double derlc[], double dernl[] );
  bool FitLoc( int n, double track_params[], int single_fit );
  int  GetTrackNumber();
  void SetTrackNumber( int value );
 
private:
  // Option set-up here
 
  /* the following 4 statements are moved to Alignment.h by wulh on 06/08/29 */
  //   static const bool m_iteration  = true;   //  */
  //   static const bool debug_mode   = false;   // More printouts */
  //   static const bool verbose_mode = false;  // Definitely more printouts (matrices,...) */
  //   static const bool verbose_reject = false;//\* added by wulh on 06/08/28 *\/ */
 
  // Max. dimensions
 
  static const int mglobl = 400; // Max. number of global parameters
  static const int mlocal = 20;  // Max. number of local parameters
  static const int mcs    = 10;  // Max. number of constraint equations
  static const int mgl    = 410; // mglobl+mlocal
 
  // Private methods
 
  bool InitUn( double cutfac );
  bool PrtGlo();
 
  double ErrPar( int i );
  double CorPar( int i, int j );
  int    SpmInv( double v[][mgl], double b[], int n, double diag[], bool flag[] );
  int    SpmInv( double v[][mlocal], double b[], int n, double diag[], bool flag[] );
  bool   SpAVAt( double v[][mlocal], double a[][mlocal], double w[][mglobl], int n, int m );
  bool   SpAX( double a[][mlocal], double x[], double y[], int n, int m );
  double chindl( int n, int nd );
 
  // Matrices
 
  double cgmat[mgl][mgl];
  double clmat[mlocal][mlocal];
  double clcmat[mglobl][mlocal];
  double corrm[mglobl][mglobl];
  double adercs[mcs][mglobl];
 
  // Vectors and useful variables
 
  double corrv[mglobl], psigm[mglobl], pparm[mglobl], dparm[mglobl];
  double scdiag[mglobl], blvec[mlocal], arhs[mcs], diag[mgl], bgvec[mgl];
 
  int indgb[mglobl], nlnpa[mglobl], indnz[mglobl], indbk[mglobl], indlc[mlocal];
 
  bool scflag[mglobl];
 
  std::vector<int>    indst;
  std::vector<double> arest;
  std::vector<double> arenl;
 
  std::vector<int>    storeind;
  std::vector<int>    storeplace;
  std::vector<double> storeare;
  std::vector<double> storenl;
 
  int store_row_size;
 
  int    m_track_number;
  double m_residual_cut_init;
  double m_residual_cut;
 
  double cfactr, cfactref;
 
  int itert, nst, nfl, ncs, nstdev;
  int loctot, locrej, nagb, nalc, nrank;
};
#endif // MILLEPEDE_H