PhotonFix Class Reference

#include <PhotonFix.h>

List of all members.

Public Member Functions
double	eta () const
double	eta () const
double	etaC () const
double	etaC () const
double	etaM () const
double	etaM () const
double	etaS () const
double	etaS () const
double	fixedEnergy () const
double	fixedEnergy () const
double	phi () const
double	phi () const
double	phiC () const
double	phiC () const
double	phiM () const
double	phiM () const
double	phiS () const
double	phiS () const
	PhotonFix (double e, double eta, double phi, double r9)
	PhotonFix (double e, double eta, double phi, double r9)
void	print () const
void	print () const
double	r9 () const
double	r9 () const
double	rawEnergy () const
double	rawEnergy () const
void	setup ()
void	setup ()
double	sigmaEnergy () const
double	sigmaEnergy () const
double	xC () const
double	xC () const
double	xM () const
double	xM () const
double	xS () const
double	xS () const
double	xZ () const
double	xZ () const
double	yC () const
double	yC () const
double	yM () const
double	yM () const
double	yS () const
double	yS () const
double	yZ () const
double	yZ () const
Static Public Member Functions
static double	asinh (double s)
static double	asinh (double s)
static void	barrelCGap (unsigned i, unsigned j, unsigned k, double c)
static void	barrelCGap (unsigned i, unsigned j, unsigned k, double c)
static void	barrelMGap (unsigned i, unsigned j, unsigned k, double c)
static void	barrelMGap (unsigned i, unsigned j, unsigned k, double c)
static void	barrelSGap (unsigned i, unsigned j, unsigned k, double c)
static void	barrelSGap (unsigned i, unsigned j, unsigned k, double c)
static void	dumpGaps (std::ostream &o)
static void	dumpGaps (std::ostream &o)
static void	dumpParameters (std::ostream &o)
static void	dumpParameters (std::ostream &o)
static void	endcapCGap (unsigned i, unsigned j, unsigned k, double c)
static void	endcapCGap (unsigned i, unsigned j, unsigned k, double c)
static void	endcapCrystal (unsigned i, unsigned j, bool c)
static void	endcapCrystal (unsigned i, unsigned j, bool c)
static void	endcapMGap (unsigned i, unsigned j, unsigned k, double c)
static void	endcapMGap (unsigned i, unsigned j, unsigned k, double c)
static void	endcapSGap (unsigned i, unsigned j, unsigned k, double c)
static void	endcapSGap (unsigned i, unsigned j, unsigned k, double c)
static double	GetaPhi (double f0, double f1)
static double	GetaPhi (double f0, double f1)
static void	getParameters (unsigned be, unsigned hl, double *p)
static void	getParameters (unsigned be, unsigned hl, double *p)
static bool	initialise (const std::string &s="Nominal")
static bool	initialise (const std::string &s="Nominal")
static bool	initialised ()
static bool	initialised ()
static bool	initialiseGeometry (const std::string &s)
static bool	initialiseGeometry (const std::string &s)
static bool	initialiseParameters (const std::string &s)
static bool	initialiseParameters (const std::string &s)
static void	printParameters (std::ostream &o)
static void	printParameters (std::ostream &o)
static void	setParameters (unsigned be, unsigned hl, const double *p)
static void	setParameters (unsigned be, unsigned hl, const double *p)
Static Private Member Functions
static double	aPhi (double f0, double f1)
static double	aPhi (double f0, double f1)
static double	dPhi (double f0, double f1)
static double	dPhi (double f0, double f1)
static double	expCorrection (double a, const double *p)
static double	expCorrection (double a, const double *p)
static double	gausCorrection (double a, const double *p)
static double	gausCorrection (double a, const double *p)
Private Attributes
double	_aC
double	_aM
double	_aS
double	_bC
unsigned	_be
double	_bM
double	_bS
double	_e
double	_eta
unsigned	_hl
double	_phi
double	_r9
Static Private Attributes
static double	_barrelCGap [169][360][2]
static double	_barrelMGap [7][18][2]
static double	_barrelSGap [33][180][2]
static double	_endcapCGap [2][7080][2]
static bool	_endcapCrystal [100][100]
static double	_endcapMGap [2][1][2]
static double	_endcapSGap [2][264][2]
static bool	_initialised = false
static double	_meanAC [2][2][4]
static double	_meanAM [2][2][4]
static double	_meanAS [2][2][4]
static double	_meanAT [2][2][4]
static double	_meanBC [2][2][4]
static double	_meanBM [2][2][4]
static double	_meanBS [2][2][4]
static double	_meanBT [2][2][4]
static double	_meanR9 [2][2][4]
static double	_meanScale [2][2][4]
static const double	_onePi
static double	_sigmaAC [2][2][4]
static double	_sigmaAM [2][2][4]
static double	_sigmaAS [2][2][4]
static double	_sigmaAT [2][2][4]
static double	_sigmaBC [2][2][4]
static double	_sigmaBM [2][2][4]
static double	_sigmaBS [2][2][4]
static double	_sigmaBT [2][2][4]
static double	_sigmaR9 [2][2][4]
static double	_sigmaScale [2][2][4]
static const double	_twoPi

---

Detailed Description

Definition at line 49 of file PhotonFix.h.

---

Constructor & Destructor Documentation

PhotonFix::PhotonFix	(	double	e,
		double	eta,
		double	phi,
		double	r9
	)

Definition at line 9 of file PhotonFix.cc.

References setup().

                                                                :
  _e(e), _eta(eta), _phi(phi), _r9(r9) {

  setup();
}

PhotonFix::PhotonFix	(	double	e,
		double	eta,
		double	phi,
		double	r9
	)

---

Member Function Documentation

double PhotonFix::aPhi	(	double	f0,
		double	f1
	)		[static, private]

Definition at line 623 of file PhotonFix.cc.

References _onePi, and dPhi().

Referenced by GetaPhi().

                                           {
  double af(0.5*(f0+f1));
  if(fabs(dPhi(af,f0))>0.5*_onePi) {
    if(af>=0.0) af-=_onePi;
    else        af+=_onePi;
  }
  
  assert(fabs(dPhi(af,f0))<0.5*_onePi);
  assert(fabs(dPhi(af,f1))<0.5*_onePi);
  
  return af;
}

static double PhotonFix::aPhi	(	double	f0,
		double	f1
	)		[static, private]

double PhotonFix::asinh

(

double

s

)

[static]

Definition at line 552 of file PhotonFix.cc.

References funct::log(), and mathSSE::sqrt().

Referenced by xZ(), and yZ().

                                {
  if(s>=0.0) return  log(sqrt(s*s+1.0)+s);
  else       return -log(sqrt(s*s+1.0)-s);
}

static double PhotonFix::asinh

(

double

s

)

[static]

static void PhotonFix::barrelCGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

void PhotonFix::barrelCGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

Definition at line 328 of file PhotonFix.cc.

References _barrelCGap, trackerHits::c, i, j, and gen::k.

Referenced by PhotonFixCMS::initialise().

                                                                      {
  _barrelCGap[i][j][k] = c;
}

static void PhotonFix::barrelMGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

void PhotonFix::barrelMGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

Definition at line 334 of file PhotonFix.cc.

References _barrelMGap, trackerHits::c, i, j, and gen::k.

Referenced by PhotonFixCMS::initialise().

                                                                      {
  _barrelMGap[i][j][k] = c;
}

static void PhotonFix::barrelSGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

void PhotonFix::barrelSGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

Definition at line 331 of file PhotonFix.cc.

References _barrelSGap, trackerHits::c, i, j, and gen::k.

Referenced by PhotonFixCMS::initialise().

                                                                      {
  _barrelSGap[i][j][k] = c;
}

double PhotonFix::dPhi	(	double	f0,
		double	f1
	)		[static, private]

Definition at line 616 of file PhotonFix.cc.

References _onePi, and _twoPi.

Referenced by aPhi(), and setup().

                                           {
  double df(f0-f1);
  if(df> _onePi) df-=_twoPi;
  if(df<-_onePi) df+=_twoPi;
  return df;
}

static double PhotonFix::dPhi	(	double	f0,
		double	f1
	)		[static, private]

static void PhotonFix::dumpGaps

(

std::ostream &

o

)

[static]

void PhotonFix::dumpGaps

(

std::ostream &

o

)

[static]

Definition at line 557 of file PhotonFix.cc.

References _barrelCGap, _barrelMGap, _barrelSGap, _endcapCGap, _endcapCrystal, _endcapMGap, _endcapSGap, i, j, and gen::k.

                                      {
  o << std::setprecision(15);
  
  for(unsigned i(0);i<169;i++) {
    for(unsigned j(0);j<360;j++) {
      for(unsigned k(0);k<2;k++) {
        o << _barrelCGap[i][j][k] << std::endl;
      }
    }
  }
  
  for(unsigned i(0);i<33;i++) {
    for(unsigned j(0);j<180;j++) {
      for(unsigned k(0);k<2;k++) {
        o << _barrelSGap[i][j][k] << std::endl;
      }
    }
  }
  
  for(unsigned i(0);i<7;i++) {
    for(unsigned j(0);j<18;j++) {
      for(unsigned k(0);k<2;k++) {
        o << _barrelMGap[i][j][k] << std::endl;
      }
    }
  }
  
  for(unsigned i(0);i<100;i++) {
    for(unsigned j(0);j<100;j++) {
      if(_endcapCrystal[i][j]) o << 0 << std::endl;
      else                     o << 1 << std::endl;
    }
  }
  
  for(unsigned i(0);i<2;i++) {
    for(unsigned j(0);j<7080;j++) {
      for(unsigned k(0);k<2;k++) {
        o << _endcapCGap[i][j][k] << std::endl;
      }
    }
  }
  
  for(unsigned i(0);i<2;i++) {
    for(unsigned j(0);j<264;j++) {
      for(unsigned k(0);k<2;k++) {
        o << _endcapSGap[i][j][k] << std::endl;
      }
    }
  }
  
  for(unsigned i(0);i<2;i++) {
    for(unsigned j(0);j<1;j++) {
      for(unsigned k(0);k<2;k++) {
        o << _endcapMGap[i][j][k] << std::endl;
      }
    }
  }
}

void PhotonFix::dumpParameters

(

std::ostream &

o

)

[static]

Definition at line 409 of file PhotonFix.cc.

