SiStripMonitorPedestals Class Reference

#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	endJob () override
void	endRun (edm::Run const &run, edm::EventSetup const &eSetup) override
	SiStripMonitorPedestals (const edm::ParameterSet &)
	~SiStripMonitorPedestals () override
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

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"

Detailed Description

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

Usage: <usage>

Definition at line 56 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(), 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 81 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 271 of file SiStripMonitorPedestals.cc.

References apvFactory_, SiStripMonitorPedestals::ModMEs::CMDistribution, SiStripMonitorPedestals::ModMEs::CMSlopeDistribution, SiStripMonitorPedestals::ModMEs::CMSubNoisePerStrip, SiStripMonitorPedestals::ModMEs::CMSubNoiseProfile, gather_cfg::cout, detcabling, digiToken_, edm::EventID::event(), HcalObjRepresent::Fill(), firstEvent, edm::EventSetup::get(), edm::Event::getByToken(), ApvAnalysisFactory::getCommonMode(), ApvAnalysisFactory::getCommonModeSlope(), SiStripDetCabling::getConnections(), ApvAnalysisFactory::getMask(), ApvAnalysisFactory::getNoise(), ApvAnalysisFactory::getPedestal(), ApvAnalysisFactory::getRawNoise(), mps_fire::i, edm::EventBase::id(), createfilelist::int, gen::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, edm::EventID::run(), RunMode2, RunMode3, runTypeFlag_, groupFilesInBlocks::temp, theEventInitNumber_, theEventIterNumber_, tmp, and ApvAnalysisFactory::update().

{

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

override

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.

References SiStripDetCabling::addActiveDetectorsRawIds(), apvFactory_, SiStripMonitorPedestals::ModMEs::BadStripsDB, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), SiStripMonitorPedestals::ModMEs::CMDistribution, SiStripMonitorPedestals::ModMEs::CMSlopeDistribution, SiStripMonitorPedestals::ModMEs::CMSubNoisePerStrip, SiStripMonitorPedestals::ModMEs::CMSubNoisePerStripDB, SiStripMonitorPedestals::ModMEs::CMSubNoiseProfile, SiStripHistoId::createHistoId(), detcabling, edm::EventSetup::get(), ApvAnalysisFactory::instantiateApvs(), SiStripDetCabling::nApvPairs(), SiStripMonitorPedestals::ModMEs::NoisyStripDistribution, SiStripMonitorPedestals::ModMEs::NoisyStrips, PedMEs, SiStripMonitorPedestals::ModMEs::PedsDistribution, SiStripMonitorPedestals::ModMEs::PedsEvolution, SiStripMonitorPedestals::ModMEs::PedsPerStrip, SiStripMonitorPedestals::ModMEs::PedsPerStripDB, pedsPSet_, edm::ESHandle< T >::product(), SiStripMonitorPedestals::ModMEs::RawNoisePerStrip, SiStripMonitorPedestals::ModMEs::RawNoiseProfile, resetMEs(), RunMode1, RunMode2, RunMode3, runTypeFlag_, SiStripFolderOrganizer::setDetectorFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and DQMStore::IBooker::tag().

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 
    ibooker.tag(local_modmes.PedsPerStripDB, detid);
    (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);
    ibooker.tag(local_modmes.CMSubNoisePerStripDB, detid);
    (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);
    ibooker.tag(local_modmes.BadStripsDB, detid);
    (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 
    ibooker.tag(local_modmes.PedsPerStrip, detid);
    (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 
    ibooker.tag(local_modmes.PedsDistribution, detid);
    (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 
    ibooker.tag(local_modmes.PedsEvolution, detid);
    (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);
    ibooker.tag(local_modmes.CMSubNoisePerStrip, detid);
    (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);
    ibooker.tag(local_modmes.RawNoisePerStrip, detid);
    (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.);
    ibooker.tag(local_modmes.CMSubNoiseProfile, detid);
    (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.);
    ibooker.tag(local_modmes.RawNoiseProfile, detid);
    (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);
    ibooker.tag(local_modmes.NoisyStrips, detid);
    (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);
    ibooker.tag(local_modmes.NoisyStripDistribution, detid);
    (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.); 
    ibooker.tag(local_modmes.CMDistribution, detid);
    (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); 
    ibooker.tag(local_modmes.CMSlopeDistribution, detid);
    (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::endJob ( void  )

override

Definition at line 467 of file SiStripMonitorPedestals.cc.

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

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

override

Definition at line 456 of file SiStripMonitorPedestals.cc.

References conf_, dqmStore_, edm::ParameterSet::getParameter(), outPutFileName, 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::fillCondDBMEs ( edm::EventSetup const &  eSetup )

private

Definition at line 501 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 474 of file SiStripMonitorPedestals.cc.

References SiStripMonitorPedestals::ModMEs::BadStripsDB, SiStripMonitorPedestals::ModMEs::CMDistribution, SiStripMonitorPedestals::ModMEs::CMSlopeDistribution, SiStripMonitorPedestals::ModMEs::CMSubNoisePerStrip, SiStripMonitorPedestals::ModMEs::CMSubNoisePerStripDB, SiStripMonitorPedestals::ModMEs::CMSubNoiseProfile, SiStripMonitorPedestals::ModMEs::NoisyStrips, PedMEs, SiStripMonitorPedestals::ModMEs::PedsDistribution, SiStripMonitorPedestals::ModMEs::PedsEvolution, SiStripMonitorPedestals::ModMEs::PedsPerStrip, SiStripMonitorPedestals::ModMEs::PedsPerStripDB, SiStripMonitorPedestals::ModMEs::RawNoisePerStrip, SiStripMonitorPedestals::ModMEs::RawNoiseProfile, MonitorElement::Reset(), RunMode1, RunMode2, RunMode3, and runTypeFlag_.

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

bool SiStripMonitorPedestals::analyzed

private

Definition at line 100 of file SiStripMonitorPedestals.h.

ApvAnalysisFactory* SiStripMonitorPedestals::apvFactory_

private

Definition at line 109 of file SiStripMonitorPedestals.h.

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

edm::ParameterSet SiStripMonitorPedestals::conf_

private

Definition at line 96 of file SiStripMonitorPedestals.h.

Referenced by endRun(), fillCondDBMEs(), and SiStripMonitorPedestals().

edm::ESHandle<SiStripDetCabling> SiStripMonitorPedestals::detcabling

private

Definition at line 98 of file SiStripMonitorPedestals.h.

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

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

private

Definition at line 67 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and SiStripMonitorPedestals().

DQMStore* SiStripMonitorPedestals::dqmStore_

private

Definition at line 95 of file SiStripMonitorPedestals.h.

Referenced by endRun().

bool SiStripMonitorPedestals::firstEvent

private

Definition at line 101 of file SiStripMonitorPedestals.h.

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

unsigned long long SiStripMonitorPedestals::m_cacheID_

private

Definition at line 115 of file SiStripMonitorPedestals.h.

Referenced by bookHistograms().

int16_t SiStripMonitorPedestals::nEvInit_

private

Definition at line 107 of file SiStripMonitorPedestals.h.

int16_t SiStripMonitorPedestals::nEvTot_

private

Definition at line 106 of file SiStripMonitorPedestals.h.

Referenced by analyze().

int16_t SiStripMonitorPedestals::nEvUpdate_

private

Definition at line 104 of file SiStripMonitorPedestals.h.

int SiStripMonitorPedestals::nIteration_

private

Definition at line 108 of file SiStripMonitorPedestals.h.

Referenced by analyze().

int SiStripMonitorPedestals::NumCMstripsInGroup_

private

Definition at line 112 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and SiStripMonitorPedestals().

std::string SiStripMonitorPedestals::outPutFileName

private

Definition at line 114 of file SiStripMonitorPedestals.h.

Referenced by endRun().

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

private

Definition at line 97 of file SiStripMonitorPedestals.h.

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

edm::ParameterSet SiStripMonitorPedestals::pedsPSet_

private

Definition at line 99 of file SiStripMonitorPedestals.h.

Referenced by createMEs(), and SiStripMonitorPedestals().

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

staticprivate

Definition at line 117 of file SiStripMonitorPedestals.h.

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

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

staticprivate

Definition at line 118 of file SiStripMonitorPedestals.h.

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

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

staticprivate

Definition at line 119 of file SiStripMonitorPedestals.h.

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

std::string SiStripMonitorPedestals::runTypeFlag_

private

Definition at line 113 of file SiStripMonitorPedestals.h.

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

int16_t SiStripMonitorPedestals::signalCutPeds_

private

Definition at line 105 of file SiStripMonitorPedestals.h.

int SiStripMonitorPedestals::theEventInitNumber_

private

Definition at line 110 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and SiStripMonitorPedestals().

int SiStripMonitorPedestals::theEventIterNumber_

private

Definition at line 111 of file SiStripMonitorPedestals.h.

Referenced by analyze(), and SiStripMonitorPedestals().