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

(

)

Definition at line 11 of file EvolutionECAL.cc.

{
}

EvolutionECAL::~EvolutionECAL ( )

virtual

Definition at line 6 of file EvolutionECAL.cc.

{
}

Member Function Documentation

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

Definition at line 364 of file EvolutionECAL.cc.

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

Referenced by EnergyResolutionVsLumi::calcampDropPhotoDetector(), and EnergyResolutionVsLumi::CalculateDegradation().

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

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

Definition at line 343 of file EvolutionECAL.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, DamageProfileEta(), DegradationMeanEM50GeV(), InducedAbsorptionEM(), and InducedAbsorptionHadronic().

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

double EvolutionECAL::DamageProfileEta

(

double

eta

)

Definition at line 183 of file EvolutionECAL.cc.

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

double EvolutionECAL::DamageProfileEtaAPD

(

double

eta

)

Definition at line 195 of file EvolutionECAL.cc.

References stringResolutionProvider_cfi::et, 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;
  }
}

double EvolutionECAL::DegradationMeanEM50GeV

(

double

mu

)

Definition at line 271 of file EvolutionECAL.cc.

References alpha, MillePedeFileConverter_cfg::e, JetChargeProducer_cfi::exp, connectstrParser::f1, connectstrParser::f2, mps_fire::i, cmsBatch::log, RPCpg::mu, and x.

Referenced by EnergyResolutionVsLumi::calcampDropTransparency(), EnergyResolutionVsLumi::CalculateDegradation(), and 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;
  

}

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

Definition at line 301 of file EvolutionECAL.cc.

References MillePedeFileConverter_cfg::e, mps_fire::i, RPCpg::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;

}

double EvolutionECAL::DoseLongitudinalProfile

(

double

z

)

Definition at line 222 of file EvolutionECAL.cc.

References alpha, beta, JetChargeProducer_cfi::exp, Gamma, PFRecoTauDiscriminationByIsolation_cfi::offset, and funct::pow().

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

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

Definition at line 241 of file EvolutionECAL.cc.

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

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

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

Definition at line 255 of file EvolutionECAL.cc.

References EquilibriumFractionColorCentersEM().

Referenced by EnergyResolutionVsLumi::calcmuEM(), EnergyResolutionVsLumi::calcmuTot(), EnergyResolutionVsLumi::CalculateDegradation(), and ChargeVPTCathode().

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

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

Definition at line 213 of file EvolutionECAL.cc.

References DamageProfileEta(), csvLumiCalc::lumi, RPCpg::mu, and funct::pow().

Referenced by EnergyResolutionVsLumi::calcmuHD(), EnergyResolutionVsLumi::calcmuTot(), EnergyResolutionVsLumi::CalculateDegradation(), and ChargeVPTCathode().

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

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

Definition at line 16 of file EvolutionECAL.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, patCaloMETCorrections_cff::C, fftjetpileupestimator_calo_uncalib_cfi::c0, alignmentValidation::c1, allConversions_cfi::d0, reco::e1, reco::e2, JetChargeProducer_cfi::exp, f, RPCpg::mu, funct::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;
}

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

Definition at line 126 of file EvolutionECAL.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, patCaloMETCorrections_cff::C, fftjetpileupestimator_calo_uncalib_cfi::c0, alignmentValidation::c1, allConversions_cfi::d0, reco::e1, reco::e2, JetChargeProducer_cfi::exp, f, RPCpg::mu, funct::pow(), dttmaxenums::R, pfBoostedDoubleSVAK8TagInfos_cfi::R0, and z.

Referenced by EnergyResolutionVsLumi::calcLightCollectionEfficiencyWeighted(), and EnergyResolutionVsLumi::calcLightCollectionEfficiencyWeighted2().

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

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

Definition at line 375 of file EvolutionECAL.cc.

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

Referenced by EnergyResolutionVsLumi::calcnoiseADC(), EnergyResolutionVsLumi::calcnoiseIncreaseADC(), and EnergyResolutionVsLumi::CalculateDegradation().

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

double EvolutionECAL::ResolutionConstantTermEM50GeV

(

double

mu

)

Definition at line 324 of file EvolutionECAL.cc.

References JetChargeProducer_cfi::exp, f, mps_fire::i, cmsBatch::log, RPCpg::mu, and x.

Referenced by EnergyResolutionVsLumi::calcresolutitonConstantTerm(), and EnergyResolutionVsLumi::CalculateDegradation().

{

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

}