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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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