MuonTrackResidualsTest.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author G. Mila - INFN Torino
 */

#include <DQMOffline/Muon/interface/MuonTrackResidualsTest.h>

// Framework
#include <FWCore/Framework/interface/Event.h>
#include "DataFormats/Common/interface/Handle.h"
#include <FWCore/Framework/interface/ESHandle.h>
#include <FWCore/Framework/interface/MakerMacros.h>
#include <FWCore/Framework/interface/EventSetup.h>
#include <FWCore/ParameterSet/interface/ParameterSet.h>

#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/Run.h"

#include <iostream>
#include <cstdio>
#include <string>
#include <cmath>
#include "TF1.h"

using namespace edm;
using namespace std;

MuonTrackResidualsTest::MuonTrackResidualsTest(const edm::ParameterSet& ps) {
  parameters = ps;

  prescaleFactor = parameters.getUntrackedParameter<int>("diagnosticPrescale", 1);

  GaussianCriterionName =
      parameters.getUntrackedParameter<string>("resDistributionTestName", "ResidualsDistributionGaussianTest");
  SigmaCriterionName = parameters.getUntrackedParameter<string>("sigmaTestName", "ResidualsSigmaInRange");
  MeanCriterionName = parameters.getUntrackedParameter<string>("meanTestName", "ResidualsMeanInRange");
}
void MuonTrackResidualsTest::dqmEndRun(DQMStore::IBooker& ibooker,
                                       DQMStore::IGetter& igetter,
                                       edm::Run const&,
                                       edm::EventSetup const&) {
  ibooker.setCurrentFolder("Muons/Tests/trackResidualsTest");

  string histName, MeanHistoName, SigmaHistoName, MeanHistoTitle, SigmaHistoTitle;
  vector<string> type;
  type.push_back("eta");
  type.push_back("theta");
  type.push_back("phi");

  for (unsigned int c = 0; c < type.size(); c++) {
    MeanHistoName = "MeanTest_" + type[c];
    SigmaHistoName = "SigmaTest_" + type[c];

    MeanHistoTitle = "Mean of the #" + type[c] + " residuals distribution";
    SigmaHistoTitle = "Sigma of the #" + type[c] + " residuals distribution";

    histName = "Res_GlbSta_" + type[c];
    histoNames[type[c]].push_back(histName);
    histName = "Res_TkGlb_" + type[c];
    histoNames[type[c]].push_back(histName);
    histName = "Res_TkSta_" + type[c];
    histoNames[type[c]].push_back(histName);

    MeanHistos[type[c]] = ibooker.book1D(MeanHistoName.c_str(), MeanHistoTitle.c_str(), 3, 0.5, 3.5);
    (MeanHistos[type[c]])->setBinLabel(1, "Res_StaGlb", 1);
    (MeanHistos[type[c]])->setBinLabel(2, "Res_TkGlb", 1);
    (MeanHistos[type[c]])->setBinLabel(3, "Res_TkSta", 1);

    SigmaHistos[type[c]] = ibooker.book1D(SigmaHistoName.c_str(), SigmaHistoTitle.c_str(), 3, 0.5, 3.5);
    (SigmaHistos[type[c]])->setBinLabel(1, "Res_StaGlb", 1);
    (SigmaHistos[type[c]])->setBinLabel(2, "Res_TkGlb", 1);
    (SigmaHistos[type[c]])->setBinLabel(3, "Res_TkSta", 1);
  }

  for (map<string, vector<string> >::const_iterator histo = histoNames.begin(); histo != histoNames.end(); histo++) {
    for (unsigned int type = 0; type < (*histo).second.size(); type++) {
      string path = "Muons/MuonRecoAnalyzer/" + (*histo).second[type];
      MonitorElement* res_histo = igetter.get(path);
      if (res_histo) {
        // gaussian test
        //  const QReport *GaussianReport = res_histo->getQReport(GaussianCriterionName);
        int BinNumber = type + 1;
        float mean = (*res_histo).getMean(1);
        float sigma = (*res_histo).getRMS(1);
        MeanHistos.find((*histo).first)->second->setBinContent(BinNumber, mean);
        SigmaHistos.find((*histo).first)->second->setBinContent(BinNumber, sigma);
      }
    }
  }
}