SiStripMonitorHLT.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiStripMonitorCluster
// Class:      SiStripMonitorHLT
//
// class SiStripMonitorHLT SiStripMonitorHLT.cc
// DQM/SiStripMonitorCluster/src/SiStripMonitorHLT.cc
#include <vector>

#include <iostream>
#include <numeric>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"

#include "DQM/SiStripMonitorCluster/interface/SiStripMonitorHLT.h"
#include "DQMServices/Core/interface/DQMStore.h"

SiStripMonitorHLT::SiStripMonitorHLT(const edm::ParameterSet& iConfig) {
  HLTDirectory = "HLTResults";
  dqmStore_ = edm::Service<DQMStore>().operator->();
  conf_ = iConfig;

  filerDecisionToken_ =
      consumes<int>(conf_.getParameter<std::string>("HLTProducer"));
  sumOfClusterToken_ =
      consumes<uint>(conf_.getParameter<std::string>("HLTProducer"));
  clusterInSubComponentsToken_ =
      consumes<std::map<uint, std::vector<SiStripCluster> > >(
          conf_.getParameter<std::string>("HLTProducer"));
}

void SiStripMonitorHLT::bookHistograms(DQMStore::IBooker& ibooker,
                                       const edm::Run& run,
                                       const edm::EventSetup& es) {
  ibooker.setCurrentFolder(HLTDirectory);
  std::string HLTProducer = conf_.getParameter<std::string>("HLTProducer");
  HLTDecision = ibooker.book1D(HLTProducer + "_HLTDecision",
                               HLTProducer + "HLTDecision", 2, -0.5, 1.5);
  // all
  SumOfClusterCharges_all = ibooker.book1D(
      "SumOfClusterCharges_all", "SumOfClusterCharges_all", 50, 0, 2000);
  ChargeOfEachClusterTIB_all =
      ibooker.book1D("ChargeOfEachClusterTIB_all", "ChargeOfEachClusterTIB_all",
                     400, -0.5, 400.5);
  ChargeOfEachClusterTOB_all =
      ibooker.book1D("ChargeOfEachClusterTOB_all", "ChargeOfEachClusterTOB_all",
                     400, -0.5, 400.5);
  ChargeOfEachClusterTEC_all =
      ibooker.book1D("ChargeOfEachClusterTEC_all", "ChargeOfEachClusterTEC_all",
                     400, -0.5, 400.5);
  NumberOfClustersAboveThreshold_all =
      ibooker.book1D("NumberOfClustersAboveThreshold_all",
                     "NumberOfClustersAboveThreshold_all", 30, 30.5, 60.5);
  // 31 = TIB2, 32 = TIB2, 33 = TIB3, 51 = TOB1, 52=TOB2, 60 = TEC
  // accepted from HLT
  SumOfClusterCharges_hlt = ibooker.book1D(
      "SumOfClusterCharges_hlt", "SumOfClusterCharges_hlt", 50, 0, 2000);
  ChargeOfEachClusterTIB_hlt =
      ibooker.book1D("ChargeOfEachClusterTIB_hlt", "ChargeOfEachClusterTIB_hlt",
                     400, -0.5, 400.5);
  ChargeOfEachClusterTOB_hlt =
      ibooker.book1D("ChargeOfEachClusterTOB_hlt", "ChargeOfEachClusterTOB_hlt",
                     400, -0.5, 400.5);
  ChargeOfEachClusterTEC_hlt =
      ibooker.book1D("ChargeOfEachClusterTEC_hlt", "ChargeOfEachClusterTEC_hlt",
                     400, -0.5, 400.5);
  NumberOfClustersAboveThreshold_hlt =
      ibooker.book1D("NumberOfClustersAboveThreshold_hlt",
                     "NumberOfClustersAboveThreshold_hlt", 30, 30.5, 60.5);
}

void SiStripMonitorHLT::analyze(const edm::Event& iEvent,
                                const edm::EventSetup& iSetup) {
  // get from event
  std::string HLTProducer = conf_.getParameter<std::string>("HLTProducer");
  edm::Handle<int> filter_decision;
  iEvent.getByToken(filerDecisionToken_, filter_decision);  // filter decision
  edm::Handle<uint> sum_of_clustch;
  iEvent.getByToken(sumOfClusterToken_,
                    sum_of_clustch);  // sum of cluster charges
  // first element of pair: layer: TIB1, ...., TEC; second element: nr of
  // clusters above threshold
  edm::Handle<std::map<uint, std::vector<SiStripCluster> > >
      clusters_in_subcomponents;
  if (HLTProducer == "ClusterMTCCFilter")
    iEvent.getByToken(clusterInSubComponentsToken_, clusters_in_subcomponents);

  // trigger decision
  HLTDecision->Fill(*filter_decision);

  // sum of charges of clusters
  SumOfClusterCharges_all->Fill(*sum_of_clustch);
  if (*filter_decision) SumOfClusterCharges_hlt->Fill(*sum_of_clustch);

  // clusters in different layers
  if (HLTProducer == "ClusterMTCCFilter") {
    // loop over layers ("subcomponents")
    for (std::map<uint, std::vector<SiStripCluster> >::const_iterator it =
             clusters_in_subcomponents->begin();
         it != clusters_in_subcomponents->end(); it++) {
      int generalized_layer = it->first;
      std::vector<SiStripCluster> theclusters = it->second;
      NumberOfClustersAboveThreshold_all->Fill(
          generalized_layer,
          theclusters.size());  // number of clusters in this generalized layer
      if (*filter_decision)
        NumberOfClustersAboveThreshold_hlt->Fill(generalized_layer,
                                                 theclusters.size());
      // loop over clusters (and detids)
      for (std::vector<SiStripCluster>::const_iterator icluster =
               theclusters.begin();
           icluster != theclusters.end(); icluster++) {
        // calculate sum of amplitudes
        unsigned int amplclus = 0;
        for (auto ia = icluster->amplitudes().begin();
             ia != icluster->amplitudes().end(); ia++) {
          if ((*ia) > 0) amplclus += (*ia);  // why should this be negative?
        }
        if (generalized_layer == 31 || generalized_layer == 32 ||
            generalized_layer ==
                33) {  // you can also ask the detid here whether is TIB
          ChargeOfEachClusterTIB_all->Fill(amplclus, 1.);
          if (*filter_decision) ChargeOfEachClusterTIB_hlt->Fill(amplclus, 1.);
        }
        if (generalized_layer == 51 || generalized_layer == 52) {
          ChargeOfEachClusterTOB_all->Fill(amplclus, 1.);
          if (*filter_decision) ChargeOfEachClusterTOB_hlt->Fill(amplclus, 1.);
        }
        if (generalized_layer == 60) {
          ChargeOfEachClusterTEC_all->Fill(amplclus, 1.);
          if (*filter_decision) ChargeOfEachClusterTEC_hlt->Fill(amplclus, 1.);
        }
      }
    }
  }
}