/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Mdc/MdcAlignAlg/MdcAlignAlg-00-01-04/MdcAlignAlg/Millepede.h

Go to the documentation of this file.

#ifndef MILLEPEDE_H 
#define MILLEPEDE_H 1

#include <vector>
#include "MdcAlignAlg/Alignment.h"

using namespace Alignment;

class Millepede{
public: 
  Millepede();
  ~Millepede( ); 

  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 

---