#include <CSCXonStrip_MatchGatti.h>

Public Member Functions
	CSCXonStrip_MatchGatti (const edm::ParameterSet &ps)

void	findXOnStrip (const CSCDetId &id, const CSCLayer *layer, const CSCStripHit &stripHit, int centralStrip, float &xWithinChamber, float &stripWidth, const float &tpeak, float &xWithinStrip, float &sigma, int &quality_flag)
	Returns fitted local x position and its estimated error. More...

void	initChamberSpecs ()
	Use specs to setup Gatti parameters. More...

void	setConditions (const CSCRecoConditions *reco)
	Cache pointer to conditions data. More...

void	setupMatrix ()
	Set matrix for XT corrections and noise. More...

	~CSCXonStrip_MatchGatti ()

Private Member Functions
double	calculateXonStripError (float stripWidth, bool ME1_1)

double	calculateXonStripPosition (float stripWidth, bool ME1_1)

	CSCXonStrip_MatchGatti (const CSCXonStrip_MatchGatti &)

double	estimated2Gatti (double Xestimated, float StripWidth, bool ME1_1)

double	estimated2GattiCorrection (double Xestimated, float StripWidth, bool ME1_1)

void	getCorrectionValues (std::string Estimator)

void	hardcodedCorrectionInitialization ()

CSCXonStrip_MatchGatti &	operator= (const CSCXonStrip_MatchGatti &)

double	xfError_Noise (double noise)

double	xfError_XTasym (double XTasym)

Private Attributes
float	a11 [3]
	Store elements of auto-correlation matrices: 0 = left, 1 = middle, 2 = right. More...

float	a12 [3]

float	a13 [3]

float	a22 [3]

float	a23 [3]

float	a33 [3]

float	chargeSignal [3][3]

float	const_syst

float	const_syst_ME12

float	const_syst_ME13

float	const_syst_ME1a

float	const_syst_ME1b

float	const_syst_ME21

float	const_syst_ME22

float	const_syst_ME31

float	const_syst_ME32

float	const_syst_ME41

double	h

double	k_1

double	k_2

double	k_3

float	noise_level

float	noise_level_ME12

float	noise_level_ME13

float	noise_level_ME1a

float	noise_level_ME1b

float	noise_level_ME21

float	noise_level_ME22

float	noise_level_ME31

float	noise_level_ME32

float	noise_level_ME41

double	norm

double	q_sum

double	q_sumC

double	q_sumL

double	q_sumR

double	r

const CSCRecoConditions *	recoConditions_

const CSCChamberSpecs *	specs_

double	sqrt_k_3

float	stripWidth

bool	use3TimeBins

bool	useCalib

float	v11 [3]
	Store elements of matrices for chi^2 computation: 0 = left, 1 = middle, 2 = right. More...

float	v12 [3]

float	v13 [3]

float	v22 [3]

float	v23 [3]

float	v33 [3]

float	x_centralVal [n_val]

float	x_correction_ME1_1 [n_SW_ME1_1][n_val]

float	x_correction_noME1_1 [n_SW_noME1_1][n_val]

float	xt_asymmetry

float	xt_asymmetry_ME12

float	xt_asymmetry_ME13

float	xt_asymmetry_ME1a

float	xt_asymmetry_ME1b

float	xt_asymmetry_ME21

float	xt_asymmetry_ME22

float	xt_asymmetry_ME31

float	xt_asymmetry_ME32

float	xt_asymmetry_ME41

float	xt_l [3][3]
	x-talks 0 = left, 1 = middle, 2 = right ; and then second [] is for time bin tmax-1, tmax, tmax+1 More...

float	xt_lr0 [3]

float	xt_lr1 [3]

float	xt_lr2 [3]

float	xt_r [3][3]

float	xtalksOffset

Static Private Attributes
static const int	n_SW_ME1_1 = 6

static const int	n_SW_noME1_1 = 11

static const int	n_val = 501

Detailed Description

When having both strip and wire hit in a layer, use Gatti "matching" to calculate position ond error of strip hit.

Author: Stoyan Stoynev - NU

Definition at line 24 of file CSCXonStrip_MatchGatti.h.

Constructor & Destructor Documentation

◆ CSCXonStrip_MatchGatti() [1/2]

CSCXonStrip_MatchGatti::CSCXonStrip_MatchGatti ( const edm::ParameterSet & ps )

explicit

Definition at line 35 of file CSCXonStrip_MatchGatti.cc.

                                                                         : recoConditions_(nullptr) {
   useCalib = ps.getParameter<bool>("CSCUseCalibrations");
   xtalksOffset = ps.getParameter<double>("CSCStripxtalksOffset");
   noise_level_ME1a = ps.getParameter<double>("NoiseLevel_ME1a");
   xt_asymmetry_ME1a = ps.getParameter<double>("XTasymmetry_ME1a");
   const_syst_ME1a = ps.getParameter<double>("ConstSyst_ME1a");
   noise_level_ME1b = ps.getParameter<double>("NoiseLevel_ME1b");
   xt_asymmetry_ME1b = ps.getParameter<double>("XTasymmetry_ME1b");
   const_syst_ME1b = ps.getParameter<double>("ConstSyst_ME1b");
   noise_level_ME12 = ps.getParameter<double>("NoiseLevel_ME12");
   xt_asymmetry_ME12 = ps.getParameter<double>("XTasymmetry_ME12");
   const_syst_ME12 = ps.getParameter<double>("ConstSyst_ME12");
   noise_level_ME13 = ps.getParameter<double>("NoiseLevel_ME13");
   xt_asymmetry_ME13 = ps.getParameter<double>("XTasymmetry_ME13");
   const_syst_ME13 = ps.getParameter<double>("ConstSyst_ME13");
   noise_level_ME21 = ps.getParameter<double>("NoiseLevel_ME21");
   xt_asymmetry_ME21 = ps.getParameter<double>("XTasymmetry_ME21");
   const_syst_ME21 = ps.getParameter<double>("ConstSyst_ME21");
   noise_level_ME22 = ps.getParameter<double>("NoiseLevel_ME22");
   xt_asymmetry_ME22 = ps.getParameter<double>("XTasymmetry_ME22");
   const_syst_ME22 = ps.getParameter<double>("ConstSyst_ME22");
   noise_level_ME31 = ps.getParameter<double>("NoiseLevel_ME31");
   xt_asymmetry_ME31 = ps.getParameter<double>("XTasymmetry_ME31");
   const_syst_ME31 = ps.getParameter<double>("ConstSyst_ME31");
   noise_level_ME32 = ps.getParameter<double>("NoiseLevel_ME32");
   xt_asymmetry_ME32 = ps.getParameter<double>("XTasymmetry_ME32");
   const_syst_ME32 = ps.getParameter<double>("ConstSyst_ME32");
   noise_level_ME41 = ps.getParameter<double>("NoiseLevel_ME41");
   xt_asymmetry_ME41 = ps.getParameter<double>("XTasymmetry_ME41");
   const_syst_ME41 = ps.getParameter<double>("ConstSyst_ME41");
  
   getCorrectionValues("StringCurrentlyNotUsed");
 }

References const_syst_ME12, const_syst_ME13, const_syst_ME1a, const_syst_ME1b, const_syst_ME21, const_syst_ME22, const_syst_ME31, const_syst_ME32, const_syst_ME41, getCorrectionValues(), edm::ParameterSet::getParameter(), noise_level_ME12, noise_level_ME13, noise_level_ME1a, noise_level_ME1b, noise_level_ME21, noise_level_ME22, noise_level_ME31, noise_level_ME32, noise_level_ME41, useCalib, xt_asymmetry_ME12, xt_asymmetry_ME13, xt_asymmetry_ME1a, xt_asymmetry_ME1b, xt_asymmetry_ME21, xt_asymmetry_ME22, xt_asymmetry_ME31, xt_asymmetry_ME32, xt_asymmetry_ME41, and xtalksOffset.

◆ ~CSCXonStrip_MatchGatti()

CSCXonStrip_MatchGatti::~CSCXonStrip_MatchGatti ( )

Definition at line 69 of file CSCXonStrip_MatchGatti.cc.

69 {}

◆ CSCXonStrip_MatchGatti() [2/2]

CSCXonStrip_MatchGatti::CSCXonStrip_MatchGatti ( const CSCXonStrip_MatchGatti & )

private

Member Function Documentation

◆ calculateXonStripError()

double CSCXonStrip_MatchGatti::calculateXonStripError	(	float	stripWidth,
		bool	ME1_1
	)

private

Definition at line 622 of file CSCXonStrip_MatchGatti.cc.

                                                                                   {
   double min;
   if (q_sumR > q_sumL) {
     min = q_sumL;
   } else {
     min = q_sumR;
   }
  
   double xf = (q_sumR - q_sumL) / (q_sumC - min) / 2;
   double xf_ErrorNoise = xfError_Noise(noise_level);
   double xf_ErrorXTasym = xfError_XTasym(xt_asymmetry);
   // x_G = x_F + correction_functon(x_F)
   // as these are correlated the error should be simply d(x_G) = |d(x_F)| + [correction_functon(x_F+|d(x_F)|) - correction_functon(x_F)]
   double d_xf = std::sqrt(std::pow(xf_ErrorNoise, 2) + std::pow(xf_ErrorXTasym, 2));
   double d_corr =
       estimated2GattiCorrection(xf + d_xf, stripWidth, ME1_1) - estimated2GattiCorrection(xf, stripWidth, ME1_1);
   double x_shift = d_xf + d_corr;
   //  double x_shift = sqrt( pow( xf_ErrorNoise, 2) + pow( xf_ErrorXTasym, 2)) *
   //(1 + (estimated2GattiCorrection(xf+0.001, stripWidth, ME1_1) -
   //  estimated2GattiCorrection(xf, stripWidth, ME1_1))*1000.);
   double x_error = std::sqrt(std::pow(std::abs(x_shift) * stripWidth, 2) + std::pow(const_syst, 2));
   return x_error;
 }

References funct::abs(), const_syst, estimated2GattiCorrection(), min(), noise_level, funct::pow(), q_sumC, q_sumL, q_sumR, mathSSE::sqrt(), stripWidth, xfError_Noise(), xfError_XTasym(), and xt_asymmetry.

Referenced by findXOnStrip().

◆ calculateXonStripPosition()

double CSCXonStrip_MatchGatti::calculateXonStripPosition	(	float	stripWidth,
		bool	ME1_1
	)

private

Definition at line 646 of file CSCXonStrip_MatchGatti.cc.

                                                                                      {
   double x_estimated = -99.;
   double min;
   if (q_sumR > q_sumL) {
     min = q_sumL;
   } else {
     min = q_sumR;
   }
   //---- This is XF ( X Florida - after the first group that used it)
   x_estimated = (q_sumR - q_sumL) / (q_sumC - min) / 2;
   double x_gatti = estimated2Gatti(x_estimated, stripWidth, ME1_1);
   return x_gatti;
 }

References estimated2Gatti(), min(), q_sumC, q_sumL, q_sumR, and stripWidth.

Referenced by findXOnStrip().

◆ estimated2Gatti()

double CSCXonStrip_MatchGatti::estimated2Gatti	(	double	Xestimated,
		float	StripWidth,
		bool	ME1_1
	)

private

Definition at line 580 of file CSCXonStrip_MatchGatti.cc.

                                                                                                {
   double x_corr = estimated2GattiCorrection(x_estimated, stripWidth, ME1_1);
   double x_gatti = x_estimated + x_corr;
  
   return x_gatti;
 }

References estimated2GattiCorrection(), and stripWidth.

Referenced by calculateXonStripPosition().

◆ estimated2GattiCorrection()

double CSCXonStrip_MatchGatti::estimated2GattiCorrection	(	double	Xestimated,
		float	StripWidth,
		bool	ME1_1
	)

private

Definition at line 502 of file CSCXonStrip_MatchGatti.cc.

                                                                                                          {
   //---- 11 "nominal" strip widths : 0.6 - 1.6 cm; for ME1_1 just 6 "nominal" strip widths : 0.3 - 0.8 cm; see HardCodedCorrectionInitialization()
   //---- Calculate corrections at specific  Xestimated (linear interpolation between points)
   int n_SW;
   int min_SW;
   if (ME1_1) {
     n_SW = n_SW_ME1_1;
     min_SW = 3;  // min SW calculated is 0.3 cm
   } else {
     n_SW = n_SW_noME1_1;
     min_SW = 6;  // min SW calculated is 0.6 cm
   }
   int stripDown = int(10. * stripWidth) - min_SW;  // 0 is at min strip width calculated
   int stripUp = stripDown + 1;
   if (stripUp > n_SW - 1) {
     //---- to be checked...
     //if(stripUp>n_SW){
     //std::cout<<" Is strip width = "<<stripWidth<<" OK?" <<std::endl;
     //}
     stripUp = n_SW - 1;
   }
  
   double half_strip_width = 0.5;
   //const int Nbins = 501;
   const int n_bins = n_val;
   double corr_2_xestim = 999.;
   if (stripDown < 0) {
     corr_2_xestim = 1;
   } else {
     //---- Parametrized x_gatti minus x_estimated differences
  
     int xestim_bin = -999;
     double delta_strip_width = 999.;
     double delta_strip_widthUpDown = 999.;
     double diff_2_strip_width = 999.;
     delta_strip_width = stripWidth - int(stripWidth * 10) / 10.;
     delta_strip_widthUpDown = 0.1;
  
     if (std::abs(x_estimated) > 0.5) {
       if (std::abs(x_estimated) > 1.) {
         corr_2_xestim = 1.;  // for now; to be investigated
       } else {
         //if(std::abs(Xestimated)>0.55){
         //std::cout<<"X position from the estimated position above 0.55 (safty margin)?! "<<std::endl;
         //CorrToXc = 999.;
         //}
         xestim_bin = int((1. - std::abs(x_estimated)) / half_strip_width * n_bins);
         if (ME1_1) {
           diff_2_strip_width = x_correction_ME1_1[stripUp][xestim_bin] - x_correction_ME1_1[stripDown][xestim_bin];
           corr_2_xestim = x_correction_ME1_1[stripDown][xestim_bin] +
                           (delta_strip_width / delta_strip_widthUpDown) * diff_2_strip_width;
         } else {
           diff_2_strip_width = x_correction_noME1_1[stripUp][xestim_bin] - x_correction_noME1_1[stripDown][xestim_bin];
           corr_2_xestim = x_correction_noME1_1[stripDown][xestim_bin] +
                           (delta_strip_width / delta_strip_widthUpDown) * diff_2_strip_width;
         }
         corr_2_xestim = -corr_2_xestim;
       }
     } else {
       xestim_bin = int((std::abs(x_estimated) / half_strip_width) * n_bins);
       if (ME1_1) {
         diff_2_strip_width = x_correction_ME1_1[stripUp][xestim_bin] - x_correction_ME1_1[stripDown][xestim_bin];
         corr_2_xestim = x_correction_ME1_1[stripDown][xestim_bin] +
                         (delta_strip_width / delta_strip_widthUpDown) * diff_2_strip_width;
       } else {
         diff_2_strip_width = x_correction_noME1_1[stripUp][xestim_bin] - x_correction_noME1_1[stripDown][xestim_bin];
         corr_2_xestim = x_correction_noME1_1[stripDown][xestim_bin] +
                         (delta_strip_width / delta_strip_widthUpDown) * diff_2_strip_width;
       }
     }
     if (x_estimated < 0.) {
       corr_2_xestim = -corr_2_xestim;
     }
   }
  
   return corr_2_xestim;
 }

References funct::abs(), createfilelist::int, n_SW_ME1_1, n_SW_noME1_1, n_val, stripWidth, x_correction_ME1_1, and x_correction_noME1_1.

Referenced by calculateXonStripError(), and estimated2Gatti().

◆ findXOnStrip()

void CSCXonStrip_MatchGatti::findXOnStrip	(	const CSCDetId &	id,
		const CSCLayer *	layer,
		const CSCStripHit &	stripHit,
		int	centralStrip,
		float &	xWithinChamber,
		float &	stripWidth,
		const float &	tpeak,
		float &	xWithinStrip,
		float &	sigma,
		int &	quality_flag
	)

Returns fitted local x position and its estimated error.

Definition at line 74 of file CSCXonStrip_MatchGatti.cc.

                                                              {
   quality_flag = 0;
   // Initialize Gatti parameters using chamberSpecs
   // Cache specs_ info for ease of access
   specs_ = layer->chamber()->specs();
   stripWidth = sWidth;
   //initChamberSpecs();
   // Initialize output parameters
   xWithinStrip = xWithinChamber;
  
   CSCStripHit::ChannelContainer const& strips = stripHit.strips();
   int nStrips = strips.size();
   int centStrip = nStrips / 2 + 1;
   std::vector<float> const& adcs = stripHit.s_adc();
   int tmax = stripHit.tmax();
  
   //float t_peak = tpeak;
   //float t_zero = 0.;
   //float adc[4];
   //
   //if ( useCalib ) {
   //
   //  for ( int t = 0; t < 4; ++t ) {
   //    int k  = t + 4 * (centStrip-1);
   //    adc[t] = adcs[k];
   //  }
   //
   //  // t_peak from peak finder is now 'absolute' i.e. in ns from start of sca time bin 0
   //  t_peak = peakTimeFinder_->peakTime( tmax, adc, t_peak );
   //  // Just for completeness, the start time of the pulse is 133 ns earlier, according to Stan :)
   //  t_zero = t_peak - 133.;
   //  // and reset tpeak since that's the way it gets passed out of this function (Argh!)
   //  tpeak = t_peak;
   //  LogTrace("CSCRecHit|CSCXonStrip_MatchGatti") << "CSCXonStrip_MatchGatti: " <<
   //     id << " strip=" << centralStrip << ", t_zero=" << t_zero << ", t_peak=" << t_peak;
   //}
  
   //---- fill the charge matrix (3x3)
   float adc[4];
   int j = 0;
   for (int i = 1; i <= nStrips; ++i) {
     if (i > (centStrip - 2) && i < (centStrip + 2)) {
       std::vector<float> adcsFit;
       for (int t = 0; t < 4; ++t) {
         int k = t + 4 * (i - 1);
         adc[t] = adcs[k];
         if (t < 3)
           chargeSignal[j][t] = adc[t];
       }
       j++;
     }
   }
  
   // Load in x-talks:
  
   if (useCalib) {
     std::vector<float> xtalks;
     recoConditions_->crossTalk(id, centralStrip, xtalks);
     float dt = 50.f * tmax - tpeak;
     // XTalks; l,r are for left, right XTalk; lr0,1,2 are for what charge "remains" in the strip
     for (int t = 0; t < 3; ++t) {
       xt_l[0][t] = xtalks[0] * (50.f * (t - 1) + dt) + xtalks[1] + xtalksOffset;
       xt_r[0][t] = xtalks[2] * (50.f * (t - 1) + dt) + xtalks[3] + xtalksOffset;
       xt_l[1][t] = xtalks[4] * (50.f * (t - 1) + dt) + xtalks[5] + xtalksOffset;
       xt_r[1][t] = xtalks[6] * (50.f * (t - 1) + dt) + xtalks[7] + xtalksOffset;
       xt_l[2][t] = xtalks[8] * (50.f * (t - 1) + dt) + xtalks[9] + xtalksOffset;
       xt_r[2][t] = xtalks[10] * (50.f * (t - 1) + dt) + xtalks[11] + xtalksOffset;
  
       xt_lr0[t] = (1. - xt_l[0][t] - xt_r[0][t]);
       xt_lr1[t] = (1. - xt_l[1][t] - xt_r[1][t]);
       xt_lr2[t] = (1. - xt_l[2][t] - xt_r[2][t]);
     }
   } else {
     for (int t = 0; t < 3; ++t) {
       xt_l[0][t] = xtalksOffset;
       xt_r[0][t] = xtalksOffset;
       xt_l[1][t] = xtalksOffset;
       xt_r[1][t] = xtalksOffset;
       xt_l[2][t] = xtalksOffset;
       xt_r[2][t] = xtalksOffset;
  
       xt_lr0[t] = (1. - xt_l[0][t] - xt_r[0][t]);
       xt_lr1[t] = (1. - xt_l[1][t] - xt_r[1][t]);
       xt_lr2[t] = (1. - xt_l[2][t] - xt_r[2][t]);
     }
   }
  
   // vector containing noise starts at tmax - 1, and tmax > 3, but....
   int tbin = tmax - 4;
  
   // .... originally, suppose to have tmax in tbin 4 or 5, but noticed in MTCC lots of
   // hits with tmax == 3, so let's allow these, and shift down noise matrix by one element...
   // This is a patch because the calibration database doesn't have elements for tbin = 2,
   // e.g. there is no element e[tmax-1,tmax+1] = e[2,4].
  
   if (tmax < 4)
     tbin = 0;  // patch
  
   // Load in auto-correlation noise matrices
   if (useCalib) {
     std::vector<float> nmatrix;
     recoConditions_->noiseMatrix(id, centralStrip, nmatrix);
     for (int istrip = 0; istrip < 3; ++istrip) {
       a11[istrip] = nmatrix[0 + tbin * 3 + istrip * 15];
       a12[istrip] = nmatrix[1 + tbin * 3 + istrip * 15];
       a13[istrip] = nmatrix[2 + tbin * 3 + istrip * 15];
       a22[istrip] = nmatrix[3 + tbin * 3 + istrip * 15];
       a23[istrip] = nmatrix[4 + tbin * 3 + istrip * 15];
       a33[istrip] = nmatrix[6 + tbin * 3 + istrip * 15];
     }
   } else {
     // FIXME:  NO HARD WIRED VALUES !!!
     for (int istrip = 0; istrip < 3; ++istrip) {
       a11[istrip] = 10.0;
       a12[istrip] = 0.0;
       a13[istrip] = 0.0;
       a22[istrip] = 10.0;
       a23[istrip] = 0.0;
       a33[istrip] = 10.0;
     }
   }
  
   //---- Set up noise, XTalk matrices
   setupMatrix();
   //---- Calculate the coordinate within the strip and associate uncertainty
  
   static const std::string ME1a("ME1/a");
   static const std::string ME1b("ME1/b");
  
   bool ME1_1 = (ME1a == specs_->chamberTypeName() || ME1b == specs_->chamberTypeName());
  
   // due to cross talks and 3 time bin sum it is in principe possible that the center strip is not the maximum strip
   // in some cases the consequences could be quite extreme
   // take some measures against the extreme cases
   bool peakMismatch = false;
   std::vector<float> charges(3);
   charges[0] = q_sumL;
   charges[1] = q_sumC;
   charges[2] = q_sumR;
   int min_index = min_element(charges.begin(), charges.end()) - charges.begin();
   int max_index = max_element(charges.begin(), charges.end()) - charges.begin();
   if (1 != max_index && (1 == min_index ||
                          // the condition below means that if the initial position estimate within strip (|xF|)
                          // is above  1.1/2 = 0.55 "strip widths" peakMismatch is set to true (so special case);
                          // in normal cases |xF|<=0.5 (0.5 is at the edge of the "central" strip)
                          (charges[max_index] - charges[min_index]) > 1.1 * (q_sumC - charges[min_index]))) {
     peakMismatch = true;
     switch (max_index) {
       case 0:
         xWithinStrip = -1;
         break;
       case 2:
         xWithinStrip = 1;
         break;
       default:
         // should be an error message here
         xWithinStrip = 0;  // in case?
         break;
     }
   }
   // we don't have the needed information (it's similar to the "edge" strip case)
   //else if(stripHit.isNearDeadStrip()){
   else if (stripHit.deadStrip() > 0) {
     xWithinStrip = 0;
   } else {
     //
     xWithinStrip = float(calculateXonStripPosition(stripWidth, ME1_1));
   }
   xWithinChamber = xWithinChamber + (xWithinStrip * stripWidth);
   if (peakMismatch) {
     sigma = stripWidth / std::sqrt(12.f);
   } else {
     //---- error estimation
     int factorStripWidth = int(std::sqrt(stripWidth / 0.38f));
     int maxConsecutiveStrips = 8;
     if (factorStripWidth) {
       maxConsecutiveStrips /= factorStripWidth;
     }
     maxConsecutiveStrips++;
  
     struct ChamberTypes {
       std::map<std::string, int> chamberTypes;
       int operator()(std::string const& s) const {
         auto p = chamberTypes.find(s);
         return p != chamberTypes.end() ? (*p).second : 0;
       }
       ChamberTypes() {
         chamberTypes["ME1/a"] = 1;
         chamberTypes["ME1/b"] = 2;
         chamberTypes["ME1/2"] = 3;
         chamberTypes["ME1/3"] = 4;
         chamberTypes["ME2/1"] = 5;
         chamberTypes["ME2/2"] = 6;
         chamberTypes["ME3/1"] = 7;
         chamberTypes["ME3/2"] = 8;
         chamberTypes["ME4/1"] = 9;
         chamberTypes["ME4/2"] = 8;
       }
     };
     static const ChamberTypes chamberTypes;
  
     switch (chamberTypes(specs_->chamberTypeName())) {
       case 1:
         noise_level = noise_level_ME1a;
         xt_asymmetry = xt_asymmetry_ME1a;
         const_syst = const_syst_ME1a;
         break;
  
       case 2:
         noise_level = noise_level_ME1b;
         xt_asymmetry = xt_asymmetry_ME1b;
         const_syst = const_syst_ME1b;
         break;
  
       case 3:
         noise_level = noise_level_ME12;
         xt_asymmetry = xt_asymmetry_ME12;
         const_syst = const_syst_ME12;
         break;
  
       case 4:
         noise_level = noise_level_ME13;
         xt_asymmetry = xt_asymmetry_ME13;
         const_syst = const_syst_ME13;
         break;
  
       case 5:
         noise_level = noise_level_ME21;
         xt_asymmetry = xt_asymmetry_ME21;
         const_syst = const_syst_ME21;
         break;
  
       case 6:
         noise_level = noise_level_ME22;
         xt_asymmetry = xt_asymmetry_ME22;
         const_syst = const_syst_ME22;
         break;
  
       case 7:
         noise_level = noise_level_ME31;
         xt_asymmetry = xt_asymmetry_ME31;
         const_syst = const_syst_ME31;
         break;
  
       case 8:
         noise_level = noise_level_ME32;
         xt_asymmetry = xt_asymmetry_ME32;
         const_syst = const_syst_ME32;
         break;
  
       case 9:
         noise_level = noise_level_ME41;
         xt_asymmetry = xt_asymmetry_ME41;
         const_syst = const_syst_ME41;
         break;
  
       default:
         noise_level = noise_level_ME22;
         xt_asymmetry = xt_asymmetry_ME22;
         const_syst = const_syst_ME22;
     }
     if (0 == stripHit.deadStrip() && stripHit.numberOfConsecutiveStrips() < maxConsecutiveStrips &&
         std::abs(stripHit.closestMaximum()) > maxConsecutiveStrips / 2) {
       sigma = float(calculateXonStripError(stripWidth, ME1_1));
     } else {  //---- near dead strip or too close maxima or too wide strip cluster
       sigma = stripWidth / std::sqrt(12.f);
     }
   }
   quality_flag = 1;
 }

