#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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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 ()

TH1 *	GetTH1FromMonitorElement (MonitorElement *me)

TH2 *	GetTH2FromMonitorElement (MonitorElement *me)

TProfile *	GetTProfileFromMonitorElement (MonitorElement *me)

void	InitStatics ()

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

void	NormalizeHistos ()

void	RecordData (MuonIterator muon)

Private Attributes
std::vector< std::string >	axis_titles

std::vector< MonitorElement * >	cd_plots

std::vector< int >	cdCompNeeded

DQMStore *	dbe

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 edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 59 of file MuIsoValidation.h.

Member Typedef Documentation

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

private

Definition at line 63 of file MuIsoValidation.h.

typedef const reco::IsoDeposit MuIsoValidation::MuIsoDepRef

private

Definition at line 64 of file MuIsoValidation.h.

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

private

Definition at line 62 of file MuIsoValidation.h.

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

private

Definition at line 61 of file MuIsoValidation.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
NOAXIS
XAXIS
YAXIS
ZAXIS

Definition at line 148 of file MuIsoValidation.h.

148 {NOAXIS,XAXIS,YAXIS,ZAXIS};

MuIsoValidation::NOAXIS

Definition: MuIsoValidation.h:148

MuIsoValidation::XAXIS

Definition: MuIsoValidation.h:148

MuIsoValidation::YAXIS

Definition: MuIsoValidation.h:148

MuIsoValidation::ZAXIS

Definition: MuIsoValidation.h:148

Constructor & Destructor Documentation

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

explicit

Definition at line 63 of file MuIsoValidation.cc.

References cd_plots, dbe, dirName, edm::ParameterSet::getUntrackedParameter(), h_1D, iConfig, InitStatics(), Muon_Tag, Muon_Token, nEvents, nIncMuons, NUM_VARS, Utilities::operator, 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
   dbe = nullptr;
   dbe = edm::Service<DQMStore>().operator->();
   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 ( )

override

Definition at line 100 of file MuIsoValidation.cc.

                                  {
   
   //Deallocate memory
   
 }

Member Function Documentation

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

overrideprivate

Definition at line 284 of file MuIsoValidation.cc.

References MonitorElement::Fill(), FillHistos(), edm::Event::getByToken(), h_nMuons, metsig::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();
   }
 }

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 355 of file MuIsoValidation.cc.

References axis_titles, create_public_lumi_plots::bin_edges, DQMStore::IBooker::book1D(), DQMStore::IBooker::bookProfile(), cd_plots, cdCompNeeded, dirName, GetTH1FromMonitorElement(), GetTProfileFromMonitorElement(), h_1D, h_nMuons, isContinuous, LOG_BINNING_ENABLED, main_titles, MakeLogBinsForProfile(), names, NUM_LOG_BINS, NUM_VARS, p_2D, param, MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), title_cd, title_cone, title_sam, JetChargeProducer_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]
                 );
     h_1D[var]->setAxisTitle(axis_titles[var],XAXIS);
     h_1D[var]->setAxisTitle("Fraction of Muons",YAXIS);
     GetTH1FromMonitorElement(h_1D[var])->Sumw2();
 
     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]
                   );
       cd_plots[var]->setAxisTitle(axis_titles[var],XAXIS);
       cd_plots[var]->setAxisTitle("Fraction of Muons",YAXIS);
       GetTH1FromMonitorElement(cd_plots[var])->Sumw2();
     }
   }//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;
       
       /*      h_2D[var1][var2] = ibooker.book2D(
           names[var1] + "_" + names[var2] + "_s",
           //title is in "y-var vs. x-var" format
           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],
           param[var2][1],
           param[var2][2]
           );
       */
       //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)
                       );
       
       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]);
     GetTProfileFromMonitorElement(p_2D[var1][var2])->SetBins(NUM_LOG_BINS, bin_edges);
     delete[] bin_edges;
       }
       /*      h_2D[var1][var2]->setAxisTitle(axis_titles[var1],XAXIS);
           h_2D[var1][var2]->setAxisTitle(axis_titles[var2],YAXIS);
           GetTH2FromMonitorElement(h_2D[var1][var2])->Sumw2();
       */
       
       p_2D[var1][var2]->setAxisTitle(axis_titles[var1],XAXIS);
       p_2D[var1][var2]->setAxisTitle(axis_titles[var2],YAXIS);
       //            GetTProfileFromMonitorElement(p_2D[var1][var2])->Sumw2();
     }
   }//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);
   
 }

void MuIsoValidation::FillHistos ( )

private

Definition at line 491 of file MuIsoValidation.cc.

References cd_plots, cdCompNeeded, h_1D, NUM_VARS, p_2D, theData, and JetChargeProducer_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
 }

TH1 * MuIsoValidation::GetTH1FromMonitorElement ( MonitorElement * me )

private

Definition at line 511 of file MuIsoValidation.cc.

References MonitorElement::getTH1().

Referenced by bookHistograms(), and NormalizeHistos().

                                                                  {
   return me->getTH1();
 }

TH2 * MuIsoValidation::GetTH2FromMonitorElement ( MonitorElement * me )

private

Definition at line 515 of file MuIsoValidation.cc.

References MonitorElement::getTH2F().

                                                                  {
   return me->getTH2F();
 }

TProfile * MuIsoValidation::GetTProfileFromMonitorElement ( MonitorElement * me )

private

Definition at line 519 of file MuIsoValidation.cc.

References DEFINE_FWK_MODULE, and MonitorElement::getTProfile().

Referenced by bookHistograms().

                                                                            {
   return me->getTProfile();
 }

void MuIsoValidation::InitStatics ( )

private

Definition at line 109 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;
 }

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

private

Definition at line 451 of file MuIsoValidation.cc.

References LOG_BINNING_RATIO, SiStripPI::max, min(), gen::n, pileupCalc::nbins, NUM_LOG_BINS, funct::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;
 }

void MuIsoValidation::NormalizeHistos ( )

private

Definition at line 469 of file MuIsoValidation.cc.

References cd_plots, cdCompNeeded, GetTH1FromMonitorElement(), h_1D, createfilelist::int, gen::n, NUM_VARS, param, and JetChargeProducer_cfi::var.

                                       {
   for(int var=0; var<NUM_VARS; var++){   
     //turn cd_plots into CDF's
     //underflow -> bin #0.  overflow -> bin #(nbins+1)
     //0th bin doesn't need changed
     
     //----normalize------
     double entries = GetTH1FromMonitorElement(h_1D[var])->GetEntries();
     if (entries==0)continue;
     GetTH1FromMonitorElement(h_1D[var])->Scale(1./entries);
 
     if (cdCompNeeded[var]) {
       int n_max = int(param[var][0])+1;
       for(int n=1; n<=n_max; ++n){
     cd_plots[var]->setBinContent(n, cd_plots[var]->getBinContent(n) + cd_plots[var]->getBinContent(n-1)); //Integrate.
       }
       //----normalize------
       GetTH1FromMonitorElement(cd_plots[var])->Scale(1./entries);
     }
   }
 }

void MuIsoValidation::RecordData ( MuonIterator muon )

private

Definition at line 311 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;
   }
 
 }