LArHVScaleCorrTool Class Reference

#include <LArHVScaleCorrTool.h>

Collaboration diagram for LArHVScaleCorrTool:

Classes
struct	HV_t
struct	HVfix_t

Public Types
typedef std::vector< HV_t >	voltageCell_t

Public Member Functions
	LArHVScaleCorrTool (const CaloCell_ID *caloCellID, MsgStream &msg, const std::vector< std::string > &fixHVStrings)
	LArHVScaleCorrTool ()=delete
virtual	~LArHVScaleCorrTool ()=default
float	getHVScale (const CaloDetDescrElement *calodde, const voltageCell_t &hvlist, MsgStream &msg) const

Private Member Functions
float	Scale_barrel (const float hv) const
float	Scale_FCAL1 (const float hv) const
float	champ_e (float hv, float d) const
float	vdrift (float e, float tempe) const
float	InvCharge (float e) const
float	Respo (float e, float e_nominal, float tempe) const
float	t_drift (float e, float e_nominal, float d, float tempe) const
float	EMEC_nominal (const float eta_r) const
float	EMEC_gap (const float eta_r, float Zeta) const
void	buildFixHVList (const std::vector< std::string > &fixHVStrings, MsgStream &msg)

Private Attributes
const LArEM_ID *	m_larem_id
const LArHEC_ID *	m_larhec_id
const LArFCAL_ID *	m_larfcal_id
const float	m_T0
std::vector< HVfix_t >	m_HVfix

Detailed Description

Definition at line 20 of file LArHVScaleCorrTool.h.

Member Typedef Documentation

voltageCell_t

typedef std::vector<HV_t> LArHVScaleCorrTool::voltageCell_t

Definition at line 30 of file LArHVScaleCorrTool.h.

Constructor & Destructor Documentation

LArHVScaleCorrTool() [1/2]

LArHVScaleCorrTool::LArHVScaleCorrTool	(	const CaloCell_ID *	caloCellID,
		MsgStream &	msg,
		const std::vector< std::string > &	fixHVStrings )

Definition at line 15 of file LArHVScaleCorrTool.cxx.

                                                                 :
  m_larem_id(calocell_id->em_idHelper()),
  m_larhec_id(calocell_id->hec_idHelper()),
  m_larfcal_id(calocell_id->fcal_idHelper()),
  m_T0(90.371)
{
  buildFixHVList(fixHVStrings, msg);
}

LArHVScaleCorrTool() [2/2]

LArHVScaleCorrTool::LArHVScaleCorrTool ( )

delete

~LArHVScaleCorrTool()

virtual LArHVScaleCorrTool::~LArHVScaleCorrTool ( )

virtualdefault

Member Function Documentation

buildFixHVList()

void LArHVScaleCorrTool::buildFixHVList ( const std::vector< std::string > & fixHVStrings, MsgStream & msg )

private

Definition at line 317 of file LArHVScaleCorrTool.cxx.

                                                              {
 
  m_HVfix.clear();
  for (const std::string& theString : fixHVStrings) {
    std::stringstream is;
    is << theString << std::endl;
 
    unsigned int iDetector, ilayer_min, ilayer_max;
    float eta_min, eta_max, phi_min, phi_max, corr;
    is >> iDetector >> ilayer_min >> ilayer_max >> eta_min >> eta_max >>
        phi_min >> phi_max >> corr;
 
    m_HVfix.push_back({.subdet = iDetector,
                       .layer_min = ilayer_min,
                       .layer_max = ilayer_max,
                       .eta_min = eta_min,
                       .eta_max = eta_max,
                       .phi_min = phi_min,
                       .phi_max = phi_max,
                       .corr = corr});
  }
 
  msg << MSG::INFO
      << "  Number of regions with overwritten HV corrections from jobOptions "
      << m_HVfix.size() << endmsg;
 
  return;
}

champ_e()

float LArHVScaleCorrTool::champ_e ( float hv, float d ) const

private

Definition at line 185 of file LArHVScaleCorrTool.cxx.

{
  float e1;
  if (hv < -3000.){ 
    return -1000.;
  }
  else
    e1 = std::abs(hv)/(d*1e3);
  if ( e1 < 0.01 ) e1 = 0.01;
  return e1;
}