References a11, a12, a13, a22, a23, a33, funct::abs(), gpuClustering::adc, calculateXonStripError(), calculateXonStripPosition(), CSCChamberSpecs::chamberTypeName(), CosmicGenFilterHelix_cfi::charges, chargeSignal, CSCStripHit::closestMaximum(), const_syst, const_syst_ME12, const_syst_ME13, const_syst_ME1a, const_syst_ME1b, const_syst_ME21, const_syst_ME22, const_syst_ME31, const_syst_ME32, const_syst_ME41, CSCRecoConditions::crossTalk(), CSCStripHit::deadStrip(), dt, f, dqmMemoryStats::float, mps_fire::i, createfilelist::int, dqmiolumiharvest::j, dqmdumpme::k, phase1PixelTopology::layer, hgcal_clustering::max_index(), HLT_FULL_cff::ME1a, noise_level, noise_level_ME12, noise_level_ME13, noise_level_ME1a, noise_level_ME1b, noise_level_ME21, noise_level_ME22, noise_level_ME31, noise_level_ME32, noise_level_ME41, CSCRecoConditions::noiseMatrix(), me0TriggerPseudoDigis_cff::nStrips, CSCStripHit::numberOfConsecutiveStrips(), AlCaHLTBitMon_ParallelJobs::p, q_sumC, q_sumL, q_sumR, recoConditions_, alignCSCRings::s, CSCStripHit::s_adc(), setupMatrix(), specs_, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, DigiDM_cff::strips, CSCStripHit::strips(), stripWidth, submitPVValidationJobs::t, tmax, CSCStripHit::tmax(), useCalib, xt_asymmetry, xt_asymmetry_ME12, xt_asymmetry_ME13, xt_asymmetry_ME1a, xt_asymmetry_ME1b, xt_asymmetry_ME21, xt_asymmetry_ME22, xt_asymmetry_ME31, xt_asymmetry_ME32, xt_asymmetry_ME41, xt_l, xt_lr0, xt_lr1, xt_lr2, xt_r, and xtalksOffset.

Referenced by CSCMake2DRecHit::hitFromStripAndWire().

◆ getCorrectionValues()

void CSCXonStrip_MatchGatti::getCorrectionValues ( std::string Estimator )

private

Definition at line 500 of file CSCXonStrip_MatchGatti.cc.

500 { hardcodedCorrectionInitialization(); }

References hardcodedCorrectionInitialization().

Referenced by CSCXonStrip_MatchGatti().

◆ hardcodedCorrectionInitialization()

void CSCXonStrip_MatchGatti::hardcodedCorrectionInitialization ( )

private

