CgemAlignment.h

Go to the documentation of this file.

#ifndef CGEMALIGNMENT_H
#define CGEMALIGNMENT_H
 
#include <string>
#include <vector>
#include <math.h>
 
namespace CgemAlignment{
 
 
    void Hello();
 
    const int LAYERNMAX = 6;
    const int GEO_METHOD = 1; 
    const int NEP = 6;  /* number of Cgem layer */
    const int NLAYER = 6;   /* number of Cgem layer */
    const int NTRKPAR = 5;  /* number of track parameters for staight line, the 3rd parameter is not used*/
    const int NTRKPARALL = 10;  /* number of track parameters and errors */
    const int NDOFALIGN = 6;    /* dx, dy, dz, rx, ry, rz */
 
    const bool g_dofs[6] = {1, 1, 1, 1, 1, 1}; // {Dx, Dy, Dz, Rx, Ry, Rz} Degrees of freedom to take into account (Boolean version)
 
 
    const double sta_fac = 1.0 / sqrt(10000.0/10000.0); // a scale factor dependent on statistics. The value is 1/sqrt(N/10000), where N is the generated MC events
    const double g_Sigm[36] = {
        /* L1                L2                        L3    */  /* L4                L5                        L6   */  
        sta_fac*0.07,    sta_fac*0.07,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,   /* dx  */
        sta_fac*0.07,    sta_fac*0.07,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,   /* dy  */
        sta_fac*0.07,    sta_fac*0.07,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,    sta_fac*0.07*0.7,   /* dz  */
        sta_fac*0.005,   sta_fac*0.005,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7,  /* rx  */
        sta_fac*0.005,   sta_fac*0.005,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7,  /* ry  */
        sta_fac*0.005,   sta_fac*0.005,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7,   sta_fac*0.005*0.7   /* rz  */
    };  /* for dx, dy, dz, rz unit ? assume they are mm and rad*/
 
    const double g_res_cut = 10.;    //Why do I set so small value : 1.2
    const double g_res_cut_init = 7.; // Larger for the first iteration (you are misaligned !!!)
    const double g_start_chi_cut = 500.; // How to optimize this value By guoaq-19-09-26
 
    const int gNsamLC = 100;
    const int gNsamGB = 100;
    ///*     const double gStepLC[5] = {0.00001, 0.000001, 0.00001, 0.0001, 0.00001};   /\* units of dr&dz are cm *\/ */
    const double gStepLC[5] = {0.05, 0.001, 0.001, 0.05, 0.002};    /* units of dr&dz are cm, probably not !! */
    const double gStepGB[36] = {
           /* L1      L2      L3   L4     L5    L6   */
        0.05, 0.05, 0.05,      0.05, 0.05, 0.05,     /* dx */
        0.05, 0.05, 0.05,      0.05, 0.05, 0.05,     /* dy */
        0.05, 0.05, 0.05,      0.05, 0.05, 0.05,     /* dz */ 
        0.001, 0.001, 0.001,   0.001, 0.001, 0.001,  /* rx */
        0.001, 0.001, 0.001,   0.001, 0.001, 0.001,  /* ry */
        0.001, 0.001, 0.001,   0.001, 0.001, 0.001   /* rz */
        };
 
}
 
#endif/* ALIGNMENT_H */