DrellYanValidation Class Reference

#include <DrellYanValidation.h>

Inheritance diagram for DrellYanValidation:

Public Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override
	DrellYanValidation (const edm::ParameterSet &)
	~DrellYanValidation () 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
	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
int	_flavor
	decay flavor More...
std::string	_name
	decay flavor name More...
MonitorElement *	cos_theta_gamma_lepton
MonitorElement *	dilep_mass
MonitorElement *	dilep_massPeak
MonitorElement *	dilep_pt
MonitorElement *	dilep_ptLog
MonitorElement *	dilep_rap
edm::ESHandle< HepPDT::ParticleDataTable >	fPDGTable
	PDT table. More...
edm::ESGetToken< HepPDT::ParticleDataTable, edm::DefaultRecord >	fPDGTableToken
MonitorElement *	gamma_energy
edm::InputTag	hepmcCollection_
edm::EDGetTokenT< edm::HepMCProduct >	hepmcCollectionToken_
MonitorElement *	leadeta
MonitorElement *	leadpt
MonitorElement *	nEvt
MonitorElement *	seceta
MonitorElement *	secpt
WeightManager	wmanager_
MonitorElement *	Zdaughters
MonitorElement *	Zmass
MonitorElement *	ZmassPeak
MonitorElement *	Zpt
MonitorElement *	ZptLog
MonitorElement *	Zrap

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

Definition at line 32 of file DrellYanValidation.h.

Constructor & Destructor Documentation

DrellYanValidation()

DrellYanValidation::DrellYanValidation ( const edm::ParameterSet &  iPSet )

explicit

Definition at line 19 of file DrellYanValidation.cc.

References fPDGTableToken, hepmcCollection_, and hepmcCollectionToken_.

    : wmanager_(iPSet, consumesCollector()),
      hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")),
      _flavor(iPSet.getParameter<int>("decaysTo")),
      _name(iPSet.getParameter<std::string>("name")) {
  hepmcCollectionToken_ = consumes<HepMCProduct>(hepmcCollection_);
  fPDGTableToken = esConsumes<edm::Transition::BeginRun>();
}

~DrellYanValidation()

DrellYanValidation::~DrellYanValidation ( )

override

Definition at line 28 of file DrellYanValidation.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Gathering the HepMCProduct information

Reimplemented from DQMEDAnalyzer.

Definition at line 84 of file DrellYanValidation.cc.

References _flavor, funct::abs(), cms::cuda::assert(), funct::cos(), cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, dqm::impl::MonitorElement::Fill(), HepMCValidationHelper::findFSRPhotons(), fPDGTable, gamma_energy, edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, mps_fire::i, iEvent, leadeta, leadpt, RazorAnalyzer::lep1, RazorAnalyzer::lep2, nEvt, LHEGenericFilter_cfi::ParticleID, HiggsValidation_cfi::pdg_id, perp(), edm::es::products(), seceta, secpt, jetUpdater_cfi::sort, HepMCValidationHelper::sortByPt(), submitPVValidationJobs::t, WeightManager::weight(), mps_merge::weight, wmanager_, x, y, z, Zdaughters, Zmass, ZmassPeak, Zpt, ZptLog, and Zrap.

                                                                                    {
  // we *DO NOT* rely on a Z entry in the particle listings!

  edm::Handle<HepMCProduct> evt;
  iEvent.getByToken(hepmcCollectionToken_, evt);

  //Get EVENT
  const HepMC::GenEvent* myGenEvent = evt->GetEvent();

  double weight = wmanager_.weight(iEvent);

  //std::cout << "weight: " << weight << std::endl;

  nEvt->Fill(0.5, weight);

  std::vector<const HepMC::GenParticle*> allproducts;

  //requires status 1 for leptons and neutrinos (except tau)
  int requiredstatus = (std::abs(_flavor) == 11 || std::abs(_flavor) == 12 || std::abs(_flavor) == 13 ||
                        std::abs(_flavor) == 14 || std::abs(_flavor) == 16)
                           ? 1
                           : 3;

  bool vetotau =
      true;  //(std::abs(_flavor) == 11 || std::abs(_flavor) == 12 || std::abs(_flavor) ==13 || std::abs(_flavor) ==14 || std::abs(_flavor) ==16) ? true : false;

  for (HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin();
       iter != myGenEvent->particles_end();
       ++iter) {
    if (vetotau) {
      if ((*iter)->status() == 3 && std::abs((*iter)->pdg_id()) == 15)
        return;
    }
    if ((*iter)->status() == requiredstatus) {
      if (std::abs((*iter)->pdg_id()) == _flavor)
        allproducts.push_back(*iter);
    }
  }

  //nothing to do if we don't have 2 particles
  if (allproducts.size() < 2)
    return;

  //sort them in pt
  std::sort(allproducts.begin(), allproducts.end(), HepMCValidationHelper::sortByPt);

  //get the first element and the first following element with opposite charge
  std::vector<const HepMC::GenParticle*> products;
  products.push_back(allproducts.front());
  const HepPDT::ParticleData* PData1 = fPDGTable->particle(HepPDT::ParticleID(allproducts.front()->pdg_id()));
  double charge1 = PData1->charge();
  for (unsigned int i = 1; i < allproducts.size(); ++i) {
    const HepPDT::ParticleData* PData2 = fPDGTable->particle(HepPDT::ParticleID(allproducts[i]->pdg_id()));
    double charge2 = PData2->charge();
    if (charge1 * charge2 < 0)
      products.push_back(allproducts[i]);
  }

  //if we did not find any opposite charge pair there is nothing to do
  if (products.size() < 2)
    return;

  assert(products[0]->momentum().perp() >= products[1]->momentum().perp());

  //leading lepton with pt > 20.
  if (products[0]->momentum().perp() < 20.)
    return;

  //assemble FourMomenta
  TLorentzVector lep1(products[0]->momentum().x(),
                      products[0]->momentum().y(),
                      products[0]->momentum().z(),
                      products[0]->momentum().t());
  TLorentzVector lep2(products[1]->momentum().x(),
                      products[1]->momentum().y(),
                      products[1]->momentum().z(),
                      products[1]->momentum().t());
  TLorentzVector dilepton_mom = lep1 + lep2;
  TLorentzVector dilepton_andphoton_mom = dilepton_mom;

  //mass > 60.
  if (dilepton_mom.M() < 60.)
    return;

  //find possible qed fsr photons
  std::vector<const HepMC::GenParticle*> fsrphotons;
  HepMCValidationHelper::findFSRPhotons(products, myGenEvent, 0.1, fsrphotons);

  Zdaughters->Fill(products[0]->pdg_id(), weight);
  Zdaughters->Fill(products[1]->pdg_id(), weight);

  std::vector<TLorentzVector> gammasMomenta;
  for (unsigned int ipho = 0; ipho < fsrphotons.size(); ++ipho) {
    TLorentzVector phomom(fsrphotons[ipho]->momentum().x(),
                          fsrphotons[ipho]->momentum().y(),
                          fsrphotons[ipho]->momentum().z(),
                          fsrphotons[ipho]->momentum().t());
    dilepton_andphoton_mom += phomom;
    Zdaughters->Fill(fsrphotons[ipho]->pdg_id(), weight);
    gammasMomenta.push_back(phomom);
  }
  //Fill Z histograms
  Zmass->Fill(dilepton_andphoton_mom.M(), weight);
  ZmassPeak->Fill(dilepton_andphoton_mom.M(), weight);
  Zpt->Fill(dilepton_andphoton_mom.Pt(), weight);
  ZptLog->Fill(log10(dilepton_andphoton_mom.Pt()), weight);
  Zrap->Fill(dilepton_andphoton_mom.Rapidity(), weight);

  //Fill dilepton histograms
  dilep_mass->Fill(dilepton_mom.M(), weight);
  dilep_massPeak->Fill(dilepton_mom.M(), weight);
  dilep_pt->Fill(dilepton_mom.Pt(), weight);
  dilep_ptLog->Fill(log10(dilepton_mom.Pt()), weight);
  dilep_rap->Fill(dilepton_mom.Rapidity(), weight);

  //Fill lepton histograms
  leadpt->Fill(lep1.Pt(), weight);
  secpt->Fill(lep2.Pt(), weight);
  leadeta->Fill(lep1.Eta(), weight);
  seceta->Fill(lep2.Eta(), weight);

  //boost everything in the Z frame
  TVector3 boost = dilepton_andphoton_mom.BoostVector();
  boost *= -1.;
  lep1.Boost(boost);
  lep2.Boost(boost);
  for (unsigned int ipho = 0; ipho < gammasMomenta.size(); ++ipho) {
    gammasMomenta[ipho].Boost(boost);
  }
  std::sort(gammasMomenta.begin(), gammasMomenta.end(), HepMCValidationHelper::GreaterByE<TLorentzVector>);

  //fill gamma histograms
  if (!gammasMomenta.empty() && dilepton_andphoton_mom.M() > 50.) {
    gamma_energy->Fill(gammasMomenta.front().E(), weight);
    double dphi = lep1.DeltaR(gammasMomenta.front()) < lep2.DeltaR(gammasMomenta.front())
                      ? lep1.DeltaPhi(gammasMomenta.front())
                      : lep2.DeltaPhi(gammasMomenta.front());
    cos_theta_gamma_lepton->Fill(cos(dphi), weight);
  }

}  //analyze

bookHistograms()

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

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 34 of file DrellYanValidation.cc.

References _name, cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, ALCARECODTCalibSynchCosmicsDQM_cff::folderName, gamma_energy, mps_fire::i, leadeta, leadpt, nEvt, seceta, secpt, AlCaHLTBitMon_QueryRunRegistry::string, Zdaughters, Zmass, ZmassPeak, Zpt, ZptLog, and Zrap.

                                                                                               {
  std::string folderName = "Generator/DrellYan";
  folderName += _name;
  DQMHelper dqm(&i);
  i.setCurrentFolder(folderName);

  // Number of analyzed events
  nEvt = dqm.book1dHisto("nEvt", "n analyzed Events", 1, 0., 1., "bin", "Number of Events");

  //Kinematics
  Zmass = dqm.book1dHisto("Zmass", "inv. Mass Z", 70, 0, 140, "M_{Z} (GeV)", "Number of Events");
  ZmassPeak = dqm.book1dHisto("ZmassPeak", "inv. Mass Z", 80, 80, 100, "M_{Z} (GeV)", "Number of Events");
  Zpt = dqm.book1dHisto("Zpt", "Z pt", 100, 0, 200, "P_{t}^{Z} (GeV)", "Number of Events");
  ZptLog =
      dqm.book1dHisto("ZptLog", "log(Z pt)", 100, 0., 5., "log_{10}(P_{t}^{Z}) (log_{10}(GeV))", "Number of Events");
  Zrap = dqm.book1dHisto("Zrap", "Z y", 100, -5, 5, "Y_{Z}", "Number of Events");
  Zdaughters = dqm.book1dHisto("Zdaughters", "Z daughters", 60, -30, 30, "Z daughters (PDG ID)", "Number of Events");

  dilep_mass = dqm.book1dHisto("dilep_mass", "inv. Mass dilepton", 70, 0, 140, "M_{ll} (GeV)", "Number of Events");
  dilep_massPeak =
      dqm.book1dHisto("dilep_massPeak", "inv. Mass dilepton", 80, 80, 100, "M_{ll} (GeV)", "Number of Events");
  dilep_pt = dqm.book1dHisto("dilep_pt", "dilepton pt", 100, 0, 200, "P_{t}^{ll} (GeV)", "Number of Events");
  dilep_ptLog = dqm.book1dHisto(
      "dilep_ptLog", "log(dilepton pt)", 100, 0., 5., "log_{10}(P_{t}^{ll}) (log_{10}(GeV))", "Number of Events");
  dilep_rap = dqm.book1dHisto("dilep_rap", "dilepton y", 100, -5, 5, "Y_{ll}", "Number of Events");

  gamma_energy = dqm.book1dHisto("gamma_energy",
                                 "photon energy in Z rest frame",
                                 200,
                                 0.,
                                 100.,
                                 "E_{#gamma}^{Z rest-frame} (GeV)",
                                 "Number of Events");
  cos_theta_gamma_lepton = dqm.book1dHisto("cos_theta_gamma_lepton",
                                           "cos_theta_gamma_lepton in Z rest frame",
                                           200,
                                           -1,
                                           1,
                                           "cos(#theta_{#gamma-lepton}^{Z rest-frame})",
                                           "Number of Events");

  leadpt = dqm.book1dHisto("leadpt", "leading lepton pt", 200, 0., 200., "P_{t}^{1st-lepton}", "Number of Events");
  secpt = dqm.book1dHisto("secpt", "second lepton pt", 200, 0., 200., "P_{t}^{2nd-lepton}", "Number of Events");
  leadeta = dqm.book1dHisto("leadeta", "leading lepton eta", 100, -5., 5., "#eta^{1st-lepton}", "Number of Events");
  seceta = dqm.book1dHisto("seceta", "second lepton eta", 100, -5., 5., "#eta^{2nd-lepton}", "Number of Events");

  return;
}