Definition at line 11 of file HardCodedCorrectionInitialization.cc.

                                                                {
   //---- Corrections depend on strip width (SW), position within the strip (XFlorida is used) and type of chamber
  
   //---- Corrections for non-ME1_1 chambers (all have practically equal Gatti parameters):
  
   //---- min SW ~ 6.6 cm (ME1/2), max SW ~ 1.6 cm (ME234/2)
   //---- Corrections are calculated for N_SW_noME1_1 strip widths (0.6 - 1.6 cm)
   //---- and for 501 values in
   float corrections_noME1_1[n_SW_noME1_1][n_val] = {
       {0,           0.00019078,  0.000383855,
        0.000579201, 0.000776796, 0.000976616,
        0.00117891,  0.00138342,  0.00159008,
        0.00179887,  0.00200977,  0.00222304,
        0.00243835,  0.00265569,  0.00287502,
        0.00309646,  0.00332001,  0.00354547,
        0.00377282,  0.00400203,  0.00423336,
        0.0044665,   0.00470141,  0.00493807,
        0.00517666,  0.005417,    0.00565901,
        0.00590266,  0.00614806,  0.00639514,
        0.00664377,  0.00689391,  0.00714568,
        0.00739901,  0.00765376,  0.0079099,
        0.00816756,  0.00842663,  0.00868701,
        0.00894865,  0.00921171,  0.00947602,
        0.0097415,   0.0100081,   0.0102761,
        0.0105451,   0.0108152,   0.0110863,
        0.0113586,   0.0116318,   0.011906,
        0.0121811,   0.0124572,   0.0127341,
        0.0130118,   0.0132904,   0.0135697,
        0.0138497,   0.0141305,   0.014412,
        0.014694,    0.0149767,   0.01526,
        0.0155438,   0.015828,    0.0161128,
        0.016398,    0.0166836,   0.0169696,
        0.0172559,   0.0175426,   0.0178295,
        0.0181166,   0.018404,    0.0186915,
        0.0189791,   0.0192669,   0.0195548,
        0.0198427,   0.0201306,   0.0204185,
        0.0207063,   0.0209941,   0.0212817,
        0.0215692,   0.0218565,   0.0221436,
        0.0224305,   0.0227171,   0.0230034,
        0.0232894,   0.023575,    0.0238602,
        0.024145,    0.0244293,   0.0247132,
        0.0249966,   0.0252793,   0.0255616,
        0.0258432,   0.0261243,   0.0264046,
        0.0266843,   0.0269633,   0.0272415,
        0.027519,    0.0277957,   0.0280716,
        0.0283466,   0.0286207,   0.028894,
        0.0291663,   0.0294377,   0.0297081,
        0.0299776,   0.030246,    0.0305133,
        0.0307796,   0.0310449,   0.0313089,
        0.0315718,   0.0318336,   0.0320942,
        0.0323536,   0.0326117,   0.0328687,
        0.0331243,   0.0333786,   0.0336316,
        0.0338833,   0.0341336,   0.0343825,
        0.0346301,   0.0348763,   0.035121,
        0.0353642,   0.035606,    0.0358463,
        0.0360851,   0.0363223,   0.036558,
        0.0367921,   0.0370248,   0.0372557,
        0.037485,    0.0377128,   0.0379389,
        0.0381633,   0.038386,    0.0386071,
        0.0388265,   0.0390442,   0.0392601,
        0.0394742,   0.0396866,   0.0398974,
        0.0401062,   0.0403132,   0.0405184,
        0.0407219,   0.0409236,   0.0411233,
        0.0413211,   0.0415171,   0.0417113,
        0.0419036,   0.042094,    0.0422824,
        0.0424689,   0.0426535,   0.0428362,
        0.043017,    0.0431956,   0.0433724,
        0.0435472,   0.04372,     0.043891,
        0.0440598,   0.0442266,   0.0443913,
        0.0445541,   0.0447149,   0.0448738,
        0.0450305,   0.0451851,   0.0453376,
        0.0454882,   0.0456367,   0.0457831,
        0.0459275,   0.0460698,   0.0462099,
        0.0463479,   0.0464839,   0.0466178,
        0.0467496,   0.0468794,   0.0470071,
        0.0471324,   0.0472557,   0.0473769,
        0.047496,    0.047613,    0.0477279,
        0.0478407,   0.0479514,   0.0480598,
        0.0481661,   0.0482702,   0.0483723,
        0.0484722,   0.0485699,   0.0486656,
        0.0487591,   0.0488505,   0.0489397,
        0.0490268,   0.0491116,   0.0491943,
        0.0492749,   0.0493533,   0.0494296,
        0.0495037,   0.0495757,   0.0496456,
        0.0497133,   0.0497789,   0.0498423,
        0.0499037,   0.0499628,   0.0500198,
        0.0500746,   0.0501273,   0.0501778,
        0.0502262,   0.0502724,   0.0503165,
        0.0503585,   0.0503984,   0.0504361,
        0.0504717,   0.0505052,   0.0505366,
        0.0505658,   0.050593,    0.050618,
        0.0506409,   0.0506617,   0.0506803,
        0.0506969,   0.0507114,   0.0507237,
        0.050734,    0.0507422,   0.0507482,
        0.0507522,   0.0507541,   0.0507539,
        0.0507517,   0.0507473,   0.0507409,
        0.0507324,   0.0507219,   0.0507093,
        0.0506946,   0.0506779,   0.0506591,
        0.0506382,   0.0506154,   0.0505905,
        0.0505635,   0.0505345,   0.0505035,
        0.0504705,   0.0504355,   0.0503984,
        0.0503594,   0.0503183,   0.0502753,
        0.0502302,   0.0501832,   0.0501342,
        0.0500832,   0.0500302,   0.0499752,
        0.0499182,   0.0498592,   0.0497983,
        0.0497354,   0.0496705,   0.0496037,
        0.049535,    0.0494644,   0.0493918,
        0.0493173,   0.0492409,   0.0491626,
        0.0490824,   0.0490003,   0.0489162,
        0.0488301,   0.0487422,   0.0486524,
        0.0485607,   0.0484671,   0.0483717,
        0.0482745,   0.0481754,   0.0480744,
        0.0479716,   0.0478668,   0.0477601,
        0.0476517,   0.0475414,   0.0474293,
        0.0473154,   0.0471997,   0.0470823,
        0.046963,    0.0468417,   0.0467186,
        0.0465938,   0.0464672,   0.0463389,
        0.0462088,   0.046077,    0.0459434,
        0.0458078,   0.0456705,   0.0455315,
        0.0453908,   0.0452484,   0.0451043,
        0.0449585,   0.0448107,   0.0446612,
        0.04451,     0.0443572,   0.0442028,
        0.0440467,   0.0438887,   0.0437289,
        0.0435675,   0.0434045,   0.0432399,
        0.0430737,   0.0429056,   0.0427358,
        0.0425643,   0.0423913,   0.0422167,
        0.0420406,   0.0418625,   0.0416828,
        0.0415016,   0.0413188,   0.0411344,
        0.0409484,   0.0407606,   0.0405712,
        0.0403803,   0.040188,    0.0399941,
        0.0397982,   0.0396009,   0.0394021,
        0.0392018,   0.039,       0.0387963,
        0.0385912,   0.0383846,   0.0381765,
        0.037967,    0.0377555,   0.0375427,
        0.0373284,   0.0371128,   0.0368954,
        0.0366764,   0.0364561,   0.0362343,
        0.0360112,   0.0357861,   0.0355597,
        0.0353319,   0.0351028,   0.034872,
        0.0346396,   0.0344059,   0.0341709,
        0.0339343,   0.0336961,   0.0334565,
        0.0332157,   0.0329735,   0.0327294,
        0.0324841,   0.0322375,   0.0319897,
        0.0317399,   0.031489,    0.0312367,
        0.0309832,   0.0307279,   0.0304713,
        0.0302135,   0.0299544,   0.0296935,
        0.0294315,   0.0291682,   0.0289035,
        0.0286372,   0.0283698,   0.0281011,
        0.0278309,   0.0275592,   0.0272864,
        0.0270125,   0.0267367,   0.0264598,
        0.0261818,   0.0259024,   0.0256214,
        0.0253393,   0.0250561,   0.0247713,
        0.0244852,   0.024198,    0.0239095,
        0.0236194,   0.0233283,   0.0230361,
        0.0227422,   0.0224472,   0.0221511,
        0.0218536,   0.0215547,   0.0212548,
        0.0209538,   0.0206511,   0.0203474,
        0.0200426,   0.0197363,   0.0194288,
        0.0191203,   0.0188104,   0.0184992,
        0.0181871,   0.0178736,   0.0175588,
        0.017243,    0.016926,    0.0166075,
        0.0162881,   0.0159677,   0.0156456,
        0.0153226,   0.0149987,   0.014673,
        0.0143465,   0.0140191,   0.01369,
        0.01336,     0.0130291,   0.0126966,
        0.0123631,   0.0120288,   0.0116929,
        0.011356,    0.0110183,   0.0106789,
        0.0103387,   0.00999765,  0.00965491,
        0.00931136,  0.00896692,  0.00862087,
        0.00827404,  0.00792623,  0.00757693,
        0.00722686,  0.00687569,  0.00652318,
        0.0061699,   0.0058154,   0.00545971,
        0.00510328,  0.00474545,  0.00438662,
        0.00402707,  0.00366593,  0.00330401,
        0.00294123,  0.00257695,  0.00221198,
        0.00184591,  0.0014786,   0.00111061,
        0.000741268, 0.000370968, 0},
  
       {0,          0.00041187, 0.000826272, 0.00124349,  0.00166337,  0.00208575, 0.00251089, 0.00293832, 0.00336842,
        0.00380078, 0.00423553, 0.00467258,  0.00511179,  0.0055533,   0.00599674, 0.00644249, 0.00689004, 0.00733967,
        0.0077911,  0.00824439, 0.00869946,  0.00915618,  0.00961462,  0.0100746,  0.0105361,  0.010999,   0.0114635,
        0.0119291,  0.0123962,  0.0128644,   0.0133338,   0.0138042,   0.0142757,  0.0147481,  0.0152215,  0.0156957,
        0.0161706,  0.0166463,  0.0171227,   0.0175996,   0.018077,    0.018555,   0.0190333,  0.019512,   0.0199909,
        0.02047,    0.0209493,  0.0214287,   0.0219081,   0.0223875,   0.0228668,  0.023346,   0.0238249,  0.0243035,
        0.0247819,  0.0252598,  0.0257373,   0.0262142,   0.0266906,   0.0271664,  0.0276415,  0.0281158,  0.0285894,
        0.0290621,  0.0295338,  0.0300047,   0.0304745,   0.0309433,   0.0314108,  0.0318773,  0.0323424,  0.0328064,
        0.033269,   0.0337302,  0.03419,     0.0346482,   0.035105,    0.0355602,  0.0360138,  0.0364656,  0.0369159,
        0.0373642,  0.0378108,  0.0382556,   0.0386984,   0.0391394,   0.0395783,  0.0400153,  0.0404501,  0.0408829,
        0.0413135,  0.041742,   0.0421683,   0.0425922,   0.043014,    0.0434333,  0.0438504,  0.044265,   0.0446772,
        0.0450871,  0.0454943,  0.0458991,   0.0463013,   0.0467009,   0.047098,   0.0474923,  0.0478841,  0.0482731,
        0.0486594,  0.0490431,  0.0494237,   0.0498018,   0.050177,    0.0505493,  0.050919,   0.0512855,  0.0516492,
        0.0520102,  0.0523678,  0.0527228,   0.0530748,   0.0534235,   0.0537695,  0.0541124,  0.0544521,  0.0547889,
        0.0551226,  0.0554531,  0.0557806,   0.056105,    0.0564261,   0.0567442,  0.0570593,  0.0573708,  0.0576794,
        0.057985,   0.058287,   0.0585859,   0.0588818,   0.0591743,   0.0594634,  0.0597496,  0.0600325,  0.0603119,
        0.0605881,  0.0608613,  0.0611311,   0.0613974,   0.0616606,   0.0619207,  0.0621773,  0.0624305,  0.0626806,
        0.0629275,  0.0631709,  0.063411,    0.0636478,   0.0638816,   0.0641118,  0.0643387,  0.0645623,  0.0647828,
        0.065,      0.0652136,  0.065424,    0.0656312,   0.0658352,   0.0660359,  0.0662331,  0.066427,   0.0666178,
        0.0668053,  0.0669897,  0.0671705,   0.067348,    0.0675224,   0.0676935,  0.0678615,  0.0680262,  0.0681875,
        0.0683455,  0.0685004,  0.0686521,   0.0688006,   0.068946,    0.0690879,  0.0692266,  0.0693622,  0.0694946,
        0.0696238,  0.06975,    0.0698729,   0.0699928,   0.0701093,   0.0702226,  0.0703328,  0.0704399,  0.0705439,
        0.0706449,  0.0707427,  0.0708375,   0.0709291,   0.0710177,   0.0711031,  0.0711853,  0.0712646,  0.0713408,
        0.071414,   0.0714842,  0.0715514,   0.0716155,   0.0716767,   0.0717349,  0.0717901,  0.0718423,  0.0718916,
        0.0719379,  0.0719812,  0.0720216,   0.072059,    0.0720935,   0.0721251,  0.0721537,  0.0721795,  0.0722025,
        0.0722225,  0.0722397,  0.0722541,   0.0722656,   0.0722743,   0.0722801,  0.0722832,  0.0722835,  0.0722809,
        0.0722756,  0.0722675,  0.0722567,   0.0722431,   0.0722268,   0.0722078,  0.072186,   0.0721615,  0.0721344,
        0.0721045,  0.072072,   0.0720369,   0.071999,    0.0719584,   0.0719152,  0.0718694,  0.0718209,  0.0717699,
        0.0717162,  0.07166,    0.0716012,   0.0715399,   0.071476,    0.0714096,  0.0713406,  0.0712692,  0.0711952,
        0.0711188,  0.0710399,  0.0709584,   0.0708743,   0.0707878,   0.0706988,  0.0706074,  0.0705137,  0.0704175,
        0.0703189,  0.070218,   0.0701147,   0.070009,    0.0699007,   0.0697901,  0.0696771,  0.0695618,  0.0694443,
        0.0693244,  0.0692023,  0.0690779,   0.0689511,   0.0688218,   0.0686903,  0.0685566,  0.0684207,  0.0682826,
        0.0681423,  0.0679999,  0.0678548,   0.0677076,   0.0675583,   0.0674068,  0.0672532,  0.0670976,  0.0669396,
        0.0667793,  0.0666169,  0.0664524,   0.0662859,   0.0661173,   0.0659466,  0.0657735,  0.0655984,  0.0654213,
        0.0652422,  0.0650612,  0.0648778,   0.0646923,   0.0645048,   0.0643155,  0.0641242,  0.0639308,  0.0637351,
        0.0635376,  0.0633381,  0.0631368,   0.0629335,   0.0627279,   0.0625205,  0.0623112,  0.0621002,  0.061887,
        0.0616717,  0.0614547,  0.0612359,   0.0610153,   0.0607924,   0.0605677,  0.0603412,  0.0601131,  0.0598829,
        0.0596506,  0.0594167,  0.0591811,   0.0589437,   0.0587041,   0.0584628,  0.0582199,  0.0579754,  0.0577285,
        0.05748,    0.05723,    0.0569783,   0.0567244,   0.0564688,   0.0562118,  0.0559531,  0.0556921,  0.0554297,
        0.0551658,  0.0549,     0.0546323,   0.0543631,   0.0540924,   0.0538197,  0.0535453,  0.0532695,  0.0529922,
        0.0527126,  0.0524317,  0.0521493,   0.0518652,   0.0515792,   0.0512918,  0.0510031,  0.0507122,  0.0504199,
        0.0501263,  0.0498308,  0.0495336,   0.0492352,   0.0489352,   0.0486333,  0.04833,    0.0480255,  0.0477189,
        0.047411,   0.0471018,  0.0467907,   0.0464781,   0.0461643,   0.0458489,  0.0455317,  0.0452134,  0.0448936,
        0.0445719,  0.0442491,  0.0439249,   0.0435988,   0.0432715,   0.042943,   0.0426125,  0.042281,   0.0419481,
        0.0416133,  0.0412775,  0.0409404,   0.0406013,   0.0402613,   0.0399199,  0.0395767,  0.0392324,  0.0388869,
        0.0385395,  0.0381912,  0.0378414,   0.03749,     0.0371377,   0.0367837,  0.0364283,  0.036072,   0.0357139,
        0.0353546,  0.0349943,  0.0346321,   0.0342689,   0.0339046,   0.0335385,  0.0331715,  0.0328031,  0.0324332,
        0.0320626,  0.03169,    0.0313165,   0.0309419,   0.0305656,   0.0301884,  0.0298098,  0.0294298,  0.029049,
        0.0286664,  0.0282829,  0.0278983,   0.027512,    0.027125,    0.0267364,  0.0263467,  0.0259561,  0.0255637,
        0.0251706,  0.0247761,  0.0243804,   0.0239839,   0.0235855,   0.0231864,  0.0227861,  0.0223844,  0.0219821,
        0.0215779,  0.021173,   0.0207668,   0.0203594,   0.0199513,   0.0195414,  0.0191309,  0.0187189,  0.018306,
        0.0178921,  0.0174766,  0.0170606,   0.0166429,   0.0162245,   0.0158049,  0.0153841,  0.0149626,  0.0145394,
        0.0141157,  0.0136904,  0.0132644,   0.0128373,   0.0124088,   0.0119799,  0.0115491,  0.0111179,  0.0106852,
        0.0102517,  0.00981723, 0.00938147,  0.00894511,  0.0085071,   0.00806868, 0.00762866, 0.00718802, 0.00674618,
        0.00630334, 0.00585968, 0.00541467,  0.00496919,  0.00452202,  0.00407449, 0.00362542, 0.00317576, 0.00272489,
        0.00227312, 0.00182046, 0.0013666,   0.000912158, 0.000456222, 0},
  
       {0,          0.000672915, 0.00134894, 0.00202799, 0.00271026, 0.00339538, 0.00408321, 0.00477386, 0.00546709,
        0.00616272, 0.00686075,  0.00756109, 0.00826349, 0.00896789, 0.0096743,  0.0103824,  0.0110922,  0.0118036,
        0.0125164,  0.0132306,   0.0139459,  0.0146623,  0.0153797,  0.0160979,  0.0168169,  0.0175365,  0.0182565,
        0.0189769,  0.0196976,   0.0204185,  0.0211393,  0.0218601,  0.0225807,  0.0233009,  0.0240207,  0.0247399,
        0.0254585,  0.0261763,   0.0268933,  0.0276092,  0.028324,   0.0290376,  0.0297499,  0.0304607,  0.03117,
        0.0318777,  0.0325836,   0.0332877,  0.03399,    0.03469,    0.035388,   0.0360838,  0.0367772,  0.0374682,
        0.0381568,  0.0388427,   0.0395259,  0.0402065,  0.0408841,  0.0415587,  0.0422304,  0.0428991,  0.0435644,
        0.0442266,  0.0448855,   0.0455409,  0.0461929,  0.0468414,  0.0474861,  0.0481274,  0.0487648,  0.0493984,
        0.0500284,  0.0506541,   0.0512761,  0.0518941,  0.0525078,  0.0531176,  0.053723,   0.0543243,  0.0549214,
        0.055514,   0.0561024,   0.0566863,  0.0572658,  0.0578408,  0.0584112,  0.0589773,  0.0595384,  0.0600953,
        0.0606472,  0.0611946,   0.0617372,  0.0622751,  0.0628083,  0.0633366,  0.0638602,  0.0643788,  0.0648927,
        0.0654015,  0.0659058,   0.0664047,  0.0668991,  0.0673882,  0.0678727,  0.0683519,  0.0688264,  0.0692957,
        0.0697602,  0.0702196,   0.0706739,  0.0711234,  0.0715676,  0.0720072,  0.0724412,  0.0728706,  0.0732947,
        0.0737139,  0.0741282,   0.0745371,  0.0749414,  0.0753403,  0.0757344,  0.0761235,  0.0765073,  0.0768865,
        0.0772604,  0.0776293,   0.0779936,  0.0783523,  0.0787064,  0.0790556,  0.0793995,  0.0797387,  0.0800731,
        0.0804022,  0.0807267,   0.0810464,  0.0813607,  0.0816704,  0.0819756,  0.0822754,  0.0825705,  0.082861,
        0.0831465,  0.0834271,   0.0837032,  0.0839746,  0.0842409,  0.0845025,  0.0847597,  0.0850123,  0.0852596,
        0.0855026,  0.085741,    0.085975,   0.0862039,  0.0864283,  0.0866483,  0.0868638,  0.0870748,  0.0872809,
        0.0874827,  0.0876801,   0.0878732,  0.0880617,  0.0882456,  0.0884251,  0.0886004,  0.0887714,  0.0889381,
        0.0891003,  0.089258,    0.0894116,  0.089561,   0.0897062,  0.0898472,  0.0899841,  0.0901165,  0.0902447,
        0.0903688,  0.0904889,   0.0906049,  0.0907168,  0.0908248,  0.0909287,  0.0910286,  0.0911242,  0.0912159,
        0.0913036,  0.0913875,   0.0914674,  0.0915435,  0.0916157,  0.0916841,  0.0917486,  0.0918093,  0.0918662,
        0.0919194,  0.0919688,   0.0920144,  0.0920562,  0.0920943,  0.0921287,  0.0921595,  0.0921866,  0.0922101,
        0.09223,    0.0922464,   0.0922592,  0.0922684,  0.0922741,  0.0922763,  0.092275,   0.0922703,  0.092262,
        0.0922504,  0.0922354,   0.0922169,  0.0921951,  0.09217,    0.0921415,  0.0921096,  0.0920745,  0.0920361,
        0.0919945,  0.0919494,   0.0919012,  0.0918497,  0.091795,   0.0917372,  0.0916762,  0.091612,   0.0915448,
        0.0914745,  0.091401,    0.0913246,  0.0912451,  0.0911625,  0.091077,   0.0909885,  0.0908967,  0.090802,
        0.0907043,  0.0906037,   0.0905003,  0.0903939,  0.0902848,  0.0901727,  0.0900579,  0.0899401,  0.0898192,
        0.0896956,  0.0895692,   0.0894402,  0.0893084,  0.0891739,  0.0890367,  0.0888966,  0.0887536,  0.088608,
        0.0884598,  0.0883091,   0.0881558,  0.0879999,  0.0878409,  0.0876794,  0.0875154,  0.0873489,  0.0871799,
        0.0870085,  0.0868341,   0.0866572,  0.086478,   0.0862963,  0.0861123,  0.0859257,  0.0857363,  0.0855446,
        0.0853506,  0.0851544,   0.0849557,  0.0847542,  0.0845504,  0.0843445,  0.0841363,  0.0839258,  0.0837125,
        0.0834971,  0.0832795,   0.0830598,  0.0828375,  0.0826128,  0.0823861,  0.0821573,  0.0819262,  0.0816925,
        0.0814567,  0.081219,    0.0809793,  0.0807367,  0.0804922,  0.0802458,  0.0799975,  0.0797463,  0.0794933,
        0.0792384,  0.0789816,   0.078722,   0.0784606,  0.0781974,  0.0779322,  0.0776644,  0.0773949,  0.0771236,
        0.0768501,  0.0765743,   0.0762968,  0.0760177,  0.0757359,  0.0754523,  0.0751671,  0.0748799,  0.0745903,
        0.0742992,  0.0740065,   0.0737111,  0.0734141,  0.0731156,  0.0728148,  0.0725121,  0.0722079,  0.0719018,
        0.0715934,  0.0712836,   0.0709722,  0.0706583,  0.0703429,  0.0700262,  0.0697069,  0.0693862,  0.069064,
        0.0687396,  0.0684135,   0.0680861,  0.0677565,  0.0674252,  0.0670926,  0.0667578,  0.0664214,  0.0660838,
        0.0657439,  0.0654025,   0.0650598,  0.0647149,  0.0643685,  0.064021,   0.063671,   0.0633198,  0.0629673,
        0.0626125,  0.0622565,   0.0618991,  0.0615395,  0.0611789,  0.0608165,  0.0604524,  0.0600872,  0.0597199,
        0.0593512,  0.0589815,   0.0586094,  0.0582363,  0.0578617,  0.0574852,  0.0571078,  0.0567284,  0.0563476,
        0.0559658,  0.0555817,   0.0551968,  0.0548102,  0.054422,   0.0540329,  0.0536416,  0.0532493,  0.0528556,
        0.0524601,  0.0520639,   0.0516655,  0.0512661,  0.0508654,  0.050463,   0.0500597,  0.0496544,  0.0492481,
        0.0488405,  0.0484312,   0.0480212,  0.047609,   0.047196,   0.0467816,  0.0463657,  0.0459489,  0.0455301,
        0.0451107,  0.0446894,   0.0442671,  0.0438436,  0.0434185,  0.0429927,  0.0425648,  0.0421363,  0.0417061,
        0.0412748,  0.0408424,   0.0404084,  0.0399738,  0.0395371,  0.0390999,  0.0386609,  0.0382211,  0.0377799,
        0.0373374,  0.0368941,   0.036449,   0.0360035,  0.0355558,  0.0351078,  0.034658,   0.0342073,  0.0337554,
        0.0333023,  0.0328482,   0.0323926,  0.0319364,  0.0314784,  0.03102,    0.0305596,  0.0300988,  0.0296363,
        0.0291731,  0.0287085,   0.028243,   0.0277763,  0.0273084,  0.0268397,  0.0263695,  0.0258987,  0.0254262,
        0.0249533,  0.0244786,   0.0240036,  0.0235267,  0.0230495,  0.0225705,  0.0220912,  0.0216101,  0.0211286,
        0.0206456,  0.0201619,   0.0196768,  0.019191,   0.018704,   0.0182161,  0.017727,   0.017237,   0.016746,
        0.016254,   0.0157609,   0.0152669,  0.0147719,  0.0142758,  0.0137788,  0.0132808,  0.0127819,  0.0122819,
        0.011781,   0.011279,    0.0107762,  0.0102724,  0.00976763, 0.00926186, 0.0087552,  0.00824756, 0.00773898,
        0.0072295,  0.006719,    0.00620769, 0.00569529, 0.00518218, 0.00466788, 0.00415298, 0.0036368,  0.00312014,
        0.00260207, 0.00208368,  0.00156373, 0.00104362, 0.00052181, 0},
  
       {0,          0.000982554, 0.00196886, 0.00295871, 0.00395189,  0.00494819, 0.00594739, 0.00694928, 0.00795364,
        0.00896026, 0.0099689,   0.0109794,  0.0119914,  0.0130048,   0.0140193,  0.0150348,  0.016051,   0.0170676,
        0.0180845,  0.0191014,   0.0201181,  0.0211344,  0.0221501,   0.023165,   0.0241788,  0.0251913,  0.0262023,
        0.0272117,  0.0282192,   0.0292245,  0.0302276,  0.0312282,   0.0322261,  0.0332211,  0.034213,   0.0352017,
        0.036187,   0.0371686,   0.0381465,  0.0391204,  0.0400901,   0.0410556,  0.0420167,  0.0429731,  0.0439248,
        0.0448715,  0.0458132,   0.0467497,  0.0476809,  0.0486067,   0.0495269,  0.0504414,  0.0513501,  0.0522529,
        0.0531496,  0.0540403,   0.0549246,  0.0558025,  0.056674,    0.057539,   0.0583974,  0.0592492,  0.0600943,
        0.0609324,  0.0617635,   0.0625877,  0.0634048,  0.064215,    0.0650181,  0.0658137,  0.066602,   0.0673832,
        0.0681571,  0.0689236,   0.0696824,  0.0704339,  0.0711781,   0.0719148,  0.0726436,  0.0733651,  0.0740791,
        0.0747853,  0.0754839,   0.0761751,  0.0768586,  0.0775342,   0.0782025,  0.0788631,  0.0795158,  0.0801612,
        0.0807988,  0.0814287,   0.0820513,  0.0826658,  0.0832729,   0.0838726,  0.0844644,  0.0850491,  0.0856256,
        0.0861951,  0.0867569,   0.0873111,  0.0878581,  0.0883973,   0.0889296,  0.0894539,  0.0899715,  0.090481,
        0.0909839,  0.0914788,   0.0919672,  0.0924476,  0.0929215,   0.0933876,  0.0938471,  0.0942991,  0.0947443,
        0.0951824,  0.0956134,   0.0960378,  0.0964547,  0.0968652,   0.0972685,  0.0976651,  0.0980553,  0.0984381,
        0.0988146,  0.0991844,   0.0995473,  0.0999041,  0.100254,    0.100597,   0.100935,   0.101265,   0.101589,
        0.101907,   0.102218,    0.102523,   0.102823,   0.103115,    0.103402,   0.103682,   0.103957,   0.104225,
        0.104487,   0.104744,    0.104995,   0.105239,   0.105477,    0.10571,    0.105938,   0.106159,   0.106375,
        0.106585,   0.106789,    0.106988,   0.107182,   0.107369,    0.107552,   0.107729,   0.1079,     0.108067,
        0.108227,   0.108383,    0.108533,   0.108679,   0.108819,    0.108954,   0.109084,   0.109208,   0.109328,
        0.109443,   0.109552,    0.109657,   0.109758,   0.109853,    0.109943,   0.110029,   0.11011,    0.110186,
        0.110257,   0.110324,    0.110386,   0.110443,   0.110497,    0.110545,   0.110589,   0.110629,   0.110664,
        0.110695,   0.110721,    0.110743,   0.110761,   0.110775,    0.110784,   0.110789,   0.11079,    0.110787,
        0.11078,    0.110768,    0.110753,   0.110733,   0.11071,     0.110682,   0.110651,   0.110615,   0.110576,
        0.110533,   0.110486,    0.110435,   0.110381,   0.110322,    0.11026,    0.110195,   0.110125,   0.110052,
        0.109976,   0.109895,    0.109811,   0.109724,   0.109633,    0.109538,   0.10944,    0.109339,   0.109234,
        0.109126,   0.109014,    0.108899,   0.10878,    0.108659,    0.108533,   0.108405,   0.108274,   0.108139,
        0.108001,   0.10786,     0.107715,   0.107568,   0.107417,    0.107263,   0.107107,   0.106947,   0.106784,
        0.106618,   0.106449,    0.106277,   0.106102,   0.105924,    0.105743,   0.105559,   0.105373,   0.105184,
        0.104992,   0.104796,    0.104598,   0.104397,   0.104194,    0.103988,   0.103778,   0.103566,   0.103352,
        0.103135,   0.102915,    0.102692,   0.102467,   0.10224,     0.10201,    0.101776,   0.10154,    0.101302,
        0.101062,   0.100819,    0.100573,   0.100325,   0.100074,    0.099821,   0.0995652,  0.0993072,  0.0990471,
        0.098784,   0.0985186,   0.0982511,  0.0979815,  0.0977087,   0.097434,   0.0971573,  0.096878,   0.0965961,
        0.0963123,  0.0960266,   0.0957379,  0.0954472,  0.0951546,   0.0948595,  0.094562,   0.0942627,  0.0939615,
        0.0936573,  0.0933514,   0.0930437,  0.0927332,  0.0924209,   0.0921068,  0.0917902,  0.0914715,  0.0911512,
        0.0908285,  0.0905036,   0.0901772,  0.0898485,  0.0895175,   0.089185,   0.0888504,  0.0885135,  0.0881751,
        0.0878346,  0.0874919,   0.0871477,  0.0868013,  0.0864529,   0.0861032,  0.0857509,  0.085397,   0.0850417,
        0.0846837,  0.0843243,   0.0839635,  0.0835999,  0.0832352,   0.0828687,  0.0824999,  0.08213,    0.0817578,
        0.0813839,  0.0810089,   0.0806311,  0.0802522,  0.0798717,   0.079489,   0.0791052,  0.0787191,  0.0783315,
        0.0779427,  0.0775514,   0.0771591,  0.0767649,  0.0763689,   0.0759719,  0.0755724,  0.0751719,  0.0747696,
        0.0743655,  0.0739605,   0.0735529,  0.0731446,  0.0727342,   0.0723224,  0.0719094,  0.0714941,  0.0710781,
        0.0706599,  0.0702405,   0.0698197,  0.069397,   0.0689735,   0.0685476,  0.068121,   0.0676924,  0.0672626,
        0.0668314,  0.0663984,   0.0659646,  0.0655285,  0.0650918,   0.064653,   0.0642132,  0.0637719,  0.0633291,
        0.0628852,  0.0624394,   0.0619931,  0.0615443,  0.061095,    0.0606437,  0.0601915,  0.0597378,  0.0592827,
        0.0588265,  0.0583686,   0.05791,    0.0574493,  0.0569882,   0.0565247,  0.0560609,  0.055595,   0.0551285,
        0.0546601,  0.054191,    0.0537202,  0.0532484,  0.0527753,   0.0523008,  0.0518254,  0.0513483,  0.0508705,
        0.0503909,  0.0499108,   0.0494287,  0.0489461,  0.0484616,   0.0479767,  0.0474897,  0.0470025,  0.0465131,
        0.0460235,  0.0455317,   0.0450398,  0.0445458,  0.0440515,   0.0435552,  0.0430586,  0.04256,    0.0420612,
        0.0415603,  0.0410592,   0.0405562,  0.0400528,  0.0395475,   0.039042,   0.0385345,  0.0380268,  0.0375171,
        0.0370073,  0.0364953,   0.0359834,  0.0354693,  0.0349552,   0.034439,   0.0339227,  0.0334045,  0.032886,
        0.0323659,  0.0318452,   0.031323,   0.0308002,  0.0302761,   0.0297512,  0.0292252,  0.028698,   0.0281702,
        0.0276409,  0.0271112,   0.0265798,  0.0260483,  0.0255148,   0.0249813,  0.0244459,  0.0239103,  0.0233731,
        0.0228354,  0.0222964,   0.0217567,  0.021216,   0.0206742,   0.0201318,  0.019588,   0.0190439,  0.018498,
        0.0179521,  0.0174044,   0.0168564,  0.0163071,  0.0157571,   0.0152062,  0.0146542,  0.0141018,  0.0135478,
        0.0129937,  0.0124378,   0.0118817,  0.0113243,  0.0107663,   0.0102073,  0.00964732, 0.00908693, 0.00852496,
        0.007963,   0.00739923,  0.00683529, 0.00627013, 0.00570424,  0.00513772, 0.00456988, 0.00400203, 0.00343225,
        0.00286246, 0.00229137,  0.00171965, 0.00114727, 0.000573636, 0},
  
       {0,          0.00135166, 0.00270752, 0.00406723, 0.00543052,  0.0067968,  0.00816576, 0.00953697, 0.01091,
        0.0122844,  0.0136597,  0.0150357,  0.0164118,  0.0177876,   0.0191627,  0.0205367,  0.0219092,  0.0232797,
        0.024648,   0.0260136,  0.027376,   0.028735,   0.0300902,   0.0314411,  0.0327875,  0.0341289,  0.0354649,
        0.0367955,  0.0381201,  0.0394383,  0.04075,    0.042055,    0.0433525,  0.0446427,  0.0459254,  0.0471998,
        0.048466,   0.0497239,  0.0509732,  0.0522132,  0.0534442,   0.0546659,  0.0558782,  0.0570807,  0.0582731,
        0.0594556,  0.0606279,  0.06179,    0.0629416,  0.0640825,   0.0652126,  0.0663318,  0.0674402,  0.0685376,
        0.0696239,  0.0706991,  0.071763,   0.0728156,  0.0738568,   0.0748866,  0.075905,   0.0769119,  0.0779072,
        0.078891,   0.0798632,  0.0808239,  0.081773,   0.0827105,   0.0836365,  0.0845509,  0.0854539,  0.0863453,
        0.0872253,  0.0880939,  0.0889509,  0.0897963,  0.0906304,   0.0914533,  0.092265,   0.0930655,  0.0938546,
        0.0946324,  0.0953993,  0.0961553,  0.0969003,  0.097634,    0.0983571,  0.0990697,  0.099771,   0.100462,
        0.101142,   0.101812,   0.102471,   0.10312,    0.103759,    0.104387,   0.105005,   0.105613,   0.106211,
        0.106799,   0.107377,   0.107945,   0.108504,   0.109053,    0.109592,   0.110122,   0.110642,   0.111153,
        0.111654,   0.112147,   0.11263,    0.113105,   0.11357,     0.114027,   0.114474,   0.114913,   0.115343,
        0.115765,   0.116178,   0.116582,   0.116979,   0.117366,    0.117746,   0.118118,   0.118481,   0.118836,
        0.119184,   0.119523,   0.119855,   0.120179,   0.120495,    0.120804,   0.121106,   0.121399,   0.121686,
        0.121965,   0.122237,   0.122501,   0.122759,   0.12301,     0.123253,   0.123489,   0.123719,   0.123942,
        0.124158,   0.124367,   0.12457,    0.124766,   0.124956,    0.125139,   0.125316,   0.125486,   0.125651,
        0.125809,   0.125961,   0.126107,   0.126246,   0.12638,     0.126507,   0.126629,   0.126745,   0.126856,
        0.126961,   0.127059,   0.127152,   0.12724,    0.127322,    0.127398,   0.127469,   0.127535,   0.127595,
        0.12765,    0.1277,     0.127745,   0.127784,   0.127819,    0.127848,   0.127872,   0.127892,   0.127906,
        0.127915,   0.12792,    0.12792,    0.127915,   0.127905,    0.12789,    0.127871,   0.127847,   0.127819,
        0.127786,   0.127749,   0.127707,   0.127661,   0.12761,     0.127555,   0.127496,   0.127432,   0.127364,
        0.127292,   0.127216,   0.127136,   0.127052,   0.126963,    0.12687,    0.126774,   0.126673,   0.126569,
        0.12646,    0.126348,   0.126232,   0.126112,   0.125989,    0.125861,   0.12573,    0.125594,   0.125456,
        0.125314,   0.125168,   0.125019,   0.124866,   0.124709,    0.124549,   0.124385,   0.124219,   0.124049,
        0.123875,   0.123697,   0.123517,   0.123333,   0.123146,    0.122956,   0.122762,   0.122565,   0.122365,
        0.122162,   0.121956,   0.121746,   0.121534,   0.121318,    0.1211,     0.120878,   0.120653,   0.120425,
        0.120195,   0.119962,   0.119725,   0.119485,   0.119243,    0.118998,   0.11875,    0.118499,   0.118245,
        0.117989,   0.11773,    0.117468,   0.117203,   0.116936,    0.116666,   0.116393,   0.116118,   0.11584,
        0.115559,   0.115276,   0.114991,   0.114702,   0.114411,    0.114118,   0.113822,   0.113523,   0.113222,
        0.112919,   0.112613,   0.112305,   0.111995,   0.111682,    0.111366,   0.111049,   0.110729,   0.110406,
        0.110082,   0.109755,   0.109425,   0.109094,   0.10876,     0.108424,   0.108086,   0.107746,   0.107403,
        0.107059,   0.106711,   0.106362,   0.106012,   0.105658,    0.105302,   0.104945,   0.104585,   0.104223,
        0.10386,    0.103494,   0.103126,   0.102756,   0.102384,    0.102011,   0.101634,   0.101257,   0.100877,
        0.100495,   0.100112,   0.0997257,  0.0993382,  0.0989492,   0.0985574,  0.0981646,  0.0977694,  0.0973724,
        0.0969742,  0.0965731,  0.096171,   0.0957667,  0.0953606,   0.0949533,  0.0945433,  0.0941323,  0.093719,
        0.0933042,  0.0928879,  0.0924693,  0.0920498,  0.0916276,   0.0912044,  0.0907793,  0.0903524,  0.0899244,
        0.0894938,  0.0890626,  0.0886288,  0.088194,   0.0877573,   0.087319,   0.0868794,  0.0864376,  0.0859952,
        0.08555,    0.0851043,  0.0846562,  0.0842071,  0.0837563,   0.0833039,  0.0828502,  0.0823946,  0.0819382,
        0.0814793,  0.08102,    0.0805581,  0.0800957,  0.079631,    0.0791655,  0.0786981,  0.0782296,  0.0777595,
        0.0772879,  0.0768152,  0.0763406,  0.0758652,  0.0753877,   0.0749096,  0.0744293,  0.0739485,  0.0734653,
        0.0729819,  0.0724959,  0.0720098,  0.0715212,  0.0710323,   0.0705411,  0.0700496,  0.0695557,  0.0690615,
        0.0685651,  0.0680684,  0.0675694,  0.06707,    0.0665685,   0.0660666,  0.0655625,  0.0650582,  0.0645516,
        0.0640448,  0.0635357,  0.0630264,  0.0625148,  0.0620032,   0.0614891,  0.0609751,  0.0604586,  0.0599421,
        0.0594234,  0.0589044,  0.0583834,  0.0578619,  0.0573387,   0.0568148,  0.0562895,  0.0557632,  0.0552357,
        0.0547069,  0.0541773,  0.0536462,  0.0531145,  0.052581,    0.0520472,  0.0515114,  0.0509755,  0.0504374,
        0.0498992,  0.0493591,  0.0488185,  0.0482765,  0.0477336,   0.0471897,  0.0466445,  0.0460987,  0.0455512,
        0.0450036,  0.0444538,  0.0439039,  0.0433523,  0.0428001,   0.0422467,  0.0416923,  0.0411372,  0.0405805,
        0.0400237,  0.0394648,  0.0389059,  0.0383452,  0.037784,    0.0372217,  0.0366583,  0.0360944,  0.0355288,
        0.0349632,  0.0343955,  0.0338277,  0.0332586,  0.0326885,   0.0321179,  0.0315457,  0.0309735,  0.0303993,
        0.0298249,  0.0292493,  0.0286727,  0.0280957,  0.027517,    0.0269383,  0.0263578,  0.025777,   0.0251952,
        0.0246122,  0.0240291,  0.023444,   0.022859,   0.0222725,   0.0216853,  0.0210977,  0.0205084,  0.0199191,
        0.0193282,  0.0187368,  0.0181447,  0.0175513,  0.0169579,   0.0163626,  0.0157671,  0.0151708,  0.0145732,
        0.0139757,  0.0133763,  0.0127767,  0.0121762,  0.0115746,   0.010973,   0.0103695,  0.00976587, 0.00916142,
        0.00855578, 0.00795015, 0.00734277, 0.00673513, 0.00612685,  0.00551723, 0.00490761, 0.00429648, 0.00368488,
        0.00307291, 0.00245935, 0.00184578, 0.00123104, 0.000615521, 0},
  
       {0,          0.00179287, 0.00359052, 0.00539216, 0.00719704, 0.00900437, 0.0108134,  0.0126233,  0.0144332,
        0.0162425,  0.0180503,  0.0198558,  0.0216583,  0.023457,   0.0252512,  0.0270402,  0.0288232,  0.0305995,
        0.0323685,  0.0341296,  0.035882,   0.037625,   0.0393583,  0.0410813,  0.0427931,  0.0444934,  0.0461819,
        0.0478578,  0.0495207,  0.0511703,  0.0528061,  0.0544276,  0.0560349,  0.0576269,  0.059204,   0.0607656,
        0.0623113,  0.0638413,  0.0653548,  0.0668521,  0.0683326,  0.0697965,  0.0712434,  0.0726733,  0.0740862,
        0.0754815,  0.0768599,  0.0782206,  0.0795639,  0.0808901,  0.0821985,  0.0834894,  0.084763,   0.0860194,
        0.0872577,  0.0884789,  0.0896829,  0.0908698,  0.0920394,  0.0931915,  0.0943266,  0.0954448,  0.0965462,
        0.0976309,  0.0986989,  0.0997503,  0.100785,   0.101804,   0.102806,   0.103793,   0.104763,   0.105717,
        0.106656,   0.10758,    0.108488,   0.10938,    0.110258,   0.111121,   0.111968,   0.112801,   0.113619,
        0.114424,   0.115214,   0.11599,    0.116751,   0.117499,   0.118233,   0.118954,   0.119661,   0.120355,
        0.121037,   0.121705,   0.12236,    0.123003,   0.123633,   0.124251,   0.124856,   0.12545,    0.126032,
        0.126601,   0.127159,   0.127705,   0.12824,    0.128764,   0.129277,   0.129778,   0.130269,   0.130749,
        0.131218,   0.131677,   0.132125,   0.132563,   0.132991,   0.133408,   0.133816,   0.134214,   0.134602,
        0.134982,   0.13535,    0.13571,    0.136061,   0.136402,   0.136734,   0.137058,   0.137372,   0.137678,
        0.137975,   0.138263,   0.138543,   0.138815,   0.139078,   0.139333,   0.13958,    0.139819,   0.140051,
        0.140274,   0.140489,   0.140697,   0.140898,   0.141091,   0.141276,   0.141454,   0.141625,   0.141789,
        0.141946,   0.142095,   0.142238,   0.142373,   0.142502,   0.142625,   0.14274,    0.14285,    0.142952,
        0.143049,   0.143138,   0.143222,   0.143299,   0.14337,    0.143435,   0.143493,   0.143546,   0.143593,
        0.143635,   0.14367,    0.1437,     0.143723,   0.143742,   0.143754,   0.143762,   0.143763,   0.14376,
        0.143751,   0.143736,   0.143717,   0.143692,   0.143662,   0.143627,   0.143587,   0.143542,   0.143492,
        0.143437,   0.143377,   0.143313,   0.143243,   0.143169,   0.143091,   0.143008,   0.142919,   0.142827,
        0.14273,    0.142628,   0.142522,   0.142412,   0.142297,   0.142179,   0.142056,   0.141928,   0.141797,
        0.141661,   0.141521,   0.141377,   0.141229,   0.141078,   0.140922,   0.140762,   0.140598,   0.14043,
        0.140259,   0.140084,   0.139905,   0.139722,   0.139536,   0.139346,   0.139153,   0.138956,   0.138754,
        0.13855,    0.138342,   0.138131,   0.137916,   0.137698,   0.137476,   0.137251,   0.137023,   0.136791,
        0.136556,   0.136318,   0.136077,   0.135832,   0.135584,   0.135333,   0.13508,    0.134822,   0.134562,
        0.134299,   0.134033,   0.133763,   0.133491,   0.133216,   0.132938,   0.132656,   0.132372,   0.132086,
        0.131796,   0.131503,   0.131208,   0.13091,    0.130609,   0.130305,   0.129999,   0.12969,    0.129378,
        0.129064,   0.128747,   0.128427,   0.128105,   0.12778,    0.127453,   0.127123,   0.126791,   0.126456,
        0.126119,   0.125779,   0.125436,   0.125092,   0.124744,   0.124395,   0.124043,   0.123689,   0.123332,
        0.122974,   0.122612,   0.122248,   0.121883,   0.121514,   0.121144,   0.120771,   0.120396,   0.12002,
        0.11964,    0.119259,   0.118876,   0.11849,    0.118102,   0.117712,   0.11732,    0.116926,   0.11653,
        0.116132,   0.115731,   0.115329,   0.114925,   0.114518,   0.114111,   0.1137,     0.113288,   0.112874,
        0.112457,   0.11204,    0.111619,   0.111198,   0.110774,   0.110348,   0.109921,   0.109491,   0.109061,
        0.108627,   0.108193,   0.107756,   0.107317,   0.106877,   0.106435,   0.105992,   0.105546,   0.105099,
        0.10465,    0.104199,   0.103747,   0.103293,   0.102837,   0.102379,   0.101921,   0.10146,    0.100998,
        0.100533,   0.100068,   0.0996003,  0.0991318,  0.098661,   0.0981892,  0.0977154,  0.0972402,  0.0967635,
        0.0962851,  0.0958054,  0.0953239,  0.0948412,  0.0943565,  0.0938709,  0.0933832,  0.0928947,  0.0924039,
        0.0919125,  0.0914187,  0.0909245,  0.0904278,  0.0899306,  0.089431,   0.0889311,  0.0884286,  0.0879258,
        0.0874206,  0.086915,   0.0864071,  0.0858986,  0.085388,   0.0848768,  0.0843635,  0.0838496,  0.0833337,
        0.082817,   0.0822986,  0.0817791,  0.0812582,  0.080736,   0.0802126,  0.0796877,  0.0791619,  0.0786344,
        0.0781062,  0.0775759,  0.0770454,  0.0765125,  0.0759796,  0.0754442,  0.0749087,  0.074371,   0.0738328,
        0.073293,   0.0727522,  0.0722102,  0.0716668,  0.0711227,  0.0705768,  0.0700306,  0.0694821,  0.0689335,
        0.0683828,  0.0678317,  0.067279,   0.0667254,  0.0661708,  0.0656147,  0.0650582,  0.0644995,  0.0639408,
        0.06338,    0.0628188,  0.0622562,  0.0616926,  0.0611282,  0.0605621,  0.059996,   0.0594274,  0.0588589,
        0.0582886,  0.0577176,  0.0571458,  0.0565724,  0.0559989,  0.0554231,  0.0548473,  0.0542699,  0.0536916,
        0.0531128,  0.0525321,  0.0519515,  0.0513688,  0.0507858,  0.0502017,  0.0496163,  0.0490308,  0.0484431,
        0.0478554,  0.0472662,  0.0466762,  0.0460857,  0.0454933,  0.044901,   0.0443069,  0.0437123,  0.043117,
        0.0425201,  0.0419231,  0.0413244,  0.0407252,  0.0401253,  0.0395238,  0.0389223,  0.0383191,  0.0377153,
        0.037111,   0.036505,   0.0358991,  0.0352915,  0.0346833,  0.0340747,  0.0334643,  0.032854,   0.0322422,
        0.0316296,  0.0310169,  0.0304022,  0.0297874,  0.0291717,  0.0285547,  0.0279378,  0.0273191,  0.0267,
        0.0260804,  0.0254592,  0.0248379,  0.0242154,  0.0235921,  0.0229687,  0.0223434,  0.0217179,  0.0210918,
        0.0204642,  0.0198366,  0.0192078,  0.0185781,  0.0179484,  0.0173169,  0.0166851,  0.016053,   0.0154192,
        0.0147853,  0.0141506,  0.0135147,  0.0128787,  0.0122414,  0.0116035,  0.0109655,  0.0103257,  0.00968566,
        0.00904539, 0.00840336, 0.00776132, 0.00711861, 0.00647456, 0.00583051, 0.00518541, 0.00453936, 0.00389331,
        0.00324586, 0.00259782, 0.00194978, 0.00130004, 0.00065002, 0},
  
       {0,          0.00232214, 0.00464938, 0.00698026, 0.00931337,  0.0116471,  0.0139801,  0.0163108,  0.0186377,
        0.0209595,  0.0232746,  0.0255817,  0.0278796,  0.0301668,   0.0324419,  0.034704,   0.036952,   0.0391842,
        0.0414,     0.0435984,  0.0457784,  0.0479391,  0.0500795,   0.0521989,  0.0542966,  0.056372,   0.0584245,
        0.0604535,  0.0624585,  0.0644391,  0.0663949,  0.0683255,   0.0702308,  0.0721103,  0.0739639,  0.0757912,
        0.0775922,  0.0793672,  0.0811159,  0.082838,   0.0845335,   0.0862031,  0.0878462,  0.0894628,  0.091054,
        0.0926184,  0.0941575,  0.0956705,  0.0971585,  0.0986207,   0.100058,   0.10147,    0.102858,   0.104222,
        0.105561,   0.106876,   0.108168,   0.109436,   0.110681,    0.111903,   0.113103,   0.11428,    0.115435,
        0.116568,   0.11768,    0.118771,   0.119841,   0.120891,    0.12192,    0.122929,   0.123918,   0.124887,
        0.125838,   0.126769,   0.127682,   0.128576,   0.129452,    0.13031,    0.131151,   0.131975,   0.13278,
        0.133569,   0.134342,   0.135099,   0.135839,   0.136562,    0.137271,   0.137964,   0.138642,   0.139304,
        0.139953,   0.140586,   0.141205,   0.14181,    0.142401,    0.142978,   0.143541,   0.144092,   0.144629,
        0.145153,   0.145665,   0.146164,   0.14665,    0.147124,    0.147586,   0.148037,   0.148474,   0.148902,
        0.149317,   0.149721,   0.150114,   0.150496,   0.150868,    0.151228,   0.151578,   0.151919,   0.152248,
        0.152568,   0.152878,   0.153177,   0.153467,   0.153748,    0.15402,    0.154281,   0.154534,   0.154777,
        0.155013,   0.155238,   0.155455,   0.155664,   0.155864,    0.156056,   0.156239,   0.156414,   0.156581,
        0.156741,   0.156893,   0.157036,   0.157172,   0.1573,      0.157421,   0.157534,   0.157641,   0.15774,
        0.157832,   0.157917,   0.157995,   0.158066,   0.15813,     0.158187,   0.158238,   0.158282,   0.15832,
        0.158351,   0.158376,   0.158395,   0.158407,   0.158414,    0.158414,   0.158409,   0.158397,   0.15838,
        0.158357,   0.158328,   0.158293,   0.158253,   0.158208,    0.158157,   0.1581,     0.158038,   0.157971,
        0.157899,   0.157821,   0.157738,   0.15765,    0.157557,    0.157459,   0.157356,   0.157249,   0.157136,
        0.15702,    0.156897,   0.15677,    0.156639,   0.156503,    0.156363,   0.156218,   0.156069,   0.155915,
        0.155757,   0.155594,   0.155427,   0.155257,   0.155082,    0.154903,   0.154719,   0.154531,   0.15434,
        0.154145,   0.153945,   0.153742,   0.153534,   0.153323,    0.153109,   0.15289,    0.152667,   0.152441,
        0.152212,   0.151979,   0.151741,   0.1515,     0.151256,    0.151009,   0.150758,   0.150503,   0.150245,
        0.149984,   0.149719,   0.149451,   0.14918,    0.148906,    0.148627,   0.148346,   0.148063,   0.147775,
        0.147484,   0.147191,   0.146894,   0.146594,   0.146291,    0.145986,   0.145677,   0.145365,   0.145051,
        0.144733,   0.144413,   0.14409,    0.143764,   0.143434,    0.143103,   0.142769,   0.142431,   0.142092,
        0.141749,   0.141404,   0.141056,   0.140705,   0.140352,    0.139997,   0.139638,   0.139277,   0.138914,
        0.138548,   0.13818,    0.137809,   0.137435,   0.13706,     0.136682,   0.136301,   0.135919,   0.135533,
        0.135146,   0.134756,   0.134363,   0.133969,   0.133572,    0.133173,   0.132772,   0.132368,   0.131963,
        0.131554,   0.131144,   0.130732,   0.130317,   0.129901,    0.129482,   0.129062,   0.128638,   0.128214,
        0.127787,   0.127358,   0.126927,   0.126493,   0.126059,    0.125621,   0.125183,   0.124742,   0.124299,
        0.123854,   0.123407,   0.122959,   0.122508,   0.122056,    0.121601,   0.121145,   0.120687,   0.120227,
        0.119766,   0.119302,   0.118837,   0.11837,    0.117902,    0.117431,   0.116959,   0.116484,   0.116009,
        0.115531,   0.115053,   0.114571,   0.114089,   0.113605,    0.113119,   0.112632,   0.112143,   0.111652,
        0.11116,    0.110666,   0.11017,    0.109673,   0.109174,    0.108674,   0.108172,   0.107669,   0.107164,
        0.106658,   0.10615,    0.10564,    0.10513,    0.104617,    0.104103,   0.103588,   0.103071,   0.102552,
        0.102033,   0.101512,   0.10099,    0.100465,   0.0999401,   0.0994127,  0.098885,   0.098355,   0.0978242,
        0.0972917,  0.0967581,  0.096223,   0.0956865,  0.095149,    0.0946096,  0.0940698,  0.0935275,  0.0929852,
        0.0924401,  0.091895,   0.0913477,  0.0907997,  0.0902502,   0.0896994,  0.0891477,  0.0885942,  0.0880403,
        0.087484,   0.0869276,  0.086369,   0.0858099,  0.0852492,   0.0846874,  0.0841247,  0.0835601,  0.0829956,
        0.0824283,  0.081861,   0.0812919,  0.0807219,  0.080151,    0.0795784,  0.0790057,  0.0784304,  0.0778551,
        0.077278,   0.0767001,  0.0761214,  0.0755408,  0.0749603,   0.0743773,  0.0737942,  0.0732097,  0.0726239,
        0.0720379,  0.0714496,  0.0708612,  0.0702712,  0.0696803,   0.0690889,  0.0684954,  0.067902,   0.0673066,
        0.0667107,  0.066114,   0.0655155,  0.064917,   0.0643166,   0.0637156,  0.063114,   0.0625105,  0.0619071,
        0.0613017,  0.0606958,  0.0600894,  0.059481,   0.0588727,   0.0582627,  0.0576519,  0.0570409,  0.0564277,
        0.0558145,  0.0552,     0.0545844,  0.0539688,  0.0533511,   0.0527331,  0.0521144,  0.051494,   0.0508737,
        0.0502518,  0.0496291,  0.0490063,  0.0483813,  0.0477562,   0.0471304,  0.046503,   0.0458757,  0.0452468,
        0.0446171,  0.0439874,  0.0433555,  0.0427235,  0.042091,    0.0414568,  0.0408225,  0.0401871,  0.0395506,
        0.0389141,  0.0382758,  0.0376371,  0.0369983,  0.0363573,   0.0357163,  0.0350748,  0.0344316,  0.0337884,
        0.0331443,  0.0324988,  0.0318534,  0.0312067,  0.0305591,   0.0299114,  0.0292622,  0.0286124,  0.0279626,
        0.0273108,  0.0266588,  0.0260068,  0.0253527,  0.0246986,   0.0240442,  0.0233879,  0.0227317,  0.0220749,
        0.0214165,  0.0207582,  0.0200992,  0.0194387,  0.0187782,   0.0181169,  0.0174543,  0.0167917,  0.0161284,
        0.0154637,  0.014799,   0.0141335,  0.0134667,  0.0128,      0.0121324,  0.0114636,  0.0107948,  0.0101252,
        0.00945437, 0.00878352, 0.00811197, 0.00743909, 0.00676621,  0.00609275, 0.00541785, 0.00474295, 0.00406763,
        0.00339073, 0.00271382, 0.00203669, 0.00135779, 0.000678897, 0},
  
       {0,          0.00295868, 0.00592246, 0.00888848, 0.0118539,   0.014816,   0.0177719,  0.020719,   0.0236545,
        0.0265758,  0.0294806,  0.0323664,  0.035231,   0.0380722,   0.0408879,  0.0436763,  0.0464357,  0.0491646,
        0.0518614,  0.0545245,  0.0571534,  0.0597465,  0.062303,    0.0648223,  0.0673034,  0.069746,   0.0721497,
        0.0745137,  0.0768385,  0.0791233,  0.081368,   0.083573,    0.0857383,  0.0878639,  0.08995,    0.0919967,
        0.0940044,  0.0959733,  0.0979038,  0.0997963,  0.101651,    0.103469,   0.105251,   0.106997,   0.108706,
        0.11038,    0.112021,   0.113626,   0.115199,   0.116738,    0.118245,   0.11972,    0.121164,   0.122576,
        0.12396,    0.125312,   0.126635,   0.127931,   0.129197,    0.130436,   0.131647,   0.132833,   0.133992,
        0.135125,   0.136233,   0.137316,   0.138375,   0.13941,     0.140421,   0.14141,    0.142376,   0.14332,
        0.144242,   0.145143,   0.146023,   0.146883,   0.147721,    0.148541,   0.149341,   0.150122,   0.150884,
        0.151627,   0.152353,   0.153061,   0.153751,   0.154424,    0.155081,   0.155721,   0.156344,   0.156952,
        0.157543,   0.15812,    0.158681,   0.159228,   0.159759,    0.160277,   0.16078,    0.161269,   0.161745,
        0.162207,   0.162655,   0.163091,   0.163514,   0.163925,    0.164322,   0.164708,   0.165082,   0.165443,
        0.165793,   0.166132,   0.166459,   0.166775,   0.167081,    0.167375,   0.167659,   0.167933,   0.168195,
        0.168448,   0.168691,   0.168925,   0.169148,   0.169361,    0.169566,   0.169761,   0.169947,   0.170124,
        0.170292,   0.170451,   0.170601,   0.170744,   0.170877,    0.171003,   0.17112,    0.171229,   0.17133,
        0.171423,   0.171508,   0.171586,   0.171656,   0.171719,    0.171775,   0.171823,   0.171864,   0.171898,
        0.171926,   0.171946,   0.171959,   0.171966,   0.171966,    0.171959,   0.171946,   0.171926,   0.171901,
        0.171868,   0.17183,    0.171786,   0.171736,   0.171679,    0.171617,   0.171549,   0.171476,   0.171396,
        0.171312,   0.171221,   0.171126,   0.171025,   0.170918,    0.170806,   0.170689,   0.170566,   0.170439,
        0.170307,   0.170169,   0.170027,   0.16988,    0.169728,    0.169571,   0.169409,   0.169243,   0.169073,
        0.168897,   0.168717,   0.168532,   0.168343,   0.16815,     0.167953,   0.167751,   0.167544,   0.167334,
        0.16712,    0.166901,   0.166678,   0.166451,   0.166221,    0.165986,   0.165747,   0.165505,   0.165258,
        0.165009,   0.164754,   0.164496,   0.164235,   0.163971,    0.163702,   0.163429,   0.163154,   0.162875,
        0.162591,   0.162305,   0.162016,   0.161723,   0.161426,    0.161127,   0.160824,   0.160517,   0.160208,
        0.159895,   0.159579,   0.15926,    0.158938,   0.158612,    0.158284,   0.157953,   0.157618,   0.157281,
        0.156941,   0.156597,   0.156251,   0.155903,   0.15555,     0.155196,   0.154838,   0.154478,   0.154115,
        0.153749,   0.15338,    0.15301,    0.152635,   0.152259,    0.15188,    0.151498,   0.151114,   0.150727,
        0.150337,   0.149946,   0.149551,   0.149154,   0.148754,    0.148352,   0.147949,   0.147541,   0.147133,
        0.146721,   0.146307,   0.145892,   0.145473,   0.145053,    0.144629,   0.144204,   0.143777,   0.143347,
        0.142916,   0.142481,   0.142046,   0.141607,   0.141166,    0.140724,   0.140279,   0.139832,   0.139383,
        0.138933,   0.138479,   0.138025,   0.137567,   0.137109,    0.136647,   0.136185,   0.13572,    0.135253,
        0.134784,   0.134313,   0.133841,   0.133366,   0.13289,     0.132412,   0.131932,   0.13145,    0.130967,
        0.130481,   0.129994,   0.129504,   0.129014,   0.128521,    0.128027,   0.12753,    0.127033,   0.126533,
        0.126033,   0.125529,   0.125025,   0.124518,   0.124011,    0.123501,   0.12299,    0.122477,   0.121963,
        0.121447,   0.12093,    0.12041,    0.11989,    0.119367,    0.118843,   0.118318,   0.117791,   0.117263,
        0.116733,   0.116201,   0.115668,   0.115134,   0.114598,    0.114061,   0.113522,   0.112982,   0.11244,
        0.111898,   0.111353,   0.110807,   0.11026,    0.109711,    0.109161,   0.108609,   0.108057,   0.107502,
        0.106947,   0.10639,    0.105832,   0.105272,   0.104712,    0.10415,    0.103586,   0.103022,   0.102455,
        0.101889,   0.10132,    0.10075,    0.100179,   0.0996065,   0.0990336,  0.0984582,  0.0978827,  0.0973051,
        0.0967267,  0.0961473,  0.0955661,  0.0949849,  0.0944009,   0.0938169,  0.0932312,  0.0926444,  0.0920571,
        0.0914674,  0.0908778,  0.0902861,  0.0896937,  0.0891006,   0.0885054,  0.0879102,  0.0873128,  0.0867149,
        0.0861161,  0.0855155,  0.0849148,  0.084312,   0.0837086,   0.0831046,  0.0824985,  0.0818925,  0.0812845,
        0.0806758,  0.0800666,  0.0794552,  0.0788439,  0.0782309,   0.077617,   0.077003,   0.0763864,  0.0757698,
        0.0751521,  0.0745329,  0.0739137,  0.0732926,  0.0726708,   0.0720487,  0.0714244,  0.0708001,  0.0701746,
        0.0695477,  0.0689209,  0.068292,   0.0676627,  0.0670331,   0.0664012,  0.0657693,  0.0651365,  0.0645021,
        0.0638678,  0.0632318,  0.062595,   0.0619582,  0.0613192,   0.0606799,  0.0600403,  0.0593986,  0.0587569,
        0.0581144,  0.0574703,  0.0568261,  0.0561808,  0.0555342,   0.0548877,  0.0542395,  0.0535906,  0.0529417,
        0.0522908,  0.0516395,  0.0509882,  0.0503347,  0.049681,    0.0490272,  0.0483713,  0.0477153,  0.0470589,
        0.0464006,  0.0457423,  0.0450835,  0.0444229,  0.0437623,   0.043101,   0.0424382,  0.0417753,  0.0411116,
        0.0404465,  0.0397813,  0.0391154,  0.038448,   0.0377806,   0.0371125,  0.0364428,  0.0357732,  0.0351028,
        0.034431,   0.0337591,  0.0330867,  0.0324126,  0.0317385,   0.031064,   0.0303878,  0.0297115,  0.029035,
        0.0283566,  0.0276781,  0.0269997,  0.0263191,  0.0256385,   0.0249579,  0.0242754,  0.0235927,  0.0229099,
        0.0222257,  0.0215408,  0.0208559,  0.0201699,  0.0194829,   0.0187959,  0.0181082,  0.0174191,  0.01673,
        0.0160406,  0.0153494,  0.0146583,  0.0139671,  0.013274,    0.0125808,  0.0118875,  0.011193,   0.0104977,
        0.00980235, 0.00910632, 0.00840896, 0.0077116,  0.00701415,  0.00631477, 0.00561538, 0.00491599, 0.00421515,
        0.00351375, 0.00281235, 0.00211021, 0.0014068,  0.000703402, 0},
  
       {0,          0.00372622, 0.00745647, 0.0111854,  0.0149078,   0.0186181,  0.0223116,  0.0259832,  0.029628,
        0.0332419,  0.0368207,  0.0403603,  0.0438576,  0.0473091,   0.0507121,  0.0540643,  0.0573628,  0.0606066,
        0.0637939,  0.0669235,  0.0699942,  0.0730053,  0.0759564,   0.078847,   0.0816773,  0.0844474,  0.0871578,
        0.0898082,  0.0923994,  0.0949326,  0.0974072,  0.0998253,   0.102187,   0.104493,   0.106745,   0.108944,
        0.11109,    0.113184,   0.115229,   0.117224,   0.11917,     0.12107,    0.122923,   0.12473,    0.126494,
        0.128214,   0.129893,   0.131529,   0.133126,   0.134684,    0.136203,   0.137684,   0.139129,   0.140538,
        0.141912,   0.143252,   0.144558,   0.145833,   0.147075,    0.148286,   0.149467,   0.150618,   0.15174,
        0.152834,   0.153901,   0.15494,    0.155953,   0.15694,     0.157902,   0.158839,   0.159751,   0.16064,
        0.161506,   0.16235,    0.163172,   0.163972,   0.16475,     0.165507,   0.166245,   0.166963,   0.16766,
        0.168339,   0.169,      0.169641,   0.170265,   0.170871,    0.171459,   0.172032,   0.172587,   0.173126,
        0.173649,   0.174157,   0.174648,   0.175126,   0.175587,    0.176035,   0.176468,   0.176888,   0.177293,
        0.177685,   0.178064,   0.178429,   0.178782,   0.179123,    0.17945,    0.179766,   0.18007,    0.180362,
        0.180643,   0.180912,   0.18117,    0.181417,   0.181653,    0.181879,   0.182094,   0.182298,   0.182493,
        0.182678,   0.182853,   0.183019,   0.183174,   0.18332,     0.183457,   0.183585,   0.183705,   0.183815,
        0.183917,   0.18401,    0.184095,   0.184171,   0.184239,    0.184299,   0.184351,   0.184396,   0.184432,
        0.184461,   0.184483,   0.184497,   0.184504,   0.184504,    0.184496,   0.184482,   0.18446,    0.184432,
        0.184397,   0.184356,   0.184308,   0.184253,   0.184192,    0.184125,   0.184052,   0.183973,   0.183887,
        0.183795,   0.183697,   0.183594,   0.183484,   0.18337,     0.183249,   0.183124,   0.182993,   0.182855,
        0.182713,   0.182565,   0.182412,   0.182255,   0.182092,    0.181924,   0.18175,    0.181572,   0.181389,
        0.181202,   0.18101,    0.180813,   0.18061,    0.180404,    0.180193,   0.179979,   0.179758,   0.179534,
        0.179305,   0.179073,   0.178835,   0.178594,   0.178348,    0.178099,   0.177845,   0.177587,   0.177325,
        0.17706,    0.17679,    0.176516,   0.176239,   0.175958,    0.175673,   0.175384,   0.175092,   0.174796,
        0.174496,   0.174194,   0.173887,   0.173577,   0.173263,    0.172947,   0.172626,   0.172302,   0.171976,
        0.171645,   0.171311,   0.170975,   0.170634,   0.170291,    0.169945,   0.169595,   0.169243,   0.168888,
        0.168528,   0.168167,   0.167802,   0.167435,   0.167065,    0.166691,   0.166315,   0.165936,   0.165554,
        0.165169,   0.164782,   0.164391,   0.163999,   0.163603,    0.163205,   0.162804,   0.1624,     0.161995,
        0.161585,   0.161175,   0.160761,   0.160344,   0.159926,    0.159504,   0.159081,   0.158654,   0.158226,
        0.157795,   0.157362,   0.156927,   0.156488,   0.156048,    0.155605,   0.155161,   0.154713,   0.154264,
        0.153813,   0.153359,   0.152903,   0.152445,   0.151985,    0.151522,   0.151058,   0.150591,   0.150123,
        0.149652,   0.149179,   0.148704,   0.148227,   0.147748,    0.147268,   0.146785,   0.1463,     0.145813,
        0.145324,   0.144834,   0.144341,   0.143847,   0.143351,    0.142853,   0.142353,   0.141851,   0.141348,
        0.140843,   0.140335,   0.139826,   0.139315,   0.138803,    0.138289,   0.137773,   0.137255,   0.136736,
        0.136215,   0.135692,   0.135168,   0.134641,   0.134114,    0.133584,   0.133054,   0.132521,   0.131987,
        0.131451,   0.130914,   0.130375,   0.129834,   0.129292,    0.128749,   0.128204,   0.127657,   0.127109,
        0.126559,   0.126009,   0.125455,   0.124902,   0.124346,    0.123789,   0.123231,   0.122671,   0.12211,
        0.121547,   0.120984,   0.120417,   0.119851,   0.119283,    0.118713,   0.118142,   0.11757,    0.116997,
        0.116421,   0.115845,   0.115267,   0.114688,   0.114109,    0.113527,   0.112945,   0.11236,    0.111775,
        0.111188,   0.1106,     0.110012,   0.109421,   0.10883,     0.108237,   0.107643,   0.107048,   0.106451,
        0.105854,   0.105255,   0.104655,   0.104055,   0.103452,    0.102849,   0.102244,   0.101638,   0.101032,
        0.100423,   0.0998144,  0.0992041,  0.0985925,  0.0979809,   0.0973666,  0.0967521,  0.0961367,  0.0955194,
        0.0949021,  0.0942828,  0.0936627,  0.0930423,  0.0924194,   0.0917966,  0.0911725,  0.0905469,  0.0899213,
        0.0892935,  0.0886652,  0.0880366,  0.0874055,  0.0867744,   0.0861423,  0.0855086,  0.0848748,  0.0842393,
        0.0836029,  0.0829665,  0.0823277,  0.0816886,  0.0810491,   0.0804074,  0.0797657,  0.0791231,  0.0784788,
        0.0778345,  0.0771887,  0.0765418,  0.0758949,  0.0752462,   0.0745967,  0.0739472,  0.0732955,  0.0726435,
        0.0719914,  0.0713369,  0.0706823,  0.0700275,  0.0693704,   0.0687133,  0.0680558,  0.0673962,  0.0667366,
        0.0660764,  0.0654143,  0.0647522,  0.0640895,  0.0634249,   0.0627604,  0.0620951,  0.0614281,  0.0607611,
        0.0600935,  0.059424,   0.0587546,  0.0580846,  0.0574128,   0.0567409,  0.0560686,  0.0553944,  0.0547202,
        0.0540456,  0.0533691,  0.0526925,  0.0520158,  0.0513369,   0.0506579,  0.049979,   0.0492979,  0.0486166,
        0.0479353,  0.0472522,  0.0465686,  0.045885,   0.0452,      0.0445141,  0.0438282,  0.0431413,  0.0424531,
        0.0417649,  0.0410762,  0.0403858,  0.0396953,  0.0390049,   0.0383122,  0.0376195,  0.0369268,  0.0362324,
        0.0355375,  0.0348426,  0.0341466,  0.0334495,  0.0327523,   0.0320548,  0.0313555,  0.0306561,  0.0299568,
        0.0292556,  0.0285541,  0.0278526,  0.02715,    0.0264463,   0.0257426,  0.0250387,  0.0243328,  0.023627,
        0.0229212,  0.0222138,  0.0215059,  0.0207979,  0.0200893,   0.0193792,  0.0186691,  0.0179591,  0.0172472,
        0.016535,   0.0158228,  0.0151099,  0.0143956,  0.0136813,   0.0129671,  0.012251,   0.0115347,  0.0108183,
        0.0101013,  0.00938289, 0.00866451, 0.00794613, 0.0072261,   0.00650568, 0.00578526, 0.00506437, 0.00434193,
        0.00361949, 0.00289705, 0.00217335, 0.0014489,  0.000724449, 0},
  
       {0,          0.00465353, 0.00930808, 0.0139535,  0.0185797,   0.0231771,  0.0277365,  0.0322493,  0.0367076,
        0.0411044,  0.045433,   0.0496882,  0.0538654,  0.0579605,   0.0619702,  0.065892,   0.0697243,  0.073466,
        0.0771168,  0.0806756,  0.0841437,  0.0875206,  0.0908083,   0.0940077,  0.0971194,  0.100146,   0.103088,
        0.105949,   0.108729,   0.111431,   0.114057,   0.116608,    0.119088,   0.121497,   0.123838,   0.126113,
        0.128323,   0.13047,    0.132557,   0.134585,   0.136556,    0.138471,   0.140333,   0.142142,   0.143901,
        0.14561,    0.147272,   0.148887,   0.150457,   0.151985,    0.153468,   0.154911,   0.156314,   0.157678,
        0.159005,   0.160293,   0.161547,   0.162766,   0.163951,    0.165103,   0.166222,   0.167311,   0.168369,
        0.169398,   0.170398,   0.17137,    0.172314,   0.173232,    0.174124,   0.17499,    0.175831,   0.176648,
        0.177442,   0.178213,   0.178961,   0.179687,   0.180391,    0.181076,   0.181738,   0.182381,   0.183006,
        0.183609,   0.184195,   0.184762,   0.185311,   0.185843,    0.186358,   0.186856,   0.187337,   0.187802,
        0.188251,   0.188685,   0.189104,   0.189507,   0.189897,    0.190271,   0.190632,   0.190979,   0.191312,
        0.191632,   0.19194,    0.192234,   0.192516,   0.192786,    0.193044,   0.193289,   0.193523,   0.193746,
        0.193958,   0.194158,   0.194347,   0.194526,   0.194694,    0.194853,   0.195001,   0.195139,   0.195267,
        0.195385,   0.195494,   0.195594,   0.195684,   0.195766,    0.195839,   0.195903,   0.195958,   0.196005,
        0.196044,   0.196074,   0.196096,   0.19611,    0.196116,    0.196115,   0.196106,   0.196089,   0.196066,
        0.196034,   0.195996,   0.19595,    0.195898,   0.195838,    0.195772,   0.195699,   0.195619,   0.195533,
        0.195441,   0.195342,   0.195237,   0.195126,   0.195009,    0.194885,   0.194756,   0.194621,   0.19448,
        0.194333,   0.194181,   0.194024,   0.19386,    0.193691,    0.193517,   0.193338,   0.193154,   0.192965,
        0.192769,   0.192569,   0.192365,   0.192156,   0.191941,    0.191722,   0.191497,   0.191269,   0.191036,
        0.190798,   0.190555,   0.190309,   0.190058,   0.189802,    0.189542,   0.189278,   0.18901,    0.188737,
        0.18846,    0.18818,    0.187895,   0.187606,   0.187313,    0.187017,   0.186716,   0.186411,   0.186104,
        0.185792,   0.185476,   0.185157,   0.184835,   0.184508,    0.184178,   0.183845,   0.183508,   0.183168,
        0.182824,   0.182476,   0.182127,   0.181773,   0.181415,    0.181056,   0.180692,   0.180325,   0.179956,
        0.179583,   0.179207,   0.178828,   0.178446,   0.178061,    0.177673,   0.177282,   0.176889,   0.176491,
        0.176092,   0.175689,   0.175284,   0.174877,   0.174466,    0.174053,   0.173636,   0.173218,   0.172796,
        0.172372,   0.171946,   0.171516,   0.171085,   0.17065,     0.170213,   0.169774,   0.169332,   0.168888,
        0.168441,   0.167993,   0.16754,    0.167087,   0.16663,     0.166172,   0.165711,   0.165248,   0.164783,
        0.164315,   0.163846,   0.163373,   0.162899,   0.162422,    0.161945,   0.161463,   0.160981,   0.160495,
        0.160009,   0.159519,   0.159029,   0.158535,   0.158041,    0.157543,   0.157044,   0.156542,   0.15604,
        0.155534,   0.155028,   0.154518,   0.154009,   0.153495,    0.152981,   0.152464,   0.151947,   0.151426,
        0.150904,   0.15038,    0.149855,   0.149328,   0.148799,    0.148268,   0.147736,   0.147201,   0.146665,
        0.146128,   0.145588,   0.145048,   0.144505,   0.143961,    0.143415,   0.142868,   0.142318,   0.141768,
        0.141215,   0.140661,   0.140106,   0.139549,   0.138991,    0.13843,    0.137869,   0.137306,   0.136741,
        0.136175,   0.135607,   0.135039,   0.134468,   0.133896,    0.133322,   0.132748,   0.132172,   0.131593,
        0.131015,   0.130434,   0.129853,   0.129269,   0.128684,    0.128099,   0.127511,   0.126923,   0.126332,
        0.125741,   0.125149,   0.124554,   0.12396,    0.123363,    0.122765,   0.122166,   0.121566,   0.120965,
        0.120361,   0.119757,   0.119152,   0.118545,   0.117938,    0.117329,   0.116719,   0.116108,   0.115495,
        0.114881,   0.114266,   0.11365,    0.113033,   0.112414,    0.111795,   0.111175,   0.110552,   0.10993,
        0.109306,   0.10868,    0.108055,   0.107427,   0.106799,    0.106169,   0.105538,   0.104907,   0.104274,
        0.10364,    0.103006,   0.102369,   0.101732,   0.101095,    0.100455,   0.0998157,  0.0991747,  0.0985322,
        0.0978897,  0.0972453,  0.0966001,  0.0959548,  0.0953071,   0.0946591,  0.0940109,  0.0933601,  0.0927093,
        0.0920581,  0.0914045,  0.090751,   0.0900967,  0.0894404,   0.0887842,  0.0881269,  0.087468,   0.086809,
        0.0861489,  0.0854873,  0.0848256,  0.0841627,  0.0834985,   0.0828342,  0.0821686,  0.0815017,  0.0808348,
        0.0801666,  0.0794971,  0.0788275,  0.0781568,  0.0774847,   0.0768126,  0.0761394,  0.0754648,  0.0747901,
        0.0741145,  0.0734373,  0.0727601,  0.0720823,  0.0714025,   0.0707228,  0.0700427,  0.0693605,  0.0686783,
        0.067996,   0.0673113,  0.0666266,  0.0659419,  0.0652552,   0.064568,   0.0638808,  0.0631921,  0.0625025,
        0.0618129,  0.0611223,  0.0604302,  0.0597382,  0.0590457,   0.0583513,  0.0576568,  0.0569624,  0.0562658,
        0.055569,   0.0548721,  0.0541739,  0.0534747,  0.0527755,   0.0520757,  0.0513742,  0.0506726,  0.049971,
        0.0492674,  0.0485635,  0.0478596,  0.0471545,  0.0464483,   0.0457421,  0.0450357,  0.0443272,  0.0436187,
        0.0429102,  0.0422002,  0.0414894,  0.0407786,  0.0400674,   0.0393543,  0.0386413,  0.0379283,  0.0372136,
        0.0364984,  0.0357831,  0.0350674,  0.0343499,  0.0336324,   0.032915,   0.032196,   0.0314763,  0.0307567,
        0.0300367,  0.0293149,  0.028593,   0.0278712,  0.0271482,   0.0264242,  0.0257002,  0.0249762,  0.0242502,
        0.023524,   0.0227979,  0.0220712,  0.0213429,  0.0206146,   0.0198863,  0.0191568,  0.0184264,  0.017696,
        0.0169656,  0.0162334,  0.0155009,  0.0147684,  0.0140356,   0.0133011,  0.0125665,  0.0118319,  0.0110966,
        0.01036,    0.00962335, 0.00888672, 0.00814893, 0.00741025,  0.00667158, 0.00593291, 0.00519273, 0.00445202,
        0.00371132, 0.00297062, 0.00222814, 0.00148543, 0.000742713, 0},
  
       {0,          0.00577596, 0.0115458,  0.0172904,  0.0229911,   0.0286301,  0.034192,   0.0396622, 0.045029,
        0.0502812,  0.0554111,  0.0604123,  0.06528,    0.0700113,   0.0746034,  0.0790576,  0.0833726, 0.087551,
        0.0915947,  0.0955061,  0.0992885,  0.102945,   0.106479,    0.109895,   0.113196,   0.116386,  0.11947,
        0.12245,    0.125331,   0.128116,   0.130809,   0.133412,    0.135931,   0.138367,   0.140724,  0.143005,
        0.145212,   0.147348,   0.149416,   0.151419,   0.153359,    0.155238,   0.157058,   0.158821,  0.160529,
        0.162185,   0.163789,   0.165344,   0.166852,   0.168314,    0.169731,   0.171105,   0.172436,  0.173728,
        0.174982,   0.176196,   0.177373,   0.178515,   0.179623,    0.180697,   0.181738,   0.182748,  0.183727,
        0.184676,   0.185596,   0.186487,   0.187352,   0.188189,    0.189001,   0.189787,   0.190548,  0.191285,
        0.192,      0.19269,    0.193359,   0.194007,   0.194631,    0.195237,   0.195822,   0.196386,  0.196933,
        0.197459,   0.197969,   0.198459,   0.198933,   0.199388,    0.199827,   0.20025,    0.200656,  0.201048,
        0.201423,   0.201784,   0.20213,    0.20246,    0.202778,    0.203081,   0.20337,    0.203646,  0.20391,
        0.20416,    0.204398,   0.204623,   0.204837,   0.205039,    0.205229,   0.205407,   0.205574,  0.205731,
        0.205876,   0.206012,   0.206136,   0.206249,   0.206353,    0.206447,   0.206532,   0.206607,  0.206673,
        0.206729,   0.206776,   0.206815,   0.206844,   0.206865,    0.206877,   0.206881,   0.206877,  0.206865,
        0.206844,   0.206816,   0.20678,    0.206736,   0.206685,    0.206627,   0.206561,   0.206488,  0.206408,
        0.206321,   0.206228,   0.206127,   0.206021,   0.205907,    0.205786,   0.205659,   0.205526,  0.205387,
        0.205243,   0.205092,   0.204935,   0.204771,   0.204602,    0.204428,   0.204248,   0.204063,  0.203871,
        0.203674,   0.203472,   0.203265,   0.203053,   0.202836,    0.202613,   0.202385,   0.202153,  0.201916,
        0.201673,   0.201425,   0.201174,   0.200918,   0.200656,    0.200391,   0.200121,   0.199847,  0.199568,
        0.199285,   0.198998,   0.198705,   0.198409,   0.19811,     0.197806,   0.197497,   0.197186,  0.19687,
        0.196549,   0.196226,   0.195899,   0.195567,   0.195232,    0.194894,   0.19455,    0.194205,  0.193855,
        0.193502,   0.193146,   0.192785,   0.192422,   0.192056,    0.191685,   0.191311,   0.190935,  0.190554,
        0.190171,   0.189785,   0.189395,   0.189003,   0.188606,    0.188208,   0.187806,   0.1874,    0.186994,
        0.186582,   0.186169,   0.185752,   0.185332,   0.18491,     0.184485,   0.184058,   0.183626,  0.183193,
        0.182757,   0.182319,   0.181878,   0.181433,   0.180988,    0.180538,   0.180087,   0.179632,  0.179176,
        0.178716,   0.178255,   0.177791,   0.177324,   0.176856,    0.176384,   0.175912,   0.175435,  0.174958,
        0.174476,   0.173995,   0.173509,   0.173023,   0.172532,    0.172042,   0.171547,   0.171052,  0.170553,
        0.170054,   0.169551,   0.169047,   0.16854,    0.168032,    0.167521,   0.16701,    0.166494,  0.165979,
        0.165459,   0.16494,    0.164417,   0.163893,   0.163366,    0.162839,   0.162308,   0.161777,  0.161243,
        0.160708,   0.160171,   0.159631,   0.159091,   0.158548,    0.158005,   0.157458,   0.15691,   0.15636,
        0.155809,   0.155256,   0.154702,   0.154146,   0.153587,    0.153029,   0.152466,   0.151904,  0.151339,
        0.150773,   0.150206,   0.149636,   0.149066,   0.148493,    0.14792,    0.147344,   0.146767,  0.146189,
        0.145609,   0.145029,   0.144445,   0.143861,   0.143275,    0.142688,   0.1421,     0.141509,  0.140918,
        0.140325,   0.13973,    0.139135,   0.138538,   0.13794,     0.13734,    0.136739,   0.136137,  0.135532,
        0.134928,   0.134321,   0.133713,   0.133105,   0.132494,    0.131883,   0.13127,    0.130656,  0.130041,
        0.129424,   0.128806,   0.128187,   0.127567,   0.126946,    0.126323,   0.125699,   0.125074,  0.124447,
        0.12382,    0.123192,   0.122562,   0.121932,   0.121299,    0.120666,   0.120032,   0.119396,  0.11876,
        0.118122,   0.117483,   0.116844,   0.116203,   0.11556,     0.114918,   0.114273,   0.113628,  0.112983,
        0.112334,   0.111686,   0.111037,   0.110386,   0.109735,    0.109083,   0.108429,   0.107775,  0.107119,
        0.106462,   0.105805,   0.105147,   0.104487,   0.103827,    0.103166,   0.102503,   0.10184,   0.101176,
        0.10051,    0.0998443,  0.0991772,  0.0985088,  0.0978403,   0.0971705,  0.0964992,  0.095828,  0.0951555,
        0.0944816,  0.0938076,  0.0931325,  0.0924558,  0.0917791,   0.0911016,  0.0904222,  0.0897428, 0.0890629,
        0.0883808,  0.0876987,  0.0870165,  0.0863318,  0.085647,    0.0849623,  0.0842752,  0.0835879, 0.0829005,
        0.0822114,  0.0815214,  0.0808314,  0.0801402,  0.0794477,   0.0787551,  0.078062,   0.0773669, 0.0766717,
        0.0759766,  0.0752791,  0.0745814,  0.0738837,  0.0731845,   0.0724843,  0.0717841,  0.0710831, 0.0703804,
        0.0696777,  0.068975,   0.06827,    0.0675648,  0.0668596,   0.0661532,  0.0654455,  0.0647379, 0.06403,
        0.0633199,  0.0626098,  0.0618997,  0.061188,   0.0604755,   0.059763,   0.0590501,  0.0583352, 0.0576203,
        0.0569054,  0.0561889,  0.0554716,  0.0547543,  0.0540368,   0.0533171,  0.0525975,  0.0518779, 0.0511569,
        0.0504349,  0.049713,   0.048991,   0.0482668,  0.0475425,   0.0468182,  0.0460932,  0.0453666, 0.04464,
        0.0439134,  0.0431853,  0.0424565,  0.0417276,  0.0409988,   0.0402677,  0.0395366,  0.0388055, 0.0380737,
        0.0373404,  0.036607,   0.0358737,  0.0351392,  0.0344036,   0.0336681,  0.0329325,  0.0321953, 0.0314576,
        0.0307198,  0.0299821,  0.0292423,  0.0285024,  0.0277624,   0.0270223,  0.0262802,  0.0255381, 0.024796,
        0.0240535,  0.0233093,  0.022565,   0.0218208,  0.0210759,   0.0203296,  0.0195832,  0.0188369, 0.0180898,
        0.0173413,  0.0165928,  0.0158444,  0.0150951,  0.0143446,   0.013594,   0.0128435,  0.0120922, 0.0113396,
        0.010587,   0.00983434, 0.00908105, 0.00832639, 0.00757173,  0.00681707, 0.00606189, 0.0053052, 0.00454851,
        0.00379181, 0.00303483, 0.00227612, 0.00151741, 0.000758707, 0.}};
  
   //---- Corrections for ME1/1 chambers (very different parameters from the other chambers):
  
   //---- min SW ~ 3.15 cm, max SW ~ 7.6 cm
   //---- Corrections are calculated for N_SW_ME1_1 strip widths (0.3 - 0.8 cm)
   //---- and for 501 values in
   float corrections_ME1_1[n_SW_ME1_1][n_val] = {
       {0,
        -0.000153178,
        -0.000304723,
        -0.000454635,
        -0.000602913,
        -0.000749555,
        -0.000894561,
        -0.00103785,
        -0.00117926,
        -0.00131904,
        -0.0014572,
        -0.00159372,
        -0.00172861,
        -0.00186188,
        -0.00199351,
        -0.00212325,
        -0.00225135,
        -0.00237783,
        -0.00250269,
        -0.00262593,
        -0.00274755,
        -0.00286755,
        -0.00298594,
        -0.00310249,
        -0.00321741,
        -0.00333073,
        -0.00344243,
        -0.00355254,
        -0.00366104,
        -0.00376795,
        -0.00387325,
        -0.00397696,
        -0.00407892,
        -0.00417925,
        -0.004278,
        -0.00437517,
        -0.00447075,
        -0.00456477,
        -0.00465721,
        -0.00474808,
        -0.00483738,
        -0.00492512,
        -0.00501128,
        -0.00509574,
        -0.00517866,
        -0.00526003,
        -0.00533986,
        -0.00541815,
        -0.0054949,
        -0.00557013,
        -0.00564383,
        -0.005716,
        -0.00578666,
        -0.0058558,
        -0.00592344,
        -0.00598956,
        -0.0060541,
        -0.00611712,
        -0.00617865,
        -0.00623871,
        -0.00629728,
        -0.00635438,
        -0.00641001,
        -0.00646418,
        -0.00651689,
        -0.00656815,
        -0.00661796,
        -0.00666633,
        -0.00671326,
        -0.00675877,
        -0.00680284,
        -0.0068455,
        -0.00688674,
        -0.00692658,
        -0.00696502,
        -0.00700205,
        -0.00703765,
        -0.00707187,
        -0.00710472,
        -0.0071362,
        -0.00716631,
        -0.00719508,
        -0.00722249,
        -0.00724857,
        -0.00727331,
        -0.00729672,
        -0.00731881,
        -0.00733959,
        -0.00735907,
        -0.00737724,
        -0.00739412,
        -0.00740972,
        -0.00742404,
        -0.00743709,
        -0.00744888,
        -0.00745941,
        -0.0074687,
        -0.00747675,
        -0.00748356,
        -0.00748915,
        -0.00749352,
        -0.00749668,
        -0.00749865,
        -0.00749942,
        -0.00749901,
        -0.00749742,
        -0.00749466,
        -0.00749074,
        -0.00748567,
        -0.00747946,
        -0.00747211,
        -0.00746363,
        -0.00745404,
        -0.00744334,
        -0.00743153,
        -0.00741864,
        -0.00740466,
        -0.00738961,
        -0.00737348,
        -0.00735631,
        -0.00733808,
        -0.00731881,
        -0.00729852,
        -0.0072772,
        -0.00725487,
        -0.00723153,
        -0.0072072,
        -0.00718189,
        -0.0071556,
        -0.00712834,
        -0.00710012,
        -0.00707096,
        -0.00704085,
        -0.00700982,
        -0.00697785,
        -0.00694495,
        -0.00691116,
        -0.00687647,
        -0.0068409,
        -0.00680445,
        -0.00676714,
        -0.00672898,
        -0.00668997,
        -0.00665013,
        -0.00660946,
        -0.00656797,
        -0.00652568,
        -0.00648259,
        -0.00643872,
        -0.00639407,
        -0.00634865,
        -0.00630247,
        -0.00625554,
        -0.00620788,
        -0.00615949,
        -0.00611038,
        -0.00606056,
        -0.00601004,
        -0.0059588,
        -0.00590688,
        -0.0058543,
        -0.00580105,
        -0.00574716,
        -0.00569262,
        -0.00563746,
        -0.00558167,
        -0.00552528,
        -0.00546828,
        -0.0054107,
        -0.00535254,
        -0.0052938,
        -0.00523451,
        -0.00517467,
        -0.00511428,
        -0.00505337,
        -0.00499193,
        -0.00492995,
        -0.00486748,
        -0.00480451,
        -0.00474107,
        -0.00467715,
        -0.00461278,
        -0.00454795,
        -0.00448268,
        -0.00441698,
        -0.00435086,
        -0.00428433,
        -0.0042174,
        -0.00415007,
        -0.00408237,
        -0.00401429,
        -0.00394583,
        -0.00387701,
        -0.00380785,
        -0.00373836,
        -0.00366854,
        -0.00359842,
        -0.00352799,
        -0.00345726,
        -0.00338625,
        -0.00331496,
        -0.00324341,
        -0.00317161,
        -0.00309955,
        -0.00302726,
        -0.00295474,
        -0.00288199,
        -0.00280904,
        -0.00273588,
        -0.00266254,
        -0.00258902,
        -0.00251532,
        -0.00244147,
        -0.00236747,
        -0.00229332,
        -0.00221904,
        -0.00214463,
        -0.00207011,
        -0.00199549,
        -0.00192076,
        -0.00184594,
        -0.00177105,
        -0.00169609,
        -0.00162107,
        -0.001546,
        -0.00147088,
        -0.00139573,
        -0.00132056,
        -0.00124538,
        -0.00117018,
        -0.00109499,
        -0.00101981,
        -0.000944646,
        -0.000869513,
        -0.000794418,
        -0.000719369,
        -0.000644374,
        -0.000569442,
        -0.000494581,
        -0.000419799,
        -0.000345105,
        -0.000270507,
        -0.000196013,
        -0.000121632,
        -4.73709e-05,
        2.67578e-05,
        0.000100743,
        0.000174583,
        0.000248269,
        0.000321794,
        0.000395148,
        0.000468326,
        0.000541318,
        0.000614116,
        0.000686715,
        0.000759104,
        0.000831278,
        0.000903227,
        0.000974946,
        0.00104642,
        0.00111763,
        0.0011886,
        0.0012593,
        0.00132973,
        0.00139988,
        0.00146975,
        0.00153933,
        0.00160862,
        0.00167759,
        0.00174626,
        0.0018146,
        0.00188262,
        0.00195031,
        0.00201765,
        0.00208463,
        0.00215125,
        0.00221752,
        0.00228341,
        0.00234893,
        0.00241407,
        0.00247883,
        0.00254319,
        0.00260715,
        0.0026707,
        0.00273384,
        0.00279656,
        0.00285885,
        0.00292071,
        0.00298214,
        0.0030431,
        0.0031036,
        0.00316365,
        0.00322324,
        0.00328236,
        0.003341,
        0.00339916,
        0.00345684,
        0.00351402,
        0.00357071,
        0.0036269,
        0.00368257,
        0.00373773,
        0.00379238,
        0.0038465,
        0.00390008,
        0.00395313,
        0.00400564,
        0.00405758,
        0.00410897,
        0.0041598,
        0.00421006,
        0.00425977,
        0.00430889,
        0.00435745,
        0.00440542,
        0.0044528,
        0.00449959,
        0.00454579,
        0.00459138,
        0.00463637,
        0.00468075,
        0.00472452,
        0.00476766,
        0.00481018,
        0.00485207,
        0.00489333,
        0.00493396,
        0.00497394,
        0.00501327,
        0.00505193,
        0.00508994,
        0.00512729,
        0.00516398,
        0.0052,
        0.00523535,
        0.00527003,
        0.00530403,
        0.00533735,
        0.00536998,
        0.00540192,
        0.00543317,
        0.00546373,
        0.00549358,
        0.00552273,
        0.00555118,
        0.00557891,
        0.00560592,
        0.00563222,
        0.0056578,
        0.00568266,
        0.00570678,
        0.00573018,
        0.00575284,
        0.00577476,
        0.00579595,
        0.00581639,
        0.00583608,
        0.00585503,
        0.00587322,
        0.00589066,
        0.00590734,
        0.00592325,
        0.00593841,
        0.0059528,
        0.00596642,
        0.00597926,
        0.00599133,
        0.00600263,
        0.00601314,
        0.00602286,
        0.00603181,
        0.00603996,
        0.00604733,
        0.00605391,
        0.00605969,
        0.00606468,
        0.00606887,
        0.00607226,
        0.00607484,
        0.00607663,
        0.0060776,
        0.00607777,
        0.00607713,
        0.00607568,
        0.00607341,
        0.00607033,
        0.00606642,
        0.0060617,
        0.00605616,
        0.00604979,
        0.0060426,
        0.00603458,
        0.00602574,
        0.00601606,
        0.00600555,
        0.00599421,
        0.00598202,
        0.005969,
        0.00595514,
        0.00594044,
        0.0059249,
        0.00590852,
        0.00589129,
        0.00587322,
        0.0058543,
        0.00583454,
        0.00581392,
        0.00579246,
        0.00577014,
        0.00574698,
        0.00572296,
        0.00569809,
        0.00567236,
        0.00564578,
        0.00561834,
        0.00559004,
        0.00556088,
        0.00553087,
        0.00549999,
        0.00546825,
        0.00543566,
        0.00540219,
        0.00536787,
        0.00533268,
        0.00529663,
        0.00525971,
        0.00522193,
        0.00518328,
        0.00514376,
        0.00510338,
        0.00506213,
        0.00502001,
        0.004977,
        0.0049331,
        0.00488833,
        0.00484269,
        0.00479618,
        0.0047488,
        0.00470054,
        0.00465142,
        0.00460142,
        0.00455055,
        0.00449881,
        0.00444619,
        0.0043927,
        0.00433834,
        0.00428311,
        0.004227,
        0.00417002,
        0.00411217,
        0.00405344,
        0.00399384,
        0.00393331,
        0.0038719,
        0.00380962,
        0.00374647,
        0.00368244,
        0.00361755,
        0.00355177,
        0.00348513,
        0.00341761,
        0.00334922,
        0.00327996,
        0.00320982,
        0.00313882,
        0.00306694,
        0.00299419,
        0.00292049,
        0.00284593,
        0.00277049,
        0.00269419,
        0.00261701,
        0.00253897,
        0.00246005,
        0.00238027,
        0.00229962,
        0.0022181,
        0.00213572,
        0.00205247,
        0.00196832,
        0.00188325,
        0.00179732,
        0.00171052,
        0.00162285,
        0.00153432,
        0.00144493,
        0.00135468,
        0.00126357,
        0.00117159,
        0.00107876,
        0.000985051,
        0.000890399,
        0.000794889,
        0.000698521,
        0.000601296,
        0.000503215,
        0.000404279,
        0.000304488,
        0.000203844,
        0.000102348,
        0},
  
       {0,           6.5726e-05,  0.000133561,
        0.000203493, 0.000275507, 0.000349591,
        0.000425731, 0.000503913, 0.000584123,
        0.000666346, 0.00075057,  0.000836778,
        0.000924956, 0.00101512,  0.00110734,
        0.00120149,  0.00129755,  0.0013955,
        0.00149533,  0.00159703,  0.00170057,
        0.00180594,  0.00191312,  0.00202214,
        0.00213306,  0.00224575,  0.00236018,
        0.00247634,  0.00259421,  0.00271377,
        0.00283502,  0.00295792,  0.00308257,
        0.00320889,  0.00333682,  0.00346633,
        0.00359741,  0.00373004,  0.0038642,
        0.00399987,  0.00413719,  0.00427598,
        0.00441622,  0.00455789,  0.00470098,
        0.00484546,  0.00499132,  0.00513867,
        0.00528737,  0.00543738,  0.00558868,
        0.00574126,  0.00589509,  0.0060502,
        0.00620661,  0.00636422,  0.006523,
        0.00668292,  0.00684398,  0.00700614,
        0.00716953,  0.00733397,  0.00749946,
        0.00766597,  0.00783349,  0.00800198,
        0.00817156,  0.00834207,  0.00851349,
        0.00868581,  0.008859,    0.00903307,
        0.00920805,  0.00938383,  0.0095604,
        0.00973774,  0.00991582,  0.0100947,
        0.0102743,   0.0104545,   0.0106355,
        0.010817,    0.0109992,   0.011182,
        0.0113654,   0.0115494,   0.0117339,
        0.0119188,   0.0121043,   0.0122903,
        0.0124768,   0.0126636,   0.0128508,
        0.0130384,   0.0132265,   0.0134148,
        0.0136034,   0.0137924,   0.0139816,
        0.014171,    0.0143607,   0.0145506,
        0.0147407,   0.014931,    0.0151214,
        0.0153119,   0.0155025,   0.0156932,
        0.0158839,   0.0160747,   0.0162655,
        0.0164563,   0.0166471,   0.0168378,
        0.0170284,   0.017219,    0.0174094,
        0.0175997,   0.0177899,   0.0179799,
        0.0181697,   0.0183593,   0.0185486,
        0.0187377,   0.0189266,   0.0191151,
        0.0193033,   0.0194912,   0.0196788,
        0.019866,    0.0200528,   0.0202392,
        0.0204251,   0.0206106,   0.0207957,
        0.0209803,   0.0211643,   0.0213479,
        0.0215309,   0.0217134,   0.0218953,
        0.0220766,   0.0222573,   0.0224374,
        0.0226168,   0.0227956,   0.0229738,
        0.0231512,   0.0233279,   0.0235038,
        0.0236791,   0.0238536,   0.0240274,
        0.0242003,   0.0243724,   0.0245437,
        0.0247142,   0.0248839,   0.0250527,
        0.0252206,   0.0253875,   0.0255536,
        0.0257188,   0.0258831,   0.0260464,
        0.0262087,   0.02637,     0.0265304,
        0.0266898,   0.0268482,   0.0270055,
        0.0271617,   0.0273169,   0.0274711,
        0.0276241,   0.0277761,   0.0279271,
        0.0280768,   0.0282254,   0.0283729,
        0.0285192,   0.0286644,   0.0288084,
        0.0289513,   0.0290929,   0.0292333,
        0.0293725,   0.0295104,   0.0296472,
        0.0297827,   0.029917,    0.03005,
        0.0301816,   0.030312,    0.0304411,
        0.030569,    0.0306955,   0.0308207,
        0.0309446,   0.0310671,   0.0311883,
        0.0313081,   0.0314265,   0.0315436,
        0.0316593,   0.0317737,   0.0318867,
        0.0319983,   0.0321083,   0.032217,
        0.0323243,   0.0324301,   0.0325346,
        0.0326376,   0.0327392,   0.0328393,
        0.032938,    0.0330352,   0.0331308,
        0.033225,    0.0333178,   0.033409,
        0.0334988,   0.033587,    0.0336738,
        0.0337591,   0.0338429,   0.0339252,
        0.0340059,   0.0340851,   0.0341627,
        0.0342388,   0.0343134,   0.0343864,
        0.0344579,   0.0345278,   0.0345962,
        0.0346631,   0.0347284,   0.0347921,
        0.0348543,   0.0349149,   0.034974,
        0.0350314,   0.0350872,   0.0351414,
        0.0351941,   0.0352452,   0.0352947,
        0.0353426,   0.0353889,   0.0354337,
        0.0354768,   0.0355183,   0.0355583,
        0.0355966,   0.0356334,   0.0356685,
        0.0357021,   0.035734,    0.0357643,
        0.035793,    0.0358201,   0.0358456,
        0.0358695,   0.0358917,   0.0359124,
        0.0359314,   0.0359488,   0.0359646,
        0.0359788,   0.0359913,   0.0360023,
        0.0360116,   0.0360193,   0.0360253,
        0.0360298,   0.0360326,   0.0360338,
        0.0360334,   0.0360314,   0.0360278,
        0.0360225,   0.0360156,   0.0360071,
        0.035997,    0.0359853,   0.0359719,
        0.035957,    0.0359404,   0.0359222,
        0.0359024,   0.035881,    0.0358579,
        0.0358333,   0.035807,    0.0357792,
        0.0357497,   0.0357186,   0.035686,
        0.0356517,   0.0356158,   0.0355783,
        0.0355393,   0.0354986,   0.0354563,
        0.0354125,   0.0353671,   0.03532,
        0.0352714,   0.0352212,   0.0351695,
        0.0351161,   0.0350612,   0.0350047,
        0.0349465,   0.0348868,   0.0348255,
        0.0347626,   0.0346982,   0.0346322,
        0.0345646,   0.0344955,   0.0344248,
        0.0343526,   0.0342788,   0.0342035,
        0.0341267,   0.0340483,   0.0339683,
        0.0338867,   0.0338036,   0.0337189,
        0.0336327,   0.033545,    0.0334558,
        0.033365,    0.0332727,   0.033179,
        0.0330837,   0.0329869,   0.0328884,
        0.0327885,   0.032687,    0.0325841,
        0.0324796,   0.0323737,   0.0322663,
        0.0321574,   0.0320471,   0.0319352,
        0.0318217,   0.0317067,   0.0315903,
        0.0314724,   0.0313531,   0.0312323,
        0.0311101,   0.0309864,   0.0308611,
        0.0307343,   0.0306061,   0.0304764,
        0.0303454,   0.0302129,   0.030079,
        0.0299436,   0.0298066,   0.0296683,
        0.0295285,   0.0293873,   0.0292448,
        0.0291008,   0.0289554,   0.0288084,
        0.02866,     0.0285102,   0.0283591,
        0.0282066,   0.0280528,   0.0278974,
        0.0277405,   0.0275823,   0.0274228,
        0.0272619,   0.0270996,   0.026936,
        0.0267707,   0.0266042,   0.0264363,
        0.0262672,   0.0260967,   0.0259247,
        0.0257513,   0.0255766,   0.0254005,
        0.0252232,   0.0250446,   0.0248644,
        0.0246829,   0.0245001,   0.024316,
        0.0241307,   0.023944,    0.0237558,
        0.0235663,   0.0233755,   0.0231836,
        0.0229903,   0.0227955,   0.0225994,
        0.0224021,   0.0222036,   0.0220039,
        0.0218025,   0.0215999,   0.0213962,
        0.0211912,   0.020985,    0.0207772,
        0.0205682,   0.0203581,   0.0201468,
        0.0199341,   0.01972,     0.0195047,
        0.0192883,   0.0190707,   0.0188516,
        0.0186313,   0.0184098,   0.0181872,
        0.0179633,   0.017738,    0.0175115,
        0.0172839,   0.0170552,   0.016825,
        0.0165935,   0.016361,    0.0161274,
        0.0158925,   0.0156562,   0.0154188,
        0.0151803,   0.0149406,   0.0146995,
        0.0144573,   0.0142141,   0.0139697,
        0.0137239,   0.0134769,   0.0132289,
        0.0129799,   0.0127293,   0.0124777,
        0.0122251,   0.0119713,   0.0117161,
        0.0114599,   0.0112026,   0.0109442,
        0.0106844,   0.0104236,   0.0101618,
        0.00989884,  0.00963451,  0.0093692,
        0.00910291,  0.0088353,   0.00856651,
        0.00829675,  0.00802603,  0.00775384,
        0.00748065,  0.0072065,   0.00693127,
        0.00665466,  0.00637711,  0.00609864,
        0.00581885,  0.00553794,  0.00525611,
        0.00497331,  0.00468906,  0.0044039,
        0.00411784,  0.00383053,  0.00354207,
        0.00325273,  0.00296242,  0.0026707,
        0.0023781,   0.00208463,  0.00178985,
        0.00149403,  0.00119734,  0.000899595,
        0.000600576, 0.000300709, 0},
  
       {0,          0.000336602, 0.000675667, 0.00101718,  0.00136154,  0.00170824, 0.00205733, 0.00240907, 0.00276303,
        0.00311932, 0.00347805,  0.00383888,  0.004202,    0.00456732,  0.00493462, 0.00530417, 0.00567568, 0.00604907,
        0.0064246,  0.00680187,  0.00718103,  0.00756203,  0.00794465,  0.00832913, 0.00871518, 0.0091028,  0.00949208,
        0.00988274, 0.010275,    0.0106685,   0.0110634,   0.0114597,   0.0118572,  0.0122559,  0.0126558,  0.0130567,
        0.0134588,  0.0138618,   0.0142658,   0.0146707,   0.0150763,   0.0154829,  0.0158901,  0.0162981,  0.0167067,
        0.0171159,  0.0175257,   0.0179359,   0.0183466,   0.0187577,   0.0191691,  0.0195808,  0.0199927,  0.0204048,
        0.0208171,  0.0212294,   0.0216418,   0.0220541,   0.0224664,   0.0228786,  0.0232907,  0.0237025,  0.024114,
        0.0245253,  0.0249362,   0.0253467,   0.0257567,   0.0261662,   0.0265752,  0.0269837,  0.0273914,  0.0277985,
        0.0282049,  0.0286104,   0.0290152,   0.0294191,   0.0298221,   0.0302242,  0.0306252,  0.0310253,  0.0314242,
        0.031822,   0.0322187,   0.0326141,   0.0330084,   0.0334013,   0.0337929,  0.0341832,  0.0345721,  0.0349596,
        0.0353455,  0.0357299,   0.0361129,   0.0364942,   0.0368739,   0.037252,   0.0376283,  0.0380031,  0.0383759,
        0.038747,   0.0391164,   0.0394838,   0.0398494,   0.0402131,   0.0405747,  0.0409346,  0.0412923,  0.041648,
        0.0420018,  0.0423533,   0.0427028,   0.0430502,   0.0433953,   0.0437383,  0.0440792,  0.0444176,  0.044754,
        0.0450881,  0.0454197,   0.0457491,   0.0460762,   0.0464007,   0.046723,   0.047043,   0.0473603,  0.0476753,
        0.047988,   0.0482978,   0.0486053,   0.0489105,   0.0492129,   0.0495127,  0.0498101,  0.0501049,  0.050397,
        0.0506865,  0.0509736,   0.0512577,   0.0515393,   0.0518183,   0.0520946,  0.0523681,  0.052639,   0.0529072,
        0.0531726,  0.0534352,   0.0536952,   0.0539526,   0.0542069,   0.0544585,  0.0547074,  0.0549536,  0.0551968,
        0.0554373,  0.055675,    0.05591,     0.056142,    0.0563711,   0.0565975,  0.0568211,  0.057042,   0.0572597,
        0.0574746,  0.0576868,   0.0578962,   0.0581027,   0.0583061,   0.0585068,  0.0587047,  0.0588998,  0.0590919,
        0.059281,   0.0594673,   0.0596508,   0.0598315,   0.0600093,   0.060184,   0.0603559,  0.0605249,  0.0606912,
        0.0608546,  0.0610152,   0.0611726,   0.0613272,   0.061479,    0.061628,   0.0617742,  0.0619175,  0.0620579,
        0.0621953,  0.0623299,   0.0624616,   0.0625906,   0.0627167,   0.06284,    0.0629605,  0.063078,   0.0631927,
        0.0633045,  0.0634136,   0.0635198,   0.0636232,   0.0637239,   0.0638218,  0.0639169,  0.0640092,  0.0640985,
        0.064185,   0.0642688,   0.0643499,   0.0644281,   0.0645037,   0.0645765,  0.0646465,  0.0647138,  0.0647784,
        0.0648402,  0.0648993,   0.0649557,   0.0650094,   0.0650602,   0.0651084,  0.0651539,  0.0651967,  0.0652368,
        0.0652742,  0.065309,    0.0653411,   0.0653706,   0.0653974,   0.0654216,  0.0654431,  0.065462,   0.0654783,
        0.065492,   0.0655031,   0.0655116,   0.0655175,   0.0655209,   0.0655216,  0.0655198,  0.0655154,  0.0655085,
        0.065499,   0.065487,    0.0654725,   0.0654555,   0.0654359,   0.0654138,  0.0653893,  0.0653623,  0.0653327,
        0.0653008,  0.0652663,   0.0652294,   0.0651901,   0.0651483,   0.0651042,  0.0650576,  0.0650086,  0.064957,
        0.0649031,  0.0648468,   0.0647881,   0.064727,    0.0646636,   0.0645979,  0.0645298,  0.0644593,  0.0643866,
        0.0643115,  0.0642342,   0.0641545,   0.0640726,   0.0639884,   0.0639017,  0.0638128,  0.0637216,  0.0636281,
        0.0635325,  0.0634346,   0.0633345,   0.0632322,   0.0631278,   0.0630211,  0.0629121,  0.0628008,  0.0626874,
        0.0625718,  0.0624541,   0.0623343,   0.0622124,   0.0620883,   0.0619621,  0.0618335,  0.0617029,  0.0615702,
        0.0614354,  0.0612986,   0.0611598,   0.0610189,   0.0608757,   0.0607304,  0.0605831,  0.0604338,  0.0602826,
        0.0601293,  0.0599741,   0.0598165,   0.059657,    0.0594955,   0.0593321,  0.0591669,  0.0589997,  0.0588301,
        0.0586587,  0.0584854,   0.0583102,   0.0581332,   0.0579542,   0.0577729,  0.0575899,  0.057405,   0.0572183,
        0.0570299,  0.0568392,   0.0566466,   0.0564522,   0.0562561,   0.0560583,  0.0558583,  0.0556563,  0.0554527,
        0.0552474,  0.0550403,   0.0548311,   0.0546201,   0.0544074,   0.054193,   0.0539769,  0.0537586,  0.0535387,
        0.0533171,  0.0530939,   0.0528688,   0.0526417,   0.0524131,   0.0521829,  0.0519509,  0.0517169,  0.0514813,
        0.0512441,  0.0510055,   0.0507646,   0.0505221,   0.0502782,   0.0500328,  0.0497852,  0.0495361,  0.0492855,
        0.0490334,  0.0487792,   0.0485235,   0.0482664,   0.0480078,   0.047747,   0.0474849,  0.0472213,  0.0469562,
        0.046689,   0.0464205,   0.0461506,   0.045879,    0.0456056,   0.0453309,  0.0450548,  0.0447767,  0.0444972,
        0.0442164,  0.043934,    0.0436498,   0.0433642,   0.0430774,   0.0427887,  0.0424985,  0.042207,   0.0419141,
        0.0416193,  0.0413233,   0.0410261,   0.0407269,   0.0404264,   0.0401247,  0.0398213,  0.0395163,  0.0392102,
        0.0389026,  0.0385933,   0.0382828,   0.0379711,   0.0376574,   0.0373427,  0.0370268,  0.0367089,  0.0363899,
        0.0360698,  0.0357479,   0.0354247,   0.0351004,   0.0347744,   0.0344471,  0.0341187,  0.0337887,  0.0334573,
        0.0331248,  0.0327908,   0.0324554,   0.032119,    0.0317809,   0.0314415,  0.0311012,  0.0307592,  0.0304159,
        0.0300717,  0.0297257,   0.0293786,   0.0290306,   0.0286807,   0.0283298,  0.0279779,  0.0276242,  0.0272695,
        0.0269138,  0.0265563,   0.026198,    0.0258384,   0.0254773,   0.0251153,  0.0247518,  0.0243872,  0.0240217,
        0.0236543,  0.0232861,   0.0229169,   0.022546,    0.0221743,   0.0218012,  0.0214268,  0.0210517,  0.0206749,
        0.0202971,  0.0199184,   0.019538,    0.019157,    0.0187745,   0.0183908,  0.0180064,  0.0176202,  0.0172333,
        0.0168452,  0.0164558,   0.0160656,   0.0156739,   0.0152812,   0.0148876,  0.0144925,  0.0140967,  0.0136994,
        0.0133012,  0.0129021,   0.0125015,   0.0121002,   0.0116975,   0.0112938,  0.0108893,  0.0104832,  0.0100766,
        0.00966846, 0.00925949,  0.00884953,  0.00843822,  0.00802638,  0.00761284, 0.00719869, 0.00678337, 0.00636692,
        0.00594983, 0.00553111,  0.0051119,   0.00469127,  0.00426981,  0.00384742, 0.00342373, 0.00299959, 0.00257369,
        0.00214735, 0.00171971,  0.00129118,  0.000861804, 0.000431112, 0},
  
       {0,          0.00067284, 0.00134879, 0.00202776, 0.00270996, 0.00339501, 0.00408276, 0.00477335, 0.00546651,
        0.00616208, 0.00686005, 0.00756033, 0.00826268, 0.00896702, 0.00967339, 0.0103815,  0.0110912,  0.0118026,
        0.0125154,  0.0132295,  0.0139448,  0.0146612,  0.0153786,  0.0160968,  0.0168158,  0.0175354,  0.0182555,
        0.0189759,  0.0196967,  0.0204176,  0.0211385,  0.0218593,  0.02258,    0.0233003,  0.0240202,  0.0247395,
        0.0254582,  0.0261762,  0.0268932,  0.0276093,  0.0283243,  0.0290381,  0.0297505,  0.0304615,  0.0311711,
        0.031879,   0.0325851,  0.0332894,  0.0339919,  0.0346922,  0.0353905,  0.0360866,  0.0367803,  0.0374715,
        0.0381604,  0.0388467,  0.0395302,  0.0402111,  0.0408891,  0.0415641,  0.0422361,  0.0429052,  0.0435709,
        0.0442335,  0.0448928,  0.0455486,  0.046201,   0.0468499,  0.0474952,  0.0481369,  0.0487748,  0.0494089,
        0.0500393,  0.0506655,  0.051288,   0.0519065,  0.0525207,  0.053131,   0.053737,   0.0543388,  0.0549365,
        0.0555296,  0.0561186,  0.056703,   0.057283,   0.0578587,  0.0584296,  0.0589963,  0.059558,   0.0601155,
        0.060668,   0.061216,   0.0617592,  0.0622977,  0.0628315,  0.0633604,  0.0638846,  0.0644038,  0.0649184,
        0.0654278,  0.0659327,  0.0664323,  0.0669273,  0.067417,   0.0679021,  0.0683819,  0.0688571,  0.069327,
        0.0697921,  0.0702521,  0.070707,   0.0711571,  0.0716019,  0.0720421,  0.0724768,  0.0729068,  0.0733315,
        0.0737513,  0.0741662,  0.0745758,  0.0749807,  0.0753801,  0.0757748,  0.0761645,  0.0765489,  0.0769287,
        0.0773031,  0.0776727,  0.0780375,  0.0783968,  0.0787515,  0.0791012,  0.0794457,  0.0797855,  0.0801204,
        0.08045,    0.080775,   0.0810952,  0.0814101,  0.0817204,  0.082026,   0.0823263,  0.0826219,  0.082913,
        0.0831989,  0.08348,    0.0837566,  0.0840285,  0.0842952,  0.0845573,  0.084815,   0.0850679,  0.0853158,
        0.0855591,  0.0857981,  0.0860324,  0.0862617,  0.0864865,  0.0867069,  0.0869229,  0.0871341,  0.0873407,
        0.0875429,  0.0877407,  0.0879342,  0.0881229,  0.0883071,  0.088487,   0.0886627,  0.088834,   0.0890011,
        0.0891634,  0.0893215,  0.0894754,  0.0896251,  0.0897706,  0.089912,   0.090049,   0.0901816,  0.0903101,
        0.0904345,  0.0905549,  0.0906711,  0.0907833,  0.0908915,  0.0909957,  0.0910955,  0.0911914,  0.0912833,
        0.0913713,  0.0914553,  0.0915355,  0.0916117,  0.0916841,  0.0917526,  0.0918173,  0.0918782,  0.0919353,
        0.0919886,  0.092038,   0.0920836,  0.0921255,  0.0921637,  0.0921982,  0.0922291,  0.0922563,  0.0922799,
        0.0922999,  0.0923163,  0.0923292,  0.0923385,  0.0923442,  0.0923465,  0.0923452,  0.0923405,  0.0923323,
        0.0923207,  0.0923056,  0.0922872,  0.0922654,  0.0922402,  0.0922117,  0.0921799,  0.0921447,  0.0921063,
        0.0920646,  0.0920197,  0.0919714,  0.0919199,  0.0918652,  0.0918073,  0.0917462,  0.091682,   0.0916146,
        0.0915442,  0.0914707,  0.0913941,  0.0913145,  0.0912319,  0.0911462,  0.0910576,  0.0909659,  0.0908711,
        0.0907733,  0.0906725,  0.0905689,  0.0904624,  0.0903531,  0.0902409,  0.0901259,  0.0900081,  0.0898871,
        0.0897633,  0.0896368,  0.0895075,  0.0893755,  0.0892408,  0.0891035,  0.0889633,  0.0888202,  0.0886744,
        0.088526,   0.088375,   0.0882214,  0.0880653,  0.0879063,  0.0877446,  0.0875803,  0.0874136,  0.0872444,
        0.0870727,  0.0868983,  0.0867211,  0.0865416,  0.0863597,  0.0861754,  0.0859887,  0.0857991,  0.0856071,
        0.0854128,  0.0852163,  0.0850175,  0.0848158,  0.0846117,  0.0844055,  0.084197,   0.0839864,  0.0837728,
        0.0835571,  0.0833392,  0.0831192,  0.0828968,  0.0826718,  0.0824447,  0.0822156,  0.0819844,  0.0817504,
        0.0815143,  0.0812763,  0.0810363,  0.0807935,  0.0805487,  0.0803019,  0.0800533,  0.079802,   0.0795486,
        0.0792933,  0.0790362,  0.0787765,  0.0785147,  0.0782512,  0.0779858,  0.0777177,  0.0774478,  0.0771762,
        0.0769024,  0.0766263,  0.0763485,  0.076069,   0.075787,   0.0755031,  0.0752175,  0.0749301,  0.0746402,
        0.0743487,  0.0740556,  0.07376,    0.0734626,  0.0731637,  0.0728628,  0.0725597,  0.0722551,  0.0719488,
        0.0716401,  0.0713298,  0.0710182,  0.070704,   0.0703882,  0.0700711,  0.0697516,  0.0694305,  0.069108,
        0.0687833,  0.0684569,  0.0681291,  0.0677993,  0.0674676,  0.0671347,  0.0667997,  0.0664629,  0.0661248,
        0.0657847,  0.0654429,  0.0650999,  0.0647547,  0.064408,   0.0640601,  0.0637099,  0.0633583,  0.0630055,
        0.0626503,  0.062294,   0.0619363,  0.0615764,  0.0612154,  0.0608528,  0.0604883,  0.0601227,  0.0597552,
        0.0593861,  0.0590161,  0.0586437,  0.0582702,  0.0578954,  0.0575185,  0.0571407,  0.056761,   0.0563799,
        0.0559978,  0.0556134,  0.0552281,  0.0548412,  0.0544526,  0.0540632,  0.0536716,  0.053279,   0.0528851,
        0.0524892,  0.0520926,  0.051694,   0.0512943,  0.0508933,  0.0504905,  0.0500869,  0.0496813,  0.0492747,
        0.0488668,  0.0484572,  0.0480468,  0.0476344,  0.047221,   0.0468063,  0.0463901,  0.0459731,  0.0455539,
        0.0451341,  0.0447126,  0.04429,    0.0438662,  0.0434407,  0.0430147,  0.0425864,  0.0421576,  0.0417272,
        0.0412956,  0.0408629,  0.0404286,  0.0399937,  0.0395567,  0.0391192,  0.03868,    0.0382398,  0.0377984,
        0.0373556,  0.036912,   0.0364666,  0.0360208,  0.0355729,  0.0351245,  0.0346744,  0.0342235,  0.0337713,
        0.0333179,  0.0328636,  0.0324077,  0.0319513,  0.0314929,  0.0310342,  0.0305736,  0.0301125,  0.0296498,
        0.0291863,  0.0287215,  0.0282556,  0.0277887,  0.0273206,  0.0268516,  0.0263811,  0.0259101,  0.0254373,
        0.0249642,  0.0244892,  0.024014,   0.0235368,  0.0230593,  0.0225801,  0.0221005,  0.0216193,  0.0211375,
        0.0206542,  0.0201703,  0.019685,   0.0191989,  0.0187117,  0.0182235,  0.0177342,  0.017244,   0.0167527,
        0.0162605,  0.0157672,  0.0152729,  0.0147777,  0.0142814,  0.0137843,  0.013286,   0.0127868,  0.0122866,
        0.0117855,  0.0112834,  0.0107803,  0.0102763,  0.00977131, 0.00926534, 0.00875847, 0.00825062, 0.00774184,
        0.00723215, 0.00672146, 0.00620995, 0.00569735, 0.00518404, 0.00466955, 0.00415446, 0.00363809, 0.00312124,
        0.00260299, 0.0020844,  0.00156428, 0.00104398, 0.00052199, 0},
  
       {0,          0.00109345, 0.00219081, 0.00329183, 0.00439625, 0.00550381, 0.00661425, 0.0077273,  0.00884269,
        0.00996015, 0.0110795,  0.0122003,  0.0133224,  0.0144455,  0.0155693,  0.0166936,  0.017818,   0.0189424,
        0.0200664,  0.0211898,  0.0223123,  0.0234336,  0.0245535,  0.0256718,  0.0267882,  0.0279024,  0.0290142,
        0.0301233,  0.0312294,  0.0323325,  0.0334321,  0.0345282,  0.0356204,  0.0367086,  0.0377926,  0.0388721,
        0.0399469,  0.0410168,  0.0420815,  0.0431411,  0.0441951,  0.0452436,  0.0462862,  0.0473229,  0.0483534,
        0.0493777,  0.0503955,  0.0514067,  0.0524112,  0.0534088,  0.0543994,  0.0553829,  0.0563592,  0.057328,
        0.0582895,  0.0592434,  0.0601896,  0.061128,   0.0620586,  0.0629813,  0.0638959,  0.0648023,  0.0657006,
        0.0665907,  0.0674726,  0.0683461,  0.0692114,  0.0700682,  0.0709164,  0.071756,   0.072587,   0.0734096,
        0.0742236,  0.0750291,  0.0758256,  0.0766134,  0.0773927,  0.0781634,  0.0789253,  0.0796782,  0.0804226,
        0.0811584,  0.0818854,  0.0826035,  0.083313,   0.0840142,  0.084706,   0.0853895,  0.0860646,  0.0867305,
        0.087388,   0.0880373,  0.0886774,  0.0893093,  0.0899329,  0.0905475,  0.091154,   0.0917519,  0.0923414,
        0.0929229,  0.0934955,  0.0940604,  0.0946164,  0.0951647,  0.0957045,  0.0962364,  0.0967601,  0.0972757,
        0.0977834,  0.098283,   0.0987748,  0.0992587,  0.0997347,  0.100203,   0.100663,   0.101116,   0.101561,
        0.101999,   0.102429,   0.102852,   0.103266,   0.103674,   0.104074,   0.104467,   0.104853,   0.105231,
        0.105602,   0.105967,   0.106324,   0.106674,   0.107017,   0.107353,   0.107682,   0.108005,   0.10832,
        0.108629,   0.108931,   0.109226,   0.109515,   0.109798,   0.110074,   0.110343,   0.110606,   0.110863,
        0.111113,   0.111357,   0.111595,   0.111827,   0.112053,   0.112272,   0.112486,   0.112693,   0.112895,
        0.113091,   0.11328,    0.113464,   0.113643,   0.113815,   0.113982,   0.114143,   0.114299,   0.114449,
        0.114594,   0.114733,   0.114867,   0.114995,   0.115118,   0.115236,   0.115348,   0.115455,   0.115558,
        0.115655,   0.115746,   0.115833,   0.115915,   0.115992,   0.116064,   0.116131,   0.116193,   0.11625,
        0.116302,   0.11635,    0.116393,   0.116431,   0.116465,   0.116494,   0.116519,   0.116539,   0.116554,
        0.116565,   0.116571,   0.116574,   0.116571,   0.116565,   0.116554,   0.116539,   0.116519,   0.116495,
        0.116468,   0.116436,   0.1164,     0.116359,   0.116315,   0.116267,   0.116215,   0.116159,   0.116099,
        0.116035,   0.115967,   0.115895,   0.115819,   0.115739,   0.115656,   0.115569,   0.115478,   0.115384,
        0.115286,   0.115185,   0.11508,    0.114971,   0.114858,   0.114742,   0.114623,   0.1145,     0.114373,
        0.114244,   0.114111,   0.113974,   0.113834,   0.113691,   0.113544,   0.113394,   0.113241,   0.113085,
        0.112926,   0.112763,   0.112596,   0.112427,   0.112255,   0.11208,    0.111902,   0.11172,    0.111535,
        0.111348,   0.111157,   0.110964,   0.110767,   0.110568,   0.110365,   0.11016,    0.109952,   0.109741,
        0.109527,   0.10931,    0.109091,   0.108869,   0.108644,   0.108416,   0.108185,   0.107953,   0.107716,
        0.107477,   0.107236,   0.106993,   0.106746,   0.106497,   0.106245,   0.105991,   0.105734,   0.105474,
        0.105213,   0.104948,   0.104681,   0.104412,   0.10414,    0.103866,   0.103589,   0.103309,   0.103028,
        0.102744,   0.102458,   0.102169,   0.101878,   0.101584,   0.101289,   0.100991,   0.100691,   0.100388,
        0.100084,   0.0997766,  0.0994672,  0.099156,   0.0988424,  0.0985263,  0.0982085,  0.0978886,  0.0975659,
        0.0972416,  0.0969154,  0.0965863,  0.0962557,  0.0959232,  0.0955879,  0.095251,   0.0949123,  0.0945709,
        0.094228,   0.0938831,  0.0935357,  0.0931869,  0.0928359,  0.0924827,  0.0921281,  0.0917711,  0.0914122,
        0.0910519,  0.0906889,  0.0903244,  0.0899584,  0.0895897,  0.0892197,  0.0888478,  0.0884737,  0.0880985,
        0.0877207,  0.0873414,  0.0869608,  0.0865774,  0.086193,   0.0858066,  0.0854182,  0.0850288,  0.0846366,
        0.0842434,  0.0838484,  0.0834513,  0.0830532,  0.0826526,  0.0822508,  0.0818473,  0.0814418,  0.0810354,
        0.0806264,  0.0802164,  0.0798046,  0.079391,   0.0789765,  0.0785593,  0.0781414,  0.0777214,  0.0773,
        0.0768773,  0.0764524,  0.0760269,  0.0755987,  0.0751698,  0.074739,   0.0743068,  0.0738733,  0.0734378,
        0.0730016,  0.0725629,  0.0721235,  0.0716821,  0.0712396,  0.0707956,  0.0703499,  0.0699033,  0.0694545,
        0.0690053,  0.0685534,  0.0681012,  0.0676468,  0.0671915,  0.0667346,  0.0662763,  0.0658168,  0.0653557,
        0.0648936,  0.0644296,  0.063965,   0.0634981,  0.0630309,  0.0625614,  0.0620914,  0.0616193,  0.0611466,
        0.060672,   0.0601966,  0.0597196,  0.0592414,  0.058762,   0.0582812,  0.0577993,  0.0573159,  0.0568316,
        0.0563456,  0.0558588,  0.0553703,  0.0548812,  0.0543902,  0.0538986,  0.0534051,  0.0529112,  0.0524153,
        0.0519189,  0.0514206,  0.0509219,  0.0504213,  0.0499202,  0.0494172,  0.0489138,  0.0484085,  0.0479028,
        0.0473952,  0.0468872,  0.0463774,  0.045867,   0.045355,   0.0448423,  0.0443281,  0.0438132,  0.0432969,
        0.0427797,  0.0422612,  0.0417417,  0.0412212,  0.0406995,  0.0401769,  0.0396529,  0.0391284,  0.0386021,
        0.0380756,  0.0375471,  0.0370186,  0.036488,   0.0359573,  0.0354247,  0.0348918,  0.0343573,  0.0338222,
        0.0332859,  0.0327486,  0.0322104,  0.031671,   0.031131,   0.0305895,  0.0300477,  0.029504,   0.0289603,
        0.0284146,  0.0278688,  0.0273214,  0.0267734,  0.0262244,  0.0256743,  0.0251236,  0.0245714,  0.0240191,
        0.0234648,  0.0229105,  0.0223545,  0.0217981,  0.0212406,  0.0206822,  0.0201231,  0.0195626,  0.019002,
        0.0184394,  0.0178768,  0.0173127,  0.0167481,  0.0161826,  0.0156159,  0.015049,   0.0144803,  0.0139115,
        0.0133412,  0.0127705,  0.0121988,  0.0116261,  0.0110531,  0.0104783,  0.00990358, 0.00932729, 0.00875055,
        0.008173,   0.00759428, 0.00701549, 0.0064348,  0.00585411, 0.00527215, 0.00468949, 0.00410634, 0.00352174,
        0.00293714, 0.00235088, 0.00176434, 0.00117696, 0.00058848, 0},
  
       {0,          0.00162336, 0.00325125, 0.00488316, 0.00651853,  0.00815664, 0.00979694, 0.0114388,  0.0130815,
        0.0147246,  0.0163673,  0.0180091,  0.0196493,  0.0212874,   0.0229227,  0.0245546,  0.0261826,  0.0278061,
        0.0294245,  0.0310375,  0.032644,   0.034244,   0.0358369,   0.0374219,  0.038999,   0.0405672,  0.0421266,
        0.0436763,  0.0452161,  0.0467457,  0.0482645,  0.0497723,   0.0512687,  0.0527534,  0.0542261,  0.0556864,
        0.0571342,  0.0585689,  0.059991,   0.0613992,  0.0627942,   0.0641754,  0.0655425,  0.0668959,  0.0682348,
        0.0695593,  0.0708697,  0.0721652,  0.0734459,  0.0747121,   0.0759637,  0.0772,     0.0784213,  0.0796279,
        0.0808195,  0.0819963,  0.0831575,  0.0843039,  0.0854353,   0.0865517,  0.0876532,  0.0887399,  0.0898116,
        0.0908685,  0.0919105,  0.0929378,  0.0939504,  0.0949484,   0.0959318,  0.0969007,  0.0978552,  0.0987953,
        0.0997213,  0.100633,   0.101531,   0.102415,   0.103285,    0.104142,   0.104985,   0.105814,   0.106629,
        0.107432,   0.108221,   0.108998,   0.109761,   0.110511,    0.111249,   0.111974,   0.112686,   0.113386,
        0.114074,   0.114749,   0.115413,   0.116065,   0.116704,    0.117332,   0.117949,   0.118553,   0.119147,
        0.119729,   0.1203,     0.120861,   0.121409,   0.121948,    0.122475,   0.122993,   0.123499,   0.123996,
        0.124481,   0.124958,   0.125423,   0.125879,   0.126325,    0.126761,   0.127188,   0.127604,   0.128012,
        0.12841,    0.128799,   0.129179,   0.12955,    0.129912,    0.130265,   0.130609,   0.130945,   0.131272,
        0.13159,    0.1319,     0.132202,   0.132495,   0.132781,    0.133059,   0.133328,   0.133589,   0.133843,
        0.134089,   0.134327,   0.134558,   0.134781,   0.134998,    0.135206,   0.135407,   0.135601,   0.135788,
        0.135969,   0.136141,   0.136307,   0.136466,   0.136619,    0.136765,   0.136904,   0.137037,   0.137163,
        0.137283,   0.137396,   0.137503,   0.137604,   0.137699,    0.137788,   0.137871,   0.137947,   0.138018,
        0.138083,   0.138142,   0.138195,   0.138242,   0.138284,    0.138321,   0.138351,   0.138377,   0.138397,
        0.138411,   0.13842,    0.138424,   0.138423,   0.138416,    0.138405,   0.138388,   0.138366,   0.13834,
        0.138308,   0.138272,   0.13823,    0.138184,   0.138133,    0.138078,   0.138018,   0.137953,   0.137883,
        0.137809,   0.13773,    0.137647,   0.137559,   0.137467,    0.137371,   0.137271,   0.137166,   0.137057,
        0.136944,   0.136826,   0.136705,   0.136579,   0.136449,    0.136316,   0.136178,   0.136037,   0.135891,
        0.135742,   0.135588,   0.135431,   0.135271,   0.135106,    0.134938,   0.134766,   0.13459,    0.134411,
        0.134228,   0.134042,   0.133852,   0.133658,   0.133461,    0.133261,   0.133057,   0.13285,    0.132639,
        0.132425,   0.132208,   0.131988,   0.131764,   0.131537,    0.131307,   0.131073,   0.130837,   0.130597,
        0.130354,   0.130108,   0.129859,   0.129607,   0.129352,    0.129094,   0.128833,   0.128569,   0.128302,
        0.128033,   0.12776,    0.127484,   0.127206,   0.126926,    0.126641,   0.126354,   0.126066,   0.125773,
        0.125478,   0.125181,   0.124881,   0.124577,   0.124272,    0.123964,   0.123653,   0.12334,    0.123024,
        0.122705,   0.122384,   0.122061,   0.121735,   0.121407,    0.121076,   0.120743,   0.120407,   0.120069,
        0.119729,   0.119386,   0.119041,   0.118693,   0.118343,    0.117992,   0.117637,   0.11728,    0.116922,
        0.11656,    0.116197,   0.115831,   0.115463,   0.115094,    0.114721,   0.114347,   0.113971,   0.113592,
        0.113212,   0.112829,   0.112444,   0.112058,   0.111668,    0.111278,   0.110885,   0.11049,    0.110093,
        0.109694,   0.109293,   0.10889,    0.108485,   0.108079,    0.10767,    0.107259,   0.106847,   0.106432,
        0.106017,   0.105598,   0.105178,   0.104756,   0.104333,    0.103908,   0.10348,    0.103052,   0.10262,
        0.102188,   0.101753,   0.101317,   0.10088,    0.10044,     0.0999993,  0.0995557,  0.0991114,  0.0986647,
        0.0982168,  0.0977671,  0.0973157,  0.0968629,  0.096408,    0.0959523,  0.0954938,  0.0950349,  0.0945733,
        0.094111,   0.0936464,  0.0931808,  0.0927133,  0.0922444,   0.0917739,  0.0913018,  0.0908285,  0.0903532,
        0.089877,   0.0893986,  0.0889195,  0.088438,   0.087956,    0.0874715,  0.0869867,  0.0864992,  0.0860115,
        0.0855211,  0.0850305,  0.0845373,  0.0840439,  0.0835479,   0.0830517,  0.0825529,  0.0820539,  0.0815524,
        0.0810507,  0.0805465,  0.080042,   0.0795351,  0.0790281,   0.0785184,  0.0780088,  0.0774965,  0.0769842,
        0.0764693,  0.0759543,  0.075437,   0.0749194,  0.0743995,   0.0738793,  0.073357,   0.0728342,  0.0723096,
        0.0717841,  0.0712572,  0.0707291,  0.0701999,  0.0696692,   0.0691378,  0.0686045,  0.0680709,  0.0675351,
        0.0669993,  0.066461,   0.0659226,  0.0653822,  0.0648413,   0.0642988,  0.0637555,  0.0632109,  0.062665,
        0.0621185,  0.0615702,  0.0610217,  0.0604709,  0.0599201,   0.0593672,  0.0588139,  0.0582592,  0.0577035,
        0.057147,   0.0565888,  0.0560305,  0.05547,    0.0549094,   0.0543469,  0.053784,   0.0532198,  0.0526544,
        0.0520886,  0.0515209,  0.0509531,  0.0503833,  0.0498132,   0.0492418,  0.0486693,  0.0480964,  0.0475216,
        0.0469468,  0.04637,    0.0457929,  0.0452147,  0.0446352,   0.0440555,  0.0434738,  0.0428921,  0.0423087,
        0.0417248,  0.04114,    0.0405538,  0.0399675,  0.0393792,   0.0387906,  0.038201,   0.0376103,  0.0370194,
        0.0364264,  0.0358334,  0.0352391,  0.0346439,  0.0340484,   0.033451,   0.0328536,  0.0322548,  0.0316552,
        0.0310553,  0.0304535,  0.0298517,  0.0292485,  0.0286446,   0.0280404,  0.0274343,  0.0268281,  0.0262208,
        0.0256125,  0.0250043,  0.0243939,  0.0237835,  0.0231721,   0.0225596,  0.0219471,  0.0213327,  0.0207181,
        0.0201029,  0.0194862,  0.0188694,  0.0182513,  0.0176325,   0.0170136,  0.0163927,  0.0157718,  0.0151501,
        0.0145272,  0.0139042,  0.0132797,  0.0126547,  0.0120296,   0.0114025,  0.0107755,  0.0101476,  0.00951858,
        0.00888951, 0.00825902, 0.00762794, 0.00699685, 0.00636378,  0.00573069, 0.00509707, 0.00446199, 0.00382691,
        0.00319078, 0.00255372, 0.00191665, 0.00127807, 0.000639035, 0}};
  
   for (int i = 0; i < n_SW_noME1_1; i++) {
     for (int j = 0; j < n_val; j++) {
       //---- maybe a vector later...
       x_correction_noME1_1[i][j] = corrections_noME1_1[i][j];
       if (i < n_SW_ME1_1) {
         x_correction_ME1_1[i][j] = corrections_ME1_1[i][j];
       }
     }
   }
 }

