HGCalSciNoiseMap Class Reference

derives from HGCalRadiation map to parse fluence parameters, provides Sci-specific functions More...

#include <HGCalSciNoiseMap.h>

Inheritance diagram for HGCalSciNoiseMap:

Public Member Functions
radiiVec	computeRadius (const HGCScintillatorDetId &)
	HGCalSciNoiseMap ()
std::pair< double, double >	scaleByDose (const HGCScintillatorDetId &, const radiiVec &)
double	scaleBySipmArea (const HGCScintillatorDetId &, const double &)
double	scaleByTileArea (const HGCScintillatorDetId &, const radiiVec &)
	returns the signal scaling and the noise More...
void	setSipmMap (const std::string &)
	~HGCalSciNoiseMap ()
Public Member Functions inherited from HGCalRadiationMap
const unsigned int &	algo ()
const HGCalDDDConstants *	ddd ()
const HGCalGeometry *	geom ()
const doseParametersMap &	getDoseMap ()
double	getDoseValue (const int, const int, const radiiVec &, bool logVal=false)
double	getFluenceValue (const int, const int, const radiiVec &, bool logVal=false)
	HGCalRadiationMap ()
void	setDoseMap (const std::string &, const unsigned int)
void	setGeometry (const CaloSubdetectorGeometry *)
const HGCalTopology *	topo ()
	~HGCalRadiationMap ()

Private Member Functions
std::unordered_map< int, float >	readSipmPars (const std::string &)

Private Attributes
const double	refEdge_
std::unordered_map< int, float >	sipmMap_

Additional Inherited Members
Public Types inherited from HGCalRadiationMap
typedef std::map< std::pair< int, int >, DoseParameters >	doseParametersMap

Detailed Description

derives from HGCalRadiation map to parse fluence parameters, provides Sci-specific functions

Definition at line 13 of file HGCalSciNoiseMap.h.

Constructor & Destructor Documentation

HGCalSciNoiseMap()

HGCalSciNoiseMap::HGCalSciNoiseMap

(

)

Definition at line 6 of file HGCalSciNoiseMap.cc.

: refEdge_(3.) {}

~HGCalSciNoiseMap()

HGCalSciNoiseMap::~HGCalSciNoiseMap ( )

inline

Definition at line 16 of file HGCalSciNoiseMap.h.

{};

Member Function Documentation

computeRadius()

radiiVec HGCalSciNoiseMap::computeRadius

(

const HGCScintillatorDetId &

cellId

)

Definition at line 87 of file HGCalSciNoiseMap.cc.

                                                                           {
  GlobalPoint global = geom()->getPosition(cellId);
 
  double radius2 = std::pow(global.x(), 2) + std::pow(global.y(), 2);  //in cm
  double radius4 = std::pow(radius2, 2);
  double radius = sqrt(radius2);
  double radius3 = radius2 * radius;
 
  double radius_m100 = radius - 100;
  double radius_m100_2 = std::pow(radius_m100, 2);
  double radius_m100_3 = radius_m100_2 * radius_m100;
  double radius_m100_4 = std::pow(radius_m100_2, 2);
 
  radiiVec radii{{radius, radius2, radius3, radius4, radius_m100, radius_m100_2, radius_m100_3, radius_m100_4}};
  return radii;
}

References HGCalRadiationMap::geom(), HGCalGeometry::getPosition(), funct::pow(), CosmicsPD_Skims::radius, mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by HGCHEbackDigitizer::runRealisticDigitizer().

readSipmPars()

std::unordered_map< int, float > HGCalSciNoiseMap::readSipmPars ( const std::string &  fullpath )

private

Definition at line 12 of file HGCalSciNoiseMap.cc.

                                                                                     {
  std::unordered_map<int, float> result;
  //no file means default sipm size
  if (fullpath.empty())
    return result;
 
  edm::FileInPath fp(fullpath);
  std::ifstream infile(fp.fullPath());
  if (!infile.is_open()) {
    throw cms::Exception("FileNotFound") << "Unable to open '" << fullpath << "'" << std::endl;
  }
  std::string line;
  while (getline(infile, line)) {
    int layer;
    float boundary;
 
    //space-separated
    std::stringstream linestream(line);
    linestream >> layer >> boundary;
 
    result[layer] = boundary;
  }
  return result;
}

References Exception, personalPlayback::fp, reco_skim_cfg_mod::fullpath, timingPdfMaker::infile, mps_splice::line, mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by setSipmMap().

