#include <SiStripMonitorHLT.h>

Inheritance diagram for SiStripMonitorHLT:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	SiStripMonitorHLT (const edm::ParameterSet &)

	~SiStripMonitorHLT () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Private Attributes
MonitorElement *	ChargeOfEachClusterTEC_all

MonitorElement *	ChargeOfEachClusterTEC_hlt

MonitorElement *	ChargeOfEachClusterTIB_all

MonitorElement *	ChargeOfEachClusterTIB_hlt

MonitorElement *	ChargeOfEachClusterTOB_all

MonitorElement *	ChargeOfEachClusterTOB_hlt

edm::EDGetTokenT< std::map< uint, std::vector< SiStripCluster > > >	clusterInSubComponentsToken_

edm::ParameterSet	conf_

edm::EDGetTokenT< int >	filerDecisionToken_

MonitorElement *	HLTDecision

std::string	HLTDirectory

MonitorElement *	NumberOfClustersAboveThreshold_all

MonitorElement *	NumberOfClustersAboveThreshold_hlt

MonitorElement *	SumOfClusterCharges_all

MonitorElement *	SumOfClusterCharges_hlt

edm::EDGetTokenT< uint >	sumOfClusterToken_

Detailed Description

Definition at line 24 of file SiStripMonitorHLT.h.

Constructor & Destructor Documentation

SiStripMonitorHLT::SiStripMonitorHLT ( const edm::ParameterSet & iConfig )

explicit

Definition at line 21 of file SiStripMonitorHLT.cc.

References clusterInSubComponentsToken_, conf_, filerDecisionToken_, edm::ParameterSet::getParameter(), HLTDirectory, AlCaHLTBitMon_QueryRunRegistry::string, and sumOfClusterToken_.

                                                                    {
   HLTDirectory = "HLTResults";
   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"));
 }

SiStripMonitorHLT::~SiStripMonitorHLT ( )

inlineoverride

Definition at line 27 of file SiStripMonitorHLT.h.

References analyze(), and bookHistograms().

27 {};

Member Function Documentation

void SiStripMonitorHLT::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 74 of file SiStripMonitorHLT.cc.

References ChargeOfEachClusterTEC_all, ChargeOfEachClusterTEC_hlt, ChargeOfEachClusterTIB_all, ChargeOfEachClusterTIB_hlt, ChargeOfEachClusterTOB_all, ChargeOfEachClusterTOB_hlt, clusterInSubComponentsToken_, conf_, filerDecisionToken_, MonitorElement::Fill(), edm::Event::getByToken(), edm::ParameterSet::getParameter(), HLTDecision, genParticles_cff::map, NumberOfClustersAboveThreshold_all, NumberOfClustersAboveThreshold_hlt, AlCaHLTBitMon_QueryRunRegistry::string, SumOfClusterCharges_all, SumOfClusterCharges_hlt, sumOfClusterToken_, and parallelization::uint().

Referenced by ~SiStripMonitorHLT().

                                                              {
   // 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.);
         }
       }
     }
   }
 }

void SiStripMonitorHLT::bookHistograms	(	DQMStore::IBooker &	,
		edm::Run const &	,
		edm::EventSetup const &
	)

override

Definition at line 34 of file SiStripMonitorHLT.cc.

References DQMStore::IBooker::book1D(), ChargeOfEachClusterTEC_all, ChargeOfEachClusterTEC_hlt, ChargeOfEachClusterTIB_all, ChargeOfEachClusterTIB_hlt, ChargeOfEachClusterTOB_all, ChargeOfEachClusterTOB_hlt, conf_, edm::ParameterSet::getParameter(), HLTDecision, HLTDirectory, NumberOfClustersAboveThreshold_all, NumberOfClustersAboveThreshold_hlt, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, SumOfClusterCharges_all, and SumOfClusterCharges_hlt.

Referenced by ~SiStripMonitorHLT().

                                                                 {
   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);
 }