CMS 3D CMS Logo

Public Member Functions | Private Member Functions | Private Attributes

SiStripNoisesDQMService Class Reference

#include <SiStripNoisesDQMService.h>

Inheritance diagram for SiStripNoisesDQMService:
SiStripBaseServiceFromDQM< SiStripNoises > SiStripCondObjBuilderBase< SiStripNoises >

List of all members.

Public Member Functions

void getObj (SiStripNoises *&obj)
 SiStripNoisesDQMService (const edm::ParameterSet &, const edm::ActivityRegistry &)
 ~SiStripNoisesDQMService ()

Private Member Functions

void readNoises ()

Private Attributes

edm::FileInPath fp_
edm::ParameterSet iConfig_

Detailed Description

Author:
M. De Mattia, S. Dutta, D. Giordano to read modules flagged by the DQM as bad and write in the database.

Definition at line 24 of file SiStripNoisesDQMService.h.


Constructor & Destructor Documentation

SiStripNoisesDQMService::SiStripNoisesDQMService ( const edm::ParameterSet iConfig,
const edm::ActivityRegistry aReg 
) [explicit]

Definition at line 13 of file SiStripNoisesDQMService.cc.

References SiStripCondObjBuilderBase< SiStripNoises >::obj_.

                                                                                                              :
  // SiStripCondObjBuilderBase<SiStripNoises>::SiStripCondObjBuilderBase(iConfig),
  SiStripBaseServiceFromDQM<SiStripNoises>::SiStripBaseServiceFromDQM(iConfig),
  iConfig_(iConfig),
  fp_(iConfig.getUntrackedParameter<edm::FileInPath>("file",edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat")))
{
  obj_ = 0;
  edm::LogInfo("SiStripNoisesDQMService") <<  "[SiStripNoisesDQMService::SiStripNoisesDQMService]";
}
SiStripNoisesDQMService::~SiStripNoisesDQMService ( )

Definition at line 23 of file SiStripNoisesDQMService.cc.

{
  edm::LogInfo("SiStripNoisesDQMService") <<  "[SiStripNoisesDQMService::~SiStripNoisesDQMService]";
}

Member Function Documentation

void SiStripNoisesDQMService::getObj ( SiStripNoises *&  obj) [inline, virtual]

Returns the CondObj

Reimplemented from SiStripCondObjBuilderBase< SiStripNoises >.

Definition at line 31 of file SiStripNoisesDQMService.h.

References SiStripCondObjBuilderBase< SiStripNoises >::obj_, and readNoises().

  {
    readNoises(); obj=obj_;
  }
void SiStripNoisesDQMService::readNoises ( ) [private]

Definition at line 28 of file SiStripNoisesDQMService.cc.

References DQMStore::cd(), gather_cfg::cout, SiStripBaseServiceFromDQM< SiStripNoises >::dqmStore_, fp_, edm::FileInPath::fullPath(), DQMStore::getAllContents(), MonitorElement::getTH1F(), edm::ParameterSet::getUntrackedParameter(), timingPdfMaker::histo, iConfig_, j, SiStripCondObjBuilderBase< SiStripNoises >::obj_, SiStripBaseServiceFromDQM< SiStripNoises >::openRequestedFile(), SiStripNoises::put(), matplotRender::reader, SiStripNoises::setData(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by getObj().

{
  std::cout << "SiStripNoisesDQMService::readNoises" << std::endl;

  openRequestedFile();

  std::cout << "[readBadComponents]: opened requested file" << std::endl;

  obj_= new SiStripNoises;

  SiStripDetInfoFileReader reader(fp_.fullPath());

  // dqmStore_->cd(iConfig_.getUntrackedParameter<std::string>("ME_DIR"));
  dqmStore_->cd();

  uint32_t stripsPerApv = 128;

  // Get the full list of monitoring elements
  // const std::vector<MonitorElement*>& MEs = dqmStore_->getAllContents(iConfig_.getUntrackedParameter<std::string>("ME_DIR","DQMData"));

  // Take a copy of the vector
  std::vector<MonitorElement*> MEs = dqmStore_->getAllContents(iConfig_.getUntrackedParameter<std::string>("ME_DIR","DQMData"));
  // Remove all but the MEs we are using
  std::vector<MonitorElement*>::iterator newEnd = remove_if(MEs.begin(), MEs.end(), StringNotMatch("CMSubNoisePerStrip__det__"));
  MEs.erase(newEnd, MEs.end());

  // The histograms are one per DetId, loop on all the DetIds and extract the corresponding histogram
  const std::map<uint32_t, SiStripDetInfoFileReader::DetInfo> DetInfos  = reader.getAllData();
  for(std::map<uint32_t, SiStripDetInfoFileReader::DetInfo>::const_iterator it = DetInfos.begin(); it != DetInfos.end(); ++it) {


    SiStripNoises::InputVector theSiStripVector;

    // Take the path for each DetId and build the complete path + histogram name


    // MonitorElement * mE = getModuleHistogram(it->first, "PedsPerStrip");


    MonitorElement * mE = 0;
    std::string MEname("CMSubNoisePerStrip__det__"+boost::lexical_cast<string>(it->first));
    for( std::vector<MonitorElement*>::const_iterator MEit = MEs.begin();
         MEit != MEs.end(); ++MEit ) {
      if( (*MEit)->getName() == MEname ) {
        mE = *MEit;
        break;
      }
    }

    // find( MEs.begin(), MEs.end(), "PedsPerStrip__det__"+boost::lexical_cast<string>(it->first), findMEbyName() );
    // MonitorElement * mE = *(find( MEs.begin(), MEs.end(), findMEbyName("PedsPerStrip__det__"+boost::lexical_cast<string>(it->first)) ));

    if( mE != 0 ) {
      TH1F* histo = mE->getTH1F();

      if( histo != 0 ) {

        // Read the noise from the histograms
        uint32_t nBinsX = histo->GetXaxis()->GetNbins();

        if( nBinsX != stripsPerApv*(it->second.nApvs) ) {
          std::cout << "ERROR: number of bin = " << nBinsX << " != number of strips = " << stripsPerApv*(it->second.nApvs) << std::endl;
        }

        // std::cout << "Bin 0 = " << histo->GetBinContent(0) << std::endl;

        // TH1 bins start from 1, 0 is the underflow, nBinsX+1 the overflow.
        for( uint32_t iBin = 1; iBin <= nBinsX; ++iBin ) {
          // encode the pedestal value and put it in the vector (push_back)
          obj_->setData( histo->GetBinContent(iBin), theSiStripVector );
        }
      }
      else {
        std::cout << "ERROR: histo = " << histo << std::endl;
      }
    }
    else {
      std::cout << "ERROR: ME = " << mE << std::endl;
    }
    // If the ME was absent fill the vector with 50 (we want a high noise to avoid these modules being considered good by mistake)
    if( theSiStripVector.empty() ) {
      for(unsigned short j=0; j<128*it->second.nApvs; ++j){
        obj_->setData(50, theSiStripVector);
      }
    }

    if ( ! obj_->put(it->first, theSiStripVector) )
      edm::LogError("SiStripNoisesFakeESSource::produce ")<<" detid already exists"<<std::endl;
  }
  dqmStore_->cd();
}

Member Data Documentation

Definition at line 41 of file SiStripNoisesDQMService.h.

Referenced by readNoises().

Reimplemented from SiStripBaseServiceFromDQM< SiStripNoises >.

Definition at line 40 of file SiStripNoisesDQMService.h.

Referenced by readNoises().