Inheritance diagram for HGCalDigiValidation:

Classes
struct	digiInfo

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override

	HGCalDigiValidation (const edm::ParameterSet &)

	~HGCalDigiValidation () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
template<class T1 , class T2 >
void	digiValidation (const T1 &detId, const T2 *geom, int layer, uint16_t adc, double charge)

void	fillDigiInfo (digiInfo &hinfo)

void	fillDigiInfo ()

void	fillOccupancyMap (std::map< int, int > &OccupancyMap, int layer)

Private Attributes
std::vector< MonitorElement * >	ADC_

std::vector< MonitorElement * >	charge_

std::vector< MonitorElement * >	DigiOccupancy_Minus_

std::vector< MonitorElement * >	DigiOccupancy_Plus_

std::vector< MonitorElement * >	DigiOccupancy_XY_

edm::EDGetToken	digiSource_

int	firstLayer_

bool	ifHCAL_

bool	ifNose_

int	layers_

MonitorElement *	MeanDigiOccupancy_Minus_

MonitorElement *	MeanDigiOccupancy_Plus_

std::string	nameDetector_

std::map< int, int >	OccupancyMap_minus_

std::map< int, int >	OccupancyMap_plus_

int	SampleIndx_

int	verbosity_

Detailed Description

Definition at line 55 of file HGCalDigiValidation.cc.

Constructor & Destructor Documentation

HGCalDigiValidation::HGCalDigiValidation ( const edm::ParameterSet & iConfig )

explicit

Definition at line 102 of file HGCalDigiValidation.cc.

References digiSource_, Exception, edm::ParameterSet::getParameter(), ifHCAL_, nameDetector_, and groupFilesInBlocks::temp.

                                                                        :
   nameDetector_(iConfig.getParameter<std::string>("DetectorName")),
   ifNose_(iConfig.getParameter<bool>("ifNose")),
   ifHCAL_(iConfig.getParameter<bool>("ifHCAL")),
   verbosity_(iConfig.getUntrackedParameter<int>("Verbosity",0)),
   SampleIndx_(iConfig.getUntrackedParameter<int>("SampleIndx",0)),
   firstLayer_(1) {
 
   auto temp = iConfig.getParameter<edm::InputTag>("DigiSource");
   if ((nameDetector_ == "HGCalEESensitive") || 
       (nameDetector_ == "HGCalHESiliconSensitive") || 
       (nameDetector_ == "HGCalHEScintillatorSensitive") ||
       (nameDetector_ == "HFNoseSensitive")) {
     digiSource_    = consumes<HGCalDigiCollection>(temp);
   } else if (nameDetector_ == "HCal") {
     if (ifHCAL_) digiSource_ = consumes<QIE11DigiCollection>(temp);
     else         digiSource_ = consumes<HGCalDigiCollection>(temp);
   } else {
     throw cms::Exception("BadHGCDigiSource")
       << "HGCal DetectorName given as " << nameDetector_ << " must be: "
       << "\"HGCalEESensitive\", \"HGCalHESiliconSensitive\", "
       << "\"HGCalHEScintillatorSensitive\", \"HFNoseSensitive\", "
       << "or \"HCal\"!"; 
   }  
 }

HGCalDigiValidation::~HGCalDigiValidation ( )

inlineoverride

Definition at line 68 of file HGCalDigiValidation.cc.

References HGCalDigiValidation::digiInfo::adc, analyze(), bookHistograms(), HGCalDigiValidation::digiInfo::charge, digiValidation(), dqmBeginRun(), fillDescriptions(), fillDigiInfo(), fillOccupancyMap(), relativeConstraints::geom, and HGCalDigiValidation::digiInfo::layer.

68 {}

Member Function Documentation

void HGCalDigiValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 139 of file HGCalDigiValidation.cc.

