EvolutionECAL Class Reference

#include <EvolutionECAL.h>

Public Member Functions
double	AgingVPT (double instLumi, double integralLumi, double eta)
double	ChargeVPTCathode (double instLumi, double eta, double integralLumi)
double	DamageProfileEta (double eta)
double	DamageProfileEtaAPD (double eta)
double	DegradationMeanEM50GeV (double mu)
double	DegradationNonLinearityEM50GeV (double mu, double ene)
double	DoseLongitudinalProfile (double z)
Double_t	EquilibriumFractionColorCentersEM (double *x, double *par)
	EvolutionECAL ()
double	InducedAbsorptionEM (double lumi, double eta)
double	InducedAbsorptionHadronic (double lumi, double eta)
double	LightCollectionEfficiency (double z, double mu)
double	LightCollectionEfficiencyWeighted (double z, double mu_ind)
double	NoiseFactorFE (double lumi, double eta)
double	ResolutionConstantTermEM50GeV (double mu)
virtual	~EvolutionECAL ()

Detailed Description

Definition at line 23 of file EvolutionECAL.h.

Constructor & Destructor Documentation

EvolutionECAL()

EvolutionECAL::EvolutionECAL

(

)

Definition at line 8 of file EvolutionECAL.cc.

{}

~EvolutionECAL()

EvolutionECAL::~EvolutionECAL ( )

virtual

Definition at line 5 of file EvolutionECAL.cc.

{}

Member Function Documentation

AgingVPT()

double EvolutionECAL::AgingVPT	(	double	instLumi,
		double	integralLumi,
		double	eta
	)

Definition at line 360 of file EvolutionECAL.cc.

                                                                               {
  if (fabs(eta) < 1.497)
    return 1.0;
  double Q = ChargeVPTCathode(instLumi, eta, integralLumi);
  double result = 0.772 + 0.228 * (3.94304e-01 * exp(-Q / 5.99232e-04) + (1 - 3.94304e-01) * exp(-Q / 1.58243e-02));
  return result;
}

References ChargeVPTCathode(), MillePedeFileConverter_cfg::e, PVValHelper::eta, JetChargeProducer_cfi::exp, muonGEMDigis_cfi::instLumi, class-composition::Q, and mps_fire::result.

ChargeVPTCathode()

double EvolutionECAL::ChargeVPTCathode	(	double	instLumi,
		double	eta,
		double	integralLumi
	)

Definition at line 345 of file EvolutionECAL.cc.

                                                                                       {
  double charge = 0.0;
  double tmpLumi = 0.0;
  double stepLumi = 1.0;
  double muEM = InducedAbsorptionEM(instLumi, eta);
  while (tmpLumi < integralLumi) {
    tmpLumi += stepLumi;
    double muHD = InducedAbsorptionHadronic(tmpLumi, eta);
    double SS0 = DegradationMeanEM50GeV(muEM + muHD);
    charge += SS0 * 0.26e-3 * DamageProfileEta(eta) * stepLumi;
  }
  return charge;
}

References ALCARECOTkAlJpsiMuMu_cff::charge, DamageProfileEta(), DegradationMeanEM50GeV(), PVValHelper::eta, InducedAbsorptionEM(), InducedAbsorptionHadronic(), and muonGEMDigis_cfi::instLumi.

Referenced by AgingVPT().

DamageProfileEta()

double EvolutionECAL::DamageProfileEta

(

double

eta

)

Definition at line 175 of file EvolutionECAL.cc.

                                                 {
  double x = fabs(eta);
  if (x < 1.497) {
    return exp(-4.11065 + 0.258478 * x);
  } else {
    return exp(-13.5112 + 7.913860 * x - 0.998649 * x * x);
  }
}

References PVValHelper::eta, JetChargeProducer_cfi::exp, and x.

Referenced by ChargeVPTCathode(), EquilibriumFractionColorCentersEM(), and InducedAbsorptionHadronic().

DamageProfileEtaAPD()

double EvolutionECAL::DamageProfileEtaAPD

(

double

eta

)

Definition at line 185 of file EvolutionECAL.cc.

                                                    {
  double x = fabs(eta);
  if (x < 1.497) {
    double et = x / 1.48 * 34.0;
    double etaprof = (9.54827 + et * 0.0379222 + et * et * (-0.00257047) + et * et * et * 0.00073546 +
                      et * et * et * et * (-1.49683e-05)) /
                     9.54827;
    return etaprof;
  } else {
    return 1.0;
  }
}

