#include <SiStripMonitorPedestals.h>

Inheritance diagram for SiStripMonitorPedestals:

Classes
struct	ModMEs

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

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

void	dqmEndRun (edm::Run const &run, edm::EventSetup const &eSetup) override

void	endJob () override

	SiStripMonitorPedestals (const edm::ParameterSet &)

	~SiStripMonitorPedestals () override

Public Member Functions inherited from DQMOneEDAnalyzer<>
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) 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)

Private Member Functions
void	createMEs (DQMStore::IBooker &, edm::EventSetup const &eSetup)

void	fillCondDBMEs (edm::EventSetup const &eSetup)

void	resetMEs (uint32_t idet)

Private Attributes
bool	analyzed

ApvAnalysisFactory *	apvFactory_

edm::ParameterSet	conf_

const SiStripDetCabling *	detcabling

edm::ESGetToken< SiStripDetCabling, SiStripDetCablingRcd >	detCablingToken_

edm::ESWatcher< SiStripDetCablingRcd >	detCablingWatcher_

edm::EDGetTokenT< edm::DetSetVector< SiStripRawDigi > >	digiToken_

DQMStore *	dqmStore_

bool	firstEvent

int16_t	nEvInit_

int16_t	nEvTot_

int16_t	nEvUpdate_

int	nIteration_

edm::ESGetToken< SiStripNoises, SiStripNoisesRcd >	noiseToken_

int	NumCMstripsInGroup_

std::string	outPutFileName

edm::ESGetToken< SiStripPedestals, SiStripPedestalsRcd >	pedestalToken_

std::map< uint32_t, ModMEs >	PedMEs

edm::ParameterSet	pedsPSet_

edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	qualityToken_

std::string	runTypeFlag_

int16_t	signalCutPeds_

int	theEventInitNumber_

int	theEventIterNumber_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

Static Private Attributes
static const std::string	RunMode1 = "ConDBPlotsOnly"

static const std::string	RunMode2 = "CalculatedPlotsOnly"

static const std::string	RunMode3 = "AllPlots"

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer<>
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from DQMOneEDAnalyzer<>
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

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)

Protected Attributes inherited from DQMOneEDAnalyzer<>
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Description: <one line="" class="" summary>="">

Usage: <usage>

Definition at line 51 of file SiStripMonitorPedestals.h.

Constructor & Destructor Documentation

◆ SiStripMonitorPedestals()

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

explicit

Definition at line 50 of file SiStripMonitorPedestals.cc.

References conf_, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::digiProducer, digiToken_, edm::ParameterSet::getParameter(), ProducerED_cfi::InputTag, noiseToken_, NumCMstripsInGroup_, pedestalToken_, pedsPSet_, qualityToken_, RunMode1, RunMode3, runTypeFlag_, AlCaHLTBitMon_QueryRunRegistry::string, theEventInitNumber_, and theEventIterNumber_.

     : dqmStore_(edm::Service<DQMStore>().operator->()),
       conf_(iConfig),
       pedsPSet_(iConfig.getParameter<edm::ParameterSet>("PedestalsPSet")),
       analyzed(false),
       firstEvent(true),
       signalCutPeds_(4),
       nEvTot_(0),
       nIteration_(0),
       apvFactory_(nullptr),
       tTopoToken_(esConsumes<edm::Transition::BeginRun>()),
       detCablingToken_(esConsumes<edm::Transition::BeginRun>()) {
   // retrieve producer name of input StripDigiCollection
   std::string digiProducer = conf_.getParameter<std::string>("DigiProducer");
   std::string digiType = "VirginRaw";
   digiToken_ = consumes<edm::DetSetVector<SiStripRawDigi>>(edm::InputTag(digiProducer, digiType));
 
   edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals  "
                                           << " Constructing....... ";
 
   theEventInitNumber_ = pedsPSet_.getParameter<int>("NumberOfEventsForInit");
   theEventIterNumber_ = pedsPSet_.getParameter<int>("NumberOfEventsForIteration");
   NumCMstripsInGroup_ = pedsPSet_.getParameter<int>("NumCMstripsInGroup");
   runTypeFlag_ = conf_.getParameter<std::string>("RunTypeFlag");
 
   if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3) {
     pedestalToken_ = esConsumes<edm::Transition::BeginRun>();
     noiseToken_ = esConsumes<edm::Transition::BeginRun>();
     qualityToken_ = esConsumes<edm::Transition::BeginRun>(
         edm::ESInputTag{"", iConfig.getParameter<std::string>("StripQualityLabel")});
   }
 }

◆ ~SiStripMonitorPedestals()

SiStripMonitorPedestals::~SiStripMonitorPedestals ( )

