#include <B2GDQM.h>

Inheritance diagram for B2GDQM:

Public Member Functions
	B2GDQM (const edm::ParameterSet &ps)

virtual	~B2GDQM ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
virtual void	analyze (edm::Event const &e, edm::EventSetup const &eSetup)

virtual void	analyzeEventInterpretation (edm::Event const &e, edm::EventSetup const &eSetup)

virtual void	beginJob ()

virtual void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)

virtual void	beginRun (edm::Run const &run, edm::EventSetup const &eSetup)

virtual void	endJob ()

virtual void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &c)

virtual void	endRun (edm::Run const &run, edm::EventSetup const &eSetup)

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)

Private Member Functions
void	bookHistos (DQMStore *bei)

Private Attributes
DQMStore *	bei_

std::vector< MonitorElement * >	boostedJet_massDrop

std::vector< MonitorElement * >	boostedJet_minMass

std::vector< MonitorElement * >	boostedJet_subjetM

std::vector< MonitorElement * >	boostedJet_subjetN

std::vector< MonitorElement * >	boostedJet_subjetPhi

std::vector< MonitorElement * >	boostedJet_subjetPt

std::vector< MonitorElement * >	boostedJet_subjetY

HLTConfigProvider	hltConfigProvider_

int	ievt

int	irun

bool	isValidHltConfig_

std::vector< edm::InputTag >	jetLabels_

std::vector< double >	jetPtMins_

std::vector< edm::EDGetTokenT < edm::View< reco::Jet > > >	jetTokens_

int	nEvents_

int	nLumiSecs_

std::vector< MonitorElement * >	pfJet_cemf

std::vector< MonitorElement * >	pfJet_chef

std::vector< MonitorElement * >	pfJet_m

std::vector< MonitorElement * >	pfJet_nemf

std::vector< MonitorElement * >	pfJet_nhef

std::vector< MonitorElement * >	pfJet_phi

std::vector< MonitorElement * >	pfJet_pt

std::vector< MonitorElement * >	pfJet_y

std::string	PFJetCorService_

MonitorElement *	pfMet_phi

MonitorElement *	pfMet_pt

edm::InputTag	PFMETLabel_

edm::EDGetTokenT< std::vector < reco::PFMET > >	PFMETToken_

edm::InputTag	theTriggerResultsCollection

edm::Handle< edm::TriggerResults >	triggerResults_

edm::EDGetTokenT < edm::TriggerResults >	triggerToken_

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 86 of file B2GDQM.h.

Constructor & Destructor Documentation

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

Definition at line 83 of file B2GDQM.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                        {
 
   edm::LogInfo("B2GDQM") <<  " Starting B2GDQM " << "\n" ;
 
 
   typedef std::vector<edm::InputTag> vtag;
 
   // Get parameters from configuration file
   // Trigger
   theTriggerResultsCollection = ps.getParameter<InputTag>("triggerResultsCollection");
   triggerToken_ = consumes< edm::TriggerResults> ( theTriggerResultsCollection );
 
   // Jets
   jetLabels_ = ps.getParameter<std::vector<edm::InputTag> >("jetLabels");
   for ( std::vector<edm::InputTag>::const_iterator jetlabel = jetLabels_.begin(),
       jetlabelEnd = jetLabels_.end(); jetlabel != jetlabelEnd; ++jetlabel ) {
     jetTokens_.push_back( consumes<edm::View<reco::Jet> >( *jetlabel ) );
   }
   jetPtMins_ = ps.getParameter<std::vector<double> > ("jetPtMins");
   PFJetCorService_ = ps.getParameter<std::string>("PFJetCorService");
 
   // MET
   PFMETLabel_         = ps.getParameter<InputTag>("pfMETCollection");  
   PFMETToken_         = consumes<std::vector<reco::PFMET> > ( PFMETLabel_ );
 
 
 }

B2GDQM::~B2GDQM ( )

virtual

Definition at line 115 of file B2GDQM.cc.

                {
   edm::LogInfo("B2GDQM") <<  " Deleting B2GDQM " << "\n" ;
 }

Member Function Documentation

void B2GDQM::analyze	(	edm::Event const &	e,
		edm::EventSetup const &	eSetup
	)

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 205 of file B2GDQM.cc.

                                                                        {
 
   analyzeEventInterpretation(iEvent, iSetup);
 }

void B2GDQM::analyzeEventInterpretation	(	edm::Event const &	e,
		edm::EventSetup const &	eSetup
	)

protectedvirtual

Definition at line 211 of file B2GDQM.cc.

References edm::View< T >::begin(), reco::PFJet::chargedEmEnergyFraction(), reco::PFJet::chargedHadronEnergyFraction(), edm::View< T >::end(), edm::Event::getByToken(), icoll, metsig::jet, diffTwoXMLs::label, reco::Candidate::mass(), reco::CATopJetProperties::minMass, reco::PFJet::neutralEmEnergyFraction(), reco::PFJet::neutralHadronEnergyFraction(), phi, reco::Candidate::phi(), RecoTauCleanerPlugins::pt, reco::Candidate::pt(), and reco::Candidate::rapidity().

                                                                                         {
 
 
   // Loop over the different types of jets, 
   //   Loop over the jets in that collection,
   //     fill PF jet information as well as substructure
   //     information for boosted jets. 
   // Utilizes the CMS top-tagging algorithm and the "mass drop" W-tagger. 
   for ( unsigned int icoll = 0; icoll < jetLabels_.size(); ++icoll ) {
 
     edm::Handle<edm::View<reco::Jet> > pfJetCollection;
     bool ValidPFJets = iEvent.getByToken(jetTokens_[icoll], pfJetCollection );
     if(!ValidPFJets) continue;
     edm::View<reco::Jet> const & pfjets = *pfJetCollection;
     
     
 
     // Jet Correction
     int countJet = 0;
     //const JetCorrector* pfcorrector = JetCorrector::getJetCorrector(PFJetCorService_,iSetup);
     
     for ( edm::View<reco::Jet>::const_iterator jet = pfjets.begin(),
         jetEnd = pfjets.end(), jetBegin = pfjets.begin();
       jet != jetEnd; ++jet ) {
       if(jet->pt()<jetPtMins_[icoll]) continue;
       pfJet_pt  [icoll]->Fill(jet->pt());
       pfJet_y   [icoll]->Fill(jet->rapidity());
       pfJet_phi [icoll]->Fill(jet->phi());
       pfJet_m   [icoll]->Fill(jet->mass());
 
       // Dynamic cast the base class (reco::Jet) to the derived class (PFJet)
       // to access the PFJet information
       reco::PFJet const * pfjet = dynamic_cast<reco::PFJet const *>( &*jet);
 
       if ( pfjet != 0 ) {
     pfJet_chef[icoll]->Fill(pfjet->chargedHadronEnergyFraction());
     pfJet_nhef[icoll]->Fill(pfjet->neutralHadronEnergyFraction());
     pfJet_cemf[icoll]->Fill(pfjet->chargedEmEnergyFraction());
     pfJet_nemf[icoll]->Fill(pfjet->neutralEmEnergyFraction());
       }
 
       // Dynamic cast the base class (reco::Jet) to the derived class (BasicJet)
       // to access the substructure information
       reco::BasicJet const * basicjet = dynamic_cast<reco::BasicJet const *>( &*jet);
       
       if ( basicjet != 0 ) {
     boostedJet_subjetN[icoll]->Fill( jet->numberOfDaughters() );
 
     for ( unsigned int ida = 0; ida < jet->numberOfDaughters(); ++ida ) {
       reco::Candidate const * subjet = jet->daughter(ida);
       boostedJet_subjetPt [icoll]->Fill ( subjet->pt() );
       boostedJet_subjetY  [icoll]->Fill ( subjet->rapidity() );
       boostedJet_subjetPhi[icoll]->Fill ( subjet->phi() );
       boostedJet_subjetM  [icoll]->Fill ( subjet->mass() );
     }
     // Check the various tagging algorithms
 
     // For top-tagging, check the minimum mass pairing
     if ( jetLabels_[icoll].label() == "cmsTopTagPFJetsCHS") {
       CATopJetHelper helper(171.2, 80.4);
       reco::CATopJetProperties properties = helper(*basicjet);
       if ( jet->numberOfDaughters() > 2 ) {
         boostedJet_minMass[icoll]->Fill ( properties.minMass );     
       } else {
         boostedJet_minMass[icoll]->Fill ( -1.0 );       
       }
       
     // For W-tagging, check the mass drop
     } else if ( jetLabels_[icoll].label() == "ca8PFJetsCHSPruned" ) {
       if ( jet->numberOfDaughters() > 1 ) {
         reco::Candidate const * da0 = jet->daughter(0);
         reco::Candidate const * da1 = jet->daughter(1);
         if ( da0->mass() > da1->mass() ) {
           boostedJet_massDrop[icoll]->Fill( da0->mass() / jet->mass() );
         } else {
           boostedJet_massDrop[icoll]->Fill( da1->mass() / jet->mass() );
         }
       } else {
         boostedJet_massDrop[icoll]->Fill( -1.0 );
       }
 
     } // end if collection is CA8 PFJets CHS Pruned
     
       } // end if basic jet != 0
       countJet++;
     }
   }
 
   
 
   // PFMETs
   edm::Handle<std::vector<reco::PFMET> > pfMETCollection;
   bool ValidPFMET = iEvent.getByToken(PFMETToken_, pfMETCollection);
   if(!ValidPFMET) return;
 
   pfMet_pt->Fill( (*pfMETCollection)[0].pt() );
   pfMet_phi->Fill( (*pfMETCollection)[0].phi() );
 
 }

void B2GDQM::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 123 of file B2GDQM.cc.

                      {
   nLumiSecs_ = 0;
   nEvents_   = 0;
 }

void B2GDQM::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	context
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 154 of file B2GDQM.cc.

                                                 { 
   //edm::LogInfo ("B2GDQM") <<"[B2GDQM]: Begin of LS transition";
 }

void B2GDQM::beginRun	(	edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 132 of file B2GDQM.cc.

References bookHistos(), cppFunctionSkipper::operator, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                  {
   edm::LogInfo ("B2GDQM") <<"[B2GDQM]: Begining of Run";
 
 
   bei_ = Service<DQMStore>().operator->();
   bei_->setCurrentFolder("Physics/B2G");
   bookHistos(bei_);
   
   // passed as parameter to HLTConfigProvider::init(), not yet used
   bool isConfigChanged = false;
   
   // isValidHltConfig_ used to short-circuit analyze() in case of problems
   //  const std::string hltProcessName( "HLT" );
   const std::string hltProcessName = theTriggerResultsCollection.process();
   isValidHltConfig_ = hltConfigProvider_.init( run, eSetup, hltProcessName, isConfigChanged );
 
 }

void B2GDQM::bookHistos ( DQMStore * bei )

private

Definition at line 163 of file B2GDQM.cc.

References DQMStore::book1D(), DQMStore::cd(), icoll, and DQMStore::setCurrentFolder().

                                     {
  
   bei->cd();
   
   //--- Event Interpretation
 
   for ( unsigned int icoll = 0; icoll < jetLabels_.size(); ++icoll ) {
     std::stringstream ss;
     ss << "Physics/B2G/" << jetLabels_[icoll].label();
     bei->setCurrentFolder(ss.str().c_str());
     pfJet_pt          .push_back( bei->book1D("pfJet_pt",     "Pt of PFJet (GeV)",      50, 0.0, 1000) );
     pfJet_y           .push_back( bei->book1D("pfJet_y",      "Rapidity of PFJet",      60, -6.0, 6.0) );
     pfJet_phi         .push_back( bei->book1D("pfJet_phi",    "#phi of PFJet (radians)",60, -3.14159, 3.14159) );
     pfJet_m           .push_back( bei->book1D("pfJet_m",      "Mass of PFJet (GeV)",    50, 0.0, 500) );
     pfJet_chef        .push_back( bei->book1D("pfJet_pfchef", "PFJetID CHEF", 50, 0.0 , 1.0));
     pfJet_nhef        .push_back( bei->book1D("pfJet_pfnhef", "PFJetID NHEF", 50, 0.0 , 1.0));
     pfJet_cemf        .push_back( bei->book1D("pfJet_pfcemf", "PFJetID CEMF", 50, 0.0 , 1.0));
     pfJet_nemf        .push_back( bei->book1D("pfJet_pfnemf", "PFJetID NEMF", 50, 0.0 , 1.0));
 
     boostedJet_subjetPt       .push_back( bei->book1D("boostedJet_subjetPt",  "Pt of subjets (GeV)", 50, 0.0 , 500));
     boostedJet_subjetY        .push_back( bei->book1D("boostedJet_subjetY",   "Rapidity of subjets", 60, -6.0, 6.0));
     boostedJet_subjetPhi      .push_back( bei->book1D("boostedJet_subjetPhi", "#phi of subjets (radians)", 60, -3.14159, 3.14159));
     boostedJet_subjetM        .push_back( bei->book1D("boostedJet_subjetM",   "Mass of subjets (GeV)", 50, 0.0 , 250.));
     boostedJet_subjetN        .push_back( bei->book1D("boostedJet_subjetN",   "Number of subjets", 10, 0, 10));
     boostedJet_massDrop       .push_back( bei->book1D("boostedJet_massDrop", "Mass drop for W-like jets", 50, 0.0 , 1.0));
     boostedJet_minMass        .push_back( bei->book1D("boostedJet_minMass",   "Minimum Mass Pairing for top-like jets", 50, 0.0 , 250.0));
 
   }
 
   bei->setCurrentFolder("Physics/B2G/MET");
   pfMet_pt                 = bei->book1D("pfMet_pt", "Pf Missing p_{T}; GeV", 50,  0.0 , 500);
   pfMet_phi                = bei->book1D("pfMet_phi", "Pf Missing p_{T} #phi;#phi (radians)", 35, -3.5, 3.5 );
 
 
   
   bei->cd();
 }

void B2GDQM::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 332 of file B2GDQM.cc.

                    {
   //edm::LogInfo("B2GDQM") <<"[B2GDQM]: endjob called!";
 }

void B2GDQM::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	c
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 313 of file B2GDQM.cc.

                                                                                               {
   //edm::LogInfo ("B2GDQM") <<"[B2GDQM]: End of LS transition, performing the DQM client operation";
   nLumiSecs_++;
   //edm::LogInfo("B2GDQM") << "============================================ " 
   //<< endl << " ===> Iteration # " << nLumiSecs_ << " " << lumiSeg.luminosityBlock() 
   //<< endl  << "============================================ " << endl;
 }

void B2GDQM::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 325 of file B2GDQM.cc.

325 {

326 }

Member Data Documentation

DQMStore* B2GDQM::bei_

private

Definition at line 118 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::boostedJet_massDrop

private

Definition at line 160 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::boostedJet_minMass

private

Definition at line 161 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetM

private

Definition at line 158 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetN

private

Definition at line 159 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetPhi

private

Definition at line 157 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetPt

private

Definition at line 155 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetY

private

Definition at line 156 of file B2GDQM.h.

HLTConfigProvider B2GDQM::hltConfigProvider_

private

Definition at line 119 of file B2GDQM.h.

int B2GDQM::ievt

private

Definition at line 116 of file B2GDQM.h.

int B2GDQM::irun

private

Definition at line 116 of file B2GDQM.h.

bool B2GDQM::isValidHltConfig_

private

Definition at line 120 of file B2GDQM.h.

std::vector<edm::InputTag> B2GDQM::jetLabels_

private

Definition at line 131 of file B2GDQM.h.

std::vector<double> B2GDQM::jetPtMins_

private

Definition at line 141 of file B2GDQM.h.

std::vector< edm::EDGetTokenT< edm::View<reco::Jet> > > B2GDQM::jetTokens_

private

Definition at line 132 of file B2GDQM.h.

int B2GDQM::nEvents_

private

Definition at line 116 of file B2GDQM.h.

int B2GDQM::nLumiSecs_

private

Definition at line 115 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::pfJet_cemf

private

Definition at line 153 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::pfJet_chef

private

Definition at line 151 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::pfJet_m

private

Definition at line 150 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::pfJet_nemf

private

Definition at line 154 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::pfJet_nhef

private

Definition at line 152 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::pfJet_phi

private

Definition at line 149 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::pfJet_pt

private

Definition at line 147 of file B2GDQM.h.

std::vector<MonitorElement*> B2GDQM::pfJet_y

private

Definition at line 148 of file B2GDQM.h.

std::string B2GDQM::PFJetCorService_

private

Definition at line 143 of file B2GDQM.h.

MonitorElement* B2GDQM::pfMet_phi

private

Definition at line 163 of file B2GDQM.h.

MonitorElement* B2GDQM::pfMet_pt

private

Definition at line 162 of file B2GDQM.h.

edm::InputTag B2GDQM::PFMETLabel_

private

Definition at line 133 of file B2GDQM.h.

edm::EDGetTokenT< std::vector<reco::PFMET> > B2GDQM::PFMETToken_

private

Definition at line 134 of file B2GDQM.h.

edm::InputTag B2GDQM::theTriggerResultsCollection

private

Definition at line 123 of file B2GDQM.h.

edm::Handle<edm::TriggerResults> B2GDQM::triggerResults_

private

Definition at line 126 of file B2GDQM.h.

edm::EDGetTokenT<edm::TriggerResults> B2GDQM::triggerToken_

private

Definition at line 124 of file B2GDQM.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation