Inheritance diagram for cms::SiPixelCondObjForHLTReader:

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

void	endJob () override

	SiPixelCondObjForHLTReader (const edm::ParameterSet &iConfig)

	~SiPixelCondObjForHLTReader () override=default

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
std::map< uint32_t, double >	_deadfrac_m

std::map< uint32_t, double >	_noisyfrac_m

TH1F *	_TH1F_Dead_all

TH1F *	_TH1F_Dead_sum

TH1F *	_TH1F_Gains_all

TH1F *	_TH1F_Gains_bpix

TH1F *	_TH1F_Gains_fpix

std::map< uint32_t, TH1F * >	_TH1F_Gains_m

TH1F *	_TH1F_Gains_sum

TH1F *	_TH1F_Noisy_all

TH1F *	_TH1F_Noisy_sum

TH1F *	_TH1F_Pedestals_all

TH1F *	_TH1F_Pedestals_bpix

TH1F *	_TH1F_Pedestals_fpix

std::map< uint32_t, TH1F * >	_TH1F_Pedestals_m

TH1F *	_TH1F_Pedestals_sum

edm::ParameterSet	conf_

std::unique_ptr< SiPixelGainCalibrationServiceBase >	SiPixelGainCalibrationService_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 50 of file SiPixelCondObjForHLTReader.cc.

Constructor & Destructor Documentation

◆ SiPixelCondObjForHLTReader()

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

explicit

Definition at line 87 of file SiPixelCondObjForHLTReader.cc.

References conf_, edm::EDConsumerBase::consumesCollector(), edm::ParameterSet::getParameter(), and SiPixelGainCalibrationService_.

       : conf_(conf), tkGeomToken_(esConsumes()) {
     if (conf_.getParameter<bool>("useSimRcd"))
       SiPixelGainCalibrationService_ =
           std::make_unique<SiPixelGainCalibrationForHLTSimService>(conf_, consumesCollector());
     else
       SiPixelGainCalibrationService_ =
           std::make_unique<SiPixelGainCalibrationForHLTService>(conf_, consumesCollector());
   }

◆ ~SiPixelCondObjForHLTReader()

cms::SiPixelCondObjForHLTReader::~SiPixelCondObjForHLTReader ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 97 of file SiPixelCondObjForHLTReader.cc.