EMEC_gap()

float LArHVScaleCorrTool::EMEC_gap ( const float eta_r, float Zeta ) const

private

Definition at line 249 of file LArHVScaleCorrTool.cxx.

{
  float EMECg;
  if (aeta<=2.5 ) EMECg = 0.9 +1.9*(  ( Zeta - 40. )/(10*std::sinh(aeta)) - 60.)/140.;
  else EMECg = 1.8 + 1.3*(  ( Zeta - 40. )/(10*std::sinh(aeta)) - 30.)/40.;
  return EMECg;
}

EMEC_nominal()

float LArHVScaleCorrTool::EMEC_nominal ( const float eta_r ) const

private

Definition at line 236 of file LArHVScaleCorrTool.cxx.

{
  if ( aeta<1.5 ) return 2500.;
  else if (aeta<1.6) return 2300.;
  else if (aeta<1.8 ) return 2100.;
  else if (aeta < 2.0 ) return 1700.;
  else if (aeta < 2.1 ) return 1500.;
  else if (aeta < 2.3 ) return 1250.;
  else if (aeta < 2.5 ) return 1000.;
  else if (aeta < 2.8 ) return 2300.;
  else return 1800.;
}

getHVScale()

float LArHVScaleCorrTool::getHVScale	(	const CaloDetDescrElement *	calodde,
		const voltageCell_t &	hvlist,
		MsgStream &	msg ) const

Definition at line 26 of file LArHVScaleCorrTool.cxx.

                                                           {
 
  double d=.2;  // nominal gap  in cm
  double nominal=2000.;  // nominal HV  in Volts
  double T=88.;  // temperature  in Kelvin
  
  bool isbarrelEM=false;
  bool isFCAL1=false;
  
  unsigned int subdet=99;
  unsigned int layer=99;
 
  const Identifier offid=calodde->identify();
  const float eta_raw = calodde->eta_raw();
  const float phi_raw = calodde->phi_raw();
 
 
  // EM calorimeter
  if (m_larem_id->is_lar_em(offid)) {
    // barrel
    if (abs(m_larem_id->barrel_ec(offid))==1) {
      subdet=0;
      layer = m_larem_id->sampling(offid);
      // EMB
      if (m_larem_id->sampling(offid)>0) {
    d=0.209;
    nominal=2000.;
    T=88.37;
    isbarrelEM=true;
      }
      // EMBPS
      else {
    nominal = 2000.;
    T = 88.37;
    const int ieta=m_larem_id->eta(offid);
    if (ieta>=0 && ieta<16)        d = 0.196; //cm
    else if (ieta>=16 && ieta<32)  d = 0.193; //cm
    else if (ieta>=32 && ieta<48)  d = 0.2; //cm
    else  d = 0.190; //cm
      }
    }
    // endcap
    else {
      subdet=1;
      layer = m_larem_id->sampling(offid);
      // End-Cap pre-sampler
      if (abs(m_larem_id->barrel_ec(offid))==2 && m_larem_id->sampling(offid)==0){
    T = 88.65; // temperature remainder Calorimeter
    nominal = -2000.;
    d = 0.2; 
      }
      //EMEC      
      else {
    T = 88.65;
    float aeta_raw = std::abs(eta_raw);
    double Zsamp;
    if ( m_larem_id->sampling(offid)==1 && m_larem_id->region(offid)>=0 )  Zsamp = 3760.; //mm
    else if ( m_larem_id->sampling(offid)==2 && m_larem_id->region(offid)<=1 )  Zsamp = 3880.; //mm
    else if(aeta_raw< 2.) Zsamp=4200.-40.*(aeta_raw-1.5);
    else Zsamp = 4180. - 160.*(aeta_raw-2.);
    nominal = EMEC_nominal(aeta_raw);
    d = EMEC_gap(aeta_raw, Zsamp)*0.1; //cm
      }//end EMEC
    }//end endcap
  }//end is-em
  // Forward Calorimeter
  else if ( m_larfcal_id->is_lar_fcal(offid)) {
    subdet=3;
    layer = m_larfcal_id->module(offid);
    T = 88.65;
    if (m_larfcal_id->module(offid)==1){ 
      d =0.0269; //cm
      nominal = 250.;
      isFCAL1=true;
    }
    else  if (m_larfcal_id->module(offid)==2){
      d =0.0376;//cm
      nominal = 375.;
    }
    else { 
      d =0.0508;//cm
      nominal = 500.;
    }
  }
  // Hadronic Calorimeter
  else if ( m_larhec_id->is_lar_hec(offid)) {
    subdet=2;
    layer = m_larhec_id->sampling(offid);
    T = 88.65;
    nominal = 1800.;
    d = 0.196;//cm
  }
  T = T - m_T0;
 
  const double E_nominal = champ_e(nominal,d);
  
  // dump values
  bool notfound=false;
  for (unsigned int i=0;i<hvlist.size();i++) {
    if (hvlist[i].hv<-10000) notfound=true;
  }
 
  if (notfound) {
    msg << MSG::WARNING << " At least one HV value not found in database for cell " << m_larem_id->show_to_string(offid) << endmsg;
  }
 
  double mycorr=0.;
  double mynorm=0.;
  for (unsigned int i=0;i<hvlist.size();i++) {
    double E = champ_e(hvlist[i].hv,d);
    
    // dont correct if E is very close to E nominal to avoid small glitches
    if (std::abs(E_nominal)>1e-3) {
      const double deviation = std::abs((E-E_nominal)/E_nominal);
      if (deviation<1e-3) E = E_nominal;
    }
 
    // barrel accordion
    if (isbarrelEM) {
      const double corr = this->Scale_barrel(hvlist[i].hv)*hvlist[i].weight;
      mycorr += corr;
    }
    //FCAL module 1
    else if (isFCAL1) {
      const double corr = this->Scale_FCAL1(hvlist[i].hv) * hvlist[i].weight;
      mycorr+=corr;
    }
    // other subdetectors
    else {
      const double corr = Respo(E,E_nominal,T)*hvlist[i].weight;
      mycorr+= corr;
    }
    mynorm += hvlist[i].weight;
  }//end loop over hvlist
  
 
  if (mycorr>1e-2) mycorr = mynorm/mycorr;
  else mycorr=1.;
  
  // Add protection against correction factor > 10
  if (mycorr>10.) {
    msg << MSG::DEBUG << "Correction factor > 10, ignore it... for cell " <<  m_larem_id->show_to_string(offid) << " " << mycorr << endmsg;;
    mycorr=1.;
  }   
 
  //Overwrite with hardcoded jobOption numbers if requested
  for (unsigned int ii=0;ii<m_HVfix.size();ii++) {
    if (subdet == m_HVfix[ii].subdet && layer >= m_HVfix[ii].layer_min && layer <= m_HVfix[ii].layer_max &&
    eta_raw >= m_HVfix[ii].eta_min && eta_raw < m_HVfix[ii].eta_max &&
    phi_raw >= m_HVfix[ii].phi_min && phi_raw < m_HVfix[ii].phi_max ) {
      msg << MSG::DEBUG << "Overwrite correction for cell " <<  m_larem_id->show_to_string(offid) << " " << m_HVfix[ii].corr << endmsg;;
      mycorr = m_HVfix[ii].corr;
    }
  }
    
  return mycorr;
}

InvCharge()

float LArHVScaleCorrTool::InvCharge ( float e ) const

private

Definition at line 206 of file LArHVScaleCorrTool.cxx.

{ 
  float q = 1.;
  if ( e > 2.) q=(1.+0.36/e);
  else if (e>1e-6) q =(1.04+0.25/e);
  return q;
}

Respo()

float LArHVScaleCorrTool::Respo ( float e, float e_nominal, float tempe ) const

private

Definition at line 215 of file LArHVScaleCorrTool.cxx.

{ 
  if (e < -999.) return 1.;
  if (e < 0.01) return 0;
  if ( e > e_nominal ) return 1;
  float den = InvCharge(e)*vdrift(e_nominal,tempe);
  float resp = den==0 ? 0 : (InvCharge(e_nominal)*vdrift(e,tempe))/den;
  return resp;
}

Scale_barrel()

float LArHVScaleCorrTool::Scale_barrel ( const float hv ) const

private