References mps_fire::i, dqmiolumiharvest::j, n_SW_ME1_1, n_SW_noME1_1, n_val, x_correction_ME1_1, and x_correction_noME1_1.

Referenced by getCorrectionValues().

◆ initChamberSpecs()

void CSCXonStrip_MatchGatti::initChamberSpecs ( )

Use specs to setup Gatti parameters.

Definition at line 474 of file CSCXonStrip_MatchGatti.cc.

                                               {
   // Not used directly but these are parameters used for extracting the correction values
   // in coordinate and error estimators
  
   // Distance between anode and cathode
   h = specs_->anodeCathodeSpacing();
   r = h / stripWidth;
  
   // Wire spacing
   double wspace = specs_->wireSpacing();
  
   // Wire radius
   double wradius = specs_->wireRadius();
  
   // Accepted parameters in Gatti function
   const double parm[5] = {.1989337e-02, -.6901542e-04, .8665786, 154.6177, -.680163e-03};
  
   k_3 = (parm[0] * wspace / h + parm[1]) *
         (parm[2] * wspace / wradius + parm[3] + parm[4] * (wspace / wradius) * (wspace / wradius));
  
   sqrt_k_3 = std::sqrt(k_3);
   norm = r * (0.5 / std::atan(sqrt_k_3));  // changed from norm to r * norm
   k_2 = M_PI_2 * (1. - sqrt_k_3 / 2.);
   k_1 = 0.25 * k_2 * sqrt_k_3 / std::atan(sqrt_k_3);
 }