References _deadfrac_m, _noisyfrac_m, _TH1F_Dead_all, _TH1F_Dead_sum, _TH1F_Gains_all, _TH1F_Gains_bpix, _TH1F_Gains_fpix, _TH1F_Gains_m, _TH1F_Gains_sum, _TH1F_Noisy_all, _TH1F_Noisy_sum, _TH1F_Pedestals_all, _TH1F_Pedestals_bpix, _TH1F_Pedestals_fpix, _TH1F_Pedestals_m, _TH1F_Pedestals_sum, conf_, TrackerGeometry::dets(), compareTotals::fs, PedestalClient_cfi::gain, edm::EventSetup::getData(), edm::ParameterSet::getUntrackedParameter(), TrackerGeometry::idToDetUnit(), TFileDirectory::make(), Skims_PA_cff::name, hgcalPlots::ncols, PixelTopology::ncolumns(), PixelTopology::nrows(), QIE10Task_cfi::ped, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, SiPixelGainCalibrationService_, PixelGeomDetUnit::specificTopology(), DetId::subdetId(), and tkGeomToken_.

                                                                                               {
     //Create Subdirectories
     edm::Service<TFileService> fs;
     TFileDirectory subDirPed = fs->mkdir("Pedestals");
     TFileDirectory subDirGain = fs->mkdir("Gains");
     char name[128];
 
     unsigned int nmodules = 0;
     uint32_t nchannels = 0;
     uint32_t ndead = 0;
     uint32_t nnoisy = 0;
 
     // Get the calibration data
     SiPixelGainCalibrationService_->setESObjects(iSetup);
     edm::LogInfo("SiPixelCondObjForHLTReader")
         << "[SiPixelCondObjForHLTReader::beginJob] End Reading CondObjForHLTects" << std::endl;
 
     // Get the Geometry
     const TrackerGeometry *tkgeom = &iSetup.getData(tkGeomToken_);
     edm::LogInfo("SiPixelCondObjForHLTReader") << " There are " << tkgeom->dets().size() << " detectors" << std::endl;
 
     // Get the list of DetId's
     std::vector<uint32_t> vdetId_ = SiPixelGainCalibrationService_->getDetIds();
 
     //Create histograms
     _TH1F_Dead_sum = fs->make<TH1F>(
         "Summary_dead", "Dead pixel fraction (0=dead, 1=alive)", vdetId_.size() + 1, 0, vdetId_.size() + 1);
     _TH1F_Dead_all = fs->make<TH1F>("DeadAll",
                                     "Dead pixel fraction (0=dead, 1=alive)",
                                     50,
                                     0.,
                                     conf_.getUntrackedParameter<double>("maxRangeDeadPixHist", 0.001));
     _TH1F_Noisy_sum = fs->make<TH1F>(
         "Summary_noisy", "Noisy pixel fraction (0=noisy, 1=alive)", vdetId_.size() + 1, 0, vdetId_.size() + 1);
     _TH1F_Noisy_all = fs->make<TH1F>("NoisyAll",
                                      "Noisy pixel fraction (0=noisy, 1=alive)",
                                      50,
                                      0.,
                                      conf_.getUntrackedParameter<double>("maxRangeDeadPixHist", 0.001));
     _TH1F_Gains_sum = fs->make<TH1F>("Summary_Gain", "Gain Summary", vdetId_.size() + 1, 0, vdetId_.size() + 1);
     _TH1F_Pedestals_sum =
         fs->make<TH1F>("Summary_Pedestal", "Pedestal Summary", vdetId_.size() + 1, 0, vdetId_.size() + 1);
     _TH1F_Pedestals_all = fs->make<TH1F>("PedestalsAll", "all Pedestals", 350, -100, 250);
     _TH1F_Pedestals_bpix = fs->make<TH1F>("PedestalsBpix", "bpix Pedestals", 350, -100, 250);
     _TH1F_Pedestals_fpix = fs->make<TH1F>("PedestalsFpix", "fpix Pedestals", 350, -100, 250);
     _TH1F_Gains_all = fs->make<TH1F>("GainsAll", "all Gains", 100, 0, 10);
     _TH1F_Gains_bpix = fs->make<TH1F>("GainsBpix", "bpix Gains", 100, 0, 10);
     _TH1F_Gains_fpix = fs->make<TH1F>("GainsFpix", "fpix Gains", 100, 0, 10);
 
     TTree *tree = new TTree("tree", "tree");
     uint32_t detid;
     double gainmeanfortree, gainrmsfortree, pedmeanfortree, pedrmsfortree;
     tree->Branch("detid", &detid, "detid/I");
     tree->Branch("ped_mean", &pedmeanfortree, "ped_mean/D");
     tree->Branch("ped_rms", &pedrmsfortree, "ped_rms/D");
     tree->Branch("gain_mean", &gainmeanfortree, "gain_mean/D");
     tree->Branch("gain_rms", &gainrmsfortree, "gain_rms/D");
 
     // Loop over DetId's
     int ibin = 1;
     for (std::vector<uint32_t>::const_iterator detid_iter = vdetId_.begin(); detid_iter != vdetId_.end();
          detid_iter++) {
       detid = *detid_iter;
 
       sprintf(name, "Pedestals_%d", detid);
       _TH1F_Pedestals_m[detid] = subDirPed.make<TH1F>(name, name, 350, -100., 250.);
       sprintf(name, "Gains_%d", detid);
       _TH1F_Gains_m[detid] = subDirGain.make<TH1F>(name, name, 100, 0., 10.);
 
       DetId detIdObject(detid);
       const PixelGeomDetUnit *_PixelGeomDetUnit =
           dynamic_cast<const PixelGeomDetUnit *>(tkgeom->idToDetUnit(DetId(detid)));
       if (_PixelGeomDetUnit == nullptr) {
         edm::LogError("SiPixelCondObjHLTDisplay") << "[SiPixelCondObjHLTReader::beginJob] the detID " << detid
                                                   << " doesn't seem to belong to Tracker" << std::endl;
         continue;
       }
 
       _deadfrac_m[detid] = 0.;
       _noisyfrac_m[detid] = 0.;
 
       nmodules++;
 
       const GeomDetUnit *geoUnit = tkgeom->idToDetUnit(detIdObject);
       const PixelGeomDetUnit *pixDet = dynamic_cast<const PixelGeomDetUnit *>(geoUnit);
       const PixelTopology &topol = pixDet->specificTopology();
 
       // Get the module sizes.
       int nrows = topol.nrows();     // rows in x
       int ncols = topol.ncolumns();  // cols in y
       float nchannelspermod = 0;
 
       for (int col_iter = 0; col_iter < ncols; col_iter++) {
         for (int row_iter = 0; row_iter < nrows; row_iter++) {
           nchannelspermod++;
           nchannels++;
 
           if (SiPixelGainCalibrationService_->isDead(detid, col_iter, row_iter)) {
             //      edm::LogPrint("SiPixelCondObjForHLTReader") << "found dead pixel " << detid << " " <<col_iter << "," << row_iter << std::endl;
             ndead++;
             _deadfrac_m[detid]++;
             continue;
           } else if (SiPixelGainCalibrationService_->isNoisy(detid, col_iter, row_iter)) {
             //      edm::LogPrint("SiPixelCondObjForHLTReader") << "found noisy pixel " << detid << " " <<col_iter << "," << row_iter << std::endl;
             nnoisy++;
             _noisyfrac_m[detid]++;
             continue;
           }
 
           float gain = SiPixelGainCalibrationService_->getGain(detid, col_iter, row_iter);
           _TH1F_Gains_m[detid]->Fill(gain);
           _TH1F_Gains_all->Fill(gain);
 
           if (detIdObject.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel))
             _TH1F_Gains_bpix->Fill(gain);
           if (detIdObject.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap))
             _TH1F_Gains_fpix->Fill(gain);
 
           float ped = SiPixelGainCalibrationService_->getPedestal(detid, col_iter, row_iter);
           _TH1F_Pedestals_m[detid]->Fill(ped);
           _TH1F_Pedestals_all->Fill(ped);
           //     edm::LogPrint("SiPixelCondObjForHLTReader")<<"detid  "<<detid<<"     ped "<<ped<<std::endl;
 
           if (detIdObject.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel))
             _TH1F_Pedestals_bpix->Fill(ped);
           if (detIdObject.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap))
             _TH1F_Pedestals_fpix->Fill(ped);
 
           //     edm::LogPrint("SiPixelCondObjForHLTReader") <<"    DetId "<<detid<<"       Col "<<col_iter<<" Row "<<row_iter<<" Ped "<<ped<<" Gain "<<gain<<std::endl;
         }
       }
 
       _deadfrac_m[detid] /= nchannelspermod;
       _noisyfrac_m[detid] /= nchannelspermod;
       _TH1F_Dead_sum->SetBinContent(ibin, _deadfrac_m[detid]);
       _TH1F_Dead_all->Fill(_deadfrac_m[detid]);
       _TH1F_Noisy_sum->SetBinContent(ibin, _noisyfrac_m[detid]);
       _TH1F_Noisy_all->Fill(_noisyfrac_m[detid]);
       _TH1F_Gains_sum->SetBinContent(ibin, _TH1F_Gains_m[detid]->GetMean());
       _TH1F_Gains_sum->SetBinError(ibin, _TH1F_Gains_m[detid]->GetRMS());
       _TH1F_Pedestals_sum->SetBinContent(ibin, _TH1F_Pedestals_m[detid]->GetMean());
       _TH1F_Pedestals_sum->SetBinError(ibin, _TH1F_Pedestals_m[detid]->GetRMS());
 
       gainmeanfortree = _TH1F_Gains_m[detid]->GetMean();
       gainrmsfortree = _TH1F_Gains_m[detid]->GetRMS();
       pedmeanfortree = _TH1F_Pedestals_m[detid]->GetMean();
       pedrmsfortree = _TH1F_Pedestals_m[detid]->GetRMS();
       //edm::LogPrint("SiPixelCondObjForHLTReader")<<"DetId "<<detid<<"       GainMean "<<gainmeanfortree<<" RMS "<<gainrmsfortree<<"      PedMean "<<pedmeanfortree<<" RMS "<<pedrmsfortree<<std::endl;
       tree->Fill();
 
       ibin++;
     }
 
     edm::LogInfo("SiPixelCondObjForHLTReader")
         << "[SiPixelCondObjForHLTReader::analyze] ---> PIXEL Modules  " << nmodules << std::endl;
     edm::LogInfo("SiPixelCondObjForHLTReader")
         << "[SiPixelCondObjForHLTReader::analyze] ---> PIXEL Channels (i.e. Number of Columns)" << nchannels
         << std::endl;
 
     edm::LogPrint("SiPixelCondObjForHLTReader") << " ---> SUMMARY :" << std::endl;
     edm::LogPrint("SiPixelCondObjForHLTReader") << "Encounted " << ndead << " dead pixels" << std::endl;
     edm::LogPrint("SiPixelCondObjForHLTReader") << "Encounted " << nnoisy << " noisy pixels" << std::endl;
     edm::LogPrint("SiPixelCondObjForHLTReader")
         << "The Gain Mean is " << _TH1F_Gains_all->GetMean() << " with rms " << _TH1F_Gains_all->GetRMS() << std::endl;
     edm::LogPrint("SiPixelCondObjForHLTReader") << "         in BPIX " << _TH1F_Gains_bpix->GetMean() << " with rms "
                                                 << _TH1F_Gains_bpix->GetRMS() << std::endl;
     edm::LogPrint("SiPixelCondObjForHLTReader") << "         in FPIX " << _TH1F_Gains_fpix->GetMean() << " with rms "
                                                 << _TH1F_Gains_fpix->GetRMS() << std::endl;
     edm::LogPrint("SiPixelCondObjForHLTReader") << "The Ped Mean is " << _TH1F_Pedestals_all->GetMean() << " with rms "
                                                 << _TH1F_Pedestals_all->GetRMS() << std::endl;
     edm::LogPrint("SiPixelCondObjForHLTReader") << "         in BPIX " << _TH1F_Pedestals_bpix->GetMean()
                                                 << " with rms " << _TH1F_Pedestals_bpix->GetRMS() << std::endl;
     edm::LogPrint("SiPixelCondObjForHLTReader") << "         in FPIX " << _TH1F_Pedestals_fpix->GetMean()
                                                 << " with rms " << _TH1F_Pedestals_fpix->GetRMS() << std::endl;
   }

◆ endJob()

void SiPixelCondObjForHLTReader::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 274 of file SiPixelCondObjForHLTReader.cc.

                                           {
     edm::LogPrint("SiPixelCondObjForHLTReader") << " ---> End job " << std::endl;
   }

◆ fillDescriptions()

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

static

Definition at line 278 of file SiPixelCondObjForHLTReader.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), and submitPVResolutionJobs::desc.

                                                                                               {
     edm::ParameterSetDescription desc;
     desc.setComment("EDAnalyzer to read per-module SiPixelGainCalibrationForHLT payloads in the EventSetup");
     desc.add<bool>("useSimRcd", false);
     desc.addUntracked<double>("maxRangeDeadPixHist", 0.001);
     descriptions.addWithDefaultLabel(desc);
   }