dqmBeginRun()

void DrellYanValidation::dqmBeginRun ( const edm::Run &  r, const edm::EventSetup &  c  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 30 of file DrellYanValidation.cc.

References c, fPDGTable, and fPDGTableToken.

                                                                            {
  fPDGTable = c.getHandle(fPDGTableToken);
}

Member Data Documentation

_flavor

int DrellYanValidation::_flavor

private

decay flavor

Definition at line 56 of file DrellYanValidation.h.

Referenced by analyze().

_name

std::string DrellYanValidation::_name

private

decay flavor name

Definition at line 58 of file DrellYanValidation.h.

Referenced by plotting.Subtract::__str__(), plotting.Transform::__str__(), plotting.FakeDuplicate::__str__(), plotting.CutEfficiency::__str__(), plotting.AggregateBins::__str__(), plotting.AggregateHistos::__str__(), plotting.ROC::__str__(), trackingPlots.TimePerEventPlot::__str__(), trackingPlots.TimePerTrackPlot::__str__(), trackingPlots.TimePerEventPlot::_create(), plotting.PlotGroup::_drawSeparate(), plotting.PlotGroup::_save(), bookHistograms(), plotting.Subtract::create(), plotting.Transform::create(), plotting.FakeDuplicate::create(), plotting.CutEfficiency::create(), plotting.AggregateBins::create(), plotting.AggregateHistos::create(), trackingPlots.TimePerTrackPlot::create(), plotting.Plot::create(), plotting.PlotGroup::draw(), plotting.Plot::getName(), plotting.PlotGroup::getName(), plotting.PlotterFolder::getName(), plotting.PlotterItem::getName(), plotting.PlotterFolder::getSelectionNameIterator(), validation.SimpleSample::name(), trackingPlots.Iteration::name(), plotting.PlotterItem::readDirs(), plotting.PlotGroup::remove(), and plotting.Plot::setName().

cos_theta_gamma_lepton

MonitorElement * DrellYanValidation::cos_theta_gamma_lepton

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

dilep_mass

MonitorElement* DrellYanValidation::dilep_mass

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

dilep_massPeak

MonitorElement * DrellYanValidation::dilep_massPeak

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

dilep_pt

MonitorElement * DrellYanValidation::dilep_pt

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

dilep_ptLog

MonitorElement * DrellYanValidation::dilep_ptLog

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

dilep_rap

MonitorElement * DrellYanValidation::dilep_rap

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

fPDGTable

edm::ESHandle<HepPDT::ParticleDataTable> DrellYanValidation::fPDGTable

private

PDT table.

Definition at line 46 of file DrellYanValidation.h.

Referenced by analyze(), and dqmBeginRun().

fPDGTableToken

edm::ESGetToken<HepPDT::ParticleDataTable, edm::DefaultRecord> DrellYanValidation::fPDGTableToken

private

Definition at line 47 of file DrellYanValidation.h.

Referenced by dqmBeginRun(), and DrellYanValidation().

gamma_energy

MonitorElement* DrellYanValidation::gamma_energy

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

hepmcCollection_

edm::InputTag DrellYanValidation::hepmcCollection_

private

Definition at line 43 of file DrellYanValidation.h.

Referenced by DrellYanValidation().

hepmcCollectionToken_

edm::EDGetTokenT<edm::HepMCProduct> DrellYanValidation::hepmcCollectionToken_

private

Definition at line 60 of file DrellYanValidation.h.

Referenced by analyze(), and DrellYanValidation().

leadeta

MonitorElement * DrellYanValidation::leadeta

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

leadpt

MonitorElement* DrellYanValidation::leadpt

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

nEvt

MonitorElement* DrellYanValidation::nEvt

private

Definition at line 49 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

seceta

MonitorElement * DrellYanValidation::seceta

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

secpt

MonitorElement * DrellYanValidation::secpt

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

wmanager_

WeightManager DrellYanValidation::wmanager_

private

Definition at line 42 of file DrellYanValidation.h.

Referenced by analyze().

Zdaughters

MonitorElement * DrellYanValidation::Zdaughters

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

Zmass

MonitorElement* DrellYanValidation::Zmass

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

ZmassPeak

MonitorElement * DrellYanValidation::ZmassPeak

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

Zpt

MonitorElement * DrellYanValidation::Zpt

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

ZptLog

MonitorElement * DrellYanValidation::ZptLog

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

Zrap

MonitorElement * DrellYanValidation::Zrap

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().