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

Definition at line 6 of file HGCalSciNoiseMap.cc.

6 : refEdge_(3.) {}

HGCalSciNoiseMap::refEdge_

const double refEdge_

Definition: HGCalSciNoiseMap.h:32

HGCalSciNoiseMap::~HGCalSciNoiseMap ( )

inline

Definition at line 16 of file HGCalSciNoiseMap.h.

References computeRadius(), readSipmPars(), scaleByDose(), scaleBySipmArea(), scaleByTileArea(), setSipmMap(), and AlCaHLTBitMon_QueryRunRegistry::string.

16 {};

Member Function Documentation

radiiVec HGCalSciNoiseMap::computeRadius ( const HGCScintillatorDetId & cellId )

Definition at line 87 of file HGCalSciNoiseMap.cc.

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

Referenced by HGCHEbackDigitizer::runRealisticDigitizer(), and ~HGCalSciNoiseMap().

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

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

private

Definition at line 12 of file HGCalSciNoiseMap.cc.

References Exception, personalPlayback::fp, edm::FileInPath::fullPath(), timingPdfMaker::infile, mps_splice::line, mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by setSipmMap(), and ~HGCalSciNoiseMap().

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

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

Definition at line 38 of file HGCalSciNoiseMap.cc.

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

Referenced by HGCHEbackDigitizer::runRealisticDigitizer(), and ~HGCalSciNoiseMap().

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

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

Definition at line 76 of file HGCalSciNoiseMap.cc.

References HGCScintillatorDetId::layer(), and sipmMap_.

Referenced by HGCHEbackDigitizer::runRealisticDigitizer(), and ~HGCalSciNoiseMap().

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

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

returns the signal scaling and the noise

Definition at line 61 of file HGCalSciNoiseMap.cc.

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

Referenced by HGCHEbackDigitizer::runRealisticDigitizer(), and ~HGCalSciNoiseMap().

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