DQMLumiMonitor.cc

Go to the documentation of this file.

/*
 * \file DQMLumiMonitor.cc
 * \author S. Dutta
 * Last Update:
 *
 * Description: Pixel Luminosity Monitoring 
 *
*/
#include "DQMServices/Components/plugins/DQMLumiMonitor.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"
#include "DataFormats/Luminosity/interface/LumiDetails.h"
#include "DataFormats/Luminosity/interface/LumiSummary.h"
#include "TPRegexp.h"

// -----------------------------
//  constructors and destructor
// -----------------------------

DQMLumiMonitor::DQMLumiMonitor(const edm::ParameterSet& ps) : parameters_(ps) {
  moduleName_ = parameters_.getParameter<std::string>("ModuleName");
  folderName_ = parameters_.getParameter<std::string>("FolderName");
  pixelClusterInputTag_ =
      consumes<edmNew::DetSetVector<SiPixelCluster> >(parameters_.getParameter<edm::InputTag>("PixelClusterInputTag"));
  lumiRecordName_ = consumes<LumiSummary, edm::InLumi>(parameters_.getParameter<std::string>("LumiRecordName"));

  nClusME_ = nullptr;
  nClusVsLSME_ = nullptr;
  intLumiVsLSME_ = nullptr;
  corrIntLumiAndClusVsLSME_ = nullptr;
}

DQMLumiMonitor::~DQMLumiMonitor() = default;

void DQMLumiMonitor::bookHistograms() {
  edm::ParameterSet ClusHistoPar = parameters_.getParameter<edm::ParameterSet>("TH1ClusPar");
  edm::ParameterSet LumiHistoPar = parameters_.getParameter<edm::ParameterSet>("TH1LumiPar");
  edm::ParameterSet LumiSecHistoPar = parameters_.getParameter<edm::ParameterSet>("TH1LSPar");

  std::string currentFolder = moduleName_ + "/" + folderName_;
  dbe_->setCurrentFolder(currentFolder);

  if (nClusME_ == nullptr)
    nClusME_ = dbe_->book1D("nPlxClus",
                            " Number of Pixel Clusters ",
                            ClusHistoPar.getParameter<int32_t>("Xbins"),
                            ClusHistoPar.getParameter<double>("Xmin"),
                            ClusHistoPar.getParameter<double>("Xmax"));
  else
    nClusME_->Reset();
  if (nClusVsLSME_ == nullptr)
    nClusVsLSME_ = dbe_->bookProfile("nClusVsLS",
                                     " Number of Pixel Cluster Vs LS number",
                                     LumiSecHistoPar.getParameter<int32_t>("Xbins"),
                                     LumiSecHistoPar.getParameter<double>("Xmin"),
                                     LumiSecHistoPar.getParameter<double>("Xmax"),
                                     0.0,
                                     0.0,
                                     "");
  else
    nClusVsLSME_->Reset();
  if (intLumiVsLSME_ == nullptr)
    intLumiVsLSME_ = dbe_->bookProfile("intLumiVsLS",
                                       " Integrated Luminosity Vs LS number",
                                       LumiSecHistoPar.getParameter<int32_t>("Xbins"),
                                       LumiSecHistoPar.getParameter<double>("Xmin"),
                                       LumiSecHistoPar.getParameter<double>("Xmax"),
                                       0.0,
                                       0.0,
                                       "");
  else
    intLumiVsLSME_->Reset();

  if (corrIntLumiAndClusVsLSME_ == nullptr)
    corrIntLumiAndClusVsLSME_ = dbe_->bookProfile2D("corrIntLumiAndClusVsLS",
                                                    " Correlation of nCluster and Integrated Luminosity Vs LS number",
                                                    LumiSecHistoPar.getParameter<int32_t>("Xbins"),
                                                    LumiSecHistoPar.getParameter<double>("Xmin"),
                                                    LumiSecHistoPar.getParameter<double>("Xmax"),
                                                    LumiHistoPar.getParameter<int32_t>("Xbins"),
                                                    LumiHistoPar.getParameter<double>("Xmin"),
                                                    LumiHistoPar.getParameter<double>("Xmax"),
                                                    0.0,
                                                    0.0);
  else
    corrIntLumiAndClusVsLSME_->Reset();
}
void DQMLumiMonitor::beginJob() {
  dbe_ = edm::Service<DQMStore>().operator->();
  intLumi_ = -1.0;
  nLumi_ = -1;
}

void DQMLumiMonitor::beginRun(edm::Run const& iRun, edm::EventSetup const& iSetup) { bookHistograms(); }
void DQMLumiMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) {
  //Access Pixel Clusters
  edm::Handle<edmNew::DetSetVector<SiPixelCluster> > siPixelClusters;
  iEvent.getByToken(pixelClusterInputTag_, siPixelClusters);

  if (!siPixelClusters.isValid()) {
    edm::LogError("PixelLumiMonotor") << "Could not find Cluster Collection ";
    return;
  }
  unsigned int nClusterPix = (*siPixelClusters).dataSize();
  nClusME_->Fill(nClusterPix);
  if (nLumi_ != -1)
    nClusVsLSME_->Fill(nLumi_, nClusterPix);
  if (intLumi_ != -1 || nLumi_ != -1)
    corrIntLumiAndClusVsLSME_->Fill(nLumi_, intLumi_, nClusterPix);
}

void DQMLumiMonitor::endLuminosityBlock(edm::LuminosityBlock const& lumiBlock, edm::EventSetup const& eSetup) {
  edm::LogInfo("PixelLumiMonotor") << " Run Number " << lumiBlock.run() << " Lumi Section Numnber "
                                   << lumiBlock.luminosityBlock();

  nLumi_ = lumiBlock.luminosityBlock();

  // Access Lumi Summary
  edm::Handle<LumiSummary> lumiSummary_;
  lumiBlock.getByToken(lumiRecordName_, lumiSummary_);
  if (lumiSummary_->isValid()) {
    intLumi_ = lumiSummary_->intgDelLumi();
    edm::LogInfo("PixelLumiMonotor") << " Luminosity in this Lumi Section " << intLumi_;
    intLumiVsLSME_->Fill(nLumi_, intLumi_);
  } else {
    edm::LogError("PixelLumiMonotor") << "No valid data found!";
  }
  /*
  // Access Lumi Details
  Handle<LumiDetails> lumiDetails;
  lumiBlock.getByLabel("expressLumiProducer", lumiDetails);
  if(lumiDetails->isValid()){
    std::cout<<"valid detail"<<std::endl;
  }else{
    std::cout << "no valid lumi detail data" <<std::endl;
  }  */
}

void DQMLumiMonitor::endRun(edm::Run const& iRun, edm::EventSetup const& iSetup) {}

void DQMLumiMonitor::endJob() {}
// Define this as a plug-in
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(DQMLumiMonitor);