DedxCalibParameters.h

Go to the documentation of this file.

#include <vector>
 
const long int RUN0 = 0;
const long int RUN1 = 0;
const long int RUN2 = 5459;
const long int RUN3 = 8093;
const long int RUN4 = 9947;
const long int RUN5 = 11397;
const long int RUN6 = 23455;
const long int RUN7 = 25244;
 
////runNO<5459, HV:2200
const double HV1_curvep0    = 500.0 * 0.1165;
const double HV1_curvep1    = 9.404;
const double HV1_curvep2    = 0.009811;
const double HV1_curvep3    = 1.621;
const double HV1_curvep4    = 1.19;
const double HV1_index_sin  = 0;
const double HV1_index_nhit = 0.9525;
const double HV1_sigmap0    = 65;
const double HV1_sigmap1    = 57.08;
const double HV1_sigmap2    = 60.03;
const double HV1_sigmap3    = -2.747;
 
////runNO:5960-8093(6201)
const double HV2_curvep0    = 550.0 * 0.0182604;
const double HV2_curvep1    = 35.0428;
const double HV2_curvep2    = 1.24185e-09;
const double HV2_curvep3    = 2.25831;
const double HV2_curvep4    = 5.81573;
const double HV2_index_sin  = 0.341267;
const double HV2_index_nhit = 0.557516;
const double HV2_sigmap0    = 33.0;
const double HV2_sigmap1    = 33.1897;
const double HV2_sigmap2    = 29.3308;
const double HV2_sigmap3    = -2.11828;
 
const double VR0CUT     = 1.0;
const double VZ0CUT     = 10.0;
const double PT0HighCut = 5.0;
const double PT0LowCut  = 0.05;
 
#define NormalMean 550 // reference value
 
#define NumHistBins 300 // set histogram range for gauss and landau distribution
#define MinHistValue 10
#define MaxHistValue 3000
#define MinHistValue1 300
#define MaxHistValue1 950
#define MinChiValue -1
#define MaxChiValue 1
 
#define MinADCValuecut 10 // choose hits again in calibration
#define MaxADCValuecut 35000
#define RPhi_PathMaxCut 3
#define Iner_RPhi_PathMinCut 0.24
#define Out_RPhi_PathMinCut 0.32
#define Iner_DriftDistCut 0.75
#define Out_DriftDistCut 1.0
 
#define NumSlices 100
#define NumSlicesDoca 40
#define NumSlicesEAng 40
#define DocaMin -1.1
#define DocaMax 1.1
#define PhiMin -0.25
#define PhiMax 0.25
 
#define NumBinCostheta 80
#define CosthetaMin -1.00
#define CosthetaMax 1.00
 
#define NumBinCalT0 35
#define NumBinTzero 135
#define TzeroMin 150
#define TzeroMax 1500
 
static const double doca_norm[43] = {
    0.585, 0.615, 0.63,  0.61, 0.59,  0.57,  0.565, 0.61,  0.8,   0.855, 0.8,
    0.855, 0.805, 0.845, 0.8,  0.845, 0.785, 0.825, 0.785, 0.815, 0.765, 0.795,
    0.82,  0.855, 0.8,   0.82, 0.86,  0.895, 0.77,  0.805, 0.83,  0.87,  0.75,
    0.78,  0.8,   0.82,  0.8,  0.815, 0.83,  0.87,  0.77,  0.795, 0.8 };
 
static const double Eangle_value_cut[15] = { -0.25, -0.06,  -0.04, -0.03, -0.02,
                                             -0.01, -0.005, 0,     0.005, 0.01,
                                             0.02,  0.03,   0.04,  0.06,  0.25 };
 
static const double Eangle_cut_bin[14] = { 1, 1, 2, 6, 6, 2, 2, 2, 2, 6, 6, 2, 1, 1 };
 
static const double tzero_bin_value[35] = {
    275,  575,  610,  630,  650,  670,  690,  710,  730,  750,  770,  790,
    810,  830,  850,  870,  890,  910,  930,  950,  970,  990,  1010, 1030,
    1050, 1070, 1090, 1110, 1130, 1150, 1170, 1190, 1225, 1275, 1400 };
 
double mylan( double* x, double* par );
double landaun( double* x, double* par );
double Landau( double* x, double* par );
double Vavilov( double* x, double* par );
double AsymGauss( double* x, double* par );
 
void dedx_pid_exp_old( int landau, int runflag, float dedx, int Nohit, float mom, float theta,
                       float t0, float lsamp, double dedx_exp[5], double ex_sigma[5],
                       double pid_prob[5], double chi_dedx[5] );
 
void dedx_pid_exp( int vflag[3], float dedx, int trkalg, int Nohit, float mom, float theta,
                   float t0, float lsamp, double dedx_exp[5], double ex_sigma[5],
                   double pid_prob[5], double chi_dedx[5], std::vector<double>& par,
                   std::vector<double>& sig_par );
 
double SpaceChargeCorrec( double m_theta, double mom, int Particle, double dEdx );
 
#define NumEvtBlocksInRun 30   // the maximum number of blocks of events in a run
#define NumEvtInBlocks 1000000 // the number of events in a block
#define MaxNumEvt 1000000000   // the maximum number of events in a run