scaleByDose()

std::pair< double, double > HGCalSciNoiseMap::scaleByDose	(	const HGCScintillatorDetId &	cellId,
		const radiiVec &	radius
	)

Definition at line 38 of file HGCalSciNoiseMap.cc.

                                                                                                                {
  if (getDoseMap().empty())
    return std::make_pair(1., 0.);
 
  //formula is: A = A0 * exp( -D^0.65 / 199.6)
  //where A0 is the response of the undamaged detector, D is the dose
  int layer = cellId.layer();
  double cellDose = getDoseValue(DetId::HGCalHSc, layer, radius);  //in kRad
  constexpr double expofactor = 1. / 199.6;
  const double dosespower = 0.65;
  double scaleFactor = std::exp(-std::pow(cellDose, dosespower) * expofactor);
 
  //formula is: N = 2.18 * sqrt(F * A / 2e13)
  //where F is the fluence and A is the SiPM area
  double cellFluence = getFluenceValue(DetId::HGCalHSc, layer, radius);  //in 1-Mev-equivalent neutrons per cm2
 
  constexpr double fluencefactor = 2. / (2 * 1e13);  //SiPM area = 2mm^2
  const double normfactor = 2.18;
  double noise = normfactor * sqrt(cellFluence * fluencefactor);
 
  return std::make_pair(scaleFactor, noise);
}

References relativeConstraints::empty, JetChargeProducer_cfi::exp, HGCalRadiationMap::getDoseMap(), HGCalRadiationMap::getDoseValue(), HGCalRadiationMap::getFluenceValue(), DetId::HGCalHSc, HGCScintillatorDetId::layer(), hgcalDigitizer_cfi::noise, funct::pow(), CosmicsPD_Skims::radius, multiplicitycorr_cfi::scaleFactor, and mathSSE::sqrt().

Referenced by HGCHEbackDigitizer::runRealisticDigitizer().

scaleBySipmArea()

double HGCalSciNoiseMap::scaleBySipmArea	(	const HGCScintillatorDetId &	cellId,
		const double &	radius
	)

Definition at line 76 of file HGCalSciNoiseMap.cc.

                                                                                                 {
  if (sipmMap_.empty())
    return 1.;
 
  int layer = cellId.layer();
  if (radius < sipmMap_[layer])
    return 2.;
  else
    return 1.;
}

References HGCScintillatorDetId::layer(), CosmicsPD_Skims::radius, and sipmMap_.

Referenced by HGCHEbackDigitizer::runRealisticDigitizer().

scaleByTileArea()

double HGCalSciNoiseMap::scaleByTileArea	(	const HGCScintillatorDetId &	cellId,
		const radiiVec &	radius
	)

returns the signal scaling and the noise

Definition at line 61 of file HGCalSciNoiseMap.cc.

                                                                                                   {
  double edge;
  if (cellId.type() == 0) {
    constexpr double factor = 2 * M_PI * 1. / 360.;
    edge = radius[0] * factor;  //1 degree
  } else {
    constexpr double factor = 2 * M_PI * 1. / 288.;
    edge = radius[0] * factor;  //1.25 degrees
  }
 
  double scaleFactor = refEdge_ / edge;  //assume reference 3cm of edge
 
  return scaleFactor;
}

References DQMScaleToClient_cfi::factor, M_PI, CosmicsPD_Skims::radius, refEdge_, multiplicitycorr_cfi::scaleFactor, and HGCScintillatorDetId::type().

Referenced by HGCHEbackDigitizer::runRealisticDigitizer().

setSipmMap()

void HGCalSciNoiseMap::setSipmMap

(

const std::string &

fullpath

)

Definition at line 9 of file HGCalSciNoiseMap.cc.

{ sipmMap_ = readSipmPars(fullpath); }

References reco_skim_cfg_mod::fullpath, readSipmPars(), and sipmMap_.

Referenced by HGCHEbackDigitizer::HGCHEbackDigitizer().

Member Data Documentation

refEdge_

const double HGCalSciNoiseMap::refEdge_

private

Definition at line 32 of file HGCalSciNoiseMap.h.

Referenced by scaleByTileArea().

sipmMap_

std::unordered_map<int, float> HGCalSciNoiseMap::sipmMap_

private

Definition at line 34 of file HGCalSciNoiseMap.h.

Referenced by scaleBySipmArea(), and setSipmMap().