Definition at line 269 of file LArHVScaleCorrTool.cxx.

{
   const float hvref[18]={1.,50.,100.,150.,200.,300.,400.,500.,600.,700.,800.,900.,1000.,1200.,1400.,1600.,1800.,2000.};
   const float hvmax = 1998.;
   const float facteur[18]={0.,0.1209,0.2135,0.2829,0.3390,0.4270,0.4961,0.5556,0.6065,0.6527,0.6906,
                0.7290,0.7626,0.8224,0.8754,0.9190,0.9606,1.};
 
// 0 HV, returns 0 response
   if (hv<-999.) {
      return 0;
   }
   else if (hv<hvref[0]) {
      float resp=facteur[0];
      return resp;
   }
// 2000 HV, returns response=1
   if (hv>hvmax) {
      float resp=facteur[17];
      return resp;
   }
 
// intermediate HV, find values between which to extrapolate
   int i1,i2;
   i1=1;
   for (int i=0;i<18;i++) {
    if (hv<hvref[i]) {
       i1=i-1;
       break;
    }
   }
   i2=i1+1;
   float resp=0.;
 
// if lowHV>=50 V, logarithmic extrapolation
   if (i1>0) {
     const float b=(std::log(facteur[i2])-std::log(facteur[i1]))/(std::log(hvref[i2])-std::log(hvref[i1]));
     const float a=std::log(facteur[i2])-b*std::log(hvref[i2]);
     resp = std::exp(a+b*std::log(hv));
   } 
// if between 0 and 50 V, scales linearly
   else {
     resp=facteur[0]*(hv/hvref[0]);
   }
   //std::cout << " hv,i1,i2,resp " << hv << " " << i1 << " " << i2 << " " << resp << std::endl;
   return resp;
}

Scale_FCAL1()

float LArHVScaleCorrTool::Scale_FCAL1 ( const float hv ) const

private

Definition at line 257 of file LArHVScaleCorrTool.cxx.

{
  if (hv<5.0) return 0;
  const float R0=-2.612;
  const float alpha=2.336;
  const float beta=0.079;
  const float scale=R0+alpha*std::pow(hv,beta);
  return scale;
}

t_drift()

float LArHVScaleCorrTool::t_drift ( float e, float e_nominal, float d, float tempe ) const

private

Definition at line 225 of file LArHVScaleCorrTool.cxx.

{
  if ( e < -999.) {
    float den = vdrift(e_nominal, tempe);
    return den==0 ? 0 : (d*1e4)/den;
  }
  if (e > e_nominal ) e = e_nominal;
  float den = vdrift(e, tempe);
  return den==0 ? 0 : (d*1e4)/den; // ns
}

vdrift()

float LArHVScaleCorrTool::vdrift ( float e, float tempe ) const

private

Definition at line 197 of file LArHVScaleCorrTool.cxx.

{       
  const float T = tempe;                
  static const float P[6] = {-0.01481,-0.0075,0.141,12.4,1.627,0.317};
  if ( e < -999.) return 0.;
  float vd = (P[0]*T+1)*( P[2]*e*std::log(1+ (P[3]/e)) + P[4]*std::pow(e,P[5])) + P[1]*T; // vdrift formula walcowialk mm/micro_s
  return vd;
}

Member Data Documentation

m_HVfix

std::vector<HVfix_t> LArHVScaleCorrTool::m_HVfix

private

Definition at line 94 of file LArHVScaleCorrTool.h.

m_larem_id

const LArEM_ID* LArHVScaleCorrTool::m_larem_id

private

Definition at line 45 of file LArHVScaleCorrTool.h.

m_larfcal_id

const LArFCAL_ID* LArHVScaleCorrTool::m_larfcal_id

private

Definition at line 47 of file LArHVScaleCorrTool.h.

m_larhec_id

const LArHEC_ID* LArHVScaleCorrTool::m_larhec_id

private

Definition at line 46 of file LArHVScaleCorrTool.h.

m_T0

const float LArHVScaleCorrTool::m_T0

private

Definition at line 48 of file LArHVScaleCorrTool.h.

---

The documentation for this class was generated from the following files:

• LArHVScaleCorrTool.h
• LArHVScaleCorrTool.cxx