References EgHLTOffHistBins_cfi::et, PVValHelper::eta, and x.

Referenced by NoiseFactorFE().

DegradationMeanEM50GeV()

double EvolutionECAL::DegradationMeanEM50GeV

(

double

mu

)

Definition at line 253 of file EvolutionECAL.cc.

                                                      {
  double retval = 1.0;
  double x = mu;
  if (x < 1e-4)
    return retval;
  if (x >= 200.0)
    x = 200.0;  // parameterization is not valid for large mu
 
  double par[11] = {1.00000e+01,
                    -4.41441e-01,
                    7.08607e-02,
                    -3.75572e-01,
                    -3.60410e-01,
                    1.30130e-01,
                    -4.72350e-01,
                    3.36315e-01,
                    -1.19872e-01,
                    1.99574e-02,
                    -1.22910e-03};
 
  double alpha = par[0];
 
  double f1 = par[1] * x + par[2] * x * x;
  double u = log(x);
  double f2 = par[10];
  for (int i = 9; i >= 3; i--)
    f2 = par[i] + f2 * u;
 
  retval = f1 / (1.0 + exp(alpha * u)) + f2 / (1.0 + exp(-alpha * u));
  retval = exp(retval);
  return retval;
}

References zMuMuMuonUserData::alpha, MillePedeFileConverter_cfg::e, JetChargeProducer_cfi::exp, DeadROC_duringRun::f1, DeadROC_duringRun::f2, mps_fire::i, dqm-mbProfile::log, amptDefaultParameters_cff::mu, and x.

Referenced by ChargeVPTCathode().

DegradationNonLinearityEM50GeV()

double EvolutionECAL::DegradationNonLinearityEM50GeV	(	double	mu,
		double	ene
	)

Definition at line 287 of file EvolutionECAL.cc.

                                                                          {
  if (ene <= 1e-3)
    return 0.0;
 
  double x = mu;
  if (mu <= 0.06)
    x = 0.06;
  if (mu >= 150.0)
    x = 150.0;
 
  double par[9] = {5.17712e-03,
                   1.97597e-02,
                   3.36596e-02,
                   2.84505e-02,
                   1.38480e-02,
                   1.11498e-02,
                   7.73634e-03,
                   -1.30767e-03,
                   -2.20628e-03};
 
  double u = log10(x);
  double slope = par[8];
  for (int i = 7; i >= 0; i--)
    slope = par[i] + slope * u;
 
  double retval = 1.0 + slope * log10(ene / 50.0);
  if (retval <= 0.0)
    retval = 0.0;
  return retval;
}

References MillePedeFileConverter_cfg::e, mps_fire::i, amptDefaultParameters_cff::mu, slope, and x.

DoseLongitudinalProfile()

double EvolutionECAL::DoseLongitudinalProfile

(

double

z

)

Definition at line 206 of file EvolutionECAL.cc.

                                                      {
  double alpha = 4.72877e+00;
  double beta = 5.91296e-01;
  double amp1 = 6.24495e+02;
  double amp2 = 1.84367e-01;
  double offset = 2.00705e+01;
  if (z >= 0.0 && z <= 22.0) {
    double term1 = (amp1 / TMath::Gamma(alpha)) * pow((beta * z), (alpha - 1)) * exp(-beta * z);
    double term2 = amp2 * (z - 11.0) * (z - 11.0) + offset;
    return (term1 + term2) / 150.44;
  } else {
    return 0;
  }
}

References zMuMuMuonUserData::alpha, zMuMuMuonUserData::beta, JetChargeProducer_cfi::exp, hltrates_dqm_sourceclient-live_cfg::offset, funct::pow(), and z.

Referenced by EquilibriumFractionColorCentersEM().

EquilibriumFractionColorCentersEM()

Double_t EvolutionECAL::EquilibriumFractionColorCentersEM	(	double *	x,
		double *	par
	)

Definition at line 222 of file EvolutionECAL.cc.

                                                                                {
  double instantLumi = par[0];
  double eta = par[1];
  double rate = DoseLongitudinalProfile(x[0]) * 5.0 * DamageProfileEta(eta) * instantLumi / 1e+34;
  if (rate <= 0.0)
    rate = 0.0;
  double alpha = par[2];
  return rate / (alpha + rate);
}

