MdcxParameters.cxx

Go to the documentation of this file.

#include "MdcxReco/MdcxParameters.h"
 
/// constants
const double MdcxParameters::c     = 29.9792458;   //cm/ns
 
/// Track attribute
const double MdcxParameters::maxTrkOmega  = 0.2;  //FIXME origin 0.08
const double MdcxParameters::maxTrkLength = 250.0;
double MdcxParameters::maxProca     = 0.6;   //OK
double MdcxParameters::minTrkProb   = 0.01;  //FIXME
 
const double MdcxParameters::maxDd0InMerge   = 2.7;    //OK
const double MdcxParameters::maxDphi0InMerge = 0.15;   //OK
const double MdcxParameters::maxPdradInMerge = 0.22;   //OK
const double MdcxParameters::maxRcsInMerge   = 18.0;   //OK
double MdcxParameters::maxRcsInAddSeg  = 50.0;   //FIXME
 
const double MdcxParameters::maxAddHitCut = 5.;   //FIXME 
double MdcxParameters::helixFitSigma= 5.;   //FIXME 3, original
 
double MdcxParameters::nSigAddHitTrk = 5.;
double MdcxParameters::addHitFactor = 10.;//yzhang 2010-05-05 
double MdcxParameters::csmax4 = 50.;   //FIXME 28, original
double MdcxParameters::csmax3 = 1.;   //FIXME 1, original
 
/// MDC Geometry
const double MdcxParameters::maxMdcRadius = 81.0;
const double MdcxParameters::firstMdcAxialRadius= 18.9;//yzhang 2010-05-05 
const double MdcxParameters::maxMdcZLen   = 170.;//115.4 original
 
//MdcxFindTracks
const double MdcxParameters::maxDlen = 150.;//FIXME
 
////N sigma of drop hits in dropHits();
float  MdcxParameters::dropHitsSigma[43] = {
  //yzhang put nSigma of layer 3,4 = 10
  10.,5.,5.,10.,// 0,1,2,3
  10.,5.,5.,10.,// 4,5,6,7
  10.,5.,5.,5.,// 8,9,10,11
  5.,5.,5.,5.,// 12,13,14,15
  5.,5.,5.,10.,// 16,17,18,19
  10.,5.,5.,5.,// 20,21,22,23
  5.,5.,5.,5.,// 24,25,26,27
  5.,5.,5.,5.,// 28,29,30,31
  5.,5.,5.,10.,// 32,33,34,35
  10.,5.,5.,5.,// 36,37,38,39
  5.,5.,10.     // 40,41,42
};
 
///-- relative to MdcxFindTracks
const int MdcxParameters::findTrkGroup[nSegCombo][3] = {
  { 4, 5, 6 },
  { 2, 1, 0 },
  { 3, 1, 0 },//yzhang add 2011-10-17 
  { 4, 1, 0 },//yzhang add 2011-10-17 
  { 9, 7, 8 },
  { 3, 5, 6 },
  { 10, 7, 8 },
  { 3, 5, 1 },
  { 4, 7, 8 },
  { 9, 5, 6 },
  { 3, 0, 1 },
  { 4, 5, 1 },
  { 2, 5, 1 }// from 2,5,1 to 3,5,1
  /*
  //yzhang bring 9,10 slayer combination forward,
  //for Pt consideration, 2010-09-14 
  { 9, 7, 8 },
  { 10, 7, 8 },
  { 9, 5, 6 },
  { 4, 5, 6 },
  { 2, 1, 0 },
  { 3, 5, 6 },
  { 3, 5, 1 },
  { 4, 7, 8 },
  { 3, 0, 1 },
  { 4, 5, 1 },
  { 2, 5, 1 }// from 2,5,1 to 3,5,1
  */
};
 
const float MdcxParameters::maxDp12[nSegCombo] = {
  //0:0.19, 1:0.22, 2:0.20, 3:0.19, 4:0.22, 5:0.19, 6:0.22, 7:0.24, 8:0.23, 9:0.18, 10:0.21
  0.5, 0.5,0.5, 0.5,0.5, 0.5,0.5, 0.5,0.5, 0.5, 0.5, 0.5, 0.5
  //0.15, 0.22, 0.20, 0.15, 0.22, 0.15, 0.22, 0.24, 0.23, 0.18, 0.21
};
 
const float MdcxParameters::maxDp13[nSegCombo] = {
  //0.33, 0.30, 0.25, 0.35, 0.25, 0.24, 0.35, 0.29, 0.30, 0.39, 0.25
  0.5, 0.5,0.5, 0.5,0.5, 0.5,0.5, 0.5,0.5, 0.5, 0.5, 0.5, 0.5
  //0.33, 0.39, 0.25, 0.35, 0.25, 0.24, 0.35, 0.29, 0.30, 0.39, 0.25
};
 
const int MdcxParameters::layerSet2AddSeg[nSegCombo][11] = {
  { 7, 8, 9, 10, 3, 2, 1, 0, 6, 5, 4 },
  { 3, 4, 5, 6, 7, 8, 9, 10, 2, 1, 0 },
  { 4, 2, 5, 6, 7, 8, 9, 10, 3, 1, 0 },
  { 3, 2, 5, 6, 7, 8, 9, 10, 4, 1, 0 },
  { 10, 6, 5, 4, 3, 2, 1, 0, 9, 8, 7 },
  { 4, 7, 8, 9, 10, 2, 1, 0, 6, 5, 3 },
  { 9, 6, 5, 4, 3, 2, 1, 0, 10, 8, 7 },
  { 4, 2, 6, 7, 8, 9, 10, 0, 5, 3, 1 },
  { 6, 5, 9, 10, 3, 2, 1, 0, 8, 7, 4 },
  { 7, 8, 10, 4, 3, 2, 1, 0, 9, 6, 5 },
  { 2, 4, 5, 6, 7, 8, 9, 10, 3, 1, 0 },
  { 3, 2, 6, 7, 8, 9, 10, 0, 5, 4, 1 },
  { 3, 4, 6, 7, 8, 9, 10, 0, 5, 2, 1 }
};
 
int MdcxParameters::debug = 0;