References _meanAC, _meanAM, _meanAS, _meanAT, _meanBC, _meanBM, _meanBS, _meanBT, _meanR9, _meanScale, _sigmaAC, _sigmaAM, _sigmaAS, _sigmaAT, _sigmaBC, _sigmaBM, _sigmaBS, _sigmaBT, _sigmaR9, _sigmaScale, and i.

                                            {
  for(unsigned be(0);be<2;be++) {
    for(unsigned hl(0);hl<2;hl++) {
      for(unsigned i(0);i<4;i++) {
        o << " _meanScale[" << be << "][" << hl << "][" << i << "]=" << _meanScale[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanAT[" << be << "][" << hl << "][" << i << "]=" << _meanAT[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanAC[" << be << "][" << hl << "][" << i << "]=" << _meanAC[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanAS[" << be << "][" << hl << "][" << i << "]=" << _meanAS[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanAM[" << be << "][" << hl << "][" << i << "]=" << _meanAM[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanBT[" << be << "][" << hl << "][" << i << "]=" << _meanBT[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanBC[" << be << "][" << hl << "][" << i << "]=" << _meanBC[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanBS[" << be << "][" << hl << "][" << i << "]=" << _meanBS[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanBM[" << be << "][" << hl << "][" << i << "]=" << _meanBM[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _meanR9[" << be << "][" << hl << "][" << i << "]=" << _meanR9[be][hl][i] << ";" << std::endl;
      }
      o << std::endl;
      
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaScale[" << be << "][" << hl << "][" << i << "]=" << _sigmaScale[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaAT[" << be << "][" << hl << "][" << i << "]=" << _sigmaAT[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaAC[" << be << "][" << hl << "][" << i << "]=" << _sigmaAC[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaAS[" << be << "][" << hl << "][" << i << "]=" << _sigmaAS[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaAM[" << be << "][" << hl << "][" << i << "]=" << _sigmaAM[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaBT[" << be << "][" << hl << "][" << i << "]=" << _sigmaBT[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaBC[" << be << "][" << hl << "][" << i << "]=" << _sigmaBC[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaBS[" << be << "][" << hl << "][" << i << "]=" << _sigmaBS[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaBM[" << be << "][" << hl << "][" << i << "]=" << _sigmaBM[be][hl][i] << ";" << std::endl;
      }
      for(unsigned i(0);i<4;i++) {
        o << " _sigmaR9[" << be << "][" << hl << "][" << i << "]=" << _sigmaR9[be][hl][i] << ";" << std::endl;
      }
      o << std::endl;
    }
  }
}

static void PhotonFix::dumpParameters

(

std::ostream &

o

)

[static]

static void PhotonFix::endcapCGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

void PhotonFix::endcapCGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

Definition at line 340 of file PhotonFix.cc.

References _endcapCGap, trackerHits::c, i, j, and gen::k.

Referenced by PhotonFixCMS::initialise().

                                                                      {
  _endcapCGap[i][j][k] = c;
}

static void PhotonFix::endcapCrystal	(	unsigned	i,
		unsigned	j,
		bool	c
	)		[static]

void PhotonFix::endcapCrystal	(	unsigned	i,
		unsigned	j,
		bool	c
	)		[static]

Definition at line 337 of file PhotonFix.cc.

References _endcapCrystal, trackerHits::c, i, and j.

Referenced by PhotonFixCMS::initialise().

                                                           {
  _endcapCrystal[i][j] = c;
}

void PhotonFix::endcapMGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

Definition at line 346 of file PhotonFix.cc.

References _endcapMGap, trackerHits::c, i, j, and gen::k.

Referenced by PhotonFixCMS::initialise().

                                                                      {
  _endcapMGap[i][j][k] = c;
}

static void PhotonFix::endcapMGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

static void PhotonFix::endcapSGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

void PhotonFix::endcapSGap	(	unsigned	i,
		unsigned	j,
		unsigned	k,
		double	c
	)		[static]

Definition at line 343 of file PhotonFix.cc.

References _endcapSGap, trackerHits::c, i, j, and gen::k.

Referenced by PhotonFixCMS::initialise().

                                                                      {
  _endcapSGap[i][j][k] = c;
}

double PhotonFix::eta

(

void

)

const

Definition at line 242 of file PhotonFix.cc.

References _eta.

                            {
  return _eta;
}

double PhotonFix::eta

(

)

const

double PhotonFix::etaC

(

)

const

Definition at line 254 of file PhotonFix.cc.

References _aC, and _be.

                             {
  assert(_be==0);
  return _aC;
}

double PhotonFix::etaC

(

)

const

double PhotonFix::etaM

(

)

const

Definition at line 264 of file PhotonFix.cc.

References _aM, and _be.

                             {
  assert(_be==0);
  return _aM;
}

double PhotonFix::etaM

(

)

const

double PhotonFix::etaS

(

)

const

Definition at line 259 of file PhotonFix.cc.

References _aS, and _be.

                             {
  assert(_be==0);
  return _aS;
}

double PhotonFix::etaS

(

)

const

static double PhotonFix::expCorrection	(	double	a,
		const double *	p
	)		[static, private]

double PhotonFix::expCorrection	(	double	a,
		const double *	p
	)		[static, private]

Definition at line 636 of file PhotonFix.cc.

References b, and funct::exp().

Referenced by fixedEnergy(), and sigmaEnergy().

                                                         {
  if(p[1]==0.0 || p[2]==0.0 || p[3]==0.0) return 0.0;
  
  double b(a-p[0]);
  if(b>=0.0) return p[1]*exp(-fabs(p[2])*b);
  else       return p[1]*exp( fabs(p[3])*b);
}

double PhotonFix::fixedEnergy

(

)

const

Definition at line 139 of file PhotonFix.cc.

References _aC, _aM, _aS, _bC, _be, _bM, _bS, _e, _eta, _hl, _meanAC, _meanAM, _meanAS, _meanAT, _meanBC, _meanBM, _meanBS, _meanBT, _meanR9, _meanScale, _r9, expCorrection(), f, xZ(), and yZ().

Referenced by PhotonFixCMS::fixedEnergy().

                                    {
  double f(0.0);
  
  // Overall scale and energy(T) dependence
  f =_meanScale[_be][_hl][0];
  f+=_meanScale[_be][_hl][1]*_e;
  f+=_meanScale[_be][_hl][2]*_e/cosh(_eta);
  f+=_meanScale[_be][_hl][3]*cosh(_eta)/_e;
  
  // General eta or zeta dependence
  if(_be==0) {
    f+=_meanAT[_be][_hl][0]*_eta*_eta;
    f+=expCorrection(_eta,_meanBT[_be][_hl]);
  } else {
    f+=_meanAT[_be][_hl][0]*xZ()*xZ();
    f+=_meanBT[_be][_hl][0]*yZ()*yZ();
  }
  
  // Eta or x crystal, submodule and module dependence
  f+=expCorrection(_aC,_meanAC[_be][_hl]);
  f+=expCorrection(_aS,_meanAS[_be][_hl]);
  f+=expCorrection(_aM,_meanAM[_be][_hl]);
  
  // Phi or y crystal, submodule and module dependence
  f+=expCorrection(_bC,_meanBC[_be][_hl]);
  f+=expCorrection(_bS,_meanBS[_be][_hl]);
  f+=expCorrection(_bM,_meanBM[_be][_hl]);
  
  // R9 dependence
  if(_hl==0) {
    f+=_meanR9[_be][_hl][1]*(_r9-_meanR9[_be][_hl][0])*(_r9-_meanR9[_be][_hl][0])
      +_meanR9[_be][_hl][2]*(_r9-_meanR9[_be][_hl][0])*(_r9-_meanR9[_be][_hl][0])*(_r9-_meanR9[_be][_hl][0]);
  } else {
    f+=_meanR9[_be][_hl][0]*_r9+_meanR9[_be][_hl][1]*_r9*_r9+_meanR9[_be][_hl][2]*_r9*_r9*_r9;
  }
  
  return _e*f;
}

double PhotonFix::fixedEnergy

(

)

const

double PhotonFix::gausCorrection	(	double	a,
		const double *	p
	)		[static, private]

Definition at line 644 of file PhotonFix.cc.

References b, and funct::exp().

                                                          {
  if(p[1]==0.0 || p[2]==0.0 || p[3]==0.0) return 0.0;
  
  double b(a-p[0]);
  if(b>=0.0) return p[1]*exp(-0.5*p[2]*p[2]*b*b);
  else       return p[1]*exp(-0.5*p[3]*p[3]*b*b);
}

static double PhotonFix::gausCorrection	(	double	a,
		const double *	p
	)		[static, private]

static double PhotonFix::GetaPhi	(	double	f0,
		double	f1
	)		[static]

double PhotonFix::GetaPhi	(	double	f0,
		double	f1
	)		[static]

Definition at line 324 of file PhotonFix.cc.

References aPhi().

Referenced by PhotonFixCMS::initialise().

                                             {
  return aPhi(f0,f1);
}

void PhotonFix::getParameters	(	unsigned	be,
		unsigned	hl,
		double *	p
	)		[static]

Definition at line 383 of file PhotonFix.cc.

References _meanAC, _meanAM, _meanAS, _meanAT, _meanBC, _meanBM, _meanBS, _meanBT, _meanR9, _meanScale, _sigmaAC, _sigmaAM, _sigmaAS, _sigmaAT, _sigmaBC, _sigmaBM, _sigmaBS, _sigmaBT, _sigmaR9, _sigmaScale, and i.

                                                                 {
  for(unsigned i(0);i<4;i++) {
    p[i+ 0*4]=_meanScale[be][hl][i];
    p[i+ 1*4]=_meanAT[be][hl][i];
    p[i+ 2*4]=_meanAC[be][hl][i];
    p[i+ 3*4]=_meanAS[be][hl][i];
    p[i+ 4*4]=_meanAM[be][hl][i];
    p[i+ 5*4]=_meanBT[be][hl][i];
    p[i+ 6*4]=_meanBC[be][hl][i];
    p[i+ 7*4]=_meanBS[be][hl][i];
    p[i+ 8*4]=_meanBM[be][hl][i];
    p[i+ 9*4]=_meanR9[be][hl][i];
    
    p[i+10*4]=_sigmaScale[be][hl][i];
    p[i+11*4]=_sigmaAT[be][hl][i];
    p[i+12*4]=_sigmaAC[be][hl][i];
    p[i+13*4]=_sigmaAS[be][hl][i];
    p[i+14*4]=_sigmaAM[be][hl][i];
    p[i+15*4]=_sigmaBT[be][hl][i];
    p[i+16*4]=_sigmaBC[be][hl][i];
    p[i+17*4]=_sigmaBS[be][hl][i];
    p[i+18*4]=_sigmaBM[be][hl][i];
    p[i+19*4]=_sigmaR9[be][hl][i];
  }
}

static void PhotonFix::getParameters	(	unsigned	be,
		unsigned	hl,
		double *	p
	)		[static]

static bool PhotonFix::initialise

(

const std::string &

s = "Nominal"

)

[static]

bool PhotonFix::initialise

(

const std::string &

s = "Nominal"

)

[static]

Definition at line 654 of file PhotonFix.cc.

References _initialised, initialiseGeometry(), and initialiseParameters().

                                             {
  if(_initialised) return false;

  
  initialiseParameters(s);
  initialiseGeometry(s);
  return true;
}

bool PhotonFix::initialised

(

)

[static]

Definition at line 651 of file PhotonFix.cc.

References _initialised.

Referenced by PhotonFixCMS::initialise().

                            {
  return _initialised;
}

static bool PhotonFix::initialised

(

)

[static]

bool PhotonFix::initialiseGeometry

(

const std::string &

s

)

[static]

Definition at line 2042 of file PhotonFix.cc.

References _barrelCGap, _barrelMGap, _barrelSGap, _endcapCGap, _endcapCrystal, _endcapMGap, _endcapSGap, gather_cfg::cout, groupFilesInBlocks::fin, i, j, and gen::k.

Referenced by initialise().

                                                     {

  std::ifstream fin("../test/PhotonFix.dat");
  assert(fin);

  std::cout << "Reading in here" << std::endl;
  for(unsigned i(0);i<169;i++) {
    for(unsigned j(0);j<360;j++) {
      for(unsigned k(0);k<2;k++) {
        fin >> _barrelCGap[i][j][k];
      }
    }
  }
 
  for(unsigned i(0);i<33;i++) {
    for(unsigned j(0);j<180;j++) {
      for(unsigned k(0);k<2;k++) {
        fin >> _barrelSGap[i][j][k];
      }
    }
  }
 
  for(unsigned i(0);i<7;i++) {
    for(unsigned j(0);j<18;j++) {
      for(unsigned k(0);k<2;k++) {
        fin >> _barrelMGap[i][j][k];
      }
    }
  }
  for(unsigned i(0);i<100;i++) {
    for(unsigned j(0);j<100;j++) {
      unsigned k;
      fin >> k;
      _endcapCrystal[i][j]=(k==0);
    }
  }
 
  for(unsigned i(0);i<2;i++) {
    for(unsigned j(0);j<7080;j++) {
      for(unsigned k(0);k<2;k++) {
        fin >> _endcapCGap[i][j][k];
      }
    }
  }
 
  for(unsigned i(0);i<2;i++) {
    for(unsigned j(0);j<264;j++) {
      for(unsigned k(0);k<2;k++) {
        fin >> _endcapSGap[i][j][k];
      }
    }
  }
 
  for(unsigned i(0);i<2;i++) {
    for(unsigned j(0);j<1;j++) {
      for(unsigned k(0);k<2;k++) {
        fin >> _endcapMGap[i][j][k];
      }
    }
  }
 
  assert(fin);
 
  return true;
}

static bool PhotonFix::initialiseGeometry

(

const std::string &

s

)

[static]

static bool PhotonFix::initialiseParameters

(

const std::string &

s

)

[static]

bool PhotonFix::initialiseParameters

(

const std::string &

s

)

[static]

Definition at line 663 of file PhotonFix.cc.

References _initialised, _meanAC, _meanAM, _meanAS, _meanAT, _meanBC, _meanBM, _meanBS, _meanBT, _meanR9, _meanScale, _sigmaAC, _sigmaAM, _sigmaAS, _sigmaAT, _sigmaBC, _sigmaBM, _sigmaBS, _sigmaBT, _sigmaR9, _sigmaScale, alignCSCRings::e, and i.

Referenced by PhotonFixCMS::initialise(), and initialise().

                                                       {
  _initialised=false;
  
  if(s=="Nominal") {
    for(unsigned be(0);be<2;be++) {
      for(unsigned hl(0);hl<2;hl++) {
        for(unsigned i(0);i<4;i++) {
          _meanScale[be][hl][i]=0;
          _meanAT[be][hl][i]=0;
          _meanAC[be][hl][i]=0;
          _meanAS[be][hl][i]=0;
          _meanAM[be][hl][i]=0;
          _meanBT[be][hl][i]=0;
          _meanBC[be][hl][i]=0;
          _meanBS[be][hl][i]=0;
          _meanBM[be][hl][i]=0;
          _meanR9[be][hl][i]=0;
          
          _sigmaScale[be][hl][i]=0;
          _sigmaAT[be][hl][i]=0;
          _sigmaAC[be][hl][i]=0;
          _sigmaAS[be][hl][i]=0;
          _sigmaAM[be][hl][i]=0;
          _sigmaBT[be][hl][i]=0;
          _sigmaBC[be][hl][i]=0;
          _sigmaBS[be][hl][i]=0;
          _sigmaBM[be][hl][i]=0;
          _sigmaR9[be][hl][i]=0;
        }
        
        _meanScale[be][hl][0]=1.0;
        if(be==0) {
          _sigmaScale[be][hl][0]=0.2;
          _sigmaScale[be][hl][1]=0.03;
          _sigmaScale[be][hl][2]=0.006;
        } else {
          _sigmaScale[be][hl][0]=0.25;
          _sigmaScale[be][hl][1]=0.05;
          _sigmaScale[be][hl][2]=0.010;
        }
      }
    }
    
    _initialised=true;
  }
  
  if(s=="3_8") {
    _meanScale[0][0][0]=0.994724;
    _meanScale[0][0][1]=1.98102e-06;
    _meanScale[0][0][2]=1.43015e-05;
    _meanScale[0][0][3]=-0.0908525;
    _meanAT[0][0][0]=0.0;
    _meanAT[0][0][1]=0.0;
    _meanAT[0][0][2]=0.0;
    _meanAT[0][0][3]=0.0;
    _meanAC[0][0][0]=-0.00352041;
    _meanAC[0][0][1]=0.00982015;
    _meanAC[0][0][2]=434.32;
    _meanAC[0][0][3]=529.508;
    _meanAS[0][0][0]=-1.1;
    _meanAS[0][0][1]=0.00135995;
    _meanAS[0][0][2]=295.712;
    _meanAS[0][0][3]=5.13202e+07;
    _meanAM[0][0][0]=-0.00140562;
    _meanAM[0][0][1]=0.156322;
    _meanAM[0][0][2]=263.097;
    _meanAM[0][0][3]=222.294;
    _meanBT[0][0][0]=0.0;
    _meanBT[0][0][1]=0.0;
    _meanBT[0][0][2]=0.0;
    _meanBT[0][0][3]=0.0;
    _meanBC[0][0][0]=-0.00294295;
    _meanBC[0][0][1]=0.011533;
    _meanBC[0][0][2]=562.905;
    _meanBC[0][0][3]=421.097;
    _meanBS[0][0][0]=-0.00204373;
    _meanBS[0][0][1]=0.00347592;
    _meanBS[0][0][2]=36.5614;
    _meanBS[0][0][3]=1265.25;
    _meanBM[0][0][0]=-0.00275381;
    _meanBM[0][0][1]=0.0812447;
    _meanBM[0][0][2]=216.885;
    _meanBM[0][0][3]=264.754;
    _meanR9[0][0][0]=0.952584;
    _meanR9[0][0][1]=22.7119;
    _meanR9[0][0][2]=402.816;
    _meanR9[0][0][3]=0;
    
    _sigmaScale[0][0][0]=0.167184;
    _sigmaScale[0][0][1]=6.14323e-11;
    _sigmaScale[0][0][2]=0.00769693;
    _sigmaScale[0][0][3]=0;
    _sigmaAT[0][0][0]=0.228255;
    _sigmaAT[0][0][1]=0;
    _sigmaAT[0][0][2]=0;
    _sigmaAT[0][0][3]=0;
    _sigmaAC[0][0][0]=-0.00411906;
    _sigmaAC[0][0][1]=0.077799;
    _sigmaAC[0][0][2]=23.1033;
    _sigmaAC[0][0][3]=-3e+17;
    _sigmaAS[0][0][0]=0;
    _sigmaAS[0][0][1]=0;
    _sigmaAS[0][0][2]=0;
    _sigmaAS[0][0][3]=0;
    _sigmaAM[0][0][0]=-0.000130695;
    _sigmaAM[0][0][1]=11.2121;
    _sigmaAM[0][0][2]=468.535;
    _sigmaAM[0][0][3]=407.652;
    _sigmaBT[0][0][0]=1.33384e-05;
    _sigmaBT[0][0][1]=8.77098;
    _sigmaBT[0][0][2]=324.048;
    _sigmaBT[0][0][3]=239.868;
    _sigmaBC[0][0][0]=-0.00281964;
    _sigmaBC[0][0][1]=0.125811;
    _sigmaBC[0][0][2]=538.949;
    _sigmaBC[0][0][3]=1358.76;
    _sigmaBS[0][0][0]=0;
    _sigmaBS[0][0][1]=0;
    _sigmaBS[0][0][2]=0;
    _sigmaBS[0][0][3]=0;
    _sigmaBM[0][0][0]=-0.00293676;
    _sigmaBM[0][0][1]=8.88276;
    _sigmaBM[0][0][2]=350.032;
    _sigmaBM[0][0][3]=580.354;
    _sigmaR9[0][0][0]=0.955876;
    _sigmaR9[0][0][1]=2254.5;
    _sigmaR9[0][0][2]=14627;
    _sigmaR9[0][0][3]=0;
    
    _meanScale[0][1][0]=0.888348;
    _meanScale[0][1][1]=1.20452e-05;
    _meanScale[0][1][2]=-1.04458e-05;
    _meanScale[0][1][3]=-0.542383;
    _meanAT[0][1][0]=0.0;
    _meanAT[0][1][1]=0.0;
    _meanAT[0][1][2]=0.0;
    _meanAT[0][1][3]=0.0;
    _meanAC[0][1][0]=-0.00320856;
    _meanAC[0][1][1]=0.0240109;
    _meanAC[0][1][2]=115.145;
    _meanAC[0][1][3]=205.859;
    _meanAS[0][1][0]=0.0349736;
    _meanAS[0][1][1]=-0.00232864;
    _meanAS[0][1][2]=318.584;
    _meanAS[0][1][3]=1.4e+09;
    _meanAM[0][1][0]=-0.00104798;
    _meanAM[0][1][1]=0.208249;
    _meanAM[0][1][2]=297.049;
    _meanAM[0][1][3]=220.609;
    _meanBT[0][1][0]=0.0;
    _meanBT[0][1][1]=0.0;
    _meanBT[0][1][2]=0.0;
    _meanBT[0][1][3]=0.0;
    _meanBC[0][1][0]=-0.00420429;
    _meanBC[0][1][1]=0.00203991;
    _meanBC[0][1][2]=172.278;
    _meanBC[0][1][3]=410.677;
    _meanBS[0][1][0]=-0.0430854;
    _meanBS[0][1][1]=0.0961883;
    _meanBS[0][1][2]=0.196958;
    _meanBS[0][1][3]=11442.2;
    _meanBM[0][1][0]=-0.00389457;
    _meanBM[0][1][1]=0.0449086;
    _meanBM[0][1][2]=78.9252;
    _meanBM[0][1][3]=103.237;
    _meanR9[0][1][0]=0.0182102;
    _meanR9[0][1][1]=-0.03752;
    _meanR9[0][1][2]=0.0198881;
    _meanR9[0][1][3]=0;
    
    _sigmaScale[0][1][0]=0.386681;
    _sigmaScale[0][1][1]=0.0913412;
    _sigmaScale[0][1][2]=0.00119232;
    _sigmaScale[0][1][3]=0;
    _sigmaAT[0][1][0]=1.36562;
    _sigmaAT[0][1][1]=0;
    _sigmaAT[0][1][2]=0;
    _sigmaAT[0][1][3]=0;
    _sigmaAC[0][1][0]=-0.00504613;
    _sigmaAC[0][1][1]=-1.09115;
    _sigmaAC[0][1][2]=8.57406;
    _sigmaAC[0][1][3]=57.1351;
    _sigmaAS[0][1][0]=0;
    _sigmaAS[0][1][1]=0;
    _sigmaAS[0][1][2]=0;
    _sigmaAS[0][1][3]=0;
    _sigmaAM[0][1][0]=-0.00014319;
    _sigmaAM[0][1][1]=5.39527;
    _sigmaAM[0][1][2]=432.566;
    _sigmaAM[0][1][3]=265.165;
    _sigmaBT[0][1][0]=-0.040161;
    _sigmaBT[0][1][1]=2.65711;
    _sigmaBT[0][1][2]=-0.398357;
    _sigmaBT[0][1][3]=-0.440649;
    _sigmaBC[0][1][0]=0.00580015;
    _sigmaBC[0][1][1]=-0.631833;
    _sigmaBC[0][1][2]=18594.3;
    _sigmaBC[0][1][3]=4.00955e+08;
    _sigmaBS[0][1][0]=0;
    _sigmaBS[0][1][1]=0;
    _sigmaBS[0][1][2]=0;
    _sigmaBS[0][1][3]=0;
    _sigmaBM[0][1][0]=-0.00376665;
    _sigmaBM[0][1][1]=3.74316;
    _sigmaBM[0][1][2]=102.72;
    _sigmaBM[0][1][3]=157.396;
    _sigmaR9[0][1][0]=-3.12696;
    _sigmaR9[0][1][1]=1.75114;
    _sigmaR9[0][1][2]=0;
    _sigmaR9[0][1][3]=0;
    
    _meanScale[1][0][0]=0.999461;
    _meanScale[1][0][1]=4.37414e-06;
    _meanScale[1][0][2]=4.92078e-06;
    _meanScale[1][0][3]=-0.121609;
    _meanAT[1][0][0]=0.0;
    _meanAT[1][0][1]=0.0;
    _meanAT[1][0][2]=0.0;
    _meanAT[1][0][3]=0.0;
    _meanAC[1][0][0]=-0.000396058;
    _meanAC[1][0][1]=0.0144837;
    _meanAC[1][0][2]=1374.93;
    _meanAC[1][0][3]=945.634;
    _meanAS[1][0][0]=-0.000871036;
    _meanAS[1][0][1]=0.0442747;
    _meanAS[1][0][2]=645.709;
    _meanAS[1][0][3]=962.845;
    _meanAM[1][0][0]=0.000434298;
    _meanAM[1][0][1]=0.0658628;
    _meanAM[1][0][2]=1928.49;
    _meanAM[1][0][3]=728.522;
    _meanBT[1][0][0]=0.0;
    _meanBT[1][0][1]=0.0;
    _meanBT[1][0][2]=0.0;
    _meanBT[1][0][3]=0.0;
    _meanBC[1][0][0]=-0.000452212;
    _meanBC[1][0][1]=0.0129968;
    _meanBC[1][0][2]=1056.08;
    _meanBC[1][0][3]=759.102;
    _meanBS[1][0][0]=-0.000786157;
    _meanBS[1][0][1]=0.0346555;
    _meanBS[1][0][2]=592.239;
    _meanBS[1][0][3]=854.285;
    _meanBM[1][0][0]=-0.0665038;
    _meanBM[1][0][1]=-0.00211713;
    _meanBM[1][0][2]=4.84395;
    _meanBM[1][0][3]=11.6644;
    _meanR9[1][0][0]=0.971355;
    _meanR9[1][0][1]=47.2751;
    _meanR9[1][0][2]=536.907;
    _meanR9[1][0][3]=0;
    
    _sigmaScale[1][0][0]=0.254641;
    _sigmaScale[1][0][1]=0.00264818;
    _sigmaScale[1][0][2]=0.0114953;
    _sigmaScale[1][0][3]=0;
    _sigmaAT[1][0][0]=0.935839;
    _sigmaAT[1][0][1]=0;
    _sigmaAT[1][0][2]=0;
    _sigmaAT[1][0][3]=0;
    _sigmaAC[1][0][0]=-0.00476475;
    _sigmaAC[1][0][1]=2.14548;
    _sigmaAC[1][0][2]=29937;
    _sigmaAC[1][0][3]=2.6e+11;
    _sigmaAS[1][0][0]=-8.17285e-05;
    _sigmaAS[1][0][1]=1.5821;
    _sigmaAS[1][0][2]=1928.83;
    _sigmaAS[1][0][3]=902.519;
    _sigmaAM[1][0][0]=0.0278577;
    _sigmaAM[1][0][1]=0.58439;
    _sigmaAM[1][0][2]=43.3575;
    _sigmaAM[1][0][3]=19.7836;
    _sigmaBT[1][0][0]=-0.456051;
    _sigmaBT[1][0][1]=0;
    _sigmaBT[1][0][2]=0;
    _sigmaBT[1][0][3]=0;
    _sigmaBC[1][0][0]=-0.00264527;
    _sigmaBC[1][0][1]=0.696043;
    _sigmaBC[1][0][2]=7.49509e+12;
    _sigmaBC[1][0][3]=96843;
    _sigmaBS[1][0][0]=0.000258933;
    _sigmaBS[1][0][1]=1.28387;
    _sigmaBS[1][0][2]=1668.71;
    _sigmaBS[1][0][3]=730.716;
    _sigmaBM[1][0][0]=0.00121506;
    _sigmaBM[1][0][1]=0.938541;
    _sigmaBM[1][0][2]=9003.57;
    _sigmaBM[1][0][3]=288.897;
    _sigmaR9[1][0][0]=1.01207;
    _sigmaR9[1][0][1]=-816.244;
    _sigmaR9[1][0][2]=-16283.8;
    _sigmaR9[1][0][3]=0;
    
    _meanScale[1][1][0]=0.324634;
    _meanScale[1][1][1]=9.48206e-05;
    _meanScale[1][1][2]=1.0e-12;
    _meanScale[1][1][3]=1.0e-12;
    _meanAT[1][1][0]=0.0;
    _meanAT[1][1][1]=0.0;
    _meanAT[1][1][2]=0.0;
    _meanAT[1][1][3]=0.0;
    _meanAC[1][1][0]=-0.00158311;
    _meanAC[1][1][1]=0.0106161;
    _meanAC[1][1][2]=338.964;
    _meanAC[1][1][3]=797.172;
    _meanAS[1][1][0]=-0.00960269;
    _meanAS[1][1][1]=-0.00496491;
    _meanAS[1][1][2]=934.472;
    _meanAS[1][1][3]=8.32667e-16;
    _meanAM[1][1][0]=-0.00219814;
    _meanAM[1][1][1]=0.653906;
    _meanAM[1][1][2]=0.0949848;
    _meanAM[1][1][3]=0.0977831;
    _meanBT[1][1][0]=0.0;
    _meanBT[1][1][1]=0.0;
    _meanBT[1][1][2]=0.0;
    _meanBT[1][1][3]=0.0;
    _meanBC[1][1][0]=-0.00423472;
    _meanBC[1][1][1]=0.0279695;
    _meanBC[1][1][2]=28073.7;
    _meanBC[1][1][3]=118612;
    _meanBS[1][1][0]=-0.0012476;
    _meanBS[1][1][1]=0.02744;
    _meanBS[1][1][2]=390.697;
    _meanBS[1][1][3]=727.861;
    _meanBM[1][1][0]=-1.36573e-05;
    _meanBM[1][1][1]=0.0667504;
    _meanBM[1][1][2]=-80154.4;
    _meanBM[1][1][3]=576.637;
    _meanR9[1][1][0]=0.113317;
    _meanR9[1][1][1]=0.0142669;
    _meanR9[1][1][2]=-0.125721;
    _meanR9[1][1][3]=0;

    _sigmaScale[1][1][0]=0.471767;
    _sigmaScale[1][1][1]=0.211196;
    _sigmaScale[1][1][2]=0.0240124;
    _sigmaScale[1][1][3]=0;
    _sigmaAT[1][1][0]=0.404395;
    _sigmaAT[1][1][1]=0;
    _sigmaAT[1][1][2]=0;
    _sigmaAT[1][1][3]=0;
    _sigmaAC[1][1][0]=0.00173151;
    _sigmaAC[1][1][1]=-0.479291;
    _sigmaAC[1][1][2]=11583.5;
    _sigmaAC[1][1][3]=-7e+09;
    _sigmaAS[1][1][0]=0.000450387;
    _sigmaAS[1][1][1]=0.662978;
    _sigmaAS[1][1][2]=924.051;
    _sigmaAS[1][1][3]=448.417;
    _sigmaAM[1][1][0]=0.00335603;
    _sigmaAM[1][1][1]=0.648407;
    _sigmaAM[1][1][2]=134.672;
    _sigmaAM[1][1][3]=27.4139;
    _sigmaBT[1][1][0]=0.602402;
    _sigmaBT[1][1][1]=0;
    _sigmaBT[1][1][2]=0;
    _sigmaBT[1][1][3]=0;
    _sigmaBC[1][1][0]=-0.00256192;
    _sigmaBC[1][1][1]=2.01276;
    _sigmaBC[1][1][2]=114558;
    _sigmaBC[1][1][3]=2.15421e+06;
    _sigmaBS[1][1][0]=0.00151576;
    _sigmaBS[1][1][1]=0.359084;
    _sigmaBS[1][1][2]=329.414;
    _sigmaBS[1][1][3]=154.509;
    _sigmaBM[1][1][0]=-0.0452587;
    _sigmaBM[1][1][1]=1.26253;
    _sigmaBM[1][1][2]=1.9e+09;
    _sigmaBM[1][1][3]=1058.76;
    _sigmaR9[1][1][0]=4.59667;
    _sigmaR9[1][1][1]=-5.14404;
    _sigmaR9[1][1][2]=0;
    _sigmaR9[1][1][3]=0;

    _initialised=true;
  }

  if(s=="3_11") {
    _meanScale[0][0][0]=0.994363;
    _meanScale[0][0][1]=4.84904e-07;
    _meanScale[0][0][2]=1.54475e-05;
    _meanScale[0][0][3]=-0.103309;
    _meanAT[0][0][0]=0.0;
    _meanAT[0][0][1]=0.0;
    _meanAT[0][0][2]=0.0;
    _meanAT[0][0][3]=0.0;
    _meanAC[0][0][0]=-0.00360057;
    _meanAC[0][0][1]=0.00970858;
    _meanAC[0][0][2]=409.406;
    _meanAC[0][0][3]=527.952;
    _meanAS[0][0][0]=-1.1;
    _meanAS[0][0][1]=0.00135995;
    _meanAS[0][0][2]=295.712;
    _meanAS[0][0][3]=5.13202e+07;
    _meanAM[0][0][0]=-0.00129854;
    _meanAM[0][0][1]=0.151466;
    _meanAM[0][0][2]=261.828;
    _meanAM[0][0][3]=214.662;
    _meanBT[0][0][0]=0.0;
    _meanBT[0][0][1]=0.0;
    _meanBT[0][0][2]=0.0;
    _meanBT[0][0][3]=0.0;
    _meanBC[0][0][0]=-0.00286864;
    _meanBC[0][0][1]=0.0114118;
    _meanBC[0][0][2]=563.962;
    _meanBC[0][0][3]=412.922;
    _meanBS[0][0][0]=-0.00210996;
    _meanBS[0][0][1]=0.00327867;
    _meanBS[0][0][2]=23.617;
    _meanBS[0][0][3]=1018.45;
    _meanBM[0][0][0]=-0.002287;
    _meanBM[0][0][1]=0.0848984;
    _meanBM[0][0][2]=235.575;
    _meanBM[0][0][3]=260.773;
    _meanR9[0][0][0]=0.951724;
    _meanR9[0][0][1]=23.7181;
    _meanR9[0][0][2]=177.34;
    _meanR9[0][0][3]=0;

    _sigmaScale[0][0][0]=0.187578;
    _sigmaScale[0][0][1]=-0.000901045;
    _sigmaScale[0][0][2]=0.00673186;
    _sigmaScale[0][0][3]=0;
    _sigmaAT[0][0][0]=0.183777;
    _sigmaAT[0][0][1]=0;
    _sigmaAT[0][0][2]=0;
    _sigmaAT[0][0][3]=0;
    _sigmaAC[0][0][0]=-0.00430202;
    _sigmaAC[0][0][1]=0.122501;
    _sigmaAC[0][0][2]=51.9772;
    _sigmaAC[0][0][3]=-3e+17;
    _sigmaAS[0][0][0]=0;
    _sigmaAS[0][0][1]=0;
    _sigmaAS[0][0][2]=0;
    _sigmaAS[0][0][3]=0;
    _sigmaAM[0][0][0]=0.00101883;
    _sigmaAM[0][0][1]=11.2009;
    _sigmaAM[0][0][2]=593.111;
    _sigmaAM[0][0][3]=345.433;
    _sigmaBT[0][0][0]=-6.02356e-05;
    _sigmaBT[0][0][1]=6.99896;
    _sigmaBT[0][0][2]=235.996;
    _sigmaBT[0][0][3]=196;
    _sigmaBC[0][0][0]=-0.00282254;
    _sigmaBC[0][0][1]=0.18764;
    _sigmaBC[0][0][2]=509.825;
    _sigmaBC[0][0][3]=1400.14;
    _sigmaBS[0][0][0]=0;
    _sigmaBS[0][0][1]=0;
    _sigmaBS[0][0][2]=0;
    _sigmaBS[0][0][3]=0;
    _sigmaBM[0][0][0]=-0.00252199;
    _sigmaBM[0][0][1]=39.1544;
    _sigmaBM[0][0][2]=612.481;
    _sigmaBM[0][0][3]=905.994;
    _sigmaR9[0][0][0]=0.95608;
    _sigmaR9[0][0][1]=2203.31;
    _sigmaR9[0][0][2]=-22454.2;
    _sigmaR9[0][0][3]=0;

    _meanScale[0][1][0]=0.889415;
    _meanScale[0][1][1]=1.21788e-05;
    _meanScale[0][1][2]=-4.3438e-06;
    _meanScale[0][1][3]=-0.629968;
    _meanAT[0][1][0]=0.0;
    _meanAT[0][1][1]=0.0;
    _meanAT[0][1][2]=0.0;
    _meanAT[0][1][3]=0.0;
    _meanAC[0][1][0]=-0.00313701;
    _meanAC[0][1][1]=0.0227998;
    _meanAC[0][1][2]=128.653;
    _meanAC[0][1][3]=234.333;
    _meanAS[0][1][0]=0.0346198;
    _meanAS[0][1][1]=-0.00261336;
    _meanAS[0][1][2]=177.983;
    _meanAS[0][1][3]=1.19839e+14;
    _meanAM[0][1][0]=-0.00100745;
    _meanAM[0][1][1]=0.264247;
    _meanAM[0][1][2]=337.255;
    _meanAM[0][1][3]=251.454;
    _meanBT[0][1][0]=0.0;
    _meanBT[0][1][1]=0.0;
    _meanBT[0][1][2]=0.0;
    _meanBT[0][1][3]=0.0;
    _meanBC[0][1][0]=-0.00397794;
    _meanBC[0][1][1]=0.00219079;
    _meanBC[0][1][2]=176.842;
    _meanBC[0][1][3]=450.29;
    _meanBS[0][1][0]=-2e+07;
    _meanBS[0][1][1]=0.0957598;
    _meanBS[0][1][2]=-8.88573e-27;
    _meanBS[0][1][3]=11442.2;
    _meanBM[0][1][0]=-0.00366315;
    _meanBM[0][1][1]=0.0622186;
    _meanBM[0][1][2]=94.5155;
    _meanBM[0][1][3]=126.404;
    _meanR9[0][1][0]=0.00636789;
    _meanR9[0][1][1]=0.000336062;
    _meanR9[0][1][2]=-0.0092699;
    _meanR9[0][1][3]=0;

    _sigmaScale[0][1][0]=0.685096;
    _sigmaScale[0][1][1]=0.129065;
    _sigmaScale[0][1][2]=-0.00212486;
    _sigmaScale[0][1][3]=0;
    _sigmaAT[0][1][0]=0.898865;
    _sigmaAT[0][1][1]=0;
    _sigmaAT[0][1][2]=0;
    _sigmaAT[0][1][3]=0;
    _sigmaAC[0][1][0]=-0.00492979;
    _sigmaAC[0][1][1]=-1.20123;
    _sigmaAC[0][1][2]=2.89231;
    _sigmaAC[0][1][3]=18.2059;
    _sigmaAS[0][1][0]=0;
    _sigmaAS[0][1][1]=0;
    _sigmaAS[0][1][2]=0;
    _sigmaAS[0][1][3]=0;
    _sigmaAM[0][1][0]=-0.000727825;
    _sigmaAM[0][1][1]=8.42395;
    _sigmaAM[0][1][2]=512.032;
    _sigmaAM[0][1][3]=415.962;
    _sigmaBT[0][1][0]=-0.0336364;
    _sigmaBT[0][1][1]=2.45182;
    _sigmaBT[0][1][2]=-0.284353;
    _sigmaBT[0][1][3]=-0.31679;
    _sigmaBC[0][1][0]=0.00510553;
    _sigmaBC[0][1][1]=-0.953869;
    _sigmaBC[0][1][2]=113872;
    _sigmaBC[0][1][3]=1.35966e+09;
    _sigmaBS[0][1][0]=0;
    _sigmaBS[0][1][1]=0;
    _sigmaBS[0][1][2]=0;
    _sigmaBS[0][1][3]=0;
    _sigmaBM[0][1][0]=-0.0034071;
    _sigmaBM[0][1][1]=4.19719;
    _sigmaBM[0][1][2]=128.952;
    _sigmaBM[0][1][3]=180.604;
    _sigmaR9[0][1][0]=-3.38988;
    _sigmaR9[0][1][1]=2.0714;
    _sigmaR9[0][1][2]=0;
    _sigmaR9[0][1][3]=0;

    _meanScale[1][0][0]=1.0009;
    _meanScale[1][0][1]=-4.79805e-06;
    _meanScale[1][0][2]=3.34625e-05;
    _meanScale[1][0][3]=-0.194267;
    _meanAT[1][0][0]=0.0;
    _meanAT[1][0][1]=0.0;
    _meanAT[1][0][2]=0.0;
    _meanAT[1][0][3]=0.0;
    _meanAC[1][0][0]=-0.000177563;
    _meanAC[1][0][1]=0.0122839;
    _meanAC[1][0][2]=1798.92;
    _meanAC[1][0][3]=776.856;
    _meanAS[1][0][0]=-0.000533039;
    _meanAS[1][0][1]=0.0642604;
    _meanAS[1][0][2]=969.596;
    _meanAS[1][0][3]=1004.15;
    _meanAM[1][0][0]=0.000163185;
    _meanAM[1][0][1]=0.085936;
    _meanAM[1][0][2]=1593.17;
    _meanAM[1][0][3]=681.623;
    _meanBT[1][0][0]=0.0;
    _meanBT[1][0][1]=0.0;
    _meanBT[1][0][2]=0.0;
    _meanBT[1][0][3]=0.0;
    _meanBC[1][0][0]=-0.000518186;
    _meanBC[1][0][1]=0.0121868;
    _meanBC[1][0][2]=1112.53;
    _meanBC[1][0][3]=933.281;
    _meanBS[1][0][0]=-0.000750734;
    _meanBS[1][0][1]=0.03859;
    _meanBS[1][0][2]=547.579;
    _meanBS[1][0][3]=775.887;
    _meanBM[1][0][0]=-0.190395;
    _meanBM[1][0][1]=-0.00362647;
    _meanBM[1][0][2]=5.25687;
    _meanBM[1][0][3]=-2.8e+08;
    _meanR9[1][0][0]=0.972346;
    _meanR9[1][0][1]=53.9185;
    _meanR9[1][0][2]=1354.5;
    _meanR9[1][0][3]=0;

    _sigmaScale[1][0][0]=0.348019;
    _sigmaScale[1][0][1]=-6.43731e-11;
    _sigmaScale[1][0][2]=0.0158647;
    _sigmaScale[1][0][3]=0;
    _sigmaAT[1][0][0]=0.215239;
    _sigmaAT[1][0][1]=0;
    _sigmaAT[1][0][2]=0;
    _sigmaAT[1][0][3]=0;
    _sigmaAC[1][0][0]=-0.00492298;
    _sigmaAC[1][0][1]=-3.40058;
    _sigmaAC[1][0][2]=17263.9;
    _sigmaAC[1][0][3]=2.6e+11;
    _sigmaAS[1][0][0]=-0.000237998;
    _sigmaAS[1][0][1]=3.0258;
    _sigmaAS[1][0][2]=1811.25;
    _sigmaAS[1][0][3]=1846.79;
    _sigmaAM[1][0][0]=0.0210134;
    _sigmaAM[1][0][1]=0.328359;
    _sigmaAM[1][0][2]=22.49;
    _sigmaAM[1][0][3]=14.5021;
    _sigmaBT[1][0][0]=-0.495072;
    _sigmaBT[1][0][1]=0;
    _sigmaBT[1][0][2]=0;
    _sigmaBT[1][0][3]=0;
    _sigmaBC[1][0][0]=-0.00265007;
    _sigmaBC[1][0][1]=0.970549;
    _sigmaBC[1][0][2]=-6.89119e+07;
    _sigmaBC[1][0][3]=180110;
    _sigmaBS[1][0][0]=0.00045833;
    _sigmaBS[1][0][1]=2.16342;
    _sigmaBS[1][0][2]=3582.4;
    _sigmaBS[1][0][3]=1100.36;
    _sigmaBM[1][0][0]=0.00188871;
    _sigmaBM[1][0][1]=1.66177;
    _sigmaBM[1][0][2]=3.2e+08;
    _sigmaBM[1][0][3]=2163.81;
    _sigmaR9[1][0][0]=-220.415;
    _sigmaR9[1][0][1]=5.19136e-08;
    _sigmaR9[1][0][2]=3.04028e-10;
    _sigmaR9[1][0][3]=0;

    _meanScale[1][1][0]=0.338011;
    _meanScale[1][1][1]=9.47815e-05;
    _meanScale[1][1][2]=-0.000238735;
    _meanScale[1][1][3]=-0.846414;
    _meanAT[1][1][0]=0.0;
    _meanAT[1][1][1]=0.0;
    _meanAT[1][1][2]=0.0;
    _meanAT[1][1][3]=0.0;
    _meanAC[1][1][0]=-0.00125367;
    _meanAC[1][1][1]=0.013324;
    _meanAC[1][1][2]=203.988;
    _meanAC[1][1][3]=431.951;
    _meanAS[1][1][0]=0.000282607;
    _meanAS[1][1][1]=0.0307431;
    _meanAS[1][1][2]=343.509;
    _meanAS[1][1][3]=274.957;
    _meanAM[1][1][0]=0.0020258;
    _meanAM[1][1][1]=0.643913;
    _meanAM[1][1][2]=0.0693877;
    _meanAM[1][1][3]=0.0816029;
    _meanBT[1][1][0]=0.0;
    _meanBT[1][1][1]=0.0;
    _meanBT[1][1][2]=0.0;
    _meanBT[1][1][3]=0.0;
    _meanBC[1][1][0]=-0.00513833;
    _meanBC[1][1][1]=5.94424e+08;
    _meanBC[1][1][2]=-62814.9;
    _meanBC[1][1][3]=118612;
    _meanBS[1][1][0]=-0.00152129;
    _meanBS[1][1][1]=0.0234694;
    _meanBS[1][1][2]=186.483;
    _meanBS[1][1][3]=754.201;
    _meanBM[1][1][0]=-0.000404987;
    _meanBM[1][1][1]=0.156384;
    _meanBM[1][1][2]=-1.7e+08;
    _meanBM[1][1][3]=1793.83;
    _meanR9[1][1][0]=0.0645278;
    _meanR9[1][1][1]=0.161614;
    _meanR9[1][1][2]=-0.215822;
    _meanR9[1][1][3]=0;

    _sigmaScale[1][1][0]=1.07376;
    _sigmaScale[1][1][1]=7.47238e-13;
    _sigmaScale[1][1][2]=0.0289594;
    _sigmaScale[1][1][3]=0;
    _sigmaAT[1][1][0]=-0.520907;
    _sigmaAT[1][1][1]=0;
    _sigmaAT[1][1][2]=0;
    _sigmaAT[1][1][3]=0;
    _sigmaAC[1][1][0]=0.00165941;
    _sigmaAC[1][1][1]=-0.351422;
    _sigmaAC[1][1][2]=8968.94;
    _sigmaAC[1][1][3]=-7e+09;
    _sigmaAS[1][1][0]=0.000490279;
    _sigmaAS[1][1][1]=0.554531;
    _sigmaAS[1][1][2]=469.111;
    _sigmaAS[1][1][3]=457.541;
    _sigmaAM[1][1][0]=0.00102079;
    _sigmaAM[1][1][1]=0.628055;
    _sigmaAM[1][1][2]=53.9452;
    _sigmaAM[1][1][3]=72.911;
    _sigmaBT[1][1][0]=-0.461542;
    _sigmaBT[1][1][1]=0;
    _sigmaBT[1][1][2]=0;
    _sigmaBT[1][1][3]=0;
    _sigmaBC[1][1][0]=-0.00219303;
    _sigmaBC[1][1][1]=0.874327;
    _sigmaBC[1][1][2]=71353.2;
    _sigmaBC[1][1][3]=2.09924e+08;
    _sigmaBS[1][1][0]=0.00104021;
    _sigmaBS[1][1][1]=0.236098;
    _sigmaBS[1][1][2]=482.954;
    _sigmaBS[1][1][3]=191.984;
    _sigmaBM[1][1][0]=-0.000116086;
    _sigmaBM[1][1][1]=2.4438;
    _sigmaBM[1][1][2]=1.9e+09;
    _sigmaBM[1][1][3]=-700.271;
    _sigmaR9[1][1][0]=4.59374;
    _sigmaR9[1][1][1]=-5.06202;
    _sigmaR9[1][1][2]=0;
    _sigmaR9[1][1][3]=0;

    _initialised=true;
  }

  if(s=="4_2") {
    _meanScale[0][0][0]=0.996757;
    _meanScale[0][0][1]=-1.02729e-05;
    _meanScale[0][0][2]=2.92397e-05;
    _meanScale[0][0][3]=-0.0892806;
    _meanAT[0][0][0]=0.000616189;
    _meanAT[0][0][1]=0;
    _meanAT[0][0][2]=0;
    _meanAT[0][0][3]=0;
    _meanAC[0][0][0]=-0.00343502;
    _meanAC[0][0][1]=0.0102123;
    _meanAC[0][0][2]=460.294;
    _meanAC[0][0][3]=510.925;
    _meanAS[0][0][0]=0;
    _meanAS[0][0][1]=0;
    _meanAS[0][0][2]=0;
    _meanAS[0][0][3]=0;
    _meanAM[0][0][0]=-0.0010363;
    _meanAM[0][0][1]=0.147364;
    _meanAM[0][0][2]=280.742;
    _meanAM[0][0][3]=200.903;
    _meanBT[0][0][0]=0;
    _meanBT[0][0][1]=0.0290516;
    _meanBT[0][0][2]=-113.185;
    _meanBT[0][0][3]=-176099;
    _meanBC[0][0][0]=-0.00262961;
    _meanBC[0][0][1]=0.0134481;
    _meanBC[0][0][2]=699.825;
    _meanBC[0][0][3]=401.452;
    _meanBS[0][0][0]=0;
    _meanBS[0][0][1]=0;
    _meanBS[0][0][2]=0;
    _meanBS[0][0][3]=0;
    _meanBM[0][0][0]=-0.00170238;
    _meanBM[0][0][1]=0.0751873;
    _meanBM[0][0][2]=232.9;
    _meanBM[0][0][3]=225.712;
    _meanR9[0][0][0]=0.946441;
    _meanR9[0][0][1]=19.5162;
    _meanR9[0][0][2]=173.503;
    _meanR9[0][0][3]=0;

    _sigmaScale[0][0][0]=0.196245;
    _sigmaScale[0][0][1]=0.0191932;
    _sigmaScale[0][0][2]=0.00660053;
    _sigmaScale[0][0][3]=0;
    _sigmaAT[0][0][0]=0.179417;
    _sigmaAT[0][0][1]=0;
    _sigmaAT[0][0][2]=0;
    _sigmaAT[0][0][3]=0;
    _sigmaAC[0][0][0]=-0.00383752;
    _sigmaAC[0][0][1]=0.0881981;
    _sigmaAC[0][0][2]=88.6384;
    _sigmaAC[0][0][3]=-3e+17;
    _sigmaAS[0][0][0]=0;
    _sigmaAS[0][0][1]=0;
    _sigmaAS[0][0][2]=0;
    _sigmaAS[0][0][3]=0;
    _sigmaAM[0][0][0]=0.000932441;
    _sigmaAM[0][0][1]=10.1637;
    _sigmaAM[0][0][2]=610.691;
    _sigmaAM[0][0][3]=290.271;
    _sigmaBT[0][0][0]=-0.00599323;
    _sigmaBT[0][0][1]=1.71943;
    _sigmaBT[0][0][2]=90.7352;
    _sigmaBT[0][0][3]=164.449;
    _sigmaBC[0][0][0]=-0.00167802;
    _sigmaBC[0][0][1]=0.303306;
    _sigmaBC[0][0][2]=1326.56;
    _sigmaBC[0][0][3]=765.552;
    _sigmaBS[0][0][0]=0;
    _sigmaBS[0][0][1]=0;
    _sigmaBS[0][0][2]=0;
    _sigmaBS[0][0][3]=0;
    _sigmaBM[0][0][0]=-0.00252598;
    _sigmaBM[0][0][1]=52.0947;
    _sigmaBM[0][0][2]=660.313;
    _sigmaBM[0][0][3]=1052.64;
    _sigmaR9[0][0][0]=0.952982;
    _sigmaR9[0][0][1]=1961.39;
    _sigmaR9[0][0][2]=22382.1;
    _sigmaR9[0][0][3]=0;

    _meanScale[0][1][0]=0.88802;
    _meanScale[0][1][1]=1.61e-05;
    _meanScale[0][1][2]=-1.08e-05;
    _meanScale[0][1][3]=-0.489;
    _meanAT[0][1][0]=-0.00165;
    _meanAT[0][1][1]=0;
    _meanAT[0][1][2]=0;
    _meanAT[0][1][3]=0;
    _meanAC[0][1][0]=-0.003195;
    _meanAC[0][1][1]=0.02012;
    _meanAC[0][1][2]=151.6;
    _meanAC[0][1][3]=255;
    _meanAS[0][1][0]=0;
    _meanAS[0][1][1]=0;
    _meanAS[0][1][2]=0;
    _meanAS[0][1][3]=0;
    _meanAM[0][1][0]=-0.00055;
    _meanAM[0][1][1]=0.1515;
    _meanAM[0][1][2]=534;
    _meanAM[0][1][3]=330;
    _meanBT[0][1][0]=0;
    _meanBT[0][1][1]=0.187;
    _meanBT[0][1][2]=299.298;
    _meanBT[0][1][3]=428.221;
    _meanBC[0][1][0]=-0.00366;
    _meanBC[0][1][1]=0.0038;
    _meanBC[0][1][2]=760;
    _meanBC[0][1][3]=760;
    _meanBS[0][1][0]=0;
    _meanBS[0][1][1]=0;
    _meanBS[0][1][2]=0;
    _meanBS[0][1][3]=0;
    _meanBM[0][1][0]=-0.00316;
    _meanBM[0][1][1]=0.0404;
    _meanBM[0][1][2]=137;
    _meanBM[0][1][3]=146;
    _meanR9[0][1][0]=0.02573;
    _meanR9[0][1][1]=-0.03722;
    _meanR9[0][1][2]=0.0144;
    _meanR9[0][1][3]=0;

    _sigmaScale[0][1][0]=0.66;
    _sigmaScale[0][1][1]=0.07465;
    _sigmaScale[0][1][2]=-0.00373;
    _sigmaScale[0][1][3]=0;
    _sigmaAT[0][1][0]=1.363;
    _sigmaAT[0][1][1]=0;
    _sigmaAT[0][1][2]=0;
    _sigmaAT[0][1][3]=0;
    _sigmaAC[0][1][0]=-0.00342;
    _sigmaAC[0][1][1]=-1.667;
    _sigmaAC[0][1][2]=6.6;
    _sigmaAC[0][1][3]=14;
    _sigmaAS[0][1][0]=0;
    _sigmaAS[0][1][1]=0;
    _sigmaAS[0][1][2]=0;
    _sigmaAS[0][1][3]=0;
    _sigmaAM[0][1][0]=-0.00094;
    _sigmaAM[0][1][1]=12.0228;
    _sigmaAM[0][1][2]=536;
    _sigmaAM[0][1][3]=464;
    _sigmaBT[0][1][0]=0.022;
    _sigmaBT[0][1][1]=3.223;
    _sigmaBT[0][1][2]=-0.361;
    _sigmaBT[0][1][3]=-0.318;
    _sigmaBC[0][1][0]=0.00567264;
    _sigmaBC[0][1][1]=-1.184;
    _sigmaBC[0][1][2]=76000;
    _sigmaBC[0][1][3]=2.14748e+08;
    _sigmaBS[0][1][0]=0;
    _sigmaBS[0][1][1]=0;
    _sigmaBS[0][1][2]=0;
    _sigmaBS[0][1][3]=0;
    _sigmaBM[0][1][0]=-0.00236;
    _sigmaBM[0][1][1]=5.78;
    _sigmaBM[0][1][2]=166;
    _sigmaBM[0][1][3]=203.5;
    _sigmaR9[0][1][0]=-2.94;
    _sigmaR9[0][1][1]=1.53445;
    _sigmaR9[0][1][2]=0;
    _sigmaR9[0][1][3]=0;

    _meanScale[1][0][0]=0.993921;
    _meanScale[1][0][1]=-8.11417e-06;
    _meanScale[1][0][2]=3.23086e-05;
    _meanScale[1][0][3]=-0.0976369;
    _meanAT[1][0][0]=0.0691979;
    _meanAT[1][0][1]=0;
    _meanAT[1][0][2]=0;
    _meanAT[1][0][3]=0;
    _meanAC[1][0][0]=-0.000310485;
    _meanAC[1][0][1]=0.0169478;
    _meanAC[1][0][2]=1301.97;
    _meanAC[1][0][3]=796.53;
    _meanAS[1][0][0]=-0.0005425;
    _meanAS[1][0][1]=0.0676613;
    _meanAS[1][0][2]=790.284;
    _meanAS[1][0][3]=970;
    _meanAM[1][0][0]=0.000134362;
    _meanAM[1][0][1]=0.364822;
    _meanAM[1][0][2]=3396.04;
    _meanAM[1][0][3]=2407.2;
    _meanBT[1][0][0]=0.0456649;
    _meanBT[1][0][1]=0;
    _meanBT[1][0][2]=0;
    _meanBT[1][0][3]=0;
    _meanBC[1][0][0]=-0.000356413;
    _meanBC[1][0][1]=0.0167018;
    _meanBC[1][0][2]=1394.52;
    _meanBC[1][0][3]=823.035;
    _meanBS[1][0][0]=-0.000484;
    _meanBS[1][0][1]=0.0542618;
    _meanBS[1][0][2]=798.07;
    _meanBS[1][0][3]=1019.88;
    _meanBM[1][0][0]=-0.175333;
    _meanBM[1][0][1]=-0.00324513;
    _meanBM[1][0][2]=2.87972;
    _meanBM[1][0][3]=7.58712e+12;
    _meanR9[1][0][0]=0.96393;
    _meanR9[1][0][1]=29.6072;
    _meanR9[1][0][2]=334.89;
    _meanR9[1][0][3]=0;

    _sigmaScale[1][0][0]=0.337911;
    _sigmaScale[1][0][1]=1.03419e-11;
    _sigmaScale[1][0][2]=0.0153846;
    _sigmaScale[1][0][3]=0;
    _sigmaAT[1][0][0]=1.12;
    _sigmaAT[1][0][1]=0;
    _sigmaAT[1][0][2]=0;
    _sigmaAT[1][0][3]=0;
    _sigmaAC[1][0][0]=48.1275;
    _sigmaAC[1][0][1]=1.50005e+08;
    _sigmaAC[1][0][2]=21231.6;
    _sigmaAC[1][0][3]=2.6e+11;
    _sigmaAS[1][0][0]=-0.00036952;
    _sigmaAS[1][0][1]=2.75194;
    _sigmaAS[1][0][2]=1045.06;
    _sigmaAS[1][0][3]=1420.87;
    _sigmaAM[1][0][0]=0.05567;
    _sigmaAM[1][0][1]=0.32;
    _sigmaAM[1][0][2]=420;
    _sigmaAM[1][0][3]=12.8233;
    _sigmaBT[1][0][0]=-0.22312;
    _sigmaBT[1][0][1]=0;
    _sigmaBT[1][0][2]=0;
    _sigmaBT[1][0][3]=0;
    _sigmaBC[1][0][0]=-0.00320087;
    _sigmaBC[1][0][1]=59.2836;
    _sigmaBC[1][0][2]=7.49509e+12;
    _sigmaBC[1][0][3]=1.32714e+07;
    _sigmaBS[1][0][0]=0.000189384;
    _sigmaBS[1][0][1]=2.89111;
    _sigmaBS[1][0][2]=2468.02;
    _sigmaBS[1][0][3]=1294.69;
    _sigmaBM[1][0][0]=0.00154629;
    _sigmaBM[1][0][1]=16.9762;
    _sigmaBM[1][0][2]=1.05e+07;
    _sigmaBM[1][0][3]=-9e+06;
    _sigmaR9[1][0][0]=99.68;
    _sigmaR9[1][0][1]=-9.37265e-08;
    _sigmaR9[1][0][2]=6.91613e-09;
    _sigmaR9[1][0][3]=0;

    _meanScale[1][1][0]=0.340125;
    _meanScale[1][1][1]=-5.02007e-05;
    _meanScale[1][1][2]=0.000208216;
    _meanScale[1][1][3]=-1.32063;
    _meanAT[1][1][0]=-0.092494;
    _meanAT[1][1][1]=0;
    _meanAT[1][1][2]=0;
    _meanAT[1][1][3]=0;
    _meanAC[1][1][0]=-0.00131984;
    _meanAC[1][1][1]=0.00950668;
    _meanAC[1][1][2]=259.804;
    _meanAC[1][1][3]=705.003;
    _meanAS[1][1][0]=-0.00705787;
    _meanAS[1][1][1]=-0.00581517;
    _meanAS[1][1][2]=7.10951e+10;
    _meanAS[1][1][3]=-68.7322;
    _meanAM[1][1][0]=-0.00496966;
    _meanAM[1][1][1]=0.655834;
    _meanAM[1][1][2]=0.0491064;
    _meanAM[1][1][3]=0.0632471;
    _meanBT[1][1][0]=-0.158344;
    _meanBT[1][1][1]=0;
    _meanBT[1][1][2]=0;
    _meanBT[1][1][3]=0;
    _meanBC[1][1][0]=-0.00498431;
    _meanBC[1][1][1]=-18315.4;
    _meanBC[1][1][2]=50760;
    _meanBC[1][1][3]=118612;
    _meanBS[1][1][0]=-0.000620577;
    _meanBS[1][1][1]=0.0247672;
    _meanBS[1][1][2]=437.005;
    _meanBS[1][1][3]=374.823;
    _meanBM[1][1][0]=0.000431046;
    _meanBM[1][1][1]=0.0224695;
    _meanBM[1][1][2]=-1.06e+07;
    _meanBM[1][1][3]=-135.039;
    _meanR9[1][1][0]=0.151826;
    _meanR9[1][1][1]=0.00348896;
    _meanR9[1][1][2]=-0.161781;
    _meanR9[1][1][3]=0;

    _sigmaScale[1][1][0]=1.19993;
    _sigmaScale[1][1][1]=1.38778e-17;
    _sigmaScale[1][1][2]=0.0251314;
    _sigmaScale[1][1][3]=0;
    _sigmaAT[1][1][0]=-0.791929;
    _sigmaAT[1][1][1]=0;
    _sigmaAT[1][1][2]=0;
    _sigmaAT[1][1][3]=0;
    _sigmaAC[1][1][0]=0.0015505;
    _sigmaAC[1][1][1]=-0.530085;
    _sigmaAC[1][1][2]=33506.7;
    _sigmaAC[1][1][3]=-7e+09;
    _sigmaAS[1][1][0]=7.98682e-05;
    _sigmaAS[1][1][1]=0.875124;
    _sigmaAS[1][1][2]=483.102;
    _sigmaAS[1][1][3]=1005.71;
    _sigmaAM[1][1][0]=-0.00510653;
    _sigmaAM[1][1][1]=1.05312;
    _sigmaAM[1][1][2]=49.6466;
    _sigmaAM[1][1][3]=7.11474e+12;
    _sigmaBT[1][1][0]=-1.30079;
    _sigmaBT[1][1][1]=0;
    _sigmaBT[1][1][2]=0;
    _sigmaBT[1][1][3]=0;
    _sigmaBC[1][1][0]=-0.00202301;
    _sigmaBC[1][1][1]=4.19368;
    _sigmaBC[1][1][2]=-1.08161e+07;
    _sigmaBC[1][1][3]=-6e+09;
    _sigmaBS[1][1][0]=0.00182147;
    _sigmaBS[1][1][1]=0.25995;
    _sigmaBS[1][1][2]=881.576;
    _sigmaBS[1][1][3]=304.891;
    _sigmaBM[1][1][0]=-0.0454995;
    _sigmaBM[1][1][1]=7.63845;
    _sigmaBM[1][1][2]=1.9e+09;
    _sigmaBM[1][1][3]=150001;
    _sigmaR9[1][1][0]=5.69741;
    _sigmaR9[1][1][1]=-6.23853;
    _sigmaR9[1][1][2]=0;
    _sigmaR9[1][1][3]=0;

    _initialised=true;
  }

  if(s=="4_2e") {
    _meanScale[0][0][0]=1.02133;
    _meanScale[0][0][1]=-0.000183073;
    _meanScale[0][0][2]=0.000252052;
    _meanScale[0][0][3]=0.300892;
    _meanAT[0][0][0]=0.016876;
    _meanAT[0][0][1]=0;
    _meanAT[0][0][2]=0;
    _meanAT[0][0][3]=0;
    _meanAC[0][0][0]=-0.00327969;
    _meanAC[0][0][1]=0.00932112;
    _meanAC[0][0][2]=430.176;
    _meanAC[0][0][3]=440.94;
    _meanAS[0][0][0]=0;
    _meanAS[0][0][1]=0;
    _meanAS[0][0][2]=0;
    _meanAS[0][0][3]=0;
    _meanAM[0][0][0]=-0.000876726;
    _meanAM[0][0][1]=0.128953;
    _meanAM[0][0][2]=262.609;
    _meanAM[0][0][3]=204.869;
    _meanBT[0][0][0]=0;
    _meanBT[0][0][1]=0;
    _meanBT[0][0][2]=0;
    _meanBT[0][0][3]=0;
    _meanBC[0][0][0]=-0.00316342;
    _meanBC[0][0][1]=0.00801235;
    _meanBC[0][0][2]=552.944;
    _meanBC[0][0][3]=507.239;
    _meanBS[0][0][0]=0;
    _meanBS[0][0][1]=0;
    _meanBS[0][0][2]=0;
    _meanBS[0][0][3]=0;
    _meanBM[0][0][0]=-0.00171652;
    _meanBM[0][0][1]=0.0746995;
    _meanBM[0][0][2]=190.606;
    _meanBM[0][0][3]=210.147;
    _meanR9[0][0][0]=0.847311;
    _meanR9[0][0][1]=-10.0049;
    _meanR9[0][0][2]=67.613;
    _meanR9[0][0][3]=0;

    _sigmaScale[0][0][0]=0.345777;
    _sigmaScale[0][0][1]=0.0375747;
    _sigmaScale[0][0][2]=-0.00516299;
    _sigmaScale[0][0][3]=0;
    _sigmaAT[0][0][0]=0.994444;
    _sigmaAT[0][0][1]=0;
    _sigmaAT[0][0][2]=0;
    _sigmaAT[0][0][3]=0;
    _sigmaAC[0][0][0]=-0.00380373;
    _sigmaAC[0][0][1]=-1.00415;
    _sigmaAC[0][0][2]=110274;
    _sigmaAC[0][0][3]=-3e+17;
    _sigmaAS[0][0][0]=0;
    _sigmaAS[0][0][1]=0;
    _sigmaAS[0][0][2]=0;
    _sigmaAS[0][0][3]=0;
    _sigmaAM[0][0][0]=-0.000881147;
    _sigmaAM[0][0][1]=8.76414;
    _sigmaAM[0][0][2]=455.387;
    _sigmaAM[0][0][3]=412.382;
    _sigmaBT[0][0][0]=0.00578517;
    _sigmaBT[0][0][1]=2.28008;
    _sigmaBT[0][0][2]=-3.36e+07;
    _sigmaBT[0][0][3]=-152.696;
    _sigmaBC[0][0][0]=-0.00161997;
    _sigmaBC[0][0][1]=50513.8;
    _sigmaBC[0][0][2]=-1.26888e+07;
    _sigmaBC[0][0][3]=3.75561e+07;
    _sigmaBS[0][0][0]=0;
    _sigmaBS[0][0][1]=0;
    _sigmaBS[0][0][2]=0;
    _sigmaBS[0][0][3]=0;
    _sigmaBM[0][0][0]=-0.00168415;
    _sigmaBM[0][0][1]=6.89324;
    _sigmaBM[0][0][2]=272.169;
    _sigmaBM[0][0][3]=315.424;
    _sigmaR9[0][0][0]=0.952571;
    _sigmaR9[0][0][1]=0;
    _sigmaR9[0][0][2]=0;
    _sigmaR9[0][0][3]=0;

    _meanScale[0][1][0]=0.870259;
    _meanScale[0][1][1]=0.000234598;
    _meanScale[0][1][2]=-0.000262422;
    _meanScale[0][1][3]=0.277914;
    _meanAT[0][1][0]=1e-06;
    _meanAT[0][1][1]=0;
    _meanAT[0][1][2]=0;
    _meanAT[0][1][3]=0;
    _meanAC[0][1][0]=-0.00324258;
    _meanAC[0][1][1]=0.0209948;
    _meanAC[0][1][2]=167.983;
    _meanAC[0][1][3]=305.464;
    _meanAS[0][1][0]=0.0333341;
    _meanAS[0][1][1]=-591.066;
    _meanAS[0][1][2]=1.48079e+11;
    _meanAS[0][1][3]=1.61987e+07;
    _meanAM[0][1][0]=-0.000991417;
    _meanAM[0][1][1]=0.121994;
    _meanAM[0][1][2]=367.004;
    _meanAM[0][1][3]=289.22;
    _meanBT[0][1][0]=0;
    _meanBT[0][1][1]=0.187;
    _meanBT[0][1][2]=299.298;
    _meanBT[0][1][3]=428.221;
    _meanBC[0][1][0]=-0.00437016;
    _meanBC[0][1][1]=0.00235015;
    _meanBC[0][1][2]=759.85;
    _meanBC[0][1][3]=-5018.79;
    _meanBS[0][1][0]=-1.2e+08;
    _meanBS[0][1][1]=0.0934286;
    _meanBS[0][1][2]=-1.226e-16;
    _meanBS[0][1][3]=11442.2;
    _meanBM[0][1][0]=-0.000153536;
    _meanBM[0][1][1]=0.0101487;
    _meanBM[0][1][2]=60.0003;
    _meanBM[0][1][3]=-4.4e+08;
    _meanR9[0][1][0]=0.0853634;
    _meanR9[0][1][1]=-0.263819;
    _meanR9[0][1][2]=0.200892;
    _meanR9[0][1][3]=0;

    _sigmaScale[0][1][0]=0.705452;
    _sigmaScale[0][1][1]=-0.0781086;
    _sigmaScale[0][1][2]=0.0039561;
    _sigmaScale[0][1][3]=0;
    _sigmaAT[0][1][0]=1.46888;
    _sigmaAT[0][1][1]=0;
    _sigmaAT[0][1][2]=0;
    _sigmaAT[0][1][3]=0;
    _sigmaAC[0][1][0]=-0.00537924;
    _sigmaAC[0][1][1]=-0.740685;
    _sigmaAC[0][1][2]=1.59714;
    _sigmaAC[0][1][3]=31.9922;
    _sigmaAS[0][1][0]=0;
    _sigmaAS[0][1][1]=0;
    _sigmaAS[0][1][2]=0;
    _sigmaAS[0][1][3]=0;
    _sigmaAM[0][1][0]=0.000163134;
    _sigmaAM[0][1][1]=7.30822;
    _sigmaAM[0][1][2]=632.013;
    _sigmaAM[0][1][3]=373.003;
    _sigmaBT[0][1][0]=0.00253725;
    _sigmaBT[0][1][1]=2.97332;
    _sigmaBT[0][1][2]=-0.516548;
    _sigmaBT[0][1][3]=-0.50994;
    _sigmaBC[0][1][0]=0.00578715;
    _sigmaBC[0][1][1]=-0.473225;
    _sigmaBC[0][1][2]=8363.08;
    _sigmaBC[0][1][3]=2.36808e+19;
    _sigmaBS[0][1][0]=0;
    _sigmaBS[0][1][1]=0;
    _sigmaBS[0][1][2]=0;
    _sigmaBS[0][1][3]=0;
    _sigmaBM[0][1][0]=0.000184199;
    _sigmaBM[0][1][1]=3.82105;
    _sigmaBM[0][1][2]=93.4062;
    _sigmaBM[0][1][3]=95.4072;
    _sigmaR9[0][1][0]=1.20664;
    _sigmaR9[0][1][1]=-3.43548;
    _sigmaR9[0][1][2]=0;
    _sigmaR9[0][1][3]=0;

    _meanScale[1][0][0]=0.966787;
    _meanScale[1][0][1]=4.77112e-05;
    _meanScale[1][0][2]=-4.47768e-05;
    _meanScale[1][0][3]=0.0854689;
    _meanAT[1][0][0]=0.0691979;
    _meanAT[1][0][1]=0;
    _meanAT[1][0][2]=0;
    _meanAT[1][0][3]=0;
    _meanAC[1][0][0]=-0.000985244;
    _meanAC[1][0][1]=0.0132327;
    _meanAC[1][0][2]=530.13;
    _meanAC[1][0][3]=1008.64;
    _meanAS[1][0][0]=-0.000571197;
    _meanAS[1][0][1]=0.0579624;
    _meanAS[1][0][2]=687.722;
    _meanAS[1][0][3]=836.21;
    _meanAM[1][0][0]=0.00052402;
    _meanAM[1][0][1]=0.402941;
    _meanAM[1][0][2]=3615.29;
    _meanAM[1][0][3]=1495.1;
    _meanBT[1][0][0]=0.0456649;
    _meanBT[1][0][1]=0;
    _meanBT[1][0][2]=0;
    _meanBT[1][0][3]=0;
    _meanBC[1][0][0]=-0.000675969;
    _meanBC[1][0][1]=0.0118283;
    _meanBC[1][0][2]=1050.1;
    _meanBC[1][0][3]=1329.8;
    _meanBS[1][0][0]=-0.000254569;
    _meanBS[1][0][1]=0.0576783;
    _meanBS[1][0][2]=853.571;
    _meanBS[1][0][3]=938.698;
    _meanBM[1][0][0]=-0.186506;
    _meanBM[1][0][1]=-0.00654349;
    _meanBM[1][0][2]=1.09286;
    _meanBM[1][0][3]=-5.8e+08;
    _meanR9[1][0][0]=0.951164;
    _meanR9[1][0][1]=-0.323083;
    _meanR9[1][0][2]=506.054;
    _meanR9[1][0][3]=0;

    _sigmaScale[1][0][0]=0.443288;
    _sigmaScale[1][0][1]=1.73472e-18;
    _sigmaScale[1][0][2]=0.0134633;
    _sigmaScale[1][0][3]=0;
    _sigmaAT[1][0][0]=6.60294;
    _sigmaAT[1][0][1]=0;
    _sigmaAT[1][0][2]=0;
    _sigmaAT[1][0][3]=0;
    _sigmaAC[1][0][0]=48.1275;
    _sigmaAC[1][0][1]=1.50005e+08;
    _sigmaAC[1][0][2]=21231.6;
    _sigmaAC[1][0][3]=2.6e+11;
    _sigmaAS[1][0][0]=0.000112234;
    _sigmaAS[1][0][1]=2.38201;
    _sigmaAS[1][0][2]=1368.36;
    _sigmaAS[1][0][3]=1014.58;
    _sigmaAM[1][0][0]=0.0227945;
    _sigmaAM[1][0][1]=1.35035;
    _sigmaAM[1][0][2]=34.8078;
    _sigmaAM[1][0][3]=41.2653;
    _sigmaBT[1][0][0]=4.6512;
    _sigmaBT[1][0][1]=0;
    _sigmaBT[1][0][2]=0;
    _sigmaBT[1][0][3]=0;
    _sigmaBC[1][0][0]=-0.00385034;
    _sigmaBC[1][0][1]=9860.11;
    _sigmaBC[1][0][2]=7.49509e+12;
    _sigmaBC[1][0][3]=5.55794e+07;
    _sigmaBS[1][0][0]=0.000287389;
    _sigmaBS[1][0][1]=2.22754;
    _sigmaBS[1][0][2]=2066.43;
    _sigmaBS[1][0][3]=926.247;
    _sigmaBM[1][0][0]=0.00110983;
    _sigmaBM[1][0][1]=4.70582;
    _sigmaBM[1][0][2]=8.47552e+10;
    _sigmaBM[1][0][3]=-4572.21;
    _sigmaR9[1][0][0]=141.283;
    _sigmaR9[1][0][1]=-6.07621e-07;
    _sigmaR9[1][0][2]=1.40907e-08;
    _sigmaR9[1][0][3]=0;

    _meanScale[1][1][0]=0.290862;
    _meanScale[1][1][1]=-5.93303e-05;
    _meanScale[1][1][2]=0.000431058;
    _meanScale[1][1][3]=0.350738;
    _meanAT[1][1][0]=-0.092494;
    _meanAT[1][1][1]=0;
    _meanAT[1][1][2]=0;
    _meanAT[1][1][3]=0;
    _meanAC[1][1][0]=-0.00105099;
    _meanAC[1][1][1]=0.0195863;
    _meanAC[1][1][2]=162.87;
    _meanAC[1][1][3]=259.479;
    _meanAS[1][1][0]=-0.00820192;
    _meanAS[1][1][1]=-0.00404515;
    _meanAS[1][1][2]=-9e+12;
    _meanAS[1][1][3]=-12.1537;
    _meanAM[1][1][0]=0.00430129;
    _meanAM[1][1][1]=0.647006;
    _meanAM[1][1][2]=0.0524654;
    _meanAM[1][1][3]=0.057336;
    _meanBT[1][1][0]=-0.158344;
    _meanBT[1][1][1]=0;
    _meanBT[1][1][2]=0;
    _meanBT[1][1][3]=0;
    _meanBC[1][1][0]=-0.00661496;
    _meanBC[1][1][1]=-0.0551321;
    _meanBC[1][1][2]=1286.49;
    _meanBC[1][1][3]=118612;
    _meanBS[1][1][0]=-0.000517633;
    _meanBS[1][1][1]=0.0294134;
    _meanBS[1][1][2]=414.964;
    _meanBS[1][1][3]=487.294;
    _meanBM[1][1][0]=0.00298075;
    _meanBM[1][1][1]=0.0146853;
    _meanBM[1][1][2]=-7.17525e+06;
    _meanBM[1][1][3]=41.9253;
    _meanR9[1][1][0]=0.213804;
    _meanR9[1][1][1]=-0.395466;
    _meanR9[1][1][2]=0.202973;
    _meanR9[1][1][3]=0;

    _sigmaScale[1][1][0]=1.29656;
    _sigmaScale[1][1][1]=9.61084e-11;
    _sigmaScale[1][1][2]=0.0270649;
    _sigmaScale[1][1][3]=0;
    _sigmaAT[1][1][0]=-0.0115339;
    _sigmaAT[1][1][1]=0;
    _sigmaAT[1][1][2]=0;
    _sigmaAT[1][1][3]=0;
    _sigmaAC[1][1][0]=0.0013797;
    _sigmaAC[1][1][1]=-0.445789;
    _sigmaAC[1][1][2]=29539.8;
    _sigmaAC[1][1][3]=-3.6e+08;
    _sigmaAS[1][1][0]=0.000213606;
    _sigmaAS[1][1][1]=0.547654;
    _sigmaAS[1][1][2]=892.47;
    _sigmaAS[1][1][3]=408.826;
    _sigmaAM[1][1][0]=-0.00121353;
    _sigmaAM[1][1][1]=1.10052;
    _sigmaAM[1][1][2]=31.303;
    _sigmaAM[1][1][3]=80.0479;
    _sigmaBT[1][1][0]=0.470361;
    _sigmaBT[1][1][1]=0;
    _sigmaBT[1][1][2]=0;
    _sigmaBT[1][1][3]=0;
    _sigmaBC[1][1][0]=-0.00222463;
    _sigmaBC[1][1][1]=18440.6;
    _sigmaBC[1][1][2]=5.80805e+07;
    _sigmaBC[1][1][3]=-4.43323e+17;
    _sigmaBS[1][1][0]=0.00156455;
    _sigmaBS[1][1][1]=0.337173;
    _sigmaBS[1][1][2]=281.318;
    _sigmaBS[1][1][3]=252.457;
    _sigmaBM[1][1][0]=-0.0455886;
    _sigmaBM[1][1][1]=1.47875;
    _sigmaBM[1][1][2]=1.9e+09;
    _sigmaBM[1][1][3]=1193.93;
    _sigmaR9[1][1][0]=8.89993;
    _sigmaR9[1][1][1]=-9.8386;
    _sigmaR9[1][1][2]=0;
    _sigmaR9[1][1][3]=0;

    _initialised=true;
  }

  assert(_initialised);
  return true;
}

double PhotonFix::phi

(

void

)

const

Definition at line 246 of file PhotonFix.cc.

References _phi.

                            {
  return _phi;
}

double PhotonFix::phi

(

)

const

double PhotonFix::phiC

(

)

const

double PhotonFix::phiC

(

)

const

Definition at line 269 of file PhotonFix.cc.

References _bC, and _be.

                             {
  assert(_be==0);
  return _bC;
}

double PhotonFix::phiM

(

)

const

Definition at line 279 of file PhotonFix.cc.

References _be, and _bM.

                             {
  assert(_be==0);
  return _bM;
}

double PhotonFix::phiM

(

)

const

double PhotonFix::phiS

(

)

const

double PhotonFix::phiS

(

)

const

Definition at line 274 of file PhotonFix.cc.

References _be, and _bS.

                             {
  assert(_be==0);
  return _bS;
}

void PhotonFix::print

(

void

)

const

Definition at line 351 of file PhotonFix.cc.

References _aC, _aM, _aS, _bC, _bM, _bS, _e, _eta, _phi, _r9, and gather_cfg::cout.

                            {
  std::cout << "PhotonFix:  e,eta,phi,r9 = " << _e << ", " << _eta << ", " << _phi << ", " << _r9 << ", gaps "
            << _aC << ", " << _aS << ", " << _aM << ", "
            << _bC << ", " << _bS << ", " << _bM << std::endl;
}

void PhotonFix::print

(

)

const

void PhotonFix::printParameters

(

std::ostream &

o

)

[static]

Definition at line 479 of file PhotonFix.cc.

References _meanAC, _meanAM, _meanAS, _meanAT, _meanBC, _meanBM, _meanBS, _meanBT, _meanR9, _meanScale, _sigmaAC, _sigmaAM, _sigmaAS, _sigmaAT, _sigmaBC, _sigmaBM, _sigmaBS, _sigmaBT, _sigmaR9, _sigmaScale, and i.

                                             {
  o << "PhotonFix::printParameters()" << std::endl;
  
  for(unsigned be(0);be<2;be++) {
    for(unsigned hl(0);hl<2;hl++) {
      o << " Parameters for " << (be==0?"barrel":"endcap")
        << ", " << (hl==0?"high":"low") << " R9" << std::endl;
      
      o << "  Mean  scaling        ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanScale[be][hl][i];
      o << std::endl;
      o << "  Mean  " << (be==0?"Eta  ":"ZetaX") << " total    ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanAT[be][hl][i];
      o << std::endl;
      o << "  Mean  " << (be==0?"Eta  ":"ZetaX") << " crystal  ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanAC[be][hl][i];
      o << std::endl;
      o << "  Mean  " << (be==0?"Eta  ":"ZetaX") << " submodule";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanAS[be][hl][i];
      o << std::endl;
      o << "  Mean  " << (be==0?"Eta  ":"ZetaX") << " module   ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanAM[be][hl][i];
      o << std::endl;
      o << "  Mean  " << (be==0?"Eta   zero     ":"ZetaY total    ");
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanBT[be][hl][i];
      o << std::endl;
      o << "  Mean  " << (be==0?"Phi  ":"ZetaY") << " crystal  ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanBC[be][hl][i];
      o << std::endl;
      o << "  Mean  " << (be==0?"Phi  ":"ZetaY") << " submodule";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanBS[be][hl][i];
      o << std::endl;
      o << "  Mean  " << (be==0?"Phi  ":"ZetaY") << " module   ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanBM[be][hl][i];
      o << std::endl;
      o << "  Mean  R9             ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _meanR9[be][hl][i];
      o << std::endl;
      
      o << "  Sigma scaling        ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaScale[be][hl][i];
      o << std::endl;
      o << "  Sigma " << (be==0?"Eta  ":"ZetaX") << " total    ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaAT[be][hl][i];
      o << std::endl;
      o << "  Sigma " << (be==0?"Eta  ":"ZetaX") << " crystal  ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaAC[be][hl][i];
      o << std::endl;
      o << "  Sigma " << (be==0?"Eta  ":"ZetaX") << " submodule";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaAS[be][hl][i];
      o << std::endl;
      o << "  Sigma " << (be==0?"Eta  ":"ZetaX") << " module   ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaAM[be][hl][i];
      o << std::endl;
      o << "  Sigma " << (be==0?"Eta  ":"ZetaY") << " total    ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaBT[be][hl][i];
      o << std::endl;
      o << "  Sigma " << (be==0?"Eta  ":"ZetaY") << " crystal  ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaBC[be][hl][i];
      o << std::endl;
      o << "  Sigma " << (be==0?"Phi  ":"ZetaY") << " submodule";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaBS[be][hl][i];
      o << std::endl;
      o << "  Sigma " << (be==0?"Phi  ":"ZetaY") << " module   ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaBM[be][hl][i];
      o << std::endl;
      o << "  Sigma R9             ";
      for(unsigned i(0);i<4;i++) o << std::setw(14) << _sigmaR9[be][hl][i];
      o << std::endl;
    }
  }
}

static void PhotonFix::printParameters

(

std::ostream &

o

)

[static]

double PhotonFix::r9

(

)

const

Definition at line 250 of file PhotonFix.cc.

References _r9.

                           {
  return _r9;
}

double PhotonFix::r9

(

)

const

double PhotonFix::rawEnergy

(

)

const

Definition at line 238 of file PhotonFix.cc.

References _e.

                                  {
  return _e;
}

double PhotonFix::rawEnergy

(

)

const

void PhotonFix::setParameters	(	unsigned	be,
		unsigned	hl,
		const double *	p
	)		[static]

Definition at line 357 of file PhotonFix.cc.

References _meanAC, _meanAM, _meanAS, _meanAT, _meanBC, _meanBM, _meanBS, _meanBT, _meanR9, _meanScale, _sigmaAC, _sigmaAM, _sigmaAS, _sigmaAT, _sigmaBC, _sigmaBM, _sigmaBS, _sigmaBT, _sigmaR9, _sigmaScale, and i.

                                                                       {
  for(unsigned i(0);i<4;i++) {
    _meanScale[be][hl][i] =p[i+ 0*4];
    _meanAT[be][hl][i]    =p[i+ 1*4];
    _meanAC[be][hl][i]    =p[i+ 2*4];
    _meanAS[be][hl][i]    =p[i+ 3*4];
    _meanAM[be][hl][i]    =p[i+ 4*4];
    _meanBT[be][hl][i]    =p[i+ 5*4];
    _meanBC[be][hl][i]    =p[i+ 6*4];
    _meanBS[be][hl][i]    =p[i+ 7*4];
    _meanBM[be][hl][i]    =p[i+ 8*4];
    _meanR9[be][hl][i]    =p[i+ 9*4];
    
    _sigmaScale[be][hl][i]=p[i+10*4];
    _sigmaAT[be][hl][i]   =p[i+11*4];
    _sigmaAC[be][hl][i]   =p[i+12*4];
    _sigmaAS[be][hl][i]   =p[i+13*4];
    _sigmaAM[be][hl][i]   =p[i+14*4];
    _sigmaBT[be][hl][i]   =p[i+15*4];
    _sigmaBC[be][hl][i]   =p[i+16*4];
    _sigmaBS[be][hl][i]   =p[i+17*4];
    _sigmaBM[be][hl][i]   =p[i+18*4];
    _sigmaR9[be][hl][i]   =p[i+19*4];
  }
}

static void PhotonFix::setParameters	(	unsigned	be,
		unsigned	hl,
		const double *	p
	)		[static]

void PhotonFix::setup

(

)

void PhotonFix::setup

(

void

)

Definition at line 15 of file PhotonFix.cc.

References _aC, _aM, _aS, _barrelCGap, _barrelMGap, _barrelSGap, _bC, _be, _bM, _bS, _endcapCGap, _endcapMGap, _endcapSGap, _eta, _hl, _initialised, _phi, _r9, dPhi(), i, j, alignCSCRings::r, diffTwoXMLs::r2, xZ(), and yZ().

Referenced by PhotonFix().

                     {
  // Check constants have been set up
  assert(_initialised);

  // Determine if EB or EE
  _be=(fabs(_eta)<1.48?0:1);
  
  // Determine if high or low R9
  if(_be==0) _hl=(_r9>=0.94?0:1);
  else       _hl=(_r9>=0.95?0:1);
  
  // Coordinates relative to cracks
  double r2Min;
  if(_be==0) {
    
    r2Min=1.0e6;
    for(unsigned i(0);i<169;i++) {
      for(unsigned j(0);j<360;j++) {
        double de(_eta-_barrelCGap[i][j][0]);
        double df(dPhi(_phi,_barrelCGap[i][j][1]));
        double r2(de*de+df*df);
        
        if(r2<r2Min) {
          r2Min=r2;
          if(i>=84) {
            _aC= de;
            _bC=-df;
          } else {
            _aC=-de;
            _bC= df;
          }
        }
      }
    }
      
    r2Min=1.0e6;
    for(unsigned i(0);i<33;i++) {
      for(unsigned j(0);j<180;j++) {
        double de(_eta-_barrelSGap[i][j][0]);
        double df(dPhi(_phi,_barrelSGap[i][j][1]));
        double r2(de*de+df*df);
          
        if(r2<r2Min) {
          r2Min=r2;
          if(i>=16) {
            _aS= de;
            _bS=-df;
          } else {
            _aS=-de;
            _bS= df;
          }
        }
      }
    }
      
    r2Min=1.0e6;
    for(unsigned i(0);i<7;i++) {
      for(unsigned j(0);j<18;j++) {
        double de(_eta-_barrelMGap[i][j][0]);
        double df(dPhi(_phi,_barrelMGap[i][j][1]));
        double r2(de*de+df*df);
          
        if(r2<r2Min) {
          r2Min=r2;
          if(i>=3) {
            _aM= de;
            _bM=-df;
          } else {
            _aM=-de;
            _bM= df;
          }
        }
      }
    }

  } else {
    unsigned iz(_eta>=0.0?0:1);
    double r[2]={xZ(),yZ()};

    r2Min=1.0e6;
    for(unsigned i(0);i<7080;i++) {
      double dx(r[0]-_endcapCGap[iz][i][0]);
      double dy(r[1]-_endcapCGap[iz][i][1]);
      double r2(dx*dx+dy*dy);

      if(r2<r2Min) {
        r2Min=r2;
        if(r[0]>0.0) _aC= dx;
        else         _aC=-dx;
        if(r[1]>0.0) _bC= dy;
        else         _bC=-dy;
      }
    }

    r2Min=1.0e6;
    for(unsigned i(0);i<264;i++) {
      double dx(r[0]-_endcapSGap[iz][i][0]);
      double dy(r[1]-_endcapSGap[iz][i][1]);
      double r2(dx*dx+dy*dy);

      if(r2<r2Min) {
        r2Min=r2;
        if(r[0]>0.0) _aS= dx;
        else         _aS=-dx;
        if(r[1]>0.0) _bS= dy;
        else         _bS=-dy;
      }
    }

    r2Min=1.0e6;
    for(unsigned i(0);i<1;i++) {
      double dx(r[0]-_endcapMGap[iz][i][0]);
      double dy(r[1]-_endcapMGap[iz][i][1]);
      double r2(dx*dx+dy*dy);

      if(r2<r2Min) {
        r2Min=r2;
        if(iz==0) {_aM= dx;_bM= dy;}
        else      {_aM=-dx;_bM=-dy;}
      }
    }
  }
}

double PhotonFix::sigmaEnergy

(

)

const

double PhotonFix::sigmaEnergy

(

)

const

Definition at line 178 of file PhotonFix.cc.

References _aC, _aM, _aS, _bC, _be, _bM, _bS, _e, _eta, _hl, _r9, _sigmaAC, _sigmaAM, _sigmaAS, _sigmaAT, _sigmaBC, _sigmaBM, _sigmaBS, _sigmaBT, _sigmaR9, _sigmaScale, expCorrection(), f, mathSSE::sqrt(), xZ(), and yZ().

Referenced by PhotonFixCMS::sigmaEnergy().

                                    {
  
  // Overall resolution scale vs energy
  double sigma;
  if(_be==0) {
    sigma =_sigmaScale[_be][_hl][0]*_sigmaScale[_be][_hl][0];
    //std::cout << "PhotonFix::sigmaEnergy 1 sigma = " << sigma << std::endl;
    sigma+=_sigmaScale[_be][_hl][1]*_sigmaScale[_be][_hl][1]*_e;
    //std::cout << "PhotonFix::sigmaEnergy 2 sigma = " << sigma << std::endl;
    sigma+=_sigmaScale[_be][_hl][2]*_sigmaScale[_be][_hl][2]*_e*_e;
    //std::cout << "PhotonFix::sigmaEnergy 3 sigma = " << sigma << std::endl;
  } else {
    sigma =_sigmaScale[_be][_hl][0]*_sigmaScale[_be][_hl][0]*cosh(_eta)*cosh(_eta);
    sigma+=_sigmaScale[_be][_hl][1]*_sigmaScale[_be][_hl][1]*_e;
    sigma+=_sigmaScale[_be][_hl][2]*_sigmaScale[_be][_hl][2]*_e*_e;
  }
  sigma=sqrt(sigma);
  
  double f(1.0);
  
  // General eta or zeta dependence
  if(_be==0) {
    f+=_sigmaAT[_be][_hl][0]*_eta*_eta;
    //std::cout << "PhotonFix::sigmaEnergy 4 f = " << f << std::endl;
    f+=expCorrection(_eta,_sigmaBT[_be][_hl]);
    //std::cout << "PhotonFix::sigmaEnergy 5 f = " << f << std::endl;
  } else {
    f+=_sigmaAT[_be][_hl][0]*xZ()*xZ();
    f+=_sigmaBT[_be][_hl][0]*yZ()*yZ();
  }
  
  // Eta or x crystal, submodule and module dependence
  f+=expCorrection(_aC,_sigmaAC[_be][_hl]);
  //std::cout << "PhotonFix::sigmaEnergy 6 f = " << f << std::endl;
  f+=expCorrection(_aS,_sigmaAS[_be][_hl]);
  //std::cout << "PhotonFix::sigmaEnergy 7 f = " << f << std::endl;
  f+=expCorrection(_aM,_sigmaAM[_be][_hl]);
  //std::cout << "PhotonFix::sigmaEnergy 8 f = " << f << std::endl;
  
  // Phi or y crystal, submodule and module dependence
  f+=expCorrection(_bC,_sigmaBC[_be][_hl]);
  //std::cout << "PhotonFix::sigmaEnergy 9 f = " << f << std::endl;
  f+=expCorrection(_bS,_sigmaBS[_be][_hl]);
  //std::cout << "PhotonFix::sigmaEnergy 10 f = " << f << std::endl;
  f+=expCorrection(_bM,_sigmaBM[_be][_hl]);
  //std::cout << "PhotonFix::sigmaEnergy 11 f = " << f << std::endl;
  
  // R9 dependence
  if(_hl==0) {
    f+=_sigmaR9[_be][_hl][1]*(_r9-_sigmaR9[_be][_hl][0])*(_r9-_sigmaR9[_be][_hl][0])
      +_sigmaR9[_be][_hl][2]*(_r9-_sigmaR9[_be][_hl][0])*(_r9-_sigmaR9[_be][_hl][0])*(_r9-_sigmaR9[_be][_hl][0]);
    //std::cout << "PhotonFix::sigmaEnergy 12 f = " << f << std::endl;
  } else {
    f+=_sigmaR9[_be][_hl][0]*_r9+_sigmaR9[_be][_hl][1]*_r9*_r9;
    //std::cout << "PhotonFix::sigmaEnergy 13 f = " << f << std::endl;
  }
  
  return sigma*f;
}

double PhotonFix::xC

(

)

const

Definition at line 289 of file PhotonFix.cc.

References _aC, and _be.

                           {
  assert(_be==1);
  return _aC;
}

double PhotonFix::xC

(

)

const

double PhotonFix::xM

(

)

const

Definition at line 299 of file PhotonFix.cc.

References _aM, and _be.

                           {
  assert(_be==1);
  return _aM;
}

double PhotonFix::xM

(

)

const

double PhotonFix::xS

(

)

const

Definition at line 294 of file PhotonFix.cc.

References _aS, and _be.

                           {
  assert(_be==1);
  return _aS;
}

double PhotonFix::xS

(

)

const

double PhotonFix::xZ

(

)

const

Definition at line 284 of file PhotonFix.cc.

References _be, _eta, _phi, asinh(), and funct::cos().

Referenced by fixedEnergy(), setup(), and sigmaEnergy().

                           {
  assert(_be==1);
  return asinh(cos(_phi)/sinh(_eta));
}

double PhotonFix::xZ

(

)

const

double PhotonFix::yC

(

)

const

double PhotonFix::yC

(

)

const

Definition at line 309 of file PhotonFix.cc.

References _bC, and _be.

                           {
  assert(_be==1);
  return _bC;
}

double PhotonFix::yM

(

)

const

double PhotonFix::yM

(

)

const

Definition at line 319 of file PhotonFix.cc.

References _be, and _bM.

                           {
  assert(_be==1);
  return _bM;
}

double PhotonFix::yS

(

)

const

Definition at line 314 of file PhotonFix.cc.

References _be, and _bS.

                           {
  assert(_be==1);
  return _bS;
}

double PhotonFix::yS

(

)

const

double PhotonFix::yZ

(

)

const

Definition at line 304 of file PhotonFix.cc.

References _be, _eta, _phi, asinh(), and funct::sin().

Referenced by fixedEnergy(), setup(), and sigmaEnergy().

                           {
  assert(_be==1);
  return asinh(sin(_phi)/sinh(_eta));
}

double PhotonFix::yZ

(

)

const

---

Member Data Documentation

double PhotonFix::_aC [private]

Definition at line 128 of file PhotonFix.h.

Referenced by etaC(), fixedEnergy(), print(), setup(), sigmaEnergy(), and xC().

double PhotonFix::_aM [private]

Definition at line 128 of file PhotonFix.h.

Referenced by etaM(), fixedEnergy(), print(), setup(), sigmaEnergy(), and xM().

double PhotonFix::_aS [private]

Definition at line 128 of file PhotonFix.h.

Referenced by etaS(), fixedEnergy(), print(), setup(), sigmaEnergy(), and xS().

static double PhotonFix::_barrelCGap [static, private]

Definition at line 162 of file PhotonFix.h.

Referenced by barrelCGap(), dumpGaps(), initialiseGeometry(), and setup().

static double PhotonFix::_barrelMGap [static, private]

Definition at line 164 of file PhotonFix.h.

Referenced by barrelMGap(), dumpGaps(), initialiseGeometry(), and setup().

static double PhotonFix::_barrelSGap [static, private]

Definition at line 163 of file PhotonFix.h.

Referenced by barrelSGap(), dumpGaps(), initialiseGeometry(), and setup().

double PhotonFix::_bC [private]

Definition at line 128 of file PhotonFix.h.

Referenced by fixedEnergy(), phiC(), print(), setup(), sigmaEnergy(), and yC().

unsigned PhotonFix::_be [private]

Definition at line 126 of file PhotonFix.h.

Referenced by etaC(), etaM(), etaS(), fixedEnergy(), phiC(), phiM(), phiS(), setup(), sigmaEnergy(), xC(), xM(), xS(), xZ(), yC(), yM(), yS(), and yZ().

double PhotonFix::_bM [private]

Definition at line 128 of file PhotonFix.h.

Referenced by fixedEnergy(), phiM(), print(), setup(), sigmaEnergy(), and yM().

double PhotonFix::_bS [private]

Definition at line 128 of file PhotonFix.h.

Referenced by fixedEnergy(), phiS(), print(), setup(), sigmaEnergy(), and yS().

double PhotonFix::_e [private]

Definition at line 127 of file PhotonFix.h.

Referenced by fixedEnergy(), print(), rawEnergy(), and sigmaEnergy().

static double PhotonFix::_endcapCGap [static, private]

Definition at line 168 of file PhotonFix.h.

Referenced by dumpGaps(), endcapCGap(), initialiseGeometry(), and setup().

static bool PhotonFix::_endcapCrystal [static, private]

Definition at line 167 of file PhotonFix.h.

Referenced by dumpGaps(), endcapCrystal(), and initialiseGeometry().

static double PhotonFix::_endcapMGap [static, private]

Definition at line 170 of file PhotonFix.h.

Referenced by dumpGaps(), endcapMGap(), initialiseGeometry(), and setup().

static double PhotonFix::_endcapSGap [static, private]

Definition at line 169 of file PhotonFix.h.

Referenced by dumpGaps(), endcapSGap(), initialiseGeometry(), and setup().

double PhotonFix::_eta [private]

Definition at line 127 of file PhotonFix.h.

Referenced by eta(), fixedEnergy(), print(), setup(), sigmaEnergy(), xZ(), and yZ().

unsigned PhotonFix::_hl [private]

Definition at line 126 of file PhotonFix.h.

Referenced by fixedEnergy(), setup(), and sigmaEnergy().

static bool PhotonFix::_initialised = false [static, private]

Definition at line 135 of file PhotonFix.h.

Referenced by initialise(), initialised(), initialiseParameters(), and setup().

static double PhotonFix::_meanAC [static, private]

Definition at line 140 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

static double PhotonFix::_meanAM [static, private]

Definition at line 142 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

static double PhotonFix::_meanAS [static, private]

Definition at line 141 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

double PhotonFix::_meanAT [static, private]

Definition at line 139 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

static double PhotonFix::_meanBC [static, private]

Definition at line 144 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

static double PhotonFix::_meanBM [static, private]

Definition at line 146 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

static double PhotonFix::_meanBS [static, private]

Definition at line 145 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

double PhotonFix::_meanBT [static, private]

Definition at line 143 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

static double PhotonFix::_meanR9 [static, private]

Definition at line 147 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

static double PhotonFix::_meanScale [static, private]

Definition at line 138 of file PhotonFix.h.

Referenced by dumpParameters(), fixedEnergy(), getParameters(), initialiseParameters(), printParameters(), and setParameters().

static const double PhotonFix::_onePi [static, private]

Definition at line 131 of file PhotonFix.h.

Referenced by aPhi(), and dPhi().

double PhotonFix::_phi [private]

Definition at line 127 of file PhotonFix.h.

Referenced by phi(), print(), setup(), xZ(), and yZ().

double PhotonFix::_r9 [private]

Definition at line 127 of file PhotonFix.h.

Referenced by fixedEnergy(), print(), r9(), setup(), and sigmaEnergy().

static double PhotonFix::_sigmaAC [static, private]

Definition at line 152 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaAM [static, private]

Definition at line 154 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaAS [static, private]

Definition at line 153 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaAT [static, private]

Definition at line 151 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaBC [static, private]

Definition at line 156 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaBM [static, private]

Definition at line 158 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaBS [static, private]

Definition at line 157 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaBT [static, private]

Definition at line 155 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaR9 [static, private]

Definition at line 159 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static double PhotonFix::_sigmaScale [static, private]

Definition at line 150 of file PhotonFix.h.

Referenced by dumpParameters(), getParameters(), initialiseParameters(), printParameters(), setParameters(), and sigmaEnergy().

static const double PhotonFix::_twoPi [static, private]

Definition at line 132 of file PhotonFix.h.

Referenced by dPhi().