References zMuMuMuonUserData::alpha, DamageProfileEta(), DoseLongitudinalProfile(), MillePedeFileConverter_cfg::e, PVValHelper::eta, RPCpg::rate(), and x.

Referenced by InducedAbsorptionEM().

InducedAbsorptionEM()

double EvolutionECAL::InducedAbsorptionEM	(	double	lumi,
		double	eta
	)

Definition at line 233 of file EvolutionECAL.cc.

                                                                 {
  double mu_max = 2.0;
  double alpha1 = 3.41488e+00;
 
  TF1 *ftmp1 = new TF1("ftmp1",
                       this,
                       &EvolutionECAL::EquilibriumFractionColorCentersEM,
                       0.0,
                       22.0,
                       3,
                       "EvolutionECAL",
                       "EquilibriumFractionColorCentersEM");
  ftmp1->SetParameters(lumi, eta, alpha1);
  double muEM = mu_max * ftmp1->Integral(0.0, 22.0) / 22.0;
 
  delete ftmp1;
  return muEM;
}

References EquilibriumFractionColorCentersEM(), PVValHelper::eta, and tools::TF1.

Referenced by ChargeVPTCathode().

InducedAbsorptionHadronic()

double EvolutionECAL::InducedAbsorptionHadronic	(	double	lumi,
		double	eta
	)

Definition at line 199 of file EvolutionECAL.cc.

                                                                       {
  double fluence = DamageProfileEta(eta) * 2.7e+13 / 500.0 * lumi;
  double mu = 2.08E-13 * pow(fluence, 1.0049);
  return mu;
}

References DamageProfileEta(), PVValHelper::eta, BXlumiParameters_cfi::lumi, amptDefaultParameters_cff::mu, and funct::pow().

Referenced by ChargeVPTCathode().

LightCollectionEfficiency()

double EvolutionECAL::LightCollectionEfficiency	(	double	z,
		double	mu
	)

Definition at line 11 of file EvolutionECAL.cc.

                                                                   {
  double f = 0;
  if (z <= 0)
    return f;
  if (z >= 0.22)
    return f;
 
  double e0 = 6.91563e-02;
  double e1 = 1.64406e+00;
  double e2 = 6.42509e-01;
  double E = e0 / (1 + exp(e1 * (log10(mu) - e2)));
 
  double d0 = 3.85334e-01;
  double d1 = -1.04647e-02;
  double D = d0 * exp(d1 * mu);
 
  double c0 = 3.77629e-01;
  double c1 = -3.23755e-01;
  double c2 = 1.50247e+00;
  double c3 = 3.03278e-01;
  double C = -1 + c0 * exp(c1 * mu) * (1 + c2 * exp(c3 * mu));
 
  double b0 = -3.33575e-01;
  double b1 = 4.44856e-01;
  double b2 = 1.91766e+00;
  double b3 = 2.69423e+00;
  double b4 = 1.06905e+00;
  double B =
      (1 / mu) * (b0 + b1 * log10(mu) + b2 * pow(log10(mu), 2) + b3 * pow(log10(mu), 3) + b4 * pow(log10(mu), 4));
 
  double a0 = 7.18248e-02;
  double a1 = 1.89016e+00;
  double a2 = 2.15651e-02;
  double a3 = 2.30786e-02;
  double A = exp(B * mu * 0.015) * (a0 / (exp(a1 * (log10(mu) + a2)) + 1) + a3);
 
  double R = 0.01 * D * (4 / (0.222 + E) / (0.222 + E) - 1 / ((0.22 - z) * (z + E)));
  f = A * exp(-B * mu * (0.22 - z)) * (1 + C * exp(R));
 
  return f;
}

References a0, testProducerWithPsetDescEmpty_cfi::a2, TtFullHadDaughter::B, b0, b1, b2, testProducerWithPsetDescEmpty_cfi::b3, gen::C, fftjetpileupestimator_calo_uncalib_cfi::c0, alignmentValidation::c1, d0, d1, StorageManager_cfg::e1, JetChargeProducer_cfi::exp, f, amptDefaultParameters_cff::mu, funct::pow(), dttmaxenums::R, and z.

LightCollectionEfficiencyWeighted()

double EvolutionECAL::LightCollectionEfficiencyWeighted	(	double	z,
		double	mu_ind
	)