References ecalMGPA::adc(), HcalCoderDb::adc2fC(), AlignmentProducer_cff::calibrations, ALCARECOTkAlJpsiMuMu_cff::charge, HGCSample::data(), HcalDetId::depth(), digiSource_, digiValidation(), edm::EventID::event(), fillDigiInfo(), muonCSCDigis_cfi::gain, relativeConstraints::geom, motionPolicyChamber::geom0, geometryDiff::geom1, HGCalTopology::geomMode(), edm::EventSetup::get(), edm::Event::getByToken(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalEndcap, HGCalGeometryMode::Hexagon8, HGCalGeometryMode::Hexagon8Full, QIE11DataFrame::id(), edm::EventBase::id(), ifHCAL_, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), HGCScintillatorDetId::layer(), HGCalDetId::layer(), HGCSiliconDetId::layer(), HFNoseDetId::layer(), ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, nameDetector_, OccupancyMap_minus_, OccupancyMap_plus_, HcalCalibrations::pedestal(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), SampleIndx_, edm::SortedCollection< T, SORT >::size(), edm::DataFrameContainer::size(), HcalDetId::subdet(), HGCSample::toa(), HGCalGeometry::topology(), HGCalGeometryMode::Trapezoid, and verbosity_.

Referenced by ~HGCalDigiValidation().

                                                { 
   OccupancyMap_plus_.clear();
   OccupancyMap_minus_.clear();
   
   const HGCalGeometry* geom0(nullptr);
   const CaloGeometry*  geom1(nullptr);
   int geomType(0);
   if (nameDetector_ == "HCal") {
     edm::ESHandle<CaloGeometry> geom;
     iSetup.get<CaloGeometryRecord>().get(geom);
     if (!geom.isValid()) 
       edm::LogVerbatim("HGCalValidation") << "HGCalDigiValidation: Cannot get "
                       << "valid Geometry Object for " 
                       << nameDetector_;
     geom1 = geom.product();
   } else {
     edm::ESHandle<HGCalGeometry> geom;
     iSetup.get<IdealGeometryRecord>().get(nameDetector_, geom);
     if (!geom.isValid()) 
       edm::LogVerbatim("HGCalValidation") << "HGCalDigiValidation: Cannot get "
                       << "valid Geometry Object for " 
                       << nameDetector_;
     geom0 = geom.product();
     HGCalGeometryMode::GeometryMode mode = geom0->topology().geomMode();
     if ((mode == HGCalGeometryMode::Hexagon8) ||
     (mode == HGCalGeometryMode::Hexagon8Full)) geomType = 1;
     else if (mode == HGCalGeometryMode::Trapezoid) geomType = 2;
     if (nameDetector_ == "HFNoseSensitive")        geomType = 3;
   }
 
   unsigned int ntot(0), nused(0);
   if (nameDetector_ == "HGCalEESensitive") {
     //HGCalEE
     edm::Handle<HGCalDigiCollection> theHGCEEDigiContainers;
     iEvent.getByToken(digiSource_, theHGCEEDigiContainers);
     if (theHGCEEDigiContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " 
                         << theHGCEEDigiContainers->size() 
                         << " element(s)";
       for (const auto & it: *(theHGCEEDigiContainers.product())) {
     ntot++; nused++;
     DetId      detId     = it.id();
     int        layer     = ((geomType == 0) ? HGCalDetId(detId).layer() :
                 (geomType == 1) ?
                 HGCSiliconDetId(detId).layer() :
                 HFNoseDetId(detId).layer());
     const HGCSample&  hgcSample = it.sample(SampleIndx_);
     uint16_t   gain      = hgcSample.toa();
     uint16_t   adc       = hgcSample.data();
     double     charge    = adc*gain;
     digiValidation(detId, geom0, layer, adc, charge);
       }
       fillDigiInfo();
     } else {
       edm::LogVerbatim("HGCalValidation") << "DigiCollection handle does not "
                       << "exist for " << nameDetector_;
     }
   } else if ((nameDetector_ == "HGCalHESiliconSensitive") || 
          (nameDetector_ == "HGCalHEScintillatorSensitive")) {
     //HGCalHE
     edm::Handle<HGCalDigiCollection> theHGCHEDigiContainers;
     iEvent.getByToken(digiSource_, theHGCHEDigiContainers);
     if (theHGCHEDigiContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " 
                         << theHGCHEDigiContainers->size()
                         << " element(s)";
       for (const auto & it: *(theHGCHEDigiContainers.product())) {
     ntot++; nused++;
     DetId      detId     = it.id();
     int        layer     = ((geomType == 0) ? HGCalDetId(detId).layer() :
                 ((geomType == 1) ? HGCSiliconDetId(detId).layer() :
                  HGCScintillatorDetId(detId).layer()));
     const HGCSample&  hgcSample = it.sample(SampleIndx_);
     uint16_t   gain      = hgcSample.toa();
     uint16_t   adc       = hgcSample.data();
     double     charge    = adc*gain;
     digiValidation(detId, geom0, layer, adc, charge);
       }
       fillDigiInfo();
     } else {
       edm::LogVerbatim("HGCalValidation") << "DigiCollection handle does not "
                       << "exist for " << nameDetector_;
     }
   } else if ((nameDetector_ == "HCal") && (!ifHCAL_)) {
     //HGCalBH
     edm::Handle<HGCalDigiCollection> theHGCBHDigiContainers;
     iEvent.getByToken(digiSource_, theHGCBHDigiContainers);
     if (theHGCBHDigiContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " 
                         << theHGCBHDigiContainers->size()
                         << " element(s)";
       for (const auto & it: *(theHGCBHDigiContainers.product())) {
     ntot++; nused++;
     HcalDetId  detId     = it.id();
     int        layer     = detId.depth();
     const HGCSample&  hgcSample = it.sample(SampleIndx_);
     uint16_t   gain      = hgcSample.toa();
     uint16_t   adc       = hgcSample.data();
     double     charge    = adc*gain;
     digiValidation(detId, geom1, layer, adc, charge);
       }
       fillDigiInfo();
     } else {
       edm::LogWarning("HGCalValidation") << "DigiCollection handle does not "
                      << "exist for " << nameDetector_;
     }
   } else if (nameDetector_ == "HCal") {
     //HE
     edm::Handle<QIE11DigiCollection>  theHEDigiContainers;
     iEvent.getByToken(digiSource_, theHEDigiContainers);
     if (theHEDigiContainers.isValid()) {
       if (verbosity_>0) 
     edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " 
                         << theHEDigiContainers->size()
                         << " element(s)";
       edm::ESHandle<HcalDbService> conditions;
       iSetup.get<HcalDbRecord > ().get(conditions);
 
       for (const auto & it: *(theHEDigiContainers.product())) {
     QIE11DataFrame df(it);
     HcalDetId detId  = (df.id());
     ntot++;
     if (detId.subdet() == HcalEndcap) {
       nused++;
       HcalCalibrations calibrations = conditions->getHcalCalibrations(detId);
       const HcalQIECoder* channelCoder = conditions->getHcalCoder(detId);
       const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
       HcalCoderDb coder(*channelCoder, *shape);
       CaloSamples tool;
       coder.adc2fC(df, tool);
       int       layer  = detId.depth();
       uint16_t  adc    = (df)[SampleIndx_].adc();
       int       capid  = (df)[SampleIndx_].capid();
       double    charge = (tool[SampleIndx_] - calibrations.pedestal(capid));
       digiValidation(detId, geom1, layer, adc, charge);
     }
       }
       fillDigiInfo();
     } else {
       edm::LogWarning("HGCalValidation") << "DigiCollection handle does not "
                      << "exist for " << nameDetector_;
     }
   } else {
     edm::LogWarning("HGCalValidation") << "invalid detector name !! " 
                        << nameDetector_;
   }
   edm::LogVerbatim("HGCalValidation") << "Event " << iEvent.id().event()
                       << " with " << ntot << " total and "
                       << nused << " used digis";
 }

