TTbarSpinCorrHepMCAnalyzer Class Reference

#include <TTbarSpinCorrHepMCAnalyzer.h>

Inheritance diagram for TTbarSpinCorrHepMCAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
	TTbarSpinCorrHepMCAnalyzer (const edm::ParameterSet &)
	~TTbarSpinCorrHepMCAnalyzer () 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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
MonitorElement *	_h_asym
MonitorElement *	_h_deltaPhi
MonitorElement *	_h_llpairM
MonitorElement *	_h_llpairPt
edm::InputTag	genEventInfoProductTag_
edm::EDGetTokenT< GenEventInfoProduct >	genEventInfoProductTagToken_
edm::InputTag	genParticlesTag_
edm::EDGetTokenT< reco::GenParticleCollection >	genParticlesTagToken_
MonitorElement *	nEvt
double	weight

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 >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 53 of file TTbarSpinCorrHepMCAnalyzer.h.

Constructor & Destructor Documentation

TTbarSpinCorrHepMCAnalyzer()

TTbarSpinCorrHepMCAnalyzer::TTbarSpinCorrHepMCAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 6 of file TTbarSpinCorrHepMCAnalyzer.cc.

    : genEventInfoProductTag_(iConfig.getParameter<edm::InputTag>("genEventInfoProductTag")),
      genParticlesTag_(iConfig.getParameter<edm::InputTag>("genParticlesTag")) {
  genEventInfoProductTagToken_ = consumes<GenEventInfoProduct>(genEventInfoProductTag_);
  genParticlesTagToken_ = consumes<reco::GenParticleCollection>(genParticlesTag_);
}

References genEventInfoProductTag_, genEventInfoProductTagToken_, genParticlesTag_, and genParticlesTagToken_.

~TTbarSpinCorrHepMCAnalyzer()

TTbarSpinCorrHepMCAnalyzer::~TTbarSpinCorrHepMCAnalyzer ( )

override

Definition at line 13 of file TTbarSpinCorrHepMCAnalyzer.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 20 of file TTbarSpinCorrHepMCAnalyzer.cc.

                                                                                            {
  using namespace edm;
 
  // --- the MC weights ---
  Handle<GenEventInfoProduct> evt_info;
  iEvent.getByToken(genEventInfoProductTagToken_, evt_info);
  if (evt_info.failedToGet())
    return;
 
  weight = evt_info->weight();
 
  // --- get genParticles ---
  Handle<reco::GenParticleCollection> genParticles;
  iEvent.getByToken(genParticlesTagToken_, genParticles);
 
  const reco::GenParticle* _lepton(nullptr);
  const reco::GenParticle* _leptonBar(nullptr);
 
  bool hasTop(false), hasTopbar(false);
  for (size_t i = 0; i < genParticles->size(); ++i) {
    const reco::GenParticle& p = (*genParticles)[i];
    if (p.pdgId() == 6)
      hasTop = true;
    if (p.pdgId() == -6)
      hasTopbar = true;
  }
 
  if (hasTop && hasTopbar) {
    // --- get status 3 leptons
    for (size_t i = 0; i < genParticles->size(); ++i) {
      const reco::GenParticle& p = (*genParticles)[i];
      if ((p.pdgId() == 11 || p.pdgId() == 13 || p.pdgId() == 15) && p.status() == 3) {
        _lepton = &p;
      }
      if ((p.pdgId() == -11 || p.pdgId() == -13 || p.pdgId() == -15) && p.status() == 3) {
        _leptonBar = &p;
      }
 
      if (_lepton && _leptonBar)
        break;
    }
 
    if (_lepton && _leptonBar) {
      const math::XYZTLorentzVector& lepton = _lepton->p4();
      const math::XYZTLorentzVector& leptonBar = _leptonBar->p4();
 
      double deltaPhi = fabs(TVector2::Phi_mpi_pi(lepton.phi() - leptonBar.phi()));
      _h_deltaPhi->Fill(deltaPhi, weight);
 
      double asym = (deltaPhi > CLHEP::halfpi) ? 0.5 : -0.5;
      _h_asym->Fill(asym, weight);
 
      math::XYZTLorentzVector llpair = lepton + leptonBar;
 
      double llpairPt = llpair.pt();
      _h_llpairPt->Fill(llpairPt, weight);
 
      double llpairM = llpair.M();
      _h_llpairM->Fill(llpairM, weight);
    }
    nEvt->Fill(0.5, weight);
  }
}

References _h_asym, _h_deltaPhi, _h_llpairM, _h_llpairPt, SiPixelRawToDigiRegional_cfi::deltaPhi, edm::HandleBase::failedToGet(), dqm::impl::MonitorElement::Fill(), genEventInfoProductTagToken_, genParticles2HepMC_cfi::genParticles, genParticlesTagToken_, mps_fire::i, iEvent, nEvt, AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::p4(), and GenEventInfoProduct::weight().

bookHistograms()

void TTbarSpinCorrHepMCAnalyzer::bookHistograms ( DQMStore::IBooker &  i, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 85 of file TTbarSpinCorrHepMCAnalyzer.cc.

                                                                                                       {
  std::string dir = "Generator/";
  dir += "TTbarSpinCorr";
  DQMHelper dqm(&i);
  i.setCurrentFolder(dir);
 
  // Number of analyzed events
  nEvt = dqm.book1dHisto("nEvt", "n analyzed Events", 1, 0., 1., "bin", "Number of Events");
  _h_asym = dqm.book1dHisto("TTbar_asym", "Asymmetr", 2, -1., 1., "Asymmetry", "Number of Events");
  _h_deltaPhi =
      dqm.book1dHisto("TTbar_deltaPhi", "#Delta#phi(ll)", 320, 0, 3.2, "#Delta#phi_{ll} (rad)", "Number of Events");
  _h_llpairPt = dqm.book1dHisto(
      "TTbar_llpairPt", "Lepton pair transverse momentum", 1000, 0, 1000, "p_{T}^{ll} (GeV)", "Number of Events");
  _h_llpairM =
      dqm.book1dHisto("TTbar_llpairM", "Lepton pair invariant mass", 1000, 0, 1000, "M_{ll} (GeV)", "Number of Events");
}

References _h_asym, _h_deltaPhi, _h_llpairM, _h_llpairPt, DeadROC_duringRun::dir, mps_fire::i, nEvt, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Data Documentation

_h_asym

MonitorElement* TTbarSpinCorrHepMCAnalyzer::_h_asym

private

Definition at line 66 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

_h_deltaPhi

MonitorElement* TTbarSpinCorrHepMCAnalyzer::_h_deltaPhi

private

Definition at line 67 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

_h_llpairM

MonitorElement* TTbarSpinCorrHepMCAnalyzer::_h_llpairM

private

Definition at line 70 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

_h_llpairPt

MonitorElement* TTbarSpinCorrHepMCAnalyzer::_h_llpairPt

private

Definition at line 69 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

genEventInfoProductTag_

edm::InputTag TTbarSpinCorrHepMCAnalyzer::genEventInfoProductTag_

private

Definition at line 72 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by TTbarSpinCorrHepMCAnalyzer().

genEventInfoProductTagToken_

edm::EDGetTokenT<GenEventInfoProduct> TTbarSpinCorrHepMCAnalyzer::genEventInfoProductTagToken_

private

Definition at line 74 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and TTbarSpinCorrHepMCAnalyzer().

genParticlesTag_

edm::InputTag TTbarSpinCorrHepMCAnalyzer::genParticlesTag_

private

Definition at line 72 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by TTbarSpinCorrHepMCAnalyzer().

genParticlesTagToken_

edm::EDGetTokenT<reco::GenParticleCollection> TTbarSpinCorrHepMCAnalyzer::genParticlesTagToken_

private

Definition at line 75 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and TTbarSpinCorrHepMCAnalyzer().

nEvt

MonitorElement* TTbarSpinCorrHepMCAnalyzer::nEvt

private

Definition at line 65 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

weight

double TTbarSpinCorrHepMCAnalyzer::weight

private

Definition at line 63 of file TTbarSpinCorrHepMCAnalyzer.h.