References CSCChamberSpecs::anodeCathodeSpacing(), k_1, k_2, k_3, M_PI_2, norm, r, specs_, mathSSE::sqrt(), sqrt_k_3, stripWidth, CSCChamberSpecs::wireRadius(), and CSCChamberSpecs::wireSpacing().

◆ operator=()

CSCXonStrip_MatchGatti& CSCXonStrip_MatchGatti::operator= ( const CSCXonStrip_MatchGatti & )

private

◆ setConditions()

void CSCXonStrip_MatchGatti::setConditions ( const CSCRecoConditions * reco )

inline

Cache pointer to conditions data.

Definition at line 51 of file CSCXonStrip_MatchGatti.h.

51 { recoConditions_ = reco; }

References dt_dqm_sourceclient_common_cff::reco, and recoConditions_.

Referenced by CSCMake2DRecHit::setConditions().

◆ setupMatrix()

void CSCXonStrip_MatchGatti::setupMatrix ( )

Set matrix for XT corrections and noise.

Definition at line 358 of file CSCXonStrip_MatchGatti.cc.

                                          {
   //---- a??? and v??[] could be skipped for now...; not used yet
  
   /*
   double dd, a11t, a12t, a13t, a22t, a23t, a33t;
   double syserr = adcSystematics;
   double xtlk_err = xtalksSystematics;
   // Left strip
   a11t = a11[0] + syserr*syserr * ChargeSignal[0][0]*ChargeSignal[0][0] + xtlk_err*xtlk_err*ChargeSignal[1][0]*ChargeSignal[1][0];
   a12t = a12[0] + syserr*syserr * ChargeSignal[0][0]*ChargeSignal[0][1];
   a13t = a13[0] + syserr*syserr * ChargeSignal[0][0]*ChargeSignal[0][2];
   a22t = a22[0] + syserr*syserr * ChargeSignal[0][1]*ChargeSignal[0][1] + xtlk_err*xtlk_err*ChargeSignal[1][1]*ChargeSignal[1][1];
   a23t = a23[0] + syserr*syserr * ChargeSignal[0][1]*ChargeSignal[0][2];
   a33t = a33[0] + syserr*syserr * ChargeSignal[0][2]*ChargeSignal[0][2] + xtlk_err*xtlk_err*ChargeSignal[1][2]*ChargeSignal[1][2];
  
   dd     = (a11t*a33t*a22t - a11t*a23t*a23t - a33t*a12t*a12t 
                        + 2.* a12t*a13t*a23t - a13t*a13t*a22t );
  
   v11[0] = (a33t*a22t - a23t*a23t)/dd;
   v12[0] =-(a33t*a12t - a13t*a23t)/dd;
   v13[0] = (a12t*a23t - a13t*a22t)/dd;
   v22[0] = (a33t*a11t - a13t*a13t)/dd;
   v23[0] =-(a23t*a11t - a12t*a13t)/dd;
   v33[0] = (a22t*a11t - a12t*a12t)/dd;
      
   // Center strip
   a11t = a11[1] + syserr*syserr * ChargeSignal[1][0]*ChargeSignal[1][0] + xtlk_err*xtlk_err*(ChargeSignal[0][0]*ChargeSignal[0][0]+ChargeSignal[2][0]*ChargeSignal[2][0]);
   a12t = a12[1] + syserr*syserr * ChargeSignal[1][0]*ChargeSignal[1][1];
   a13t = a13[1] + syserr*syserr * ChargeSignal[1][0]*ChargeSignal[1][2];
   a22t = a22[1] + syserr*syserr * ChargeSignal[1][1]*ChargeSignal[1][1] + xtlk_err*xtlk_err*(ChargeSignal[0][1]*ChargeSignal[0][1]+ChargeSignal[2][1]*ChargeSignal[2][1]);
   a23t = a23[1] + syserr*syserr * ChargeSignal[1][1]*ChargeSignal[1][2];
   a33t = a33[1] + syserr*syserr * ChargeSignal[1][2]*ChargeSignal[1][2] + xtlk_err*xtlk_err*(ChargeSignal[0][2]*ChargeSignal[0][2]+ChargeSignal[2][2]*ChargeSignal[2][2]);
  
   dd     = (a11t*a33t*a22t - a11t*a23t*a23t - a33t*a12t*a12t
                        + 2.* a12t*a13t*a23t - a13t*a13t*a22t );
  
   v11[1] = (a33t*a22t - a23t*a23t)/dd;
   v12[1] =-(a33t*a12t - a13t*a23t)/dd;
   v13[1] = (a12t*a23t - a13t*a22t)/dd;
   v22[1] = (a33t*a11t - a13t*a13t)/dd;
   v23[1] =-(a23t*a11t - a12t*a13t)/dd;
   v33[1] = (a22t*a11t - a12t*a12t)/dd;
  
   // Right strip
   a11t = a11[2] + syserr*syserr * ChargeSignal[2][0]*ChargeSignal[2][0] + xtlk_err*xtlk_err*ChargeSignal[1][0]*ChargeSignal[1][0];
   a12t = a12[2] + syserr*syserr * ChargeSignal[2][0]*ChargeSignal[2][1];
   a13t = a13[2] + syserr*syserr * ChargeSignal[2][0]*ChargeSignal[2][2];
   a22t = a22[2] + syserr*syserr * ChargeSignal[2][1]*ChargeSignal[2][1] + xtlk_err*xtlk_err*ChargeSignal[1][1]*ChargeSignal[1][1];
   a23t = a23[2] + syserr*syserr * ChargeSignal[2][1]*ChargeSignal[2][2];
   a33t = a33[2] + syserr*syserr * ChargeSignal[2][2]*ChargeSignal[2][2] + xtlk_err*xtlk_err*ChargeSignal[1][2]*ChargeSignal[1][2];
  
   dd     = (a11t*a33t*a22t - a11t*a23t*a23t - a33t*a12t*a12t
                         +2.* a12t*a13t*a23t - a13t*a13t*a22t );
  
   v11[2] = (a33t*a22t - a23t*a23t)/dd;
   v12[2] =-(a33t*a12t - a13t*a23t)/dd;
   v13[2] = (a12t*a23t - a13t*a22t)/dd;
   v22[2] = (a33t*a11t - a13t*a13t)/dd;
   v23[2] =-(a23t*a11t - a12t*a13t)/dd;
   v33[2] = (a22t*a11t - a12t*a12t)/dd;
 */
   //---- Find the inverted XTalk matrix and apply it to the charge (3x3)
   //---- Thus the charge before the XTalk is obtained
   CLHEP::HepMatrix cross_talks_inv(3, 3);
   int err = 0;
   //---- q_sum is 3 time bins summed; L, C, R - left, central, right strips
   q_sum = q_sumL = q_sumC = q_sumR = 0.;
   double charge = 0.;
   for (int iTime = 0; iTime < 3; iTime++) {
     cross_talks_inv(1, 1) = xt_lr0[iTime];
     cross_talks_inv(1, 2) = xt_l[1][iTime];
     cross_talks_inv(1, 3) = 0.;
     cross_talks_inv(2, 1) = xt_r[0][iTime];
     cross_talks_inv(2, 2) = xt_lr1[iTime];
     cross_talks_inv(2, 3) = xt_l[2][iTime];
     cross_talks_inv(3, 1) = 0.;
     cross_talks_inv(3, 2) = xt_r[1][iTime];
     cross_talks_inv(3, 3) = xt_lr2[iTime];
     cross_talks_inv.invert(err);
     if (err != 0) {
       edm::LogWarning("FailedXTalkiInversionNoCrosstalkCorrection")
           << "Failed to invert XTalks matrix. No cross-talk correction for this rechit.";
       //edm::LogError("CSCRecHit") << "Failed to invert XTalks matrix. No cross-talk correction for this rechit.";
       return;
     }
     //---- "charge" is XT-corrected charge
     charge = chargeSignal[0][iTime] * cross_talks_inv(1, 1) + chargeSignal[1][iTime] * cross_talks_inv(1, 2) +
              chargeSignal[2][iTime] * cross_talks_inv(1, 3);
     //---- Negative charge? According to studies (and logic) - better use 0 charge
     //----- Later studies suggest that this only do harm. I am still worried about
     // charges of -50 ADC and below (0.5% of the cases) but let see
     //if(charge<0.){
     //charge = 0.;
     //}
     q_sum += charge;
     q_sumL += charge;
     charge = chargeSignal[0][iTime] * cross_talks_inv(2, 1) + chargeSignal[1][iTime] * cross_talks_inv(2, 2) +
              chargeSignal[2][iTime] * cross_talks_inv(2, 3);
     //if(charge<0.){
     //charge = 0.;
     //}
     q_sum += charge;
     q_sumC += charge;
     charge = chargeSignal[0][iTime] * cross_talks_inv(3, 1) + chargeSignal[1][iTime] * cross_talks_inv(3, 2) +
              chargeSignal[2][iTime] * cross_talks_inv(3, 3);
     //if(charge<0.){
     //charge = 0.;
     //}
     q_sum += charge;
     q_sumR += charge;
   }
 }

