Cgem/CgemAlignAlg/CgemAlignAlg-00-00-08/CgemAlignAlg/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 <vector>
#include <iostream>
#include <fstream>
#include <fstream>
/* #include "MdcAlignAlg/Alignment.h" */
 
/* using namespace Alignment; */
 
class Millepede{
public: 
  /// Standard constructor
  Millepede();
  ~Millepede( ); ///< Destructor
 
  /// Initialization
  bool initialize();
 
  bool InitMille(bool DOF[], double Sigm[], int nlay, int nglo_on_lay
         , 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:
 
  //std::ofstream fmill;
 
// Option set-up here 
 
  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;
 
// Max. dimensions
 
  static const int mglobl       = 40;  // Max. number of global parameters
  static const int mlocal       = 10;  // Max. number of local parameters
  static const int mcs          = 36;  // Max. number of constraint equations, default value is 10
  //static const int mgl            = 410; // mglobl+mlocal
  static const int mgl          = 100; // mglobl+mlocal
 
 
  bool record_mill;   // More printouts */
  bool record_err;  // Definitely more printouts (matrices,...) */
 
 
// 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;
 
  // cfactref : Reference value for Chi^2/ndof cut 
  double cfactr, cfactref; // what is the first variables ? it is the initial cut factor !!
 
  int itert, nst, nfl, ncs, nstdev;
  int loctot, locrej, nagb, nglolay, nalc, nrank; 
 
};
#endif // MILLEPEDE_H