#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 DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

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

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

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::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

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

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)

	~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

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, 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

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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_

edm::ESHandle< SiStripDetCabling >	detcabling

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

DQMStore *	dqmStore_

bool	firstEvent

unsigned long long	m_cacheID_

int16_t	nEvInit_

int16_t	nEvTot_

int16_t	nEvUpdate_

int	nIteration_

int	NumCMstripsInGroup_

std::string	outPutFileName

std::map< uint32_t, ModMEs >	PedMEs

edm::ParameterSet	pedsPSet_

std::string	runTypeFlag_

int16_t	signalCutPeds_

int	theEventInitNumber_

int	theEventIterNumber_

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 DQMEDAnalyzer
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
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 edm::ProducerBase
ProducesCollector	producesCollector ()

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

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

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

Usage: <usage>

Definition at line 54 of file SiStripMonitorPedestals.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 53 of file SiStripMonitorPedestals.cc.

References conf_, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::digiProducer, digiToken_, edm::ParameterSet::getParameter(), HLT_2018_cff::InputTag, NumCMstripsInGroup_, pedsPSet_, 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),
       m_cacheID_(0) {
   // 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");
 }

SiStripMonitorPedestals::~SiStripMonitorPedestals ( )

override

Definition at line 80 of file SiStripMonitorPedestals.cc.

References apvFactory_.

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

Member Function Documentation

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

override

Definition at line 273 of file SiStripMonitorPedestals.cc.

                                                                                          {
   edm::LogInfo("SiStripMonitorPedestals")
       << "SiStripMonitorPedestals::analyze: Run " << iEvent.id().run() << " Event " << iEvent.id().event();
 
   eSetup.get<SiStripDetCablingRcd>().get(detcabling);
 
   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;
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 91 of file SiStripMonitorPedestals.cc.

References createMEs(), detcabling, fillCondDBMEs(), edm::EventSetup::get(), m_cacheID_, RunMode1, RunMode3, and runTypeFlag_.

                                                                           {
   unsigned long long cacheID = eSetup.get<SiStripDetCablingRcd>().cacheIdentifier();
   if (m_cacheID_ != cacheID) {
     m_cacheID_ = cacheID;
     eSetup.get<SiStripDetCablingRcd>().get(detcabling);
     edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals::bookHistograms: "
                                             << " Creating MEs for new Cabling ";
     createMEs(ibooker, eSetup);
   }
 
   if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3)
     fillCondDBMEs(eSetup);
 }

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

private

Definition at line 110 of file SiStripMonitorPedestals.cc.

Referenced by bookHistograms().

                                                                                          {
   // Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology *const tTopo = tTopoHandle.product();
 
   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();
 }

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 451 of file SiStripMonitorPedestals.cc.

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

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

void SiStripMonitorPedestals::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 462 of file SiStripMonitorPedestals.cc.

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

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

private

Definition at line 496 of file SiStripMonitorPedestals.cc.

References SiStripMonitorPedestals::ModMEs::BadStripsDB, SiStripMonitorPedestals::ModMEs::CMSubNoisePerStripDB, conf_, DEFINE_FWK_MODULE, detcabling, MillePedeFileConverter_cfg::e, HcalObjRepresent::Fill(), edm::EventSetup::get(), SiStripNoises::getNoise(), edm::ParameterSet::getParameter(), SiStripPedestals::getPed(), SiStripNoises::getRange(), SiStripPedestals::getRange(), SiStripBadStrip::getRange(), mps_fire::i, SiStripQuality::IsStripBad(), SiStripDetCabling::nApvPairs(), PedMEs, SiStripMonitorPedestals::ModMEs::PedsPerStripDB, AlCaHLTBitMon_QueryRunRegistry::string, and cms::Exception::what().

Referenced by bookHistograms().

                                                                        {
   // get Pedestal and Noise  ES handle
   edm::ESHandle<SiStripPedestals> pedestalHandle;
   edm::ESHandle<SiStripNoises> noiseHandle;
   edm::ESHandle<SiStripQuality> qualityHandle;
 
   eSetup.get<SiStripPedestalsRcd>().get(pedestalHandle);
   eSetup.get<SiStripNoisesRcd>().get(noiseHandle);
   std::string quality_label = conf_.getParameter<std::string>("StripQualityLabel");
   eSetup.get<SiStripQualityRcd>().get(quality_label, qualityHandle);
 
   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 = noiseHandle->getRange(detid);
     SiStripPedestals::Range pedRange = pedestalHandle->getRange(detid);
     SiStripQuality::Range qualityRange = qualityHandle->getRange(detid);
 
     for (int istrip = 0; istrip < nStrip; ++istrip) {
       try {
         // Fill Pedestals
         (local_modmes.PedsPerStripDB)->Fill(istrip + 1, pedestalHandle->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, noiseHandle->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, qualityHandle->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
   }
 }

void SiStripMonitorPedestals::resetMEs ( uint32_t idet )

private

Definition at line 469 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();
     }
   }
 }