References ALCARECOTkAlJpsiMuMu_cff::charge, chargeSignal, submitPVResolutionJobs::err, q_sum, q_sumC, q_sumL, q_sumR, xt_l, xt_lr0, xt_lr1, xt_lr2, and xt_r.

Referenced by findXOnStrip().

◆ xfError_Noise()

double CSCXonStrip_MatchGatti::xfError_Noise ( double noise )

private

Definition at line 587 of file CSCXonStrip_MatchGatti.cc.

                                                          {
   double min, max;
   if (q_sumR > q_sumL) {
     min = q_sumL;
     max = q_sumR;
   } else {
     min = q_sumR;
     max = q_sumL;
   }
   //---- Error propagation...
   //---- Names here are fake! Due to technical features
   double dr_L2 = pow(q_sumR - q_sumL, 2);
   double dr_C2 = pow(q_sumC - min, 2);
   double dr_R2 = pow(q_sumC - max, 2);
   double error = std::sqrt(dr_L2 + dr_C2 + dr_R2) * noise / std::pow(q_sumC - min, 2) / 2;
  
   return error;
 }

References relativeConstraints::error, SiStripPI::max, min(), gpuVertexFinder::noise, funct::pow(), q_sumC, q_sumL, q_sumR, and mathSSE::sqrt().