Definition at line 120 of file EvolutionECAL.cc.

                                                                               {
  if (z <= 0)
    return 0;
  if (z >= 1)
    return 0;
  if (mu_ind < 0)
    return 1;
 
  double mu = mu_ind + 0.7;
  double lmu = log10(mu);
 
  double e0 = 6.91563e-02;
  double e1 = 1.64406e+00;
  double e2 = 6.42509e-01;
  double E = e0 / (1 + exp(e1 * (lmu - e2)));
 
  double d0 = 3.85334e-01;
  double d1 = -1.04647e-02;
  double D = d0 * exp(d1 * mu);
 
  double c0 = 3.77629e-01;
  double c1 = -3.23755e-01;
  double c2 = 1.50247e+00;
  double c3 = 3.03278e-01;
  double C = -1 + c0 * exp(c1 * mu) * (1 + c2 * exp(c3 * mu));
 
  double b0 = -3.33575e-01;
  double b1 = 4.44856e-01;
  double b2 = 1.91766e+00;
  double b3 = 2.69423e+00;
  double b4 = 1.06905e+00;
  double B = (1 / mu) * (b0 + b1 * lmu + b2 * pow(lmu, 2) + b3 * pow(lmu, 3) + b4 * pow(lmu, 4));
 
  double a0 = 7.18248e-02;
  double a1 = 1.89016e+00;
  double a2 = 2.15651e-02;
  double a3 = 2.30786e-02;
  double A = exp(B * mu * 0.015) * (a0 / (exp(a1 * (lmu + a2)) + 1) + a3);
 
  double R = 0.01 * D * (4 / (0.222 + E) / (0.222 + E) - 1 / ((0.22 * 0.22) * (1. - z) * (z + E / 0.22)));
 
  // for undamaged crystal, mu0 = 0.7
  double A0 = 0.0631452;
  double B0 = -0.52267;
  double C0 = -0.139646;
  double D0 = 0.382522;
  double E0 = 0.054473;
  double R0 = 0.01 * D0 * (4 / (0.222 + E0) / (0.222 + E0) - 1 / ((0.22 * 0.22) * (1. - z) * (z + E0 / 0.22)));
 
  double f = A / A0 * exp(-(B * mu - B0 * 0.7) * 0.22 * (1. - z)) * (1 + C * exp(R)) / (1 + C0 * exp(R0));
 
  return f;
}

References a0, testProducerWithPsetDescEmpty_cfi::a2, TtFullHadDaughter::B, b0, b1, b2, testProducerWithPsetDescEmpty_cfi::b3, gen::C, fftjetpileupestimator_calo_uncalib_cfi::c0, alignmentValidation::c1, d0, d1, StorageManager_cfg::e1, JetChargeProducer_cfi::exp, f, amptDefaultParameters_cff::mu, funct::pow(), dttmaxenums::R, HLT_FULL_cff::R0, and z.

NoiseFactorFE()

double EvolutionECAL::NoiseFactorFE	(	double	lumi,
		double	eta
	)

Definition at line 369 of file EvolutionECAL.cc.

                                                           {
  double x = fabs(eta);
  if (x < 1.497) {
    return sqrt(1.0 + 0.495 * 0.03512 * lumi * DamageProfileEtaAPD(eta));
  } else {
    return 1.0;
  }
}

References DamageProfileEtaAPD(), PVValHelper::eta, mathSSE::sqrt(), and x.

ResolutionConstantTermEM50GeV()

double EvolutionECAL::ResolutionConstantTermEM50GeV

(

double

mu

)

Definition at line 319 of file EvolutionECAL.cc.

                                                             {
  double x = mu;
  if (mu <= 0.01)
    x = 0.01;
  if (mu >= 200.0)
    x = 200.0;
 
  double par[10] = {-6.21503e+00,
                    1.59759e+00,
                    -4.75221e-02,
                    -3.90299e-02,
                    3.97269e-03,
                    2.29574e-03,
                    -1.05280e-04,
                    -9.60963e-05,
                    -1.29594e-06,
                    1.70850e-06};
 
  double u = log(x);
  double f = par[9];
  for (int i = 8; i >= 0; i--)
    f = par[i] + f * u;
  return exp(f);
}

References JetChargeProducer_cfi::exp, f, mps_fire::i, dqm-mbProfile::log, amptDefaultParameters_cff::mu, and x.