override

Definition at line 85 of file SiStripMonitorPedestals.cc.

References apvFactory_.

                                                   {
   edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals  "
                                           << " Destructing...... ";
   if (apvFactory_) {
     delete apvFactory_;
   }
 }

Member Function Documentation

◆ analyze()

void SiStripMonitorPedestals::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	eSetup
	)

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 273 of file SiStripMonitorPedestals.cc.

References apvFactory_, SiStripMonitorPedestals::ModMEs::CMDistribution, SiStripMonitorPedestals::ModMEs::CMSlopeDistribution, SiStripMonitorPedestals::ModMEs::CMSubNoisePerStrip, SiStripMonitorPedestals::ModMEs::CMSubNoiseProfile, gather_cfg::cout, detcabling, detCablingToken_, digiToken_, HcalObjRepresent::Fill(), firstEvent, ApvAnalysisFactory::getCommonMode(), ApvAnalysisFactory::getCommonModeSlope(), SiStripDetCabling::getConnections(), edm::EventSetup::getData(), ApvAnalysisFactory::getMask(), ApvAnalysisFactory::getNoise(), ApvAnalysisFactory::getPedestal(), ApvAnalysisFactory::getRawNoise(), mps_fire::i, iEvent, createfilelist::int, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmdumpme::k, nEvTot_, nIteration_, SiStripMonitorPedestals::ModMEs::NoisyStripDistribution, SiStripMonitorPedestals::ModMEs::NoisyStrips, NumCMstripsInGroup_, PedMEs, SiStripMonitorPedestals::ModMEs::PedsDistribution, SiStripMonitorPedestals::ModMEs::PedsEvolution, SiStripMonitorPedestals::ModMEs::PedsPerStrip, SiStripMonitorPedestals::ModMEs::RawNoisePerStrip, SiStripMonitorPedestals::ModMEs::RawNoiseProfile, RunMode2, RunMode3, runTypeFlag_, groupFilesInBlocks::temp, theEventInitNumber_, theEventIterNumber_, createJobs::tmp, and ApvAnalysisFactory::update().

                                                                                          {
   edm::LogInfo("SiStripMonitorPedestals")
       << "SiStripMonitorPedestals::analyze: Run " << iEvent.id().run() << " Event " << iEvent.id().event();
 
   detcabling = &eSetup.getData(detCablingToken_);
 
   if (runTypeFlag_ != RunMode2 && runTypeFlag_ != RunMode3)
     return;
 
   // Increment # of Events
   nEvTot_++;
 
   // get DigiCollection object from Event
   edm::Handle<edm::DetSetVector<SiStripRawDigi>> digi_collection;
   // you have a collection as there are all the digis for the event for every
   // detector
   iEvent.getByToken(digiToken_, digi_collection);
 
   // Increase the number of iterations ...
   if ((nEvTot_ - theEventInitNumber_) % theEventIterNumber_ == 1)
     nIteration_++;
 
   // loop over all MEs
   for (std::map<uint32_t, ModMEs>::const_iterator i = PedMEs.begin(); i != PedMEs.end(); i++) {
     uint32_t detid = i->first;
     ModMEs local_modmes = i->second;
     // get iterators for digis belonging to one DetId, it is an iterator, i.e.
     // one element of the vector
     std::vector<edm::DetSet<SiStripRawDigi>>::const_iterator digis = digi_collection->find(detid);
     if (digis == digi_collection->end() || digis->data.empty() || digis->data.size() > 768) {
       if (digis == digi_collection->end()) {
         edm::LogError("SiStripMonitorPedestals") << " SiStripMonitorPedestals::analyze: Event " << nEvTot_ << " DetId "
                                                  << detid << " at the end of Digi Collection!!!";
       } else {
         edm::LogError("SiStripMonitorPedestals") << " [SiStripMonitorPedestals::analyze: Event " << nEvTot_ << " DetId "
                                                  << detid << " # of Digis " << digis->data.size();
       }
       std::vector<const FedChannelConnection *> fed_conns = detcabling->getConnections(detid);
       bool firstchannel(true);
       for (unsigned int k = 0; k < fed_conns.size(); k++) {
         if (fed_conns[k] && fed_conns[k]->isConnected()) {
           if (firstchannel) {
             edm::LogError("SiStripMonitorPedestals") << " SiStripMonitorPedestals::analyze: Fed Id "
                                                      << fed_conns[k]->fedId() << " Channel " << fed_conns[k]->fedCh();
             firstchannel = false;
           } else
             edm::LogError("SiStripMonitorPedestals")
                 << "  SiStripMonitorPedestals::analyze: Channel " << fed_conns[k]->fedCh();
         }
       }
       std::cout << std::endl;
       continue;
     }
 
     if (digis->data.empty()) {
       edm::LogError("MonitorDigi_tmp") << "[SiStripRawDigiToRaw::createFedBuffers] Zero digis found!";
     }
     uint32_t id = detid;
     //    cout <<"Data size "<<digis->data.size()<<endl;
     apvFactory_->update(id, (*digis));
 
     if (nEvTot_ > theEventInitNumber_) {
       if (local_modmes.CMDistribution != nullptr) {
         std::vector<float> tmp;
         tmp.clear();
         apvFactory_->getCommonMode(id, tmp);
         // unpacking the info looking at the right topology
         int numberCMBlocks = int(128. / NumCMstripsInGroup_);
         int ibin = 0;
         for (std::vector<float>::const_iterator iped = tmp.begin(); iped != tmp.end(); iped++) {
           int iapv = int(ibin / numberCMBlocks);
           (local_modmes.CMDistribution)->Fill(iapv, static_cast<float>(*iped));
           ibin++;
         }
       }
       if (local_modmes.CMSlopeDistribution != nullptr) {
         std::vector<float> tmp;
         tmp.clear();
         int iapv = 0;
         apvFactory_->getCommonModeSlope(id, tmp);
         for (std::vector<float>::const_iterator it = tmp.begin(); it != tmp.end(); it++) {
           (local_modmes.CMSlopeDistribution)->Fill(iapv, static_cast<float>(*it));
           iapv++;
         }
       }
     }
 
     // asking for the status
     if ((nEvTot_ - theEventInitNumber_ - theEventIterNumber_) % theEventIterNumber_ == 1) {
       std::vector<float> tmp;
       tmp.clear();
       apvFactory_->getPedestal(id, tmp);
       if (local_modmes.PedsPerStrip != nullptr) {
         int numberOfApvs = int(tmp.size() / 128.);
         for (int i = 0; i < numberOfApvs; i++) {
           std::vector<float> myPedPerApv;
           apvFactory_->getPedestal(id, i, myPedPerApv);
           float avarage = 0;
           avarage = std::accumulate(myPedPerApv.begin(), myPedPerApv.end(), avarage);
           avarage = avarage / 128.;
           (local_modmes.PedsEvolution)->setBinContent(i + 1, nIteration_, avarage);
         }
         int ibin = 0;
 
         for (std::vector<float>::const_iterator iped = tmp.begin(); iped != tmp.end(); iped++) {
           int napv = int(ibin / 128.);
           ibin++;
           float last_value = (local_modmes.PedsPerStrip)->getBinContent(ibin);
           if (last_value != 0.) {
             (local_modmes.PedsPerStrip)->setBinContent(ibin, (static_cast<float>(*iped) + last_value) / 2.);
           } else {
             (local_modmes.PedsPerStrip)->setBinContent(ibin, static_cast<float>(*iped));
           }
           (local_modmes.PedsDistribution)->Fill(napv, static_cast<float>(*iped));
         }
       }
 
       if (local_modmes.CMSubNoisePerStrip != nullptr && local_modmes.CMSubNoiseProfile != nullptr) {
         tmp.clear();
         apvFactory_->getNoise(id, tmp);
         int ibin = 0;
         for (std::vector<float>::const_iterator iped = tmp.begin(); iped != tmp.end(); iped++) {
           ibin++;
           (local_modmes.CMSubNoiseProfile)->Fill(static_cast<double>(ibin * 1.), static_cast<float>(*iped));
 
           float last_value = (local_modmes.CMSubNoisePerStrip)->getBinContent(ibin);
           if (last_value != 0.) {
             (local_modmes.CMSubNoisePerStrip)->setBinContent(ibin, (static_cast<float>(*iped) + last_value) / 2.);
           } else {
             (local_modmes.CMSubNoisePerStrip)->setBinContent(ibin, static_cast<float>(*iped));
           }
         }
       }
 
       if (local_modmes.RawNoisePerStrip != nullptr && local_modmes.RawNoiseProfile != nullptr) {
         tmp.clear();
         apvFactory_->getRawNoise(id, tmp);
         int ibin = 0;
         for (std::vector<float>::const_iterator iped = tmp.begin(); iped != tmp.end(); iped++) {
           ibin++;
           (local_modmes.RawNoiseProfile)->Fill(static_cast<double>(ibin * 1.), static_cast<float>(*iped));
           float last_value = (local_modmes.RawNoisePerStrip)->getBinContent(ibin);
           if (last_value != 0.) {
             (local_modmes.RawNoisePerStrip)->setBinContent(ibin, (static_cast<float>(*iped) + last_value) / 2.);
           } else {
             (local_modmes.RawNoisePerStrip)->setBinContent(ibin, static_cast<float>(*iped));
           }
         }
       }
 
       if (local_modmes.NoisyStrips != nullptr) {
         TkApvMask::MaskType temp;
         apvFactory_->getMask(id, temp);
         int ibin = 0;
         for (TkApvMask::MaskType::const_iterator iped = temp.begin(); iped != temp.end(); iped++) {
           ibin++;
 
           if (nIteration_ < 2) {
             if (*iped == 1)
               (local_modmes.NoisyStrips)->Fill(ibin, 3.);
             if (*iped == 2)
               (local_modmes.NoisyStrips)->Fill(ibin, 4.);
             if (*iped == 0)
               (local_modmes.NoisyStrips)->Fill(ibin, 0.);
           } else {
             (local_modmes.NoisyStrips)->Fill(ibin, static_cast<float>(*iped));
             (local_modmes.NoisyStripDistribution)->Fill(static_cast<float>(*iped));
           }
         }
       }
     }
   }
   if (firstEvent)
     firstEvent = false;
 }