Referenced by calculateXonStripError().

◆ xfError_XTasym()

double CSCXonStrip_MatchGatti::xfError_XTasym ( double XTasym )

private

Definition at line 606 of file CSCXonStrip_MatchGatti.cc.

                                                             {
   double min;
   if (q_sumR > q_sumL) {
     min = q_sumL;
   } else {
     min = q_sumR;
   }
   //---- Error propagation
   double dXTL = (std::pow(q_sumC, 2) + std::pow(q_sumR, 2) - q_sumL * q_sumR - q_sumR * q_sumC);
   double dXTR = (std::pow(q_sumC, 2) + std::pow(q_sumL, 2) - q_sumL * q_sumR - q_sumL * q_sumC);
   double dXT = std::sqrt(std::pow(dXTL, 2) + std::pow(dXTR, 2)) / std::pow((q_sumC - min), 2) / 2;
   double error = dXT * xt_asym;
  
   return error;
 }

References relativeConstraints::error, min(), funct::pow(), q_sumC, q_sumL, q_sumR, and mathSSE::sqrt().

Referenced by calculateXonStripError().

Member Data Documentation

◆ a11

float CSCXonStrip_MatchGatti::a11[3]

private

Store elements of auto-correlation matrices: 0 = left, 1 = middle, 2 = right.

