#include <TTbar_Kinematics.h>

Inheritance diagram for TTbar_Kinematics:

Public Member Functions
	TTbar_Kinematics (const edm::ParameterSet &)
	~TTbar_Kinematics ()
Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	beginRun (edm::Run const &, edm::EventSetup const &)
virtual void	endJob ()
virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	endRun (edm::Run const &, edm::EventSetup const &)
Private Attributes
DQMStore *	dbe
	ME's "container".
edm::InputTag	genEventInfoProductTag_
MonitorElement *	hBottomE
MonitorElement *	hBottomEta
MonitorElement *	hBottomEtaPt
MonitorElement *	hBottomEtaPz
MonitorElement *	hBottomMass
MonitorElement *	hBottomMassDeltaY
MonitorElement *	hBottomMassEta
MonitorElement *	hBottomMassPz
MonitorElement *	hBottomMassY
MonitorElement *	hBottomPt
MonitorElement *	hBottomPtPz
MonitorElement *	hBottomPz
MonitorElement *	hBottomY
MonitorElement *	hBottomYPz
edm::InputTag	hepmcCollection_
MonitorElement *	hTopMass
MonitorElement *	hTopPt
MonitorElement *	hTopY
MonitorElement *	hTTbarMass
MonitorElement *	hTTbarPt
MonitorElement *	hTTbarY
MonitorElement *	hWminPz
MonitorElement *	hWplusPz
MonitorElement *	nEvt
TLorentzVector	tlv_Bottom
TLorentzVector	tlv_BottomBar
TLorentzVector	tlv_Top
TLorentzVector	tlv_TopBar
TLorentzVector	tlv_TTbar
TLorentzVector	tlv_Wmin
TLorentzVector	tlv_Wplus
double	weight

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 59 of file TTbar_Kinematics.h.

Constructor & Destructor Documentation

TTbar_Kinematics::TTbar_Kinematics ( const edm::ParameterSet & iConfig ) [explicit]

Definition at line 8 of file TTbar_Kinematics.cc.

References dbe, and cppFunctionSkipper::operator.

                                                                 :
  genEventInfoProductTag_(iConfig.getParameter<edm::InputTag>("genEventInfoProductTag"))
{
  dbe = 0;
  dbe = edm::Service<DQMStore>().operator->();
}

TTbar_Kinematics::~TTbar_Kinematics ( )

Definition at line 16 of file TTbar_Kinematics.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void TTbar_Kinematics::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		`[private, virtual]`

Implements edm::EDAnalyzer.

Definition at line 31 of file TTbar_Kinematics.cc.

References MonitorElement::Fill(), genEventInfoProductTag_, TtGenEvtProducer_cfi::genEvt, edm::Event::getByLabel(), hBottomE, hBottomEta, hBottomEtaPt, hBottomEtaPz, hBottomMass, hBottomMassDeltaY, hBottomMassEta, hBottomMassPz, hBottomMassY, hBottomPt, hBottomPtPz, hBottomPz, hBottomY, hBottomYPz, hTopMass, hTopPt, hTopY, hTTbarMass, hTTbarPt, hTTbarY, hWminPz, hWplusPz, edm::HandleBase::isValid(), nEvt, reco::LeafCandidate::p4(), tlv_Bottom, tlv_BottomBar, tlv_Top, tlv_TopBar, tlv_TTbar, tlv_Wmin, tlv_Wplus, and weight.

{

  // --- the MC weights ---
  Handle<GenEventInfoProduct> evt_info;
  iEvent.getByLabel(genEventInfoProductTag_, evt_info);
  if(!evt_info.isValid()) return;
  weight = evt_info->weight() ;

  // --- get TopQuarkAnalysis TtGenEvent
  Handle<TtGenEvent> genEvt;
  iEvent.getByLabel("genEvt", genEvt);

  if(!genEvt.isValid())return;

  const reco::GenParticle*        top = 0;
  const reco::GenParticle*    antitop = 0;
  const reco::GenParticle*     bottom = 0;
  const reco::GenParticle* antibottom = 0;
  const reco::GenParticle*      Wplus = 0;
  const reco::GenParticle*      Wmin  = 0;

  top        = genEvt->top();
  antitop    = genEvt->topBar();
  bottom     = genEvt->b();
  antibottom = genEvt->bBar();
  Wplus      = genEvt->wPlus();
  Wmin       = genEvt->wMinus();

  tlv_Top        = TLorentzVector(0,0,0,0) ;
  tlv_TopBar     = TLorentzVector(0,0,0,0) ;
  tlv_Bottom     = TLorentzVector(0,0,0,0) ;
  tlv_BottomBar  = TLorentzVector(0,0,0,0) ;
  tlv_Wplus      = TLorentzVector(0,0,0,0) ;
  tlv_Wmin       = TLorentzVector(0,0,0,0) ;
  tlv_TTbar      = TLorentzVector(0,0,0,0);

  if(top)        tlv_Top.SetPxPyPzE(top->p4().px(),top->p4().py(),top->p4().pz(),top->p4().e());
  if(antitop)    tlv_TopBar.SetPxPyPzE(antitop->p4().px(),antitop->p4().py(),antitop->p4().pz(),antitop->p4().e());
  if(bottom)     tlv_Bottom.SetPxPyPzE(bottom->p4().px(),bottom->p4().py(),bottom->p4().pz(),bottom->p4().e());
  if(antibottom) tlv_BottomBar.SetPxPyPzE(antibottom->p4().px(),antibottom->p4().py(),antibottom->p4().pz(),antibottom->p4().e());
  if(Wplus)      tlv_Wplus.SetPxPyPzE(Wplus->p4().px(),Wplus->p4().py(),Wplus->p4().pz(),Wplus->p4().e());
  if(Wmin)       tlv_Wmin.SetPxPyPzE(Wmin->p4().px(),Wmin->p4().py(),Wmin->p4().pz(),Wmin->p4().e());
  tlv_TTbar = tlv_Top + tlv_TopBar ;

  //---topquarkquantities---
  nEvt->Fill(0.5,weight);
  if(top && antitop){
    hTopPt->Fill(tlv_Top.Pt(),weight);
    hTopPt->Fill(tlv_TopBar.Pt(),weight);
    
    hTopY->Fill(tlv_Top.Rapidity(),weight);
    hTopY->Fill(tlv_TopBar.Rapidity(),weight);
    
    hTopMass->Fill(tlv_Top.M(),weight);
    hTopMass->Fill(tlv_TopBar.M(),weight);
    
    //---ttbarpairquantities---
    hTTbarPt->Fill(tlv_TTbar.Pt(),weight);
    hTTbarPt->Fill(tlv_TTbar.Pt(),weight);
    
    hTTbarY->Fill(tlv_TTbar.Rapidity(),weight);
    hTTbarY->Fill(tlv_TTbar.Rapidity(),weight);
    
    hTTbarMass->Fill(tlv_TTbar.M(),weight);
    hTTbarMass->Fill(tlv_TTbar.M(),weight);
  }
  if(bottom && antibottom){
    hBottomPt->Fill(tlv_Bottom.Pt(),weight);
    hBottomPt->Fill(tlv_BottomBar.Pt(),weight);
    
    hBottomEta->Fill(tlv_Bottom.Eta(),weight);
    hBottomEta->Fill(tlv_BottomBar.Eta(),weight);
    
    //hBottomY->Fill(math::XYZTLorentzVector(bottom->momentum()).Rapidity(),weight);
    //hBottomY->Fill(math::XYZTLorentzVector(antibottom->momentum()).Rapidity(),weight);
    
    hBottomY->Fill(tlv_Bottom.Rapidity(),weight);
    hBottomY->Fill(tlv_BottomBar.Rapidity(),weight);
    
    hBottomPz->Fill(tlv_Bottom.Pz(),weight);
    hBottomPz->Fill(tlv_BottomBar.Pz(),weight);
    
    hBottomE->Fill(tlv_Bottom.E(),weight);
    hBottomE->Fill(tlv_BottomBar.E(),weight);
    
    hBottomMass->Fill(tlv_Bottom.M(),weight);
    hBottomMass->Fill(tlv_BottomBar.M(),weight);
    
    hBottomPtPz->Fill(tlv_Bottom.Pt(),tlv_Bottom.Pz(),weight);
    hBottomPtPz->Fill(tlv_BottomBar.Pt(),tlv_BottomBar.Pz(),weight);
    
    hBottomEtaPz->Fill(tlv_Bottom.Eta(),tlv_Bottom.Pz(),weight);
    hBottomEtaPz->Fill(tlv_BottomBar.Eta(),tlv_BottomBar.Pz(),weight);
    
    hBottomEtaPt->Fill(tlv_Bottom.Eta(),tlv_Bottom.Pt(),weight);
    hBottomEtaPt->Fill(tlv_BottomBar.Eta(),tlv_BottomBar.Pt(),weight);
    
    hBottomYPz->Fill(tlv_Bottom.Rapidity(),tlv_Bottom.Pz(),weight);
    hBottomYPz->Fill(tlv_BottomBar.Rapidity(),tlv_BottomBar.Pz(),weight);
    
    hBottomMassPz->Fill(tlv_Bottom.M(),tlv_Bottom.Pz(),weight);
    hBottomMassPz->Fill(tlv_BottomBar.M(),tlv_BottomBar.Pz(),weight);
    
    hBottomMassEta->Fill(tlv_Bottom.M(),tlv_Bottom.Eta(),weight);
    hBottomMassEta->Fill(tlv_BottomBar.M(),tlv_BottomBar.Eta(),weight);
    
    hBottomMassY->Fill(tlv_Bottom.M(),tlv_Bottom.Rapidity(),weight);
    hBottomMassY->Fill(tlv_BottomBar.M(),tlv_BottomBar.Rapidity(),weight);
    
    hBottomMassDeltaY->Fill(tlv_Bottom.M(),tlv_Bottom.Eta()-tlv_Bottom.Rapidity(),weight);
    hBottomMassDeltaY->Fill(tlv_BottomBar.M(),tlv_BottomBar.Eta()-tlv_BottomBar.Rapidity(),weight);
  }
  if(Wplus && Wmin){
    hWplusPz->Fill(tlv_Wplus.Pz(),weight);
    hWminPz->Fill(tlv_Wmin.Pz(),weight);
  }
}

void TTbar_Kinematics::beginJob ( void ) [private, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 153 of file TTbar_Kinematics.cc.

References DQMStore::book1D(), DQMStore::book2D(), dbe, hBottomE, hBottomEta, hBottomEtaPt, hBottomEtaPz, hBottomMass, hBottomMassDeltaY, hBottomMassEta, hBottomMassPz, hBottomMassY, hBottomPt, hBottomPtPz, hBottomPz, hBottomY, hBottomYPz, hTopMass, hTopPt, hTopY, hTTbarMass, hTTbarPt, hTTbarY, hWminPz, hWplusPz, nEvt, MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().

{
  if(!dbe) return;
  dbe->setCurrentFolder("Generator/TTbar");

  nEvt = dbe->book1D("nEvt", "n analyzed Events", 1, 0., 1.);

  hTopPt         = dbe->book1D("TTbar_TopPt","t quark transverse momentum",1000,0.,1000.); hTopPt->setAxisTitle("t quark transverse momentum");
  hTopY          = dbe->book1D("TTbar_TopY","t quark rapidity",200,-5.,5.);                hTopY->setAxisTitle("t quark rapidity");
  hTopMass       = dbe->book1D("TTbar_TopMass","t quark mass",500,0.,500.);                hTopMass->setAxisTitle("t quark mass");
 
  hTTbarPt       = dbe->book1D("TTbar_TTbarPt","tt pair transverse momentum",1000,0.,1000.); hTTbarPt->setAxisTitle("tt pair transverse momentum");
  hTTbarY        = dbe->book1D("TTbar_TTbarY","tt pair rapidity",200,-5.,5.);                hTTbarY->setAxisTitle("tt pair rapidity");
  hTTbarMass     = dbe->book1D("TTbar_TTbarMass","tt pair mass",1000,0.,1000.);              hTTbarMass->setAxisTitle("tt pair mass");

  hBottomPt      = dbe->book1D("TTbar_BottomPt","b quark transverse momentum",1000,0.,1000.);     hBottomPt->setAxisTitle("b quark transverse momentum");
  hBottomEta     = dbe->book1D("TTbar_BottomEta","b quark pseudo-rapidity",200,-5.,5.);           hBottomEta->setAxisTitle("b quark pseudo-rapidity");
  hBottomY       = dbe->book1D("TTbar_BottomY","b quark rapidity",200,-5.,5.);                    hBottomY->setAxisTitle("b quark rapidity");
  hBottomPz      = dbe->book1D("TTbar_BottomPz","b quark longitudinal momentum",200,-100.,100.);  hBottomPz->setAxisTitle("b quark longitudinal momentum");
  hBottomE       = dbe->book1D("TTbar_BottomE","b quark energy",1000,0.,1000.);                   hBottomE->setAxisTitle("b quark energy");
  hBottomMass    = dbe->book1D("TTbar_BottomMass","b quark mass",50,0.,5.);                       hBottomMass->setAxisTitle("b quark mass");
  
  hBottomPtPz    = dbe->book2D("TTbar_BottomPtPz","b quark longitudinal vs transverse momentum",1000,0.,1000.,200,-100.,100.);     hBottomPtPz->setAxisTitle("P_{z} (GeV)",1);hBottomPtPz->setAxisTitle("P_{t} (GeV)",2);
  hBottomEtaPz   = dbe->book2D("TTbar_BottomEtaPz","b quark longitudinal momentum vs pseudorapidity",200,-5.,5.,200,-100.,100.);   hBottomEtaPz->setAxisTitle("#eta",1);hBottomEtaPz->setAxisTitle("P_{z} (GeV)",1);
  hBottomEtaPt   = dbe->book2D("TTbar_BottomEtaPt"," quark transveral   momentum vs pseudorapidity",200,-5.,5.,1000,0.,1000.);   hBottomEtaPt->setAxisTitle("#eta");hBottomEtaPt->setAxisTitle("P_{t} (GeV)");
  hBottomYPz     = dbe->book2D("TTbar_BottomYPz","b quark longitudinal momentum vs rapidity",200,-5.,5.,200,-100.,100.);           hBottomYPz->setAxisTitle("Y");hBottomYPz->setAxisTitle("P_{z} (GeV)");
  hBottomMassPz  = dbe->book2D("TTbar_BottomMassPz","b quark longitudinal momentum vs mass",50,0.,5.,200,-100.,100.);              hBottomMassPz->setAxisTitle("M (GeV)");hBottomMassPz->setAxisTitle("P_{z} (GeV)");
  hBottomMassEta = dbe->book2D("TTbar_BottomMassEta","b quark pseudorapidity vs mass",50,0.,5.,200,-5.,5.);                        hBottomMassEta->setAxisTitle("M (GeV)");hBottomMassEta->setAxisTitle("#eta");
  hBottomMassY   = dbe->book2D("TTbar_BottomMassY","b quark rapidity vs mass",50,0.,5.,200,-5.,5.);                                hBottomMassY->setAxisTitle("M (GeV)"); hBottomMassY->setAxisTitle("Y");
  hBottomMassDeltaY = dbe->book2D("TTbar_BottomMassDeltaY","b quark pseudorapidity - rapidity vs mass",50,0.,50.,2000,-5.,5.);     hBottomMassDeltaY->setAxisTitle("M (GeV)");hBottomMassDeltaY->setAxisTitle("Y");

  hWplusPz       = dbe->book1D("TTbar_WplusPz","W+ boson longitudinal momentum",200,-100.,100.);  hWplusPz->setAxisTitle("W+ boson longitudinal momentum");
  hWminPz        = dbe->book1D("TTbar_WminPz","W- boson longitudinal momentum",200,-100.,100.);   hWminPz->setAxisTitle("W- boson longitudinal momentum");

}

void TTbar_Kinematics::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)		`[private, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 209 of file TTbar_Kinematics.cc.

{
}

void TTbar_Kinematics::beginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)		`[private, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 197 of file TTbar_Kinematics.cc.

{
}

void TTbar_Kinematics::endJob ( void ) [private, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 191 of file TTbar_Kinematics.cc.

{
}

void TTbar_Kinematics::endLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)		`[private, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 215 of file TTbar_Kinematics.cc.

{
}

void TTbar_Kinematics::endRun	(	edm::Run const &	,
		edm::EventSetup const &
	)		`[private, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 203 of file TTbar_Kinematics.cc.

{
}

void TTbar_Kinematics::fillDescriptions ( edm::ConfigurationDescriptions & descriptions ) [static]

Reimplemented from edm::EDAnalyzer.

Definition at line 221 of file TTbar_Kinematics.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                             {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}