SiStripBaselineValidator.cc

Go to the documentation of this file.

// Original Author:  Ivan Amos Cali
//         Created:  Mon Jul 28 14:10:52 CEST 2008
//
//
 
// system include files
#include <memory>
#include <iostream>
 
// user include files
#include "DQM/SiStripMonitorDigi/interface/SiStripBaselineValidator.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DataFormats/Common/interface/DetSet.h"
#include "DQMServices/Core/interface/DQMStore.h"
 
//ROOT inclusion
#include "TH1F.h"
#include "TH2F.h"
#include "TString.h"
#include <cassert>
#include <fstream>
 
class TFile;
 
using namespace edm;
using namespace std;
 
SiStripBaselineValidator::SiStripBaselineValidator(const edm::ParameterSet& conf) {
  srcProcessedRawDigi_ = conf.getParameter<edm::InputTag>("srcProcessedRawDigi");
  moduleRawDigiToken_ =
      consumes<edm::DetSetVector<SiStripDigi> >(conf.getParameter<edm::InputTag>("srcProcessedRawDigi"));
}
 
SiStripBaselineValidator::~SiStripBaselineValidator() {}
 
void SiStripBaselineValidator::bookHistograms(DQMStore::IBooker& ibooker,
                                              const edm::Run& run,
                                              const edm::EventSetup& es) {
  ibooker.setCurrentFolder("SiStrip/BaselineValidator");
 
  h1NumbadAPVsRes_ = ibooker.book1D("ResAPVs", ";#ResAPVs", 100, 1.0, 10001);
 
  h1ADC_vs_strip_ = ibooker.book2D("ADCvsAPVs", ";ADCvsAPVs", 768, -0.5, 767.5, 1023, -0.5, 1022.5);
 
  return;
}
 
void SiStripBaselineValidator::analyze(const edm::Event& e, const edm::EventSetup& es) {
  edm::Handle<edm::DetSetVector<SiStripDigi> > moduleRawDigi;
  e.getByToken(moduleRawDigiToken_, moduleRawDigi);
  edm::DetSetVector<SiStripDigi>::const_iterator itRawDigis = moduleRawDigi->begin();
 
  int NumResAPVs = 0;
  for (; itRawDigis != moduleRawDigi->end(); ++itRawDigis) {  
 
    edm::DetSet<SiStripDigi>::const_iterator itRaw = itRawDigis->begin();
    int strip = 0, totStripAPV = 0, apv = 0, prevapv = itRaw->strip() / 128;
 
    for (; itRaw != itRawDigis->end(); ++itRaw) {  
 
      strip = itRaw->strip();
      apv = strip / 128;
      float adc = itRaw->adc();
      h1ADC_vs_strip_->Fill(strip, adc);  
 
      if (prevapv != apv) {
        if (totStripAPV > 64) {
          NumResAPVs++;
        }
        prevapv = apv;
        totStripAPV = 0;
      }
      if (adc > 0)
        ++totStripAPV;
 
    }  
    if (totStripAPV > 64) {
      NumResAPVs++;
    }
 
  }  
 
 
  h1NumbadAPVsRes_->Fill(NumResAPVs);  
 
}  
 
//define this as a plug-in
DEFINE_FWK_MODULE(SiStripBaselineValidator);