Inheritance diagram for RecoTauPlotDiscriminator:

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

	RecoTauPlotDiscriminator (const edm::ParameterSet &pset)

	~RecoTauPlotDiscriminator () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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
typedef std::map< std::string, HistoMap >	DiscMap

typedef std::map< std::string, TH1 * >	HistoMap

typedef std::vector< edm::InputTag >	VInputTag

Private Attributes
VInputTag	discs_

DiscMap	histos_

edm::InputTag	pileupInfoSrc_

double	pileupPtCut_

edm::InputTag	pileupVerticesSrc_

bool	plotPU_

edm::InputTag	src_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, ESInputTag const &tag)

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)

Detailed Description

Definition at line 34 of file RecoTauPlotDiscriminator.cc.

Member Typedef Documentation

typedef std::map<std::string, HistoMap> RecoTauPlotDiscriminator::DiscMap

private

Definition at line 42 of file RecoTauPlotDiscriminator.cc.

typedef std::map<std::string, TH1*> RecoTauPlotDiscriminator::HistoMap

private

Definition at line 41 of file RecoTauPlotDiscriminator.cc.

typedef std::vector<edm::InputTag> RecoTauPlotDiscriminator::VInputTag

private

Definition at line 43 of file RecoTauPlotDiscriminator.cc.

Constructor & Destructor Documentation

RecoTauPlotDiscriminator::RecoTauPlotDiscriminator ( const edm::ParameterSet & pset )

Definition at line 52 of file RecoTauPlotDiscriminator.cc.

References discs_, edm::ParameterSet::getParameter(), histos_, edm::InputTag::label(), SiStripPI::max, min(), pileupCalc::nbins, pileupInfoSrc_, pileupPtCut_, pileupVerticesSrc_, plotPU_, AlCaHLTBitMon_QueryRunRegistry::string, and GlobalPosition_Frontier_DevDB_cff::tag.

   :src_(pset.getParameter<edm::InputTag>("src")) {
   uint32_t nbins = pset.getParameter<uint32_t>("nbins");
   double min = pset.getParameter<double>("min");
   double max = pset.getParameter<double>("max");
   edm::Service<TFileService> fs;
   // Get the discriminators
   discs_ = pset.getParameter<std::vector<edm::InputTag> >("discriminators");
 
   plotPU_ = pset.getParameter<bool>("plotPU");
   if (plotPU_) {
     pileupPtCut_ = pset.getParameter<double>("pileupTauPtCut");
     pileupInfoSrc_ = pset.getParameter<edm::InputTag>("pileupInfo");
     pileupVerticesSrc_ = pset.getParameter<edm::InputTag>("pileupVertices");
   }
 
   BOOST_FOREACH(const edm::InputTag &tag, discs_) {
     HistoMap discMap;
     discMap["plain"] =
         fs->make<TH1F>(tag.label().c_str(), tag.label().c_str(),
                        nbins, min, max);
 
     // Make correlation plots w.r.t tau pt
     std::string vs_pt_name = tag.label()+"_pt";
     discMap["vs_pt"] =
         fs->make<TH2F>(vs_pt_name.c_str(), vs_pt_name.c_str(),
                        nbins, min, max, 100, 0, 200);
 
     // W.r.t. jet pt
     std::string vs_jetpt_name = tag.label()+"_jetPt";
     discMap["vs_jetPt"] =
         fs->make<TH2F>(vs_jetpt_name.c_str(), vs_jetpt_name.c_str(),
                        nbins, min, max, 100, 0, 200);
 
     // W.r.t. embedded pt in alternat lorentz vector (used to hold gen tau pt)
     std::string vs_embedpt_name = tag.label()+"_embedPt";
     discMap["vs_embedPt"] =
         fs->make<TH2F>(vs_embedpt_name.c_str(), vs_embedpt_name.c_str(),
                        nbins, min, max, 100, 0, 200);
 
     // 3D histogram with tau pt & jet pt
     std::string vs_pt_jetPt_name = tag.label()+"_pt_jetPt";
     discMap["vs_pt_jetPt"] =
         fs->make<TH3F>(vs_pt_jetPt_name.c_str(), vs_pt_jetPt_name.c_str(),
                        nbins, min, max, 100, 0, 200, 100, 0, 200);
 
     std::string vs_pt_embedPt_name = tag.label()+"_pt_embedPt";
     discMap["vs_pt_embedPt"] =
         fs->make<TH3F>(vs_pt_embedPt_name.c_str(), vs_pt_embedPt_name.c_str(),
                        nbins, min, max, 100, 0, 200, 100, 0, 200);
 
 
     std::string vs_eta_name = tag.label()+"_eta";
     discMap["vs_eta"] =
         fs->make<TH2F>(vs_eta_name.c_str(), vs_eta_name.c_str(),
                        nbins, min, max, 100, -2.5, 2.5);
 
     std::string vs_dm_name = tag.label()+"_dm";
     discMap["vs_dm"] =
         fs->make<TH2F>(vs_dm_name.c_str(), vs_dm_name.c_str(),
                        nbins, min, max, 15, -0.5, 14.5);
 
     if (plotPU_) {
       std::string vs_truePU_name = tag.label()+"_truePU";
       discMap["vs_truePU"] = fs->make<TH2F>(vs_truePU_name.c_str(),
           vs_truePU_name.c_str(), nbins, min, max, 15, -0.5, 14.5);
       std::string vs_recoPU_name = tag.label()+"_recoPU";
       discMap["vs_recoPU"] = fs->make<TH2F>(vs_recoPU_name.c_str(),
           vs_recoPU_name.c_str(), nbins, min, max, 15, -0.5, 14.5);
     }
 
     histos_[tag.label()] = discMap;
   }
 }

