MuIsoValidation Class Reference

#include <MuIsoValidation.h>

Inheritance diagram for MuIsoValidation:

Public Member Functions
	MuIsoValidation (const edm::ParameterSet &)
	~MuIsoValidation () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Types
enum	{ NOAXIS, XAXIS, YAXIS, ZAXIS }
typedef edm::Handle< reco::IsoDepositMap >	MuIsoDepHandle
typedef const reco::IsoDeposit	MuIsoDepRef
typedef edm::RefToBase< reco::Muon >	MuonBaseRef
typedef edm::View< reco::Muon >::const_iterator	MuonIterator

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	FillHistos ()
void	InitStatics ()
void	MakeLogBinsForProfile (Double_t *bin_edges, const double min, const double max)
void	RecordData (MuonIterator muon)

Private Attributes
std::vector< std::string >	axis_titles
std::vector< MonitorElement * >	cd_plots
std::vector< int >	cdCompNeeded
std::string	dirName
edm::InputTag	ecalIsoDeposit_Tag
std::vector< MonitorElement * >	h_1D
MonitorElement *	h_nMuons
edm::InputTag	hcalIsoDeposit_Tag
edm::InputTag	hoIsoDeposit_Tag
edm::ParameterSet	iConfig
std::vector< int >	isContinuous
double	L_BIN_WIDTH
int	LOG_BINNING_ENABLED
double	LOG_BINNING_RATIO
std::vector< std::string >	main_titles
edm::InputTag	Muon_Tag
edm::EDGetTokenT< edm::View< reco::Muon > >	Muon_Token
std::vector< std::string >	names
int	nEvents
int	nIncMuons
int	NUM_LOG_BINS
std::vector< std::vector< MonitorElement * > >	p_2D
std::vector< std::vector< double > >	param
bool	requireCombinedMuon
std::string	rootfilename
double	S_BIN_WIDTH
std::string	subDirName
std::string	subsystemname_
double	theData [NUM_VARS]
int	theMuonData
std::string	title_cd
std::string	title_cone
std::string	title_sam
edm::InputTag	tkIsoDeposit_Tag

Static Private Attributes
static const int	NUM_VARS = 21

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 55 of file MuIsoValidation.h.

Member Typedef Documentation

MuIsoDepHandle

typedef edm::Handle<reco::IsoDepositMap> MuIsoValidation::MuIsoDepHandle

private

Definition at line 59 of file MuIsoValidation.h.

MuIsoDepRef

typedef const reco::IsoDeposit MuIsoValidation::MuIsoDepRef

private

Definition at line 60 of file MuIsoValidation.h.

MuonBaseRef

typedef edm::RefToBase<reco::Muon> MuIsoValidation::MuonBaseRef

private

Definition at line 58 of file MuIsoValidation.h.

MuonIterator

typedef edm::View<reco::Muon>::const_iterator MuIsoValidation::MuonIterator

private

Definition at line 57 of file MuIsoValidation.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
NOAXIS
XAXIS
YAXIS
ZAXIS

Definition at line 137 of file MuIsoValidation.h.

{ NOAXIS, XAXIS, YAXIS, ZAXIS };

Constructor & Destructor Documentation

MuIsoValidation()

MuIsoValidation::MuIsoValidation ( const edm::ParameterSet &  ps )

explicit

Definition at line 61 of file MuIsoValidation.cc.

