#ifndef _constants_H
#define _constants_H

//
// detector specifications
//

const unsigned int nStation = 4;// Number of detector planes
const unsigned int nSensors = 8;// Number of sensors on one side plane
const unsigned int nSides = 1; // Number of sides plane 1=double-sided sensor!
const unsigned int Stations[nStation*nSensors] = {0,1,2,3,4,5,6,7,
						  8,9,10,11,12,13,14,15,
						  16,17,18,19,20,21,22,23,
						  24,25,26,27,28,29,30,31};// stations from 0-31
/* const unsigned int Stations[nStation*nSensors] = {0,0,0,0,4,0,0,0, */
/* 						  0,0,0,0,12,0,0,0, */
/* 						  0,0,0,19,0,0,0,0, */
/* 						  0,0,0,0,28,0,0,0};// stations from 0-31 */
/* // */
/* // Misalignment scales */
/* // */
/* const double offsetX = 0.05; */
/* const double offsetY = 0.05; */
/* const double offsetZ = 0.05; */

/* const double offsetA = 0.005; */
/* const double offsetB = 0.005; */
/* const double offsetC = 0.005; */


//
// MILLEPEDE tuning
//

const bool m_dofs[6] = {1,1,1,1,1,1}; // Degrees of freedom to take into account (Boolean version)
//const bool m_dofs[6] = {1,1,0,0,0,1}; // Degrees of freedom to take into account (Boolean version)
//const bool m_dofs[6] = {1,1,0,0,0,0}; // Degrees of freedom to take into account (Boolean version)
//const bool m_dofs[6] = {0,0,0,1,1,1}; // Degrees of freedom to take into account (Boolean version)

// Individual constraint on alignment constants (set via ParSig in Millepede)
//const double m_Sigm[6] = {0.0010,0.0010,0.0010,0.00001,0.00001,0.00001}; // for Dx< 100 mkm!
//const double m_Sigm[6] = {0.010,0.010,0.010,0.00001,0.00001,0.00001};// for Dx ~ 200 mkm!
//const double m_Sigm[6] = {0.00010,0.00010,0.00010,0.001,0.001,0.001};// for Dx ~ 0, Dgamma < 1 rad
const double m_Sigm[6] = {0.0010,0.0010,0.0010,0.001,0.001,0.001};// !!! it's work for transl only
//const double m_Sigm[6] = {0.0000010,0.0000010,0.0000010,0.0001,0.0001,0.0001};// !!! does it work for rot only?
// Cuts on residuals (when fitting the tracks) (in mm)
// !!! Adapt them to your situation !!!

const double m_res_cut = 1.;
const double m_res_cut_init = 1.; // Larger for the first iteration (you are misaligned !!!)
const double m_start_chi_cut = 100.;

#endif