RecoTauPlotDiscriminator::~RecoTauPlotDiscriminator ( )

inlineoverride

Definition at line 37 of file RecoTauPlotDiscriminator.cc.

References analyze().

37 {}

Member Function Documentation

void RecoTauPlotDiscriminator::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 128 of file RecoTauPlotDiscriminator.cc.

References DEFINE_FWK_MODULE, discs_, HcalObjRepresent::Fill(), edm::Event::getByLabel(), histos_, input, edm::InputTag::label(), pileupInfoSrc_, pileupPtCut_, pileupVerticesSrc_, plotPU_, mps_fire::result, src_, GlobalPosition_Frontier_DevDB_cff::tag, and metsig::tau.

Referenced by ~RecoTauPlotDiscriminator().

                                                            {
   // Get the input collection to clean
   edm::Handle<reco::CandidateView> input;
   evt.getByLabel(src_, input);
 
   // Cast the input candidates to Refs to real taus
   reco::PFTauRefVector inputRefs =
       reco::tau::castView<reco::PFTauRefVector>(input);
 
   edm::Handle<PileupSummaryInfo> puInfo;
   edm::Handle<reco::VertexCollection> puVertices;
   if (plotPU_) {
     evt.getByLabel(pileupInfoSrc_, puInfo);
     evt.getByLabel(pileupVerticesSrc_, puVertices);
   }
 
   // Plot the discriminator output for each of our taus
   BOOST_FOREACH(const reco::PFTauRef& tau, inputRefs) {
     // Plot each discriminator
     BOOST_FOREACH(const edm::InputTag &tag, discs_) {
       edm::Handle<reco::PFTauDiscriminator> discHandle;
       evt.getByLabel(tag, discHandle);
       //const HistoMap &discHistos = disc.second;
       double result = (*discHandle)[tau];
       HistoMap& mymap = histos_[tag.label()];
       mymap["plain"]->Fill(result);
       mymap["vs_pt"]->Fill(result, tau->pt());
       mymap["vs_jetPt"]->Fill(result, tau->jetRef()->pt());
       mymap["vs_embedPt"]->Fill(result, tau->alternatLorentzVect().pt());
       dynamic_cast<TH3F*>(mymap["vs_pt_jetPt"])->Fill(
           result, tau->pt(), tau->jetRef()->pt());
       dynamic_cast<TH3F*>(mymap["vs_pt_embedPt"])->Fill(
           result, tau->pt(), tau->alternatLorentzVect().pt());
       mymap["vs_eta"]->Fill(result, tau->eta());
       mymap["vs_dm"]->Fill(result, tau->decayMode());
       if (plotPU_ && tau->pt() > pileupPtCut_) {
         if (puInfo.isValid())
           mymap["vs_truePU"]->Fill(result, puInfo->getPU_NumInteractions());
         mymap["vs_recoPU"]->Fill(result, puVertices->size());
       }
     }
   }
 }