References cd_plots, dirName, edm::ParameterSet::getUntrackedParameter(), h_1D, iConfig, InitStatics(), Muon_Tag, Muon_Token, nEvents, nIncMuons, NUM_VARS, p_2D, requireCombinedMuon, AlCaHLTBitMon_QueryRunRegistry::string, and subsystemname_.

                                                          {
  iConfig = ps;

  //  rootfilename = iConfig.getUntrackedParameter<string>("rootfilename"); // comment out for inclusion
  requireCombinedMuon = iConfig.getUntrackedParameter<bool>("requireCombinedMuon");
  dirName = iConfig.getUntrackedParameter<std::string>("directory");
  //subDirName = iConfig.getParameter<std::string>("@module_label");

  //dirName += subDirName;

  //--------Initialize tags-------
  Muon_Tag = iConfig.getUntrackedParameter<edm::InputTag>("Global_Muon_Label");
  Muon_Token = consumes<edm::View<reco::Muon> >(Muon_Tag);

  //-------Initialize counters----------------
  nEvents = 0;
  nIncMuons = 0;
  //  nCombinedMuons = 0;

  InitStatics();

  //Set up DAQ
  subsystemname_ = iConfig.getUntrackedParameter<std::string>("subSystemFolder", "YourSubsystem");

  //------"allocate" space for the data vectors-------
  h_1D.resize(NUM_VARS);
  cd_plots.resize(NUM_VARS);
  p_2D.resize(NUM_VARS, vector<MonitorElement*>(NUM_VARS));
}

~MuIsoValidation()

MuIsoValidation::~MuIsoValidation ( )

override

Definition at line 94 of file MuIsoValidation.cc.

                                  {
  //Deallocate memory
}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 315 of file MuIsoValidation.cc.

References dqm::impl::MonitorElement::Fill(), FillHistos(), h_nMuons, iEvent, HLT_2023v12_cff::muon, Muon_Token, nEvents, nIncMuons, RecordData(), requireCombinedMuon, theMuonData, and parallelization::uint.

                                                                                 {
  ++nEvents;
  edm::LogInfo("Tutorial") << "\nInvestigating event #" << nEvents << "\n";

  // Get Muon Collection
  edm::Handle<edm::View<reco::Muon> > muonsHandle;  //
  iEvent.getByToken(Muon_Token, muonsHandle);

  //Fill event entry in histogram of number of muons
  edm::LogInfo("Tutorial") << "Number of Muons: " << muonsHandle->size();
  theMuonData = muonsHandle->size();
  h_nMuons->Fill(theMuonData);

  //Fill historgams concerning muon isolation
  uint iMuon = 0;
  for (MuonIterator muon = muonsHandle->begin(); muon != muonsHandle->end(); ++muon, ++iMuon) {
    ++nIncMuons;
    if (requireCombinedMuon) {
      if (muon->combinedMuon().isNull())
        continue;
    }
    RecordData(muon);
    FillHistos();
  }
}

bookHistograms()

void MuIsoValidation::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &    )

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 384 of file MuIsoValidation.cc.

References axis_titles, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), cd_plots, cdCompNeeded, dirName, h_1D, h_nMuons, isContinuous, LOG_BINNING_ENABLED, main_titles, MakeLogBinsForProfile(), names, NUM_LOG_BINS, NUM_VARS, p_2D, param, dqm::impl::MonitorElement::setAxisTitle(), dqm::implementation::NavigatorBase::setCurrentFolder(), title_cd, title_cone, title_sam, trigObjTnPSource_cfi::var, XAXIS, and YAXIS.

                                                                                                       {
  ibooker.setCurrentFolder(dirName);
  //---initialize number of muons histogram---
  h_nMuons = ibooker.book1D("nMuons", title_sam + "Number of Muons", 20, 0., 20.);
  h_nMuons->setAxisTitle("Number of Muons", XAXIS);
  h_nMuons->setAxisTitle("Fraction of Events", YAXIS);

  //---Initialize 1D Histograms---
  for (int var = 0; var < NUM_VARS; var++) {
    h_1D[var] = ibooker.book1D(names[var],
                               title_sam + main_titles[var] + title_cone,
                               (int)param[var][0],
                               param[var][1],
                               param[var][2],
                               [](TH1* th1) { th1->Sumw2(); });
    h_1D[var]->setAxisTitle(axis_titles[var], XAXIS);
    h_1D[var]->setAxisTitle("Fraction of Muons", YAXIS);

    if (cdCompNeeded[var]) {
      cd_plots[var] = ibooker.book1D(names[var] + "_cd",
                                     title_sam + title_cd + main_titles[var] + title_cone,
                                     (int)param[var][0],
                                     param[var][1],
                                     param[var][2],
                                     [](TH1* th1) { th1->Sumw2(); });
      cd_plots[var]->setAxisTitle(axis_titles[var], XAXIS);
      cd_plots[var]->setAxisTitle("Fraction of Muons", YAXIS);
    }
  }  //Finish 1D

  //---Initialize 2D Histograms---
  for (int var1 = 0; var1 < NUM_VARS; var1++) {
    for (int var2 = 0; var2 < NUM_VARS; var2++) {
      if (var1 == var2)
        continue;

      //Monitor elements is weird and takes y axis parameters as well
      //as x axis parameters for a 1D profile plot
      p_2D[var1][var2] = ibooker.bookProfile(names[var1] + "_" + names[var2],
                                             title_sam + main_titles[var2] + " <vs> " + main_titles[var1] + title_cone,
                                             (int)param[var1][0],
                                             param[var1][1],
                                             param[var1][2],
                                             (int)param[var2][0],  //documentation says this is disregarded
                                             param[var2][1],       //does this do anything?
                                             param[var2][2],       //does this do anything?
                                             " ",                  //profile errors = spread/sqrt(num_datums)
                                             [&](TProfile* tprof) {
                                               if (LOG_BINNING_ENABLED && isContinuous[var1]) {
                                                 Double_t* bin_edges = new Double_t[NUM_LOG_BINS + 1];
                                                 // nbins+1 because there is one more edge than there are bins
                                                 MakeLogBinsForProfile(bin_edges, param[var1][1], param[var1][2]);
                                                 tprof->SetBins(NUM_LOG_BINS, bin_edges);
                                                 delete[] bin_edges;
                                               }
                                             });

      p_2D[var1][var2]->setAxisTitle(axis_titles[var1], XAXIS);
      p_2D[var1][var2]->setAxisTitle(axis_titles[var2], YAXIS);
    }
  }  //Finish 2D

  //avg pT not defined for zero tracks.
  //MonitorElement is inflxible and won't let me change the
  //number of bins!  I guess all I'm doing here is changing
  //range of the x axis when it is printed, not the actual
  //bins that are filled
  p_2D[4][9]->setAxisRange(0.5, 15.5, XAXIS);
}

FillHistos()

void MuIsoValidation::FillHistos ( )

private

Definition at line 470 of file MuIsoValidation.cc.

References cd_plots, cdCompNeeded, h_1D, NUM_VARS, p_2D, theData, and trigObjTnPSource_cfi::var.

Referenced by analyze().

                                 {
  //----------Fill 1D histograms---------------
  for (int var = 0; var < NUM_VARS; var++) {
    h_1D[var]->Fill(theData[var]);
    if (cdCompNeeded[var])
      cd_plots[var]->Fill(theData[var]);  //right now, this is a regular PDF (just like h_1D)
  }                                       //Finish 1D

  //----------Fill 2D histograms---------------
  for (int var1 = 0; var1 < NUM_VARS; ++var1) {
    for (int var2 = 0; var2 < NUM_VARS; ++var2) {
      if (var1 == var2)
        continue;
      //change below to regular int interating!
      //      h_2D[var1][var2]->Fill(theData[var1], theData[var2]);
      p_2D[var1][var2]->Fill(theData[var1], theData[var2]);
    }
  }  //Finish 2D
}

InitStatics()

void MuIsoValidation::InitStatics ( )

private

Definition at line 101 of file MuIsoValidation.cc.

References axis_titles, cdCompNeeded, createfilelist::int, isContinuous, L_BIN_WIDTH, LOG_BINNING_ENABLED, LOG_BINNING_RATIO, main_titles, names, NUM_LOG_BINS, NUM_VARS, param, S_BIN_WIDTH, title_cd, title_cone, and title_sam.

Referenced by MuIsoValidation().

                                  {
  //-----------Initialize primatives-----------
  S_BIN_WIDTH = 1.0;  //in GeV
  L_BIN_WIDTH = 2.0;  //in GeV
  LOG_BINNING_ENABLED = 1;
  NUM_LOG_BINS = 15;
  LOG_BINNING_RATIO = 1.1;  //ratio by which each bin is wider than the last for log binning
  //i.e.  bin widths are (x), (r*x), (r^2*x), ..., (r^(nbins)*x)

  //-------Initialize Titles---------
  title_sam = "";   //"[Sample b-jet events] ";
  title_cone = "";  //" [in R=0.3 IsoDeposit Cone]";
  //The above two pieces of info will be printed on the title of the whole page,
  //not for each individual histogram
  title_cd = "C.D. of ";

  //-------"Allocate" memory for vectors
  main_titles.resize(NUM_VARS);
  axis_titles.resize(NUM_VARS);
  names.resize(NUM_VARS);
  param.resize(NUM_VARS, vector<double>(3));
  isContinuous.resize(NUM_VARS);
  cdCompNeeded.resize(NUM_VARS);

  //-----Titles of the plots-----------
  main_titles[0] = "Total Tracker Momentum";
  main_titles[1] = "Total EM Cal Energy";
  main_titles[2] = "Total Had Cal Energy";
  main_titles[3] = "Total HO Cal Energy";
  main_titles[4] = "Number of Tracker Tracks";
  main_titles[5] = "Number of Jets around Muon";
  main_titles[6] = "Tracker p_{T} within veto cone";
  main_titles[7] = "EM E_{T} within veto cone";
  main_titles[8] = "Had E_{T} within veto cone";
  main_titles[9] = "HO E_{T} within veto cone";
  main_titles[10] = "Muon p_{T}";
  main_titles[11] = "Muon #eta";
  main_titles[12] = "Muon #phi";
  main_titles[13] = "Average Momentum per Track ";
  main_titles[14] = "Weighted Energy";
  main_titles[15] = "PF Sum of Charged Hadron Pt";
  main_titles[16] = "PF Sum of Total Hadron Pt";
  main_titles[17] = "PF Sum of E,Mu Pt";
  main_titles[18] = "PF Sum of Neutral Hadron Et";
  main_titles[19] = "PF Sum of Photon Et";
  main_titles[20] = "PF Sum of Pt from non-PV";

  //------Titles on the X or Y axis------------
  axis_titles[0] = "#Sigma p_{T}   (GeV)";
  axis_titles[1] = "#Sigma E_{T}^{EM}   (GeV)";
  axis_titles[2] = "#Sigma E_{T}^{Had}   (GeV)";
  axis_titles[3] = "#Sigma E_{T}^{HO}   (GeV)";
  axis_titles[4] = "N_{Tracks}";
  axis_titles[5] = "N_{Jets}";
  axis_titles[6] = "#Sigma p_{T,veto} (GeV)";
  axis_titles[7] = "#Sigma E_{T,veto}^{EM}   (GeV)";
  axis_titles[8] = "#Sigma E_{T,veto}^{Had}   (GeV)";
  axis_titles[9] = "#Sigma E_{T,veto}^{HO}   (GeV)";
  axis_titles[10] = "p_{T,#mu} (GeV)";
  axis_titles[11] = "#eta_{#mu}";
  axis_titles[12] = "#phi_{#mu}";
  axis_titles[13] = "#Sigma p_{T} / N_{Tracks} (GeV)";
  axis_titles[14] = "(1.5) X #Sigma E_{T}^{EM} + #Sigma E_{T}^{Had}";
  axis_titles[15] = "#Sigma p_{T}^{PFHadCha}   (GeV)";
  axis_titles[16] = "#Sigma p_{T}^{PFTotCha}   (GeV)";
  axis_titles[17] = "#Sigma p_{T}^{PFEMu}   (GeV)";
  axis_titles[18] = "#Sigma E_{T}^{PFHadNeu}   (GeV)";
  axis_titles[19] = "#Sigma E_{T}^{PFPhot}   (GeV)";
  axis_titles[20] = "#Sigma p_{T}^{PFPU}   (GeV)";

  //-----------Names given for the root file----------
  names[0] = "sumPt";
  names[1] = "emEt";
  names[2] = "hadEt";
  names[3] = "hoEt";
  names[4] = "nTracks";
  names[5] = "nJets";
  names[6] = "trackerVetoPt";
  names[7] = "emVetoEt";
  names[8] = "hadVetoEt";
  names[9] = "hoVetoEt";
  names[10] = "muonPt";
  names[11] = "muonEta";
  names[12] = "muonPhi";
  names[13] = "avgPt";
  names[14] = "weightedEt";
  names[15] = "PFsumChargedHadronPt";
  names[16] = "PFsumChargedTotalPt";
  names[17] = "PFsumEMuPt";
  names[18] = "PFsumNeutralHadronEt";
  names[19] = "PFsumPhotonEt";
  names[20] = "PFsumPUPt";

  //----------Parameters for binning of histograms---------
  //param[var][0] is the number of bins
  //param[var][1] is the low edge of the low bin
  //param[var][2] is the high edge of the high bin
  //
  // maximum value------,
  //                    |
  //                    V
  param[0][0] = (int)(20.0 / S_BIN_WIDTH);
  param[0][1] = 0.0;
  param[0][2] = param[0][0] * S_BIN_WIDTH;
  param[1][0] = (int)(20.0 / S_BIN_WIDTH);
  param[1][1] = 0.0;
  param[1][2] = param[1][0] * S_BIN_WIDTH;
  param[2][0] = (int)(20.0 / S_BIN_WIDTH);
  param[2][1] = 0.0;
  param[2][2] = param[2][0] * S_BIN_WIDTH;
  param[3][0] = 20;
  param[3][1] = 0.0;
  param[3][2] = 2.0;
  param[4][0] = 16;
  param[4][1] = -0.5;
  param[4][2] = param[4][0] - 0.5;
  param[5][0] = 4;
  param[5][1] = -0.5;
  param[5][2] = param[5][0] - 0.5;
  param[6][0] = (int)(40.0 / S_BIN_WIDTH);
  param[6][1] = 0.0;
  param[6][2] = param[6][0] * S_BIN_WIDTH;
  param[7][0] = 20;
  param[7][1] = 0.0;
  param[7][2] = 10.0;
  param[8][0] = (int)(20.0 / S_BIN_WIDTH);
  param[8][1] = 0.0;
  param[8][2] = param[8][0] * S_BIN_WIDTH;
  param[9][0] = 20;
  param[9][1] = 0.0;
  param[9][2] = 5.0;
  param[10][0] = (int)(40.0 / S_BIN_WIDTH);
  param[10][1] = 0.0;
  param[10][2] = param[10][0] * S_BIN_WIDTH;
  param[11][0] = 24;
  param[11][1] = -2.4;
  param[11][2] = 2.4;
  param[12][0] = 32;
  param[12][1] = -3.2;
  param[12][2] = 3.2;
  param[13][0] = (int)(15.0 / S_BIN_WIDTH);
  param[13][1] = 0.0;
  param[13][2] = param[13][0] * S_BIN_WIDTH;
  param[14][0] = (int)(20.0 / S_BIN_WIDTH);
  param[14][1] = 0.0;
  param[14][2] = param[14][0] * S_BIN_WIDTH;
  param[15][0] = (int)(20.0 / S_BIN_WIDTH);
  param[15][1] = 0.0;
  param[15][2] = param[15][0] * S_BIN_WIDTH;
  param[16][0] = (int)(20.0 / S_BIN_WIDTH);
  param[15][1] = 0.0;
  param[16][2] = param[16][0] * S_BIN_WIDTH;
  param[17][0] = (int)(20.0 / S_BIN_WIDTH) + 1;
  param[17][1] = -S_BIN_WIDTH;
  param[17][2] = param[17][0] * S_BIN_WIDTH;
  param[18][0] = (int)(20.0 / S_BIN_WIDTH);
  param[18][1] = 0.0;
  param[18][2] = param[18][0] * S_BIN_WIDTH;
  param[19][0] = (int)(20.0 / S_BIN_WIDTH);
  param[19][1] = 0.0;
  param[19][2] = param[19][0] * S_BIN_WIDTH;
  param[20][0] = (int)(20.0 / S_BIN_WIDTH);
  param[20][1] = 0.0;
  param[20][2] = param[20][0] * S_BIN_WIDTH;

  //--------------Is the variable continuous (i.e. non-integer)?-------------
  //---------(Log binning will only be used for continuous variables)--------
  isContinuous[0] = 1;
  isContinuous[1] = 1;
  isContinuous[2] = 1;
  isContinuous[3] = 1;
  isContinuous[4] = 0;
  isContinuous[5] = 0;
  isContinuous[6] = 1;
  isContinuous[7] = 1;
  isContinuous[8] = 1;
  isContinuous[9] = 1;
  isContinuous[10] = 1;
  isContinuous[11] = 1;
  isContinuous[12] = 1;
  isContinuous[13] = 1;
  isContinuous[14] = 1;
  isContinuous[15] = 1;
  isContinuous[16] = 1;
  isContinuous[17] = 1;
  isContinuous[18] = 1;
  isContinuous[19] = 1;
  isContinuous[20] = 1;

  //----Should the cumulative distribution be calculated for this variable?-----
  cdCompNeeded[0] = 1;
  cdCompNeeded[1] = 1;
  cdCompNeeded[2] = 1;
  cdCompNeeded[3] = 1;
  cdCompNeeded[4] = 1;
  cdCompNeeded[5] = 1;
  cdCompNeeded[6] = 1;
  cdCompNeeded[7] = 1;
  cdCompNeeded[8] = 1;
  cdCompNeeded[9] = 1;
  cdCompNeeded[10] = 0;
  cdCompNeeded[11] = 0;
  cdCompNeeded[12] = 0;
  cdCompNeeded[13] = 1;
  cdCompNeeded[14] = 1;
  cdCompNeeded[15] = 1;
  cdCompNeeded[16] = 1;
  cdCompNeeded[17] = 1;
  cdCompNeeded[18] = 1;
  cdCompNeeded[19] = 1;
  cdCompNeeded[20] = 1;
}

MakeLogBinsForProfile()

void MuIsoValidation::MakeLogBinsForProfile ( Double_t *  bin_edges, const double  min, const double  max  )

private

Definition at line 454 of file MuIsoValidation.cc.

References LOG_BINNING_RATIO, SiStripPI::max, SiStripPI::min, dqmiodumpmetadata::n, LaserClient_cfi::nbins, NUM_LOG_BINS, conifer::pow(), and alignCSCRings::r.

Referenced by bookHistograms().

                                                                                                   {
  const double& r = LOG_BINNING_RATIO;
  const int& nbins = NUM_LOG_BINS;

  const double first_bin_width = (r > 1.0) ?  //so we don't divide by zero
                                     (max - min) * (1 - r) / (1 - pow(r, nbins))
                                           : (max - min) / nbins;

  bin_edges[0] = min;
  bin_edges[1] = min + first_bin_width;
  for (int n = 2; n < nbins; ++n) {
    bin_edges[n] = bin_edges[n - 1] + (bin_edges[n - 1] - bin_edges[n - 2]) * r;
  }
  bin_edges[nbins] = max;
}

RecordData()

void MuIsoValidation::RecordData ( MuonIterator  muon )

private

Definition at line 342 of file MuIsoValidation.cc.

References theData.

Referenced by analyze().

                                                  {
  theData[0] = muon->isolationR03().sumPt;
  theData[1] = muon->isolationR03().emEt;
  theData[2] = muon->isolationR03().hadEt;
  theData[3] = muon->isolationR03().hoEt;

  theData[4] = muon->isolationR03().nTracks;
  theData[5] = muon->isolationR03().nJets;
  theData[6] = muon->isolationR03().trackerVetoPt;
  theData[7] = muon->isolationR03().emVetoEt;
  theData[8] = muon->isolationR03().hadVetoEt;
  theData[9] = muon->isolationR03().hoVetoEt;

  theData[10] = muon->pt();
  theData[11] = muon->eta();
  theData[12] = muon->phi();

  // make sure nTracks != 0 before filling this one
  if (theData[4] != 0)
    theData[13] = (double)theData[0] / (double)theData[4];
  else
    theData[13] = -99;

  theData[14] = 1.5 * theData[1] + theData[2];

  // Now PF isolation
  theData[15] = -99.;
  theData[16] = -99.;
  theData[17] = -99.;
  theData[18] = -99.;
  theData[19] = -99.;
  theData[20] = -99.;
  if (muon->isPFMuon() && muon->isPFIsolationValid()) {
    theData[15] = muon->pfIsolationR03().sumChargedHadronPt;
    theData[16] = muon->pfIsolationR03().sumChargedParticlePt;
    theData[17] = muon->pfIsolationR03().sumChargedParticlePt - muon->pfIsolationR03().sumChargedHadronPt;
    theData[18] = muon->pfIsolationR03().sumNeutralHadronEt;
    theData[19] = muon->pfIsolationR03().sumPhotonEt;
    theData[20] = muon->pfIsolationR03().sumPUPt;
  }
}

Member Data Documentation

axis_titles

std::vector<std::string> MuIsoValidation::axis_titles

private

Definition at line 109 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

cd_plots

std::vector<MonitorElement*> MuIsoValidation::cd_plots

private

Definition at line 127 of file MuIsoValidation.h.

Referenced by bookHistograms(), FillHistos(), and MuIsoValidation().

cdCompNeeded

std::vector<int> MuIsoValidation::cdCompNeeded

private

Definition at line 113 of file MuIsoValidation.h.

Referenced by bookHistograms(), FillHistos(), and InitStatics().

dirName

std::string MuIsoValidation::dirName

private

Definition at line 90 of file MuIsoValidation.h.

Referenced by bookHistograms(), MuIsoValidation(), and core.JSONAnalyzer.JSONAnalyzer::write().

ecalIsoDeposit_Tag

edm::InputTag MuIsoValidation::ecalIsoDeposit_Tag

private

Definition at line 84 of file MuIsoValidation.h.

h_1D

std::vector<MonitorElement*> MuIsoValidation::h_1D

private

Definition at line 126 of file MuIsoValidation.h.

Referenced by bookHistograms(), FillHistos(), and MuIsoValidation().

h_nMuons

MonitorElement* MuIsoValidation::h_nMuons

private

Definition at line 125 of file MuIsoValidation.h.

Referenced by analyze(), and bookHistograms().

hcalIsoDeposit_Tag

edm::InputTag MuIsoValidation::hcalIsoDeposit_Tag

private

Definition at line 83 of file MuIsoValidation.h.

hoIsoDeposit_Tag

edm::InputTag MuIsoValidation::hoIsoDeposit_Tag

private

Definition at line 85 of file MuIsoValidation.h.

iConfig

edm::ParameterSet MuIsoValidation::iConfig

private

Definition at line 119 of file MuIsoValidation.h.

Referenced by MuIsoValidation().

isContinuous

std::vector<int> MuIsoValidation::isContinuous

private

Definition at line 112 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

L_BIN_WIDTH

double MuIsoValidation::L_BIN_WIDTH

private

Definition at line 97 of file MuIsoValidation.h.

Referenced by InitStatics().

LOG_BINNING_ENABLED

int MuIsoValidation::LOG_BINNING_ENABLED

private

Definition at line 99 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

LOG_BINNING_RATIO

double MuIsoValidation::LOG_BINNING_RATIO

private

Definition at line 101 of file MuIsoValidation.h.

Referenced by InitStatics(), and MakeLogBinsForProfile().

main_titles

std::vector<std::string> MuIsoValidation::main_titles

private

Definition at line 108 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

Muon_Tag

edm::InputTag MuIsoValidation::Muon_Tag

private

Definition at line 80 of file MuIsoValidation.h.

Referenced by MuIsoValidation().

Muon_Token

edm::EDGetTokenT<edm::View<reco::Muon> > MuIsoValidation::Muon_Token

private

Definition at line 81 of file MuIsoValidation.h.

Referenced by analyze(), and MuIsoValidation().

names

std::vector<std::string> MuIsoValidation::names

private

Definition at line 110 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

nEvents

int MuIsoValidation::nEvents

private

Definition at line 132 of file MuIsoValidation.h.

Referenced by analyze(), and MuIsoValidation().

nIncMuons

int MuIsoValidation::nIncMuons

private

Definition at line 133 of file MuIsoValidation.h.

Referenced by analyze(), and MuIsoValidation().

NUM_LOG_BINS

int MuIsoValidation::NUM_LOG_BINS

private

Definition at line 100 of file MuIsoValidation.h.

Referenced by bookHistograms(), InitStatics(), and MakeLogBinsForProfile().

NUM_VARS

const int MuIsoValidation::NUM_VARS = 21

staticprivate

Definition at line 96 of file MuIsoValidation.h.

Referenced by bookHistograms(), FillHistos(), InitStatics(), and MuIsoValidation().

p_2D

std::vector<std::vector<MonitorElement*> > MuIsoValidation::p_2D

private

Definition at line 129 of file MuIsoValidation.h.

Referenced by bookHistograms(), FillHistos(), and MuIsoValidation().

param

std::vector<std::vector<double> > MuIsoValidation::param

private

Definition at line 111 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

requireCombinedMuon

bool MuIsoValidation::requireCombinedMuon

private

Definition at line 102 of file MuIsoValidation.h.

Referenced by analyze(), and MuIsoValidation().

rootfilename

std::string MuIsoValidation::rootfilename

private

Definition at line 88 of file MuIsoValidation.h.

S_BIN_WIDTH

double MuIsoValidation::S_BIN_WIDTH

private

Definition at line 98 of file MuIsoValidation.h.

Referenced by InitStatics().

subDirName

std::string MuIsoValidation::subDirName

private

Definition at line 91 of file MuIsoValidation.h.

subsystemname_

std::string MuIsoValidation::subsystemname_

private

Definition at line 93 of file MuIsoValidation.h.

Referenced by MuIsoValidation().

theData

double MuIsoValidation::theData[NUM_VARS]

private

Definition at line 122 of file MuIsoValidation.h.

Referenced by FillHistos(), and RecordData().

theMuonData

int MuIsoValidation::theMuonData

private

Definition at line 121 of file MuIsoValidation.h.

Referenced by analyze().

title_cd

std::string MuIsoValidation::title_cd

private

Definition at line 106 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

title_cone

std::string MuIsoValidation::title_cone

private

Definition at line 105 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

title_sam

std::string MuIsoValidation::title_sam

private

Definition at line 104 of file MuIsoValidation.h.

Referenced by bookHistograms(), and InitStatics().

tkIsoDeposit_Tag

edm::InputTag MuIsoValidation::tkIsoDeposit_Tag

private

Definition at line 82 of file MuIsoValidation.h.