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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

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

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 55 of file TTbarSpinCorrHepMCAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file TTbarSpinCorrHepMCAnalyzer.cc.

References genEventInfoProductTag_, genEventInfoProductTagToken_, genParticlesTag_, and genParticlesTagToken_.

                                                                                     :
  genEventInfoProductTag_(iConfig.getParameter<edm::InputTag>("genEventInfoProductTag")),
  genParticlesTag_(iConfig.getParameter<edm::InputTag>("genParticlesTag"))
{

  genEventInfoProductTagToken_=consumes<GenEventInfoProduct>(genEventInfoProductTag_);
  genParticlesTagToken_=consumes<reco::GenParticleCollection>(genParticlesTag_);

}

TTbarSpinCorrHepMCAnalyzer::~TTbarSpinCorrHepMCAnalyzer ( )

override

Definition at line 17 of file TTbarSpinCorrHepMCAnalyzer.cc.

{}

Member Function Documentation

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

override

Definition at line 24 of file TTbarSpinCorrHepMCAnalyzer.cc.

References _h_asym, _h_deltaPhi, _h_llpairM, _h_llpairPt, hiPixelPairStep_cff::deltaPhi, edm::HandleBase::failedToGet(), MonitorElement::Fill(), genEventInfoProductTagToken_, GenHFHadronMatcher_cfi::genParticles, genParticlesTagToken_, edm::Event::getByToken(), mps_fire::i, nEvt, AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::p4(), reco::LeafCandidate::pdgId(), reco::LeafCandidate::status(), and GenEventInfoProduct::weight().

                                                                                           {
  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) ;
  }
}

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

override

Setting the DQM top directories

Definition at line 89 of file TTbarSpinCorrHepMCAnalyzer.cc.

References _h_asym, _h_deltaPhi, _h_llpairM, _h_llpairPt, DQMHelper::book1dHisto(), dir, nEvt, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                        {
  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");
}

Member Data Documentation

MonitorElement* TTbarSpinCorrHepMCAnalyzer::_h_asym

private

Definition at line 68 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TTbarSpinCorrHepMCAnalyzer::_h_deltaPhi

private

Definition at line 69 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TTbarSpinCorrHepMCAnalyzer::_h_llpairM

private

Definition at line 72 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TTbarSpinCorrHepMCAnalyzer::_h_llpairPt

private

Definition at line 71 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag TTbarSpinCorrHepMCAnalyzer::genEventInfoProductTag_

private

Definition at line 74 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by TTbarSpinCorrHepMCAnalyzer().

edm::EDGetTokenT<GenEventInfoProduct> TTbarSpinCorrHepMCAnalyzer::genEventInfoProductTagToken_

private

Definition at line 76 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and TTbarSpinCorrHepMCAnalyzer().

edm::InputTag TTbarSpinCorrHepMCAnalyzer::genParticlesTag_

private

Definition at line 74 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by TTbarSpinCorrHepMCAnalyzer().

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

private

Definition at line 77 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and TTbarSpinCorrHepMCAnalyzer().

MonitorElement* TTbarSpinCorrHepMCAnalyzer::nEvt

private

Definition at line 67 of file TTbarSpinCorrHepMCAnalyzer.h.

Referenced by analyze(), and bookHistograms().

double TTbarSpinCorrHepMCAnalyzer::weight

private

Definition at line 65 of file TTbarSpinCorrHepMCAnalyzer.h.