◆ bookHistograms()

void SiStripMonitorPedestals::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

overridevirtual

Implements DQMOneEDAnalyzer<>.

Definition at line 96 of file SiStripMonitorPedestals.cc.

References edm::ESWatcher< T >::check(), createMEs(), detcabling, detCablingToken_, detCablingWatcher_, fillCondDBMEs(), edm::EventSetup::getData(), RunMode1, RunMode3, and runTypeFlag_.

                                                                           {
   if (detCablingWatcher_.check(eSetup)) {
     detcabling = &eSetup.getData(detCablingToken_);
     edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals::bookHistograms: "
                                             << " Creating MEs for new Cabling ";
     createMEs(ibooker, eSetup);
   }
 
   if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3)
     fillCondDBMEs(eSetup);
 }

◆ createMEs()

void SiStripMonitorPedestals::createMEs	(	DQMStore::IBooker &	ibooker,
		edm::EventSetup const &	eSetup
	)

private

Definition at line 113 of file SiStripMonitorPedestals.cc.

Referenced by bookHistograms().

                                                                                          {
   const auto tTopo = &es.getData(tTopoToken_);
 
   std::vector<uint32_t> SelectedDetIds;
 
   // ApvAnalysisFactory
   if (apvFactory_) {
     delete apvFactory_;
   }
   apvFactory_ = new ApvAnalysisFactory(pedsPSet_);
 
   detcabling->addActiveDetectorsRawIds(SelectedDetIds);
 
   // use SistripHistoId for producing histogram id (and title)
   SiStripHistoId hidmanager;
   // create SiStripFolderOrganizer
   SiStripFolderOrganizer folder_organizer;
 
   edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals::createMEs: "
                                           << "Number of Detector Present in cabling " << SelectedDetIds.size();
 
   for (std::vector<uint32_t>::const_iterator idetid = SelectedDetIds.begin(), iEnd = SelectedDetIds.end();
        idetid != iEnd;
        ++idetid) {
     uint32_t detid = *idetid;
 
     // Check consistency in DetId
     if (detid == 0 || detid == 0xFFFFFFFF) {
       edm::LogError("SiStripMonitorPedestals") << "SiStripMonitorPedestals::createMEs: "
                                                << "Wrong DetId !!!!!! " << detid << " Neglecting !!!!!! ";
       continue;
     }
 
     unsigned int apv_pairs = detcabling->nApvPairs(detid);
 
     // Check consistency in Apv numbers
     if (apv_pairs < 1 || apv_pairs > 3) {
       edm::LogError("SiStripMonitorPedestals") << "SiStripMonitorPedestals::createMEs: Wrong APV Pairs  => detId "
                                                << detid << " APV pairs " << apv_pairs << " Neglecting !!!!!! ";
       continue;
     }
     unsigned int napvs = apv_pairs * 2;
     unsigned int nStrip = napvs * 128;
 
     bool newDetId = apvFactory_->instantiateApvs(detid, napvs);
 
     if (newDetId) {
       ModMEs local_modmes;
       local_modmes.PedsPerStrip = nullptr;
       local_modmes.PedsDistribution = nullptr;
       local_modmes.PedsEvolution = nullptr;
       local_modmes.CMSubNoisePerStrip = nullptr;
       local_modmes.RawNoisePerStrip = nullptr;
       local_modmes.CMSubNoiseProfile = nullptr;
       local_modmes.RawNoiseProfile = nullptr;
       local_modmes.NoisyStrips = nullptr;
       local_modmes.NoisyStripDistribution = nullptr;
       local_modmes.CMDistribution = nullptr;
       local_modmes.CMSlopeDistribution = nullptr;
       local_modmes.PedsPerStripDB = nullptr;
       local_modmes.CMSubNoisePerStripDB = nullptr;
       local_modmes.BadStripsDB = nullptr;
 
       std::string hid;
       // set appropriate folder using SiStripFolderOrganizer
       folder_organizer.setDetectorFolder(detid, tTopo);  // pass the detid to this method
 
       // if the deid already exists in the map, then reset MEs otherwise create
       // them
       resetMEs(detid);
 
       if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3) {
         // Pedestals histos
         hid = hidmanager.createHistoId("PedestalFromCondDB", "det", detid);
         local_modmes.PedsPerStripDB =
             ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);  // to modify the size binning
         (local_modmes.PedsPerStripDB)->setAxisTitle("Pedestal from CondDB(ADC) vs Strip Number", 1);
 
         hid = hidmanager.createHistoId("NoiseFromCondDB", "det", detid);
         local_modmes.CMSubNoisePerStripDB = ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);
         (local_modmes.CMSubNoisePerStripDB)->setAxisTitle("CMSubNoise from CondDB(ADC) vs Strip Number", 1);
 
         hid = hidmanager.createHistoId("BadStripFlagCondDB", "det", detid);
         local_modmes.BadStripsDB = ibooker.book2D(hid, hid, nStrip, 0.5, nStrip + 0.5, 6, -0.5, 5.5);
         (local_modmes.BadStripsDB)->setAxisTitle("Strip Flag from CondDB(ADC) vs Strip Number", 1);
       }
       if (runTypeFlag_ == RunMode2 || runTypeFlag_ == RunMode3) {
         // Pedestals histos
         hid = hidmanager.createHistoId("PedsPerStrip", "det", detid);
         local_modmes.PedsPerStrip = ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);  // to modify the size binning
         (local_modmes.PedsPerStrip)->setAxisTitle("Pedestal (ADC)  vs Strip Number ", 1);
 
         hid = hidmanager.createHistoId("PedsDistribution", "det", detid);
         local_modmes.PedsDistribution = ibooker.book2D(hid,
                                                        hid,
                                                        napvs,
                                                        -0.5,
                                                        napvs - 0.5,
                                                        300,
                                                        200,
                                                        500);  // to modify the size binning
         (local_modmes.PedsDistribution)->setAxisTitle("Apv Number", 1);
         (local_modmes.PedsDistribution)->setAxisTitle("Mean Pedestal Value (ADC)", 2);
 
         hid = hidmanager.createHistoId("PedsEvolution", "det", detid);
         local_modmes.PedsEvolution = ibooker.book2D(hid,
                                                     hid,
                                                     napvs,
                                                     -0.5,
                                                     napvs - 0.5,
                                                     50,
                                                     0.,
                                                     50.);  // to modify the size binning
         (local_modmes.PedsEvolution)->setAxisTitle("Apv Number", 1);
         (local_modmes.PedsEvolution)->setAxisTitle("Iteration Number", 2);
 
         // Noise histos
         hid = hidmanager.createHistoId("CMSubNoisePerStrip", "det", detid);
         local_modmes.CMSubNoisePerStrip = ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);
         (local_modmes.CMSubNoisePerStrip)->setAxisTitle("CMSubNoise (ADC) vs Strip Number", 1);
 
         hid = hidmanager.createHistoId("RawNoisePerStrip", "det", detid);
         local_modmes.RawNoisePerStrip = ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);
         (local_modmes.RawNoisePerStrip)->setAxisTitle("RawNoise(ADC) vs Strip Number", 1);
 
         hid = hidmanager.createHistoId("CMSubNoiseProfile", "det", detid);
         local_modmes.CMSubNoiseProfile = ibooker.bookProfile(hid, hid, nStrip, 0.5, nStrip + 0.5, 100, 0., 100.);
         (local_modmes.CMSubNoiseProfile)->setAxisTitle("Mean of CMSubNoise (ADC) vs Strip Number", 1);
 
         hid = hidmanager.createHistoId("RawNoiseProfile", "det", detid);
         local_modmes.RawNoiseProfile = ibooker.bookProfile(hid, hid, nStrip, 0.5, nStrip + 0.5, 100, 0., 100.);
         (local_modmes.RawNoiseProfile)->setAxisTitle("Mean of RawNoise (ADC) vs Strip Number", 1);
 
         hid = hidmanager.createHistoId("NoisyStrips", "det", detid);
         local_modmes.NoisyStrips = ibooker.book2D(hid, hid, nStrip, 0.5, nStrip + 0.5, 6, -0.5, 5.5);
         (local_modmes.NoisyStrips)->setAxisTitle("Strip Number", 1);
         (local_modmes.NoisyStrips)->setAxisTitle("Flag Value", 2);
 
         hid = hidmanager.createHistoId("NoisyStripDistribution", "det", detid);
         local_modmes.NoisyStripDistribution = ibooker.book1D(hid, hid, 11, -0.5, 10.5);
         (local_modmes.NoisyStripDistribution)->setAxisTitle("Flag Value", 1);
 
         // Common Mode histos
         hid = hidmanager.createHistoId("CMDistribution", "det", detid);
         local_modmes.CMDistribution = ibooker.book2D(hid, hid, napvs, -0.5, napvs - 0.5, 150, -15., 15.);
         (local_modmes.CMDistribution)->setAxisTitle("Common Mode (ADC) vs APV Number", 1);
 
         hid = hidmanager.createHistoId("CMSlopeDistribution", "det", detid);
         local_modmes.CMSlopeDistribution = ibooker.book2D(hid, hid, napvs, -0.5, napvs - 0.5, 100, -0.05, 0.05);
         (local_modmes.CMSlopeDistribution)->setAxisTitle("Common Mode Slope vs APV Number", 1);
       }
       // data from CondDB
       // append to PedMEs
       PedMEs.insert(std::make_pair(detid, local_modmes));
     }  // newDetId
   }
   edm::LogInfo("SiStripMonitorPedestals")
       << "SiStripMonitorPedestals::createMEs: Number of DETS used " << PedMEs.size();
 }

◆ dqmEndRun()

void SiStripMonitorPedestals::dqmEndRun	(	edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 451 of file SiStripMonitorPedestals.cc.

References conf_, dqmStore_, edm::ParameterSet::getParameter(), outPutFileName, dqm::implementation::DQMStore::save(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                       {
   bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
   if (outputMEsInRootFile) {
     std::string outPutFileName = conf_.getParameter<std::string>("OutPutFileName");
     dqmStore_->save(outPutFileName);
   }
 }

◆ endJob()

void SiStripMonitorPedestals::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 461 of file SiStripMonitorPedestals.cc.

                                          {
   edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals::EndJob: "
                                           << " Finishing!! ";
 }

◆ fillCondDBMEs()

void SiStripMonitorPedestals::fillCondDBMEs ( edm::EventSetup const & eSetup )

private

Definition at line 495 of file SiStripMonitorPedestals.cc.

References SiStripMonitorPedestals::ModMEs::BadStripsDB, SiStripMonitorPedestals::ModMEs::CMSubNoisePerStripDB, detcabling, MillePedeFileConverter_cfg::e, HcalObjRepresent::Fill(), edm::EventSetup::getData(), mps_fire::i, SiStripDetCabling::nApvPairs(), hgcalLayerClustersEE_cfi::noises, noiseToken_, pedestalToken_, PedMEs, SiStripMonitorPedestals::ModMEs::PedsPerStripDB, quality, and qualityToken_.

Referenced by bookHistograms().

                                                                        {
   const auto &pedestals = eSetup.getData(pedestalToken_);
   const auto &noises = eSetup.getData(noiseToken_);
   const auto &quality = eSetup.getData(qualityToken_);
 
   for (std::map<uint32_t, ModMEs>::const_iterator i = PedMEs.begin(); i != PedMEs.end(); i++) {
     uint32_t detid = i->first;
     ModMEs local_modmes = i->second;
     edm::LogInfo("SiStripMonitorPedestals") << " SiStripMonitorPedestals::analyze: "
                                             << " Get Ped/Noise/Bad Strips from CondDb for DetId " << detid;
     int nStrip = detcabling->nApvPairs(detid) * 256;
     // Get range of pedestal and noise for the detid
     SiStripNoises::Range noiseRange = noises.getRange(detid);
     SiStripPedestals::Range pedRange = pedestals.getRange(detid);
     SiStripQuality::Range qualityRange = quality.getRange(detid);
 
     for (int istrip = 0; istrip < nStrip; ++istrip) {
       try {
         // Fill Pedestals
         (local_modmes.PedsPerStripDB)->Fill(istrip + 1, pedestals.getPed(istrip, pedRange));
       } catch (cms::Exception &e) {
         edm::LogError("SiStripMonitorPedestals") << "[SiStripMonitorPedestals::analyze]  cms::Exception accessing "
                                                     "SiStripPedestalsService_.getPedestal("
                                                  << detid << "," << istrip << ") :  "
                                                  << " " << e.what();
       }
       try {
         // Fill Noises
         (local_modmes.CMSubNoisePerStripDB)->Fill(istrip + 1, noises.getNoise(istrip, noiseRange));
 
       } catch (cms::Exception &e) {
         edm::LogError("SiStripMonitorPedestals") << "[SiStripMonitorPedestals::analyze]  cms::Exception accessing "
                                                     "SiStripNoiseService_.getNoise("
                                                  << detid << "," << istrip << ") :  "
                                                  << " " << e.what();
       }
       try {
         // Fill BadStripsNoise
         (local_modmes.BadStripsDB)->Fill(istrip + 1, quality.IsStripBad(qualityRange, istrip) ? 1. : 0.);
 
       } catch (cms::Exception &e) {
         edm::LogError("SiStripMonitorPedestals") << "[SiStripMonitorPedestals::analyze]  cms::Exception accessing "
                                                     "SiStripNoiseService_.getDisable("
                                                  << detid << "," << istrip << ") :  "
                                                  << " " << e.what();
       }
     }  // close istrip loop
   }
 }

◆ resetMEs()

void SiStripMonitorPedestals::resetMEs ( uint32_t idet )

private

Definition at line 468 of file SiStripMonitorPedestals.cc.

Referenced by createMEs().

                                                     {
   std::map<uint32_t, ModMEs>::iterator pos = PedMEs.find(idet);
   if (pos != PedMEs.end()) {
     ModMEs mod_me = pos->second;
 
     if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3) {
       mod_me.PedsPerStripDB->Reset();
       mod_me.CMSubNoisePerStripDB->Reset();
       mod_me.BadStripsDB->Reset();
     }
     if (runTypeFlag_ == RunMode2 || runTypeFlag_ == RunMode3) {
       mod_me.PedsPerStrip->Reset();
       mod_me.PedsDistribution->Reset();
       mod_me.PedsEvolution->Reset();
       mod_me.CMSubNoisePerStrip->Reset();
       mod_me.RawNoisePerStrip->Reset();
       mod_me.CMSubNoiseProfile->Reset();
       mod_me.RawNoiseProfile->Reset();
       mod_me.NoisyStrips->Reset();
       mod_me.CMDistribution->Reset();
       mod_me.CMSlopeDistribution->Reset();
     }
   }
 }

Member Data Documentation

◆ analyzed

bool SiStripMonitorPedestals::analyzed

private

Definition at line 94 of file SiStripMonitorPedestals.h.

◆ apvFactory_

ApvAnalysisFactory* SiStripMonitorPedestals::apvFactory_

private

Definition at line 103 of file SiStripMonitorPedestals.h.

Referenced by analyze(), createMEs(), and ~SiStripMonitorPedestals().

◆ conf_

edm::ParameterSet SiStripMonitorPedestals::conf_

private

Definition at line 90 of file SiStripMonitorPedestals.h.

Referenced by dqmEndRun(), and SiStripMonitorPedestals().

◆ detcabling

const SiStripDetCabling* SiStripMonitorPedestals::detcabling

private

Definition at line 92 of file SiStripMonitorPedestals.h.

Referenced by analyze(), bookHistograms(), createMEs(), and fillCondDBMEs().

◆ detCablingToken_

edm::ESGetToken<SiStripDetCabling, SiStripDetCablingRcd> SiStripMonitorPedestals::detCablingToken_

private

Definition at line 111 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and bookHistograms().

◆ detCablingWatcher_

edm::ESWatcher<SiStripDetCablingRcd> SiStripMonitorPedestals::detCablingWatcher_

private

Definition at line 112 of file SiStripMonitorPedestals.h.

Referenced by bookHistograms().

◆ digiToken_

edm::EDGetTokenT<edm::DetSetVector<SiStripRawDigi> > SiStripMonitorPedestals::digiToken_

private

Definition at line 62 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and SiStripMonitorPedestals().

◆ dqmStore_

DQMStore* SiStripMonitorPedestals::dqmStore_

private

Definition at line 89 of file SiStripMonitorPedestals.h.

Referenced by dqmEndRun().

◆ firstEvent

bool SiStripMonitorPedestals::firstEvent

private

Definition at line 95 of file SiStripMonitorPedestals.h.

Referenced by analyze(), looper.Looper::loop(), and Printer.Printer::process().

◆ nEvInit_

int16_t SiStripMonitorPedestals::nEvInit_

private

Definition at line 101 of file SiStripMonitorPedestals.h.

◆ nEvTot_

int16_t SiStripMonitorPedestals::nEvTot_

private

Definition at line 100 of file SiStripMonitorPedestals.h.

Referenced by analyze().

◆ nEvUpdate_

int16_t SiStripMonitorPedestals::nEvUpdate_

private

Definition at line 98 of file SiStripMonitorPedestals.h.

◆ nIteration_

int SiStripMonitorPedestals::nIteration_

private

Definition at line 102 of file SiStripMonitorPedestals.h.

Referenced by analyze().

◆ noiseToken_

edm::ESGetToken<SiStripNoises, SiStripNoisesRcd> SiStripMonitorPedestals::noiseToken_

private

Definition at line 114 of file SiStripMonitorPedestals.h.

Referenced by fillCondDBMEs(), and SiStripMonitorPedestals().

◆ NumCMstripsInGroup_

int SiStripMonitorPedestals::NumCMstripsInGroup_

private

Definition at line 106 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and SiStripMonitorPedestals().

◆ outPutFileName

std::string SiStripMonitorPedestals::outPutFileName

private

Definition at line 108 of file SiStripMonitorPedestals.h.

Referenced by dqmEndRun().

◆ pedestalToken_

edm::ESGetToken<SiStripPedestals, SiStripPedestalsRcd> SiStripMonitorPedestals::pedestalToken_

private

Definition at line 113 of file SiStripMonitorPedestals.h.

Referenced by fillCondDBMEs(), and SiStripMonitorPedestals().

◆ PedMEs

std::map<uint32_t, ModMEs> SiStripMonitorPedestals::PedMEs

private

Definition at line 91 of file SiStripMonitorPedestals.h.

Referenced by analyze(), createMEs(), fillCondDBMEs(), and resetMEs().

◆ pedsPSet_

edm::ParameterSet SiStripMonitorPedestals::pedsPSet_

private

Definition at line 93 of file SiStripMonitorPedestals.h.

Referenced by createMEs(), and SiStripMonitorPedestals().

◆ qualityToken_

edm::ESGetToken<SiStripQuality, SiStripQualityRcd> SiStripMonitorPedestals::qualityToken_

private

Definition at line 115 of file SiStripMonitorPedestals.h.

Referenced by fillCondDBMEs(), and SiStripMonitorPedestals().

◆ RunMode1

const std::string SiStripMonitorPedestals::RunMode1 = "ConDBPlotsOnly"

staticprivate

Definition at line 117 of file SiStripMonitorPedestals.h.

Referenced by bookHistograms(), createMEs(), resetMEs(), and SiStripMonitorPedestals().

◆ RunMode2

const std::string SiStripMonitorPedestals::RunMode2 = "CalculatedPlotsOnly"

staticprivate

Definition at line 118 of file SiStripMonitorPedestals.h.

Referenced by analyze(), createMEs(), and resetMEs().

◆ RunMode3

const std::string SiStripMonitorPedestals::RunMode3 = "AllPlots"

staticprivate

Definition at line 119 of file SiStripMonitorPedestals.h.

Referenced by analyze(), bookHistograms(), createMEs(), resetMEs(), and SiStripMonitorPedestals().

◆ runTypeFlag_

std::string SiStripMonitorPedestals::runTypeFlag_

private

Definition at line 107 of file SiStripMonitorPedestals.h.

Referenced by analyze(), bookHistograms(), createMEs(), resetMEs(), and SiStripMonitorPedestals().

◆ signalCutPeds_

int16_t SiStripMonitorPedestals::signalCutPeds_

private

Definition at line 99 of file SiStripMonitorPedestals.h.

◆ theEventInitNumber_

int SiStripMonitorPedestals::theEventInitNumber_

private

Definition at line 104 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and SiStripMonitorPedestals().

◆ theEventIterNumber_

int SiStripMonitorPedestals::theEventIterNumber_

private

Definition at line 105 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and SiStripMonitorPedestals().

◆ tTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> SiStripMonitorPedestals::tTopoToken_

private

Definition at line 110 of file SiStripMonitorPedestals.h.

Referenced by createMEs().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SiStripMonitorPedestals()

◆ ~SiStripMonitorPedestals()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ createMEs()

◆ dqmEndRun()

◆ endJob()

◆ fillCondDBMEs()

◆ resetMEs()

Member Data Documentation

◆ analyzed

◆ apvFactory_

◆ conf_

◆ detcabling

◆ detCablingToken_

◆ detCablingWatcher_

◆ digiToken_

◆ dqmStore_

◆ firstEvent

◆ nEvInit_

◆ nEvTot_

◆ nEvUpdate_

◆ nIteration_

◆ noiseToken_

◆ NumCMstripsInGroup_

◆ outPutFileName

◆ pedestalToken_

◆ PedMEs

◆ pedsPSet_

◆ qualityToken_

◆ RunMode1

◆ RunMode2

◆ RunMode3

◆ runTypeFlag_

◆ signalCutPeds_

◆ theEventInitNumber_

◆ theEventIterNumber_

◆ tTopoToken_