DrellYanValidation Class Reference

#include <DrellYanValidation.h>

Inheritance diagram for DrellYanValidation:

Public Member Functions
virtual void	analyze (edm::Event const &, edm::EventSetup const &) override
virtual void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
virtual void	dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override
	DrellYanValidation (const edm::ParameterSet &)
virtual	~DrellYanValidation ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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...
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 edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 33 of file DrellYanValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 19 of file DrellYanValidation.cc.

References 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_);
}

DrellYanValidation::~DrellYanValidation ( )

virtual

Definition at line 28 of file DrellYanValidation.cc.

{}

Member Function Documentation

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

overridevirtual

Gathering the HepMCProduct information

Implements edm::stream::EDAnalyzerBase.

Definition at line 69 of file DrellYanValidation.cc.

References _flavor, funct::abs(), assert(), funct::cos(), cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, MonitorElement::Fill(), HepMCValidationHelper::findFSRPhotons(), fPDGTable, gamma_energy, edm::Event::getByToken(), hepmcCollectionToken_, i, leadeta, leadpt, RazorAnalyzer::lep1, RazorAnalyzer::lep2, nEvt, perp(), edm::es::products(), seceta, secpt, HepMCValidationHelper::sortByPt(), lumiQTWidget::t, WeightManager::weight(), puppiForMET_cff::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 = (abs(_flavor) == 11 || abs(_flavor) == 12 || abs(_flavor) ==13 || abs(_flavor) ==14 || abs(_flavor) ==16) ? 1 : 3;

  bool vetotau = true; //(abs(_flavor) == 11 || abs(_flavor) == 12 || abs(_flavor) ==13 || abs(_flavor) ==14 || 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 && abs((*iter)->pdg_id() == 15) ) return; 
    }
    if((*iter)->status()==requiredstatus) {
      if(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.size() != 0 && 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

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, DQMHelper::book1dHisto(), cos_theta_gamma_lepton, dilep_mass, dilep_massPeak, dilep_pt, dilep_ptLog, dilep_rap, gamma_energy, leadeta, leadpt, nEvt, seceta, secpt, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, Zdaughters, Zmass, ZmassPeak, Zpt, ZptLog, and Zrap.

                                                                                                {

    std::string folderName = "Generator/DrellYan";
    folderName+=_name;
    DQMHelper dqm(&i); i.setCurrentFolder(folderName.c_str());
    
    // 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;
}

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 30 of file DrellYanValidation.cc.

References fPDGTable, and edm::EventSetup::getData().

                                                                            {
  c.getData( fPDGTable );
}

Member Data Documentation

int DrellYanValidation::_flavor

private

decay flavor

Definition at line 57 of file DrellYanValidation.h.

Referenced by analyze().

std::string DrellYanValidation::_name

private

decay flavor name

Definition at line 59 of file DrellYanValidation.h.

Referenced by plotting.Subtract::__str__(), plotting.FakeDuplicate::__str__(), plotting.AggregateBins::__str__(), plotting.AggregateHistos::__str__(), plotting.ROC::__str__(), plotting.PlotGroup::_drawSeparate(), plotting.PlotGroup::_save(), bookHistograms(), plotting.Subtract::create(), plotting.FakeDuplicate::create(), plotting.AggregateBins::create(), plotting.AggregateHistos::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(), and plotting.PlotGroup::remove().

MonitorElement * DrellYanValidation::cos_theta_gamma_lepton

private

Definition at line 54 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* DrellYanValidation::dilep_mass

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::dilep_massPeak

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::dilep_pt

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::dilep_ptLog

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::dilep_rap

private

Definition at line 52 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

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

private

PDT table.

Definition at line 48 of file DrellYanValidation.h.

Referenced by analyze(), and dqmBeginRun().

MonitorElement* DrellYanValidation::gamma_energy

private

Definition at line 54 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag DrellYanValidation::hepmcCollection_

private

Definition at line 45 of file DrellYanValidation.h.

Referenced by DrellYanValidation().

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

private

Definition at line 61 of file DrellYanValidation.h.

Referenced by analyze(), and DrellYanValidation().

MonitorElement * DrellYanValidation::leadeta

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* DrellYanValidation::leadpt

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* DrellYanValidation::nEvt

private

Definition at line 50 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::seceta

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::secpt

private

Definition at line 53 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

WeightManager DrellYanValidation::wmanager_

private

Definition at line 44 of file DrellYanValidation.h.

Referenced by analyze().

MonitorElement * DrellYanValidation::Zdaughters

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* DrellYanValidation::Zmass

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::ZmassPeak

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::Zpt

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::ZptLog

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * DrellYanValidation::Zrap

private

Definition at line 51 of file DrellYanValidation.h.

Referenced by analyze(), and bookHistograms().