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"

/*#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackCandidate/interface/TrackCandidateCollection.h" 
#include "DataFormats/TrackCandidate/interface/TrackCandidate.h" 
#include "DataFormats/TrackReco/interface/Track.h" 
#include "DataFormats/TrackReco/interface/TrackFwd.h" 
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.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" );
 // hiSelectedTracks =  conf.getParameter<edm::InputTag>( "hiSelectedTracks" );
  createOutputFile_ = conf.getUntrackedParameter<bool>("saveFile",false);
  outputFile_   = conf.getParameter<std::string>("outputFile");
  dbe = &*edm::Service<DQMStore>();
  moduleRawDigiToken_ = consumes<edm::DetSetVector<SiStripRawDigi> >(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);
  //dbe->tag(h1NumbadAPVsRes_->getFullname(),1);
  ibooker.tag(h1NumbadAPVsRes_,1);
  
  h1ADC_vs_strip_ = ibooker.book2D("ADCvsAPVs",";ADCvsAPVs", 768,-0.5,767.5,  1023, -0.5, 1022.5);
  ibooker.tag(h1ADC_vs_strip_,2);
  
  return;
}

// ------------ method called once each job just before starting event loop  ------------
void SiStripBaselineValidator::beginJob()
{

}

void SiStripBaselineValidator::analyze(const edm::Event& e, const edm::EventSetup& es)
{
  edm::Handle< edm::DetSetVector<SiStripRawDigi> > moduleRawDigi;
  e.getByToken( moduleRawDigiToken_, moduleRawDigi );
  edm::DetSetVector<SiStripRawDigi>::const_iterator itRawDigis = moduleRawDigi->begin();
 
   int NumResAPVs=0;
   for (; itRawDigis != moduleRawDigi->end(); ++itRawDigis) {   
     

     edm::DetSet<SiStripRawDigi>::const_iterator itRaw = itRawDigis->begin(); 
     int strip =0, totADC=0;

     for(;itRaw != itRawDigis->end(); ++itRaw, ++strip){  
       
       float adc = itRaw->adc();
       h1ADC_vs_strip_->Fill(strip,adc); 



       totADC+= adc;
       
       if(strip%127 ==0){
     if(totADC!= 0){
       totADC =0;
       
       NumResAPVs++;

     }
       }
       
     } 
   
     
   }  


   h1NumbadAPVsRes_->Fill(NumResAPVs); 



} 



// ------------ method called once each job just after ending the event loop  ------------
void 
SiStripBaselineValidator::endJob() {

    if (!outputFile_.empty() && createOutputFile_) {
       dbe->save(outputFile_);
    }  
}

//define this as a plug-in
DEFINE_FWK_MODULE(SiStripBaselineValidator);