CMS 3D CMS Logo
Public Member Functions

EvolutionECAL Class Reference

#include <EvolutionECAL.h>

List of all members.

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(), alignCSCRings::e, create_public_lumi_plots::exp, and query::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 DeDxDiscriminatorTools::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 create_public_lumi_plots::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 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, alignCSCRings::e, create_public_lumi_plots::exp, python::connectstrParser::f1, python::connectstrParser::f2, i, create_public_lumi_plots::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 alignCSCRings::e, 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, create_public_lumi_plots::exp, Gamma, evf::evtn::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(), alignCSCRings::e, 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(), fjr2json::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 funct::A, funct::C, alignmentValidation::c1, funct::D, create_public_lumi_plots::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 funct::A, funct::C, alignmentValidation::c1, funct::D, create_public_lumi_plots::exp, f, RPCpg::mu, funct::pow(), dttmaxenums::R, 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 create_public_lumi_plots::exp, f, i, create_public_lumi_plots::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);

}
Generated for CMSSW Reference Manual by doxygen 1.7.3