Definition at line 75 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip().

◆ a12

float CSCXonStrip_MatchGatti::a12[3]

private

Definition at line 75 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip().

◆ a13

float CSCXonStrip_MatchGatti::a13[3]

private

Definition at line 75 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip().

◆ a22

float CSCXonStrip_MatchGatti::a22[3]

private

Definition at line 75 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip().

◆ a23

float CSCXonStrip_MatchGatti::a23[3]

private

Definition at line 75 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip().

◆ a33

float CSCXonStrip_MatchGatti::a33[3]

private

Definition at line 75 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip().

◆ chargeSignal

float CSCXonStrip_MatchGatti::chargeSignal[3][3]

private

Definition at line 65 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip(), and setupMatrix().

◆ const_syst

float CSCXonStrip_MatchGatti::const_syst

private

Definition at line 116 of file CSCXonStrip_MatchGatti.h.

Referenced by calculateXonStripError(), and findXOnStrip().

◆ const_syst_ME12

float CSCXonStrip_MatchGatti::const_syst_ME12

private

Definition at line 126 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ const_syst_ME13

float CSCXonStrip_MatchGatti::const_syst_ME13

private

Definition at line 129 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ const_syst_ME1a

float CSCXonStrip_MatchGatti::const_syst_ME1a

private

Definition at line 120 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ const_syst_ME1b

float CSCXonStrip_MatchGatti::const_syst_ME1b

private

Definition at line 123 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ const_syst_ME21

float CSCXonStrip_MatchGatti::const_syst_ME21

private

Definition at line 132 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ const_syst_ME22

float CSCXonStrip_MatchGatti::const_syst_ME22

private

Definition at line 135 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ const_syst_ME31

float CSCXonStrip_MatchGatti::const_syst_ME31

private

Definition at line 138 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ const_syst_ME32

float CSCXonStrip_MatchGatti::const_syst_ME32

private

Definition at line 141 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ const_syst_ME41

float CSCXonStrip_MatchGatti::const_syst_ME41

private

Definition at line 144 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ h

double CSCXonStrip_MatchGatti::h

private

Definition at line 58 of file CSCXonStrip_MatchGatti.h.

◆ k_1

double CSCXonStrip_MatchGatti::k_1

private

Definition at line 62 of file CSCXonStrip_MatchGatti.h.

Referenced by initChamberSpecs().

◆ k_2

double CSCXonStrip_MatchGatti::k_2

private

Definition at line 62 of file CSCXonStrip_MatchGatti.h.

Referenced by initChamberSpecs().

◆ k_3

double CSCXonStrip_MatchGatti::k_3

private

Definition at line 62 of file CSCXonStrip_MatchGatti.h.

Referenced by initChamberSpecs().

◆ n_SW_ME1_1

const int CSCXonStrip_MatchGatti::n_SW_ME1_1 = 6

staticprivate

Definition at line 107 of file CSCXonStrip_MatchGatti.h.

Referenced by estimated2GattiCorrection(), and hardcodedCorrectionInitialization().

◆ n_SW_noME1_1

const int CSCXonStrip_MatchGatti::n_SW_noME1_1 = 11

staticprivate

Definition at line 106 of file CSCXonStrip_MatchGatti.h.

Referenced by estimated2GattiCorrection(), and hardcodedCorrectionInitialization().

◆ n_val

const int CSCXonStrip_MatchGatti::n_val = 501

staticprivate

Definition at line 108 of file CSCXonStrip_MatchGatti.h.

Referenced by estimated2GattiCorrection(), and hardcodedCorrectionInitialization().

◆ noise_level

float CSCXonStrip_MatchGatti::noise_level

private

Definition at line 114 of file CSCXonStrip_MatchGatti.h.

Referenced by calculateXonStripError(), and findXOnStrip().

◆ noise_level_ME12

float CSCXonStrip_MatchGatti::noise_level_ME12

private

Definition at line 124 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ noise_level_ME13

float CSCXonStrip_MatchGatti::noise_level_ME13

private

Definition at line 127 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ noise_level_ME1a

float CSCXonStrip_MatchGatti::noise_level_ME1a

private

Definition at line 118 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ noise_level_ME1b

float CSCXonStrip_MatchGatti::noise_level_ME1b

private

Definition at line 121 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ noise_level_ME21

float CSCXonStrip_MatchGatti::noise_level_ME21

private

Definition at line 130 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ noise_level_ME22

float CSCXonStrip_MatchGatti::noise_level_ME22

private

Definition at line 133 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ noise_level_ME31

float CSCXonStrip_MatchGatti::noise_level_ME31

private

Definition at line 136 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ noise_level_ME32

float CSCXonStrip_MatchGatti::noise_level_ME32

private

Definition at line 139 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ noise_level_ME41

float CSCXonStrip_MatchGatti::noise_level_ME41

private

Definition at line 142 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ norm

double CSCXonStrip_MatchGatti::norm

private

Definition at line 62 of file CSCXonStrip_MatchGatti.h.

Referenced by initChamberSpecs().

◆ q_sum

double CSCXonStrip_MatchGatti::q_sum

private

Definition at line 82 of file CSCXonStrip_MatchGatti.h.

Referenced by setupMatrix().

◆ q_sumC

double CSCXonStrip_MatchGatti::q_sumC

private

Definition at line 82 of file CSCXonStrip_MatchGatti.h.

Referenced by calculateXonStripError(), calculateXonStripPosition(), findXOnStrip(), setupMatrix(), xfError_Noise(), and xfError_XTasym().

◆ q_sumL

double CSCXonStrip_MatchGatti::q_sumL

private

Definition at line 82 of file CSCXonStrip_MatchGatti.h.

Referenced by calculateXonStripError(), calculateXonStripPosition(), findXOnStrip(), setupMatrix(), xfError_Noise(), and xfError_XTasym().

◆ q_sumR

double CSCXonStrip_MatchGatti::q_sumR

private

Definition at line 82 of file CSCXonStrip_MatchGatti.h.

Referenced by calculateXonStripError(), calculateXonStripPosition(), findXOnStrip(), setupMatrix(), xfError_Noise(), and xfError_XTasym().

◆ r

double CSCXonStrip_MatchGatti::r

private

Definition at line 60 of file CSCXonStrip_MatchGatti.h.

Referenced by initChamberSpecs().

◆ recoConditions_

const CSCRecoConditions* CSCXonStrip_MatchGatti::recoConditions_

private

Definition at line 90 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip(), and setConditions().

◆ specs_

const CSCChamberSpecs* CSCXonStrip_MatchGatti::specs_

private

Definition at line 78 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip(), and initChamberSpecs().

◆ sqrt_k_3

double CSCXonStrip_MatchGatti::sqrt_k_3

private

Definition at line 62 of file CSCXonStrip_MatchGatti.h.

Referenced by initChamberSpecs().

◆ stripWidth

float CSCXonStrip_MatchGatti::stripWidth

private

Definition at line 59 of file CSCXonStrip_MatchGatti.h.

Referenced by calculateXonStripError(), calculateXonStripPosition(), estimated2Gatti(), estimated2GattiCorrection(), findXOnStrip(), and initChamberSpecs().

◆ use3TimeBins

bool CSCXonStrip_MatchGatti::use3TimeBins

private

Definition at line 86 of file CSCXonStrip_MatchGatti.h.

◆ useCalib

bool CSCXonStrip_MatchGatti::useCalib

private

Definition at line 85 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ v11

float CSCXonStrip_MatchGatti::v11[3]

private

Store elements of matrices for chi^2 computation: 0 = left, 1 = middle, 2 = right.

Definition at line 72 of file CSCXonStrip_MatchGatti.h.

◆ v12

float CSCXonStrip_MatchGatti::v12[3]

private

Definition at line 72 of file CSCXonStrip_MatchGatti.h.

◆ v13

float CSCXonStrip_MatchGatti::v13[3]

private

Definition at line 72 of file CSCXonStrip_MatchGatti.h.

◆ v22

float CSCXonStrip_MatchGatti::v22[3]

private

Definition at line 72 of file CSCXonStrip_MatchGatti.h.

◆ v23

float CSCXonStrip_MatchGatti::v23[3]

private

Definition at line 72 of file CSCXonStrip_MatchGatti.h.

◆ v33

float CSCXonStrip_MatchGatti::v33[3]

private

Definition at line 72 of file CSCXonStrip_MatchGatti.h.

◆ x_centralVal

float CSCXonStrip_MatchGatti::x_centralVal[n_val]

private

Definition at line 112 of file CSCXonStrip_MatchGatti.h.

◆ x_correction_ME1_1

float CSCXonStrip_MatchGatti::x_correction_ME1_1[n_SW_ME1_1][n_val]

private

Definition at line 111 of file CSCXonStrip_MatchGatti.h.

Referenced by estimated2GattiCorrection(), and hardcodedCorrectionInitialization().

◆ x_correction_noME1_1

float CSCXonStrip_MatchGatti::x_correction_noME1_1[n_SW_noME1_1][n_val]

private

Definition at line 110 of file CSCXonStrip_MatchGatti.h.

Referenced by estimated2GattiCorrection(), and hardcodedCorrectionInitialization().

◆ xt_asymmetry

float CSCXonStrip_MatchGatti::xt_asymmetry

private

Definition at line 115 of file CSCXonStrip_MatchGatti.h.

Referenced by calculateXonStripError(), and findXOnStrip().

◆ xt_asymmetry_ME12

float CSCXonStrip_MatchGatti::xt_asymmetry_ME12

private

Definition at line 125 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_asymmetry_ME13

float CSCXonStrip_MatchGatti::xt_asymmetry_ME13

private

Definition at line 128 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_asymmetry_ME1a

float CSCXonStrip_MatchGatti::xt_asymmetry_ME1a

private

Definition at line 119 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_asymmetry_ME1b

float CSCXonStrip_MatchGatti::xt_asymmetry_ME1b

private

Definition at line 122 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_asymmetry_ME21

float CSCXonStrip_MatchGatti::xt_asymmetry_ME21

private

Definition at line 131 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_asymmetry_ME22

float CSCXonStrip_MatchGatti::xt_asymmetry_ME22

private

Definition at line 134 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_asymmetry_ME31

float CSCXonStrip_MatchGatti::xt_asymmetry_ME31

private

Definition at line 137 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_asymmetry_ME32

float CSCXonStrip_MatchGatti::xt_asymmetry_ME32

private

Definition at line 140 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_asymmetry_ME41

float CSCXonStrip_MatchGatti::xt_asymmetry_ME41

private

Definition at line 143 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

◆ xt_l

float CSCXonStrip_MatchGatti::xt_l[3][3]

private

x-talks 0 = left, 1 = middle, 2 = right ; and then second [] is for time bin tmax-1, tmax, tmax+1

Definition at line 68 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip(), and setupMatrix().

◆ xt_lr0

float CSCXonStrip_MatchGatti::xt_lr0[3]

private

Definition at line 69 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip(), and setupMatrix().

◆ xt_lr1

float CSCXonStrip_MatchGatti::xt_lr1[3]

private

Definition at line 69 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip(), and setupMatrix().

◆ xt_lr2

float CSCXonStrip_MatchGatti::xt_lr2[3]

private

Definition at line 69 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip(), and setupMatrix().

◆ xt_r

float CSCXonStrip_MatchGatti::xt_r[3][3]

private

Definition at line 68 of file CSCXonStrip_MatchGatti.h.

Referenced by findXOnStrip(), and setupMatrix().

◆ xtalksOffset

float CSCXonStrip_MatchGatti::xtalksOffset

private

Definition at line 87 of file CSCXonStrip_MatchGatti.h.

Referenced by CSCXonStrip_MatchGatti(), and findXOnStrip().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CSCXonStrip_MatchGatti() [1/2]

◆ ~CSCXonStrip_MatchGatti()

◆ CSCXonStrip_MatchGatti() [2/2]

Member Function Documentation

◆ calculateXonStripError()

◆ calculateXonStripPosition()

◆ estimated2Gatti()

◆ estimated2GattiCorrection()

◆ findXOnStrip()

◆ getCorrectionValues()

◆ hardcodedCorrectionInitialization()

◆ initChamberSpecs()

◆ operator=()

◆ setConditions()

◆ setupMatrix()

◆ xfError_Noise()

◆ xfError_XTasym()

Member Data Documentation

◆ a11

◆ a12

◆ a13

◆ a22

◆ a23

◆ a33

◆ chargeSignal

◆ const_syst

◆ const_syst_ME12

◆ const_syst_ME13

◆ const_syst_ME1a

◆ const_syst_ME1b

◆ const_syst_ME21

◆ const_syst_ME22

◆ const_syst_ME31

◆ const_syst_ME32

◆ const_syst_ME41

◆ h

◆ k_1

◆ k_2

◆ k_3

◆ n_SW_ME1_1

◆ n_SW_noME1_1

◆ n_val

◆ noise_level

◆ noise_level_ME12

◆ noise_level_ME13

◆ noise_level_ME1a

◆ noise_level_ME1b

◆ noise_level_ME21

◆ noise_level_ME22

◆ noise_level_ME31

◆ noise_level_ME32

◆ noise_level_ME41

◆ norm

◆ q_sum

◆ q_sumC

◆ q_sumL

◆ q_sumR

◆ r

◆ recoConditions_

◆ specs_

◆ sqrt_k_3

◆ stripWidth

◆ use3TimeBins

◆ useCalib

◆ v11

◆ v12

◆ v13

◆ v22

◆ v23

◆ v33

◆ x_centralVal

◆ x_correction_ME1_1

◆ x_correction_noME1_1

◆ xt_asymmetry