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

#include <HGCalSciNoiseMap.h>

Inheritance diagram for HGCalSciNoiseMap:

Public Member Functions
	HGCalSciNoiseMap ()

std::pair< double, double >	scaleByDose (const HGCScintillatorDetId &, const double)

double	scaleBySipmArea (const HGCScintillatorDetId &, const double)

double	scaleByTileArea (const HGCScintillatorDetId &, const double)
	returns the signal scaling and the noise More...

void	setSipmMap (const std::string &)

	~HGCalSciNoiseMap ()

Public Member Functions inherited from HGCalRadiationMap
const unsigned int &	algo ()

double	computeRadius (const HGCScintillatorDetId &)

const HGCalDDDConstants *	ddd ()

const HGCalGeometry *	geom ()

const doseParametersMap &	getDoseMap ()

double	getDoseValue (const int, const int, const double, bool logVal=false)

double	getFluenceValue (const int, const int, const double, bool logVal=false)

	HGCalRadiationMap ()

void	setDoseMap (const std::string &, const unsigned int)

void	setFluenceScaleFactor (double val)

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.

6 : refEdge_(3.) {}

◆ ~HGCalSciNoiseMap()

HGCalSciNoiseMap::~HGCalSciNoiseMap ( )

inline

Definition at line 16 of file HGCalSciNoiseMap.h.

16 {};

Member Function Documentation

◆ 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 double	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 double	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 * factor;  //1 degree
   } else {
     constexpr double factor = 2 * M_PI * 1. / 288.;
     edge = radius * factor;  //1.25 degrees
   }
  
   double scaleFactor = refEdge_ / edge;  //assume reference 3cm of edge
  
   return scaleFactor;
 }