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)
	~EvolutionECAL ()=default

Detailed Description

Definition at line 23 of file EvolutionECAL.h.

Constructor & Destructor Documentation

EvolutionECAL()

EvolutionECAL::EvolutionECAL ( )

inline

Definition at line 25 of file EvolutionECAL.h.

{};

~EvolutionECAL()

EvolutionECAL::~EvolutionECAL ( )

default

Member Function Documentation

AgingVPT()

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

Definition at line 353 of file EvolutionECAL.cc.

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

                                                                               {
  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;
}

ChargeVPTCathode()

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

Definition at line 338 of file EvolutionECAL.cc.

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

Referenced by AgingVPT().

                                                                                       {
  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;
}

DamageProfileEta()

double EvolutionECAL::DamageProfileEta

(

double

eta

)

Definition at line 169 of file EvolutionECAL.cc.

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

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

                                                 {
  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);
  }
}

DamageProfileEtaAPD()

double EvolutionECAL::DamageProfileEtaAPD

(

double

eta

)

Definition at line 179 of file EvolutionECAL.cc.

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

Referenced by NoiseFactorFE().

                                                    {
  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;
  }
}

DegradationMeanEM50GeV()

double EvolutionECAL::DegradationMeanEM50GeV

(

double

mu

)

Definition at line 246 of file EvolutionECAL.cc.

References simBeamSpotPI::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().

                                                      {
  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;
}

DegradationNonLinearityEM50GeV()

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

Definition at line 280 of file EvolutionECAL.cc.

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

                                                                          {
  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;
}

DoseLongitudinalProfile()

double EvolutionECAL::DoseLongitudinalProfile

(

double

z

)

Definition at line 200 of file EvolutionECAL.cc.

References simBeamSpotPI::alpha, HLT_2023v12_cff::beta, JetChargeProducer_cfi::exp, hltrates_dqm_sourceclient-live_cfg::offset, conifer::pow(), and z.

Referenced by EquilibriumFractionColorCentersEM().

                                                      {
  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;
  }
}

EquilibriumFractionColorCentersEM()

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

Definition at line 216 of file EvolutionECAL.cc.

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

Referenced by InducedAbsorptionEM().

                                                                                {
  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);
}

InducedAbsorptionEM()

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

Definition at line 227 of file EvolutionECAL.cc.

References EquilibriumFractionColorCentersEM(), and PVValHelper::eta.

Referenced by ChargeVPTCathode().

                                                                 {
  double mu_max = 2.0;
  double alpha1 = 3.41488e+00;

  std::unique_ptr<TF1> ftmp1 = std::make_unique<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;

  return muEM;
}

InducedAbsorptionHadronic()

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

Definition at line 193 of file EvolutionECAL.cc.

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

Referenced by ChargeVPTCathode().

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

LightCollectionEfficiency()

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

Definition at line 5 of file EvolutionECAL.cc.

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

                                                                   {
  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;
}

LightCollectionEfficiencyWeighted()

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

Definition at line 114 of file EvolutionECAL.cc.

References A, a0, testProducerWithPsetDescEmpty_cfi::a2, B, b0, b1, b2, b3, b4, correctionTermsCaloMet_cff::C, fftjetpileupestimator_calo_uncalib_cfi::c0, alignmentValidation::c1, d0, d1, StorageManager_cfg::e1, JetChargeProducer_cfi::exp, f, amptDefaultParameters_cff::mu, conifer::pow(), dttmaxenums::R, HLT_2023v12_cff::R0, and z.

                                                                               {
  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;
}

NoiseFactorFE()

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

Definition at line 362 of file EvolutionECAL.cc.

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

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

ResolutionConstantTermEM50GeV()

double EvolutionECAL::ResolutionConstantTermEM50GeV

(

double

mu

)

Definition at line 312 of file EvolutionECAL.cc.

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

                                                             {
  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);
}