void HGCalDigiValidation::bookHistograms	(	DQMStore::IBooker &	iB,
		edm::Run const &	,
		edm::EventSetup const &
	)

override

Definition at line 378 of file HGCalDigiValidation.cc.

References ADC_, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), charge_, DEFINE_FWK_MODULE, DigiOccupancy_Minus_, DigiOccupancy_Plus_, DigiOccupancy_XY_, firstLayer_, layers_, MeanDigiOccupancy_Minus_, MeanDigiOccupancy_Plus_, nameDetector_, and DQMStore::IBooker::setCurrentFolder().

Referenced by ~HGCalDigiValidation().

                                            {
   
   iB.setCurrentFolder("HGCAL/HGCalDigisV/"+nameDetector_);
 
   std::ostringstream histoname;
   for (int il = 0; il < layers_; ++il) {
     int ilayer = firstLayer_ + il;
     histoname.str(""); histoname << "charge_"<< "layer_" << ilayer;
     charge_.push_back(iB.book1D(histoname.str().c_str(),"charge_",100,-25,25));
       
     histoname.str(""); histoname << "ADC_" << "layer_" << ilayer;
     ADC_.push_back(iB.book1D(histoname.str().c_str(), "DigiOccupancy",200,0,1000));
       
     histoname.str(""); histoname << "DigiOccupancy_XY_" << "layer_" << ilayer;
     DigiOccupancy_XY_.push_back(iB.book2D(histoname.str().c_str(), "DigiOccupancy", 50, -500, 500, 50, -500, 500));
       
     histoname.str(""); histoname << "DigiOccupancy_Plus_" << "layer_" << ilayer;
     DigiOccupancy_Plus_.push_back(iB.book1D(histoname.str().c_str(), "DigiOccupancy +z", 100, 0, 1000));
     histoname.str(""); histoname << "DigiOccupancy_Minus_" << "layer_" << ilayer;
     DigiOccupancy_Minus_.push_back(iB.book1D(histoname.str().c_str(), "DigiOccupancy -z", 100, 0, 1000));
   }
 
   histoname.str(""); histoname << "SUMOfDigiOccupancy_Plus";
   MeanDigiOccupancy_Plus_ = iB.book1D( histoname.str().c_str(), "SUMOfDigiOccupancy_Plus", layers_, -0.5, layers_-0.5);
   histoname.str(""); histoname << "SUMOfRecDigiOccupancy_Minus";
   MeanDigiOccupancy_Minus_ = iB.book1D( histoname.str().c_str(), "SUMOfDigiOccupancy_Minus", layers_, -0.5,layers_-0.5);
 }

