B2GDQM Class Reference

#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
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) 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();
      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, DQMStore::setCurrentFolder(), and contentValuesCheck::ss.

                                    {
 
  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.

                                                                 {
}

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.