template<class T1 , class T2 >

void HGCalDigiValidation::digiValidation	(	const T1 &	detId,
		const T2 *	geom,
		int	layer,
		uint16_t	adc,
		double	charge
	)

private

Definition at line 295 of file HGCalDigiValidation.cc.

References ecalMGPA::adc(), HGCalDigiValidation::digiInfo::adc, ALCARECOTkAlJpsiMuMu_cff::charge, HGCalDigiValidation::digiInfo::charge, TauDecayModes::dec, fillDigiInfo(), fillOccupancyMap(), firstLayer_, globals_cff::id1, HGCalDigiValidation::digiInfo::layer, OccupancyMap_minus_, OccupancyMap_plus_, verbosity_, HGCalDigiValidation::digiInfo::x, HGCalDigiValidation::digiInfo::y, and HGCalDigiValidation::digiInfo::z.

Referenced by analyze(), and ~HGCalDigiValidation().

                                     {
   
   if (verbosity_>1) edm::LogVerbatim("HGCalValidation") << std::hex 
                             << detId.rawId()
                             << std::dec;
   DetId id1 = DetId(detId.rawId());
   const GlobalPoint& global1 = geom->getPosition(id1);
   
   if (verbosity_>1) 
     edm::LogVerbatim("HGCalValidation") << " adc = "         <<  adc
                     << " charge = "      <<  charge;
   
   digiInfo   hinfo; 
   hinfo.x       =  global1.x();
   hinfo.y       =  global1.y();
   hinfo.z       =  global1.z();
   hinfo.adc     =  adc;
   hinfo.charge  =  charge;
   hinfo.layer   =  layer-firstLayer_;
   
   if (verbosity_>1) 
     edm::LogVerbatim("HGCalValidation") << "gx =  "  << hinfo.x
                     << " gy = "  << hinfo.y
                     << " gz = "  << hinfo.z;
   
   fillDigiInfo(hinfo);
 
   if (global1.eta() > 0)  fillOccupancyMap(OccupancyMap_plus_,  hinfo.layer);
   else                    fillOccupancyMap(OccupancyMap_minus_, hinfo.layer);
   
 }

void HGCalDigiValidation::dqmBeginRun	(	const edm::Run &	,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 355 of file HGCalDigiValidation.cc.

References HGCalDDDConstants::firstLayer(), firstLayer_, edm::EventSetup::get(), HcalDDDRecConstants::getMaxDepth(), HGCalDDDConstants::layers(), layers_, nameDetector_, and verbosity_.

Referenced by ~HGCalDigiValidation().

                                                    {
   if (nameDetector_ == "HCal") {
     edm::ESHandle<HcalDDDRecConstants> pHRNDC;
     iSetup.get<HcalRecNumberingRecord>().get( pHRNDC );
     const HcalDDDRecConstants *hcons   = &(*pHRNDC);
     layers_ = hcons->getMaxDepth(1);
   } else {
     edm::ESHandle<HGCalDDDConstants>  pHGDC;
     iSetup.get<IdealGeometryRecord>().get(nameDetector_, pHGDC);
     const HGCalDDDConstants & hgcons_ = (*pHGDC);
     layers_     = hgcons_.layers(true);
     firstLayer_ = hgcons_.firstLayer();
   }
   
   if (verbosity_>0) 
     edm::LogVerbatim("HGCalValidation") << "current DQM directory:  "
                     << "HGCAL/HGCalDigisV/" 
                     << nameDetector_ << "  layer = "
                     << layers_ << " with the first one at "
                     << firstLayer_;
 }  

void HGCalDigiValidation::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 128 of file HGCalDigiValidation.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~HGCalDigiValidation().

                                                                                      {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("DetectorName","HGCalEESensitive");
   desc.add<edm::InputTag>("DigiSource",edm::InputTag("hgcalDigis","EE"));
   desc.add<bool>("ifNose",false);
   desc.add<bool>("ifHCAL",false);
   desc.addUntracked<int>("Verbosity",0);
   desc.addUntracked<int>("SampleIndx",0);
   descriptions.add("hgcalDigiValidationEEDefault",desc);
 }

void HGCalDigiValidation::fillDigiInfo ( digiInfo & hinfo )

private

Definition at line 335 of file HGCalDigiValidation.cc.

References HGCalDigiValidation::digiInfo::adc, ADC_, HGCalDigiValidation::digiInfo::charge, charge_, DigiOccupancy_XY_, HGCalDigiValidation::digiInfo::layer, HGCalDigiValidation::digiInfo::x, and HGCalDigiValidation::digiInfo::y.

                                                       {
   int ilayer = hinfo.layer;
   charge_.at(ilayer)->Fill(hinfo.charge);
   DigiOccupancy_XY_.at(ilayer)->Fill(hinfo.x, hinfo.y);
   ADC_.at(ilayer)->Fill(hinfo.adc);
 }

void HGCalDigiValidation::fillDigiInfo ( )

private

Definition at line 342 of file HGCalDigiValidation.cc.

References DigiOccupancy_Minus_, DigiOccupancy_Plus_, OccupancyMap_minus_, and OccupancyMap_plus_.

Referenced by analyze(), digiValidation(), and ~HGCalDigiValidation().

                                        {
   for (const auto & itr : OccupancyMap_plus_) {
     int layer     = itr.first;
     int occupancy = itr.second;
     DigiOccupancy_Plus_.at(layer)->Fill(occupancy);
   }
   for (const auto & itr : OccupancyMap_minus_) {
     int layer     = itr.first;
     int occupancy = itr.second;
     DigiOccupancy_Minus_.at(layer)->Fill(occupancy);
   }
 }

void HGCalDigiValidation::fillOccupancyMap	(	std::map< int, int > &	OccupancyMap,
		int	layer
	)

private

Definition at line 329 of file HGCalDigiValidation.cc.

Referenced by digiValidation(), and ~HGCalDigiValidation().

                                   {
   if (OccupancyMap.find(layer) != OccupancyMap.end()) OccupancyMap[layer] ++;
   else                                                OccupancyMap[layer] = 1;
 }