#include <DQMExample_Step1.h>

Inheritance diagram for DQMExample_Step1:

Public Member Functions
	DQMExample_Step1 (const edm::ParameterSet &ps)

	~DQMExample_Step1 () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override

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

void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

void	endRun (edm::Run const &run, edm::EventSetup const &eSetup) override

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)

Private Member Functions
void	bookHistos (DQMStore::IBooker &)

double	calcDeltaPhi (double phi1, double phi2)

double	Distance (const reco::Candidate &c1, const reco::Candidate &c2)

double	DistancePhi (const reco::Candidate &c1, const reco::Candidate &c2)

bool	MediumEle (const edm::Event &iEvent, const edm::EventSetup &iESetup, const reco::GsfElectron &electron)

Private Attributes
MonitorElement *	h_eEta_leading

MonitorElement *	h_eEta_leading_HLT

MonitorElement *	h_eEta_leading_HLT_matched

MonitorElement *	h_eEta_leading_matched

MonitorElement *	h_eMultiplicity

MonitorElement *	h_eMultiplicity_HLT

MonitorElement *	h_ePhi_leading

MonitorElement *	h_ePhi_leading_HLT

MonitorElement *	h_ePhi_leading_HLT_matched

MonitorElement *	h_ePhi_leading_matched

MonitorElement *	h_ePt_diff

MonitorElement *	h_ePt_leading

MonitorElement *	h_ePt_leading_HLT

MonitorElement *	h_ePt_leading_HLT_matched

MonitorElement *	h_ePt_leading_matched

MonitorElement *	h_jEta_leading

MonitorElement *	h_jMultiplicity

MonitorElement *	h_jPhi_leading

MonitorElement *	h_jPt_leading

MonitorElement *	h_pfMet

MonitorElement *	h_vertex_number

int	nBJets

int	nElectrons

double	ptThrJet_

double	ptThrL1_

double	ptThrL2_

double	ptThrMet_

math::XYZPoint	PVPoint_

edm::EDGetTokenT< reco::BeamSpot >	theBSCollection_

edm::EDGetTokenT< reco::CaloJetCollection >	theCaloJetCollection_

edm::EDGetTokenT< reco::ConversionCollection >	theConversionCollection_

edm::EDGetTokenT< reco::GsfElectronCollection >	theElectronCollection_

edm::EDGetTokenT< reco::PFMETCollection >	thePfMETCollection_

edm::EDGetTokenT< reco::VertexCollection >	thePVCollection_

edm::EDGetTokenT< trigger::TriggerEvent >	triggerEvent_

edm::InputTag	triggerFilter_

std::string	triggerPath_

edm::EDGetTokenT< edm::TriggerResults >	triggerResults_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 48 of file DQMExample_Step1.h.

Constructor & Destructor Documentation

DQMExample_Step1::DQMExample_Step1 ( const edm::ParameterSet & ps )

Definition at line 22 of file DQMExample_Step1.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ptThrJet_, ptThrL1_, ptThrL2_, ptThrMet_, AlCaHLTBitMon_QueryRunRegistry::string, theBSCollection_, theCaloJetCollection_, theConversionCollection_, theElectronCollection_, thePfMETCollection_, thePVCollection_, triggerEvent_, triggerFilter_, triggerPath_, and triggerResults_.

 {
   edm::LogInfo("DQMExample_Step1") <<  "Constructor  DQMExample_Step1::DQMExample_Step1 " << std::endl;
   
   // Get parameters from configuration file
   theElectronCollection_   = consumes<reco::GsfElectronCollection>(ps.getParameter<edm::InputTag>("electronCollection"));
   theCaloJetCollection_    = consumes<reco::CaloJetCollection>(ps.getParameter<edm::InputTag>("caloJetCollection"));
   thePfMETCollection_      = consumes<reco::PFMETCollection>(ps.getParameter<edm::InputTag>("pfMETCollection"));
   theConversionCollection_ = consumes<reco::ConversionCollection>(ps.getParameter<edm::InputTag>("conversionsCollection"));
   thePVCollection_         = consumes<reco::VertexCollection>(ps.getParameter<edm::InputTag>("PVCollection"));
   theBSCollection_         = consumes<reco::BeamSpot>(ps.getParameter<edm::InputTag>("beamSpotCollection"));
   triggerEvent_            = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("TriggerEvent"));
   triggerResults_          = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
   triggerFilter_           = ps.getParameter<edm::InputTag>("TriggerFilter");
   triggerPath_             = ps.getParameter<std::string>("TriggerPath");
 
 
   // cuts:
   ptThrL1_ = ps.getUntrackedParameter<double>("PtThrL1");
   ptThrL2_ = ps.getUntrackedParameter<double>("PtThrL2");
   ptThrJet_ = ps.getUntrackedParameter<double>("PtThrJet");
   ptThrMet_ = ps.getUntrackedParameter<double>("PtThrMet");
  
 }

DQMExample_Step1::~DQMExample_Step1 ( )

override

Definition at line 50 of file DQMExample_Step1.cc.

 {
   edm::LogInfo("DQMExample_Step1") <<  "Destructor DQMExample_Step1::~DQMExample_Step1 " << std::endl;
 }

Member Function Documentation

void DQMExample_Step1::analyze	(	edm::Event const &	e,
		edm::EventSetup const &	eSetup
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 75 of file DQMExample_Step1.cc.

References funct::abs(), edm::HLTGlobalStatus::accept(), deltaR(), Distance(), ElectronMcFakeValidator_cfi::electronCollection, reco::Particle::eta(), reco::LeafCandidate::eta(), edm::HandleBase::failedToGet(), MonitorElement::Fill(), HcalObjRepresent::Fill(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), h_eEta_leading, h_eEta_leading_HLT, h_eEta_leading_HLT_matched, h_eEta_leading_matched, h_eMultiplicity, h_eMultiplicity_HLT, h_ePhi_leading, h_ePhi_leading_HLT, h_ePhi_leading_HLT_matched, h_ePhi_leading_matched, h_ePt_diff, h_ePt_leading, h_ePt_leading_HLT, h_ePt_leading_HLT_matched, h_ePt_leading_matched, h_jEta_leading, h_jMultiplicity, h_jPhi_leading, h_jPt_leading, h_pfMet, h_vertex_number, edm::HandleBase::isValid(), relativeConstraints::keys, MediumEle(), trigger::TriggerObject::particle(), reco::Particle::pdgId(), reco::Particle::phi(), reco::LeafCandidate::phi(), reco::Particle::pt(), reco::LeafCandidate::pt(), ptThrJet_, ptThrL1_, ptThrL2_, PVPoint_, edm::TriggerNames::size(), trigger::TriggerEvent::sizeFilters(), theCaloJetCollection_, theElectronCollection_, thePfMETCollection_, thePVCollection_, triggerEvent_, triggerFilter_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), TriggerAnalyzer::triggerObjects, triggerPath_, triggerResults_, trigNames, funct::true, findQualityFiles::v, and edm::HLTGlobalStatus::wasrun().

 {
   edm::LogInfo("DQMExample_Step1") <<  "DQMExample_Step1::analyze" << std::endl;
 
 
   //-------------------------------
   //--- Vertex Info
   //-------------------------------
   edm::Handle<reco::VertexCollection> vertexHandle;
   e.getByToken(thePVCollection_, vertexHandle);
   if ( !vertexHandle.isValid() ) 
     {
       edm::LogError ("DQMClientExample") << "invalid collection: vertex" << "\n";
       return;
     }
   
   int vertex_number = vertexHandle->size();
   reco::VertexCollection::const_iterator v = vertexHandle->begin();
 
   math::XYZPoint PVPoint(-999, -999, -999);
   if(vertex_number != 0)
     PVPoint = math::XYZPoint(v->position().x(), v->position().y(), v->position().z());
   
   PVPoint_=PVPoint;
 
   //-------------------------------
   //--- MET
   //-------------------------------
   edm::Handle<reco::PFMETCollection> pfMETCollection;
   e.getByToken(thePfMETCollection_, pfMETCollection);
   if ( !pfMETCollection.isValid() )    
     {
       edm::LogError ("DQMClientExample") << "invalid collection: MET" << "\n";
       return;
     }
   //-------------------------------
   //--- Electrons
   //-------------------------------
   edm::Handle<reco::GsfElectronCollection> electronCollection;
   e.getByToken(theElectronCollection_, electronCollection);
   if ( !electronCollection.isValid() )
     {
       edm::LogError ("DQMClientExample") << "invalid collection: electrons" << "\n";
       return;
     }
 
   float nEle=0;
   int posEle=0, negEle=0;
   const reco::GsfElectron* ele1 = nullptr;
   const reco::GsfElectron* ele2 = nullptr;
   for (reco::GsfElectronCollection::const_iterator recoElectron=electronCollection->begin(); recoElectron!=electronCollection->end(); ++recoElectron)
     {
       //decreasing pT
       if( MediumEle(e,eSetup,*recoElectron) )
     {
       if(!ele1 && recoElectron->pt() > ptThrL1_)
         ele1 = &(*recoElectron);
       
       else if(!ele2 && recoElectron->pt() > ptThrL2_)
         ele2 = &(*recoElectron);
 
     }
       
       if(recoElectron->charge()==1)
     posEle++;
       else if(recoElectron->charge()==-1)
     negEle++;
 
     } // end of loop over electrons
   
   nEle = posEle+negEle;
   
   //-------------------------------
   //--- Jets
   //-------------------------------
   edm::Handle<reco::CaloJetCollection> caloJetCollection;
   e.getByToken (theCaloJetCollection_,caloJetCollection);
   if ( !caloJetCollection.isValid() ) 
     {
       edm::LogError ("DQMClientExample") << "invalid collection: jets" << "\n";
       return;
     }
 
   int   nJet = 0;
   const reco::CaloJet* jet1 = nullptr;
   const reco::CaloJet* jet2 = nullptr;
   
   for (reco::CaloJetCollection::const_iterator i_calojet = caloJetCollection->begin(); i_calojet != caloJetCollection->end(); ++i_calojet) 
     {
       //remove jet-ele matching
       if(ele1)
     if (Distance(*i_calojet,*ele1) < 0.3) continue;
       
       if(ele2)
     if (Distance(*i_calojet,*ele2) < 0.3) continue;
       
       if (i_calojet->pt() < ptThrJet_) continue;
 
       nJet++;
       
       if (!jet1) 
     jet1 = &(*i_calojet);
       
       else if (!jet2)
     jet2 = &(*i_calojet);
     }
   
   // ---------------------------
   // ---- Analyze Trigger Event
   // ---------------------------
 
   //check what is in the menu
   edm::Handle<edm::TriggerResults> hltresults;
   e.getByToken(triggerResults_,hltresults);
   
   if(!hltresults.isValid())
     {
       edm::LogError ("DQMClientExample") << "invalid collection: TriggerResults" << "\n";
       return;
     }
   
   bool hasFired = false;
   const edm::TriggerNames& trigNames = e.triggerNames(*hltresults);
   unsigned int numTriggers = trigNames.size();
   
   for( unsigned int hltIndex=0; hltIndex<numTriggers; ++hltIndex )
     {
       if (trigNames.triggerName(hltIndex)==triggerPath_ &&  hltresults->wasrun(hltIndex) &&  hltresults->accept(hltIndex))
     hasFired = true;
     }
   
 
 
   //access the trigger event
   edm::Handle<trigger::TriggerEvent> triggerEvent;
   e.getByToken(triggerEvent_, triggerEvent);
   if( triggerEvent.failedToGet() )
     {
       edm::LogError ("DQMClientExample") << "invalid collection: TriggerEvent" << "\n";
       return;
     }
 
 
   reco::Particle* ele1_HLT = nullptr;
   int nEle_HLT = 0;
 
   size_t filterIndex = triggerEvent->filterIndex( triggerFilter_ );
   trigger::TriggerObjectCollection triggerObjects = triggerEvent->getObjects();
   if( !(filterIndex >= triggerEvent->sizeFilters()) )
     {
       const trigger::Keys& keys = triggerEvent->filterKeys( filterIndex );
       std::vector<reco::Particle> triggeredEle;
       
       for( size_t j = 0; j < keys.size(); ++j ) 
     {
       trigger::TriggerObject foundObject = triggerObjects[keys[j]];
       if( abs( foundObject.particle().pdgId() ) != 11 )  continue; //make sure that it is an electron
       
       triggeredEle.push_back( foundObject.particle() );
       ++nEle_HLT;
     }
       
       if( !triggeredEle.empty() ) 
     ele1_HLT = &(triggeredEle.at(0));
     }
 
   //-------------------------------
   //--- Fill the histos
   //-------------------------------
 
   //vertex
   h_vertex_number -> Fill( vertex_number );
 
   //met
   h_pfMet -> Fill( pfMETCollection->begin()->et() );
 
   //multiplicities
   h_eMultiplicity->Fill(nEle);       
   h_jMultiplicity->Fill(nJet);
   h_eMultiplicity_HLT->Fill(nEle_HLT);
 
   //leading not matched
   if(ele1)
     {
       h_ePt_leading->Fill(ele1->pt());
       h_eEta_leading->Fill(ele1->eta());
       h_ePhi_leading->Fill(ele1->phi());
     }
   if(ele1_HLT)
     {
       h_ePt_leading_HLT->Fill(ele1_HLT->pt());
       h_eEta_leading_HLT->Fill(ele1_HLT->eta());
       h_ePhi_leading_HLT->Fill(ele1_HLT->phi());
     }
   //leading Jet
   if(jet1)
     {
       h_jPt_leading->Fill(jet1->pt());
       h_jEta_leading->Fill(jet1->eta());
       h_jPhi_leading->Fill(jet1->phi());
     }
 
 
   //fill only when the trigger candidate mathes with the reco one
   if( ele1 && ele1_HLT && deltaR(*ele1_HLT,*ele1) < 0.3 && hasFired==true )
     {
       h_ePt_leading_matched->Fill(ele1->pt());
       h_eEta_leading_matched->Fill(ele1->eta());
       h_ePhi_leading_matched->Fill(ele1->phi());
       
       h_ePt_leading_HLT_matched->Fill(ele1_HLT->pt());
       h_eEta_leading_HLT_matched->Fill(ele1_HLT->eta());
       h_ePhi_leading_HLT_matched->Fill(ele1_HLT->phi());
 
       h_ePt_diff->Fill(ele1->pt()-ele1_HLT->pt());
     }
 }

void DQMExample_Step1::bookHistograms	(	DQMStore::IBooker &	ibooker_,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 65 of file DQMExample_Step1.cc.

References bookHistos().

 {
   edm::LogInfo("DQMExample_Step1") <<  "DQMExample_Step1::bookHistograms" << std::endl;
   
   //book at beginRun
   bookHistos(ibooker_);
 }

void DQMExample_Step1::bookHistos ( DQMStore::IBooker & ibooker_ )

private

Definition at line 305 of file DQMExample_Step1.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), h_eEta_leading, h_eEta_leading_HLT, h_eEta_leading_HLT_matched, h_eEta_leading_matched, h_eMultiplicity, h_eMultiplicity_HLT, h_ePhi_leading, h_ePhi_leading_HLT, h_ePhi_leading_HLT_matched, h_ePhi_leading_matched, h_ePt_diff, h_ePt_leading, h_ePt_leading_HLT, h_ePt_leading_HLT_matched, h_ePt_leading_matched, h_jEta_leading, h_jMultiplicity, h_jPhi_leading, h_jPt_leading, h_pfMet, h_vertex_number, and DQMStore::IBooker::setCurrentFolder().

Referenced by bookHistograms().

 {
   ibooker_.cd();
   ibooker_.setCurrentFolder("Physics/TopTest");
 
   h_vertex_number = ibooker_.book1D("Vertex_number", "Number of event vertices in collection", 40,-0.5,   39.5 );
 
   h_pfMet        = ibooker_.book1D("pfMet",        "Pf Missing E_{T}; GeV"          , 20,  0.0 , 100);
 
   h_eMultiplicity = ibooker_.book1D("NElectrons","# of electrons per event",10,0.,10.);
   h_ePt_leading_matched = ibooker_.book1D("ElePt_leading_matched","Pt of leading electron",50,0.,100.);
   h_eEta_leading_matched = ibooker_.book1D("EleEta_leading_matched","Eta of leading electron",50,-5.,5.);
   h_ePhi_leading_matched = ibooker_.book1D("ElePhi_leading_matched","Phi of leading electron",50,-3.5,3.5);
 
   h_ePt_leading = ibooker_.book1D("ElePt_leading","Pt of leading electron",50,0.,100.);
   h_eEta_leading = ibooker_.book1D("EleEta_leading","Eta of leading electron",50,-5.,5.);
   h_ePhi_leading = ibooker_.book1D("ElePhi_leading","Phi of leading electron",50,-3.5,3.5);
 
   h_jMultiplicity = ibooker_.book1D("NJets","# of electrons per event",10,0.,10.);
   h_jPt_leading = ibooker_.book1D("JetPt_leading","Pt of leading Jet",150,0.,300.);
   h_jEta_leading = ibooker_.book1D("JetEta_leading","Eta of leading Jet",50,-5.,5.);
   h_jPhi_leading = ibooker_.book1D("JetPhi_leading","Phi of leading Jet",50,-3.5,3.5);
 
   h_eMultiplicity_HLT = ibooker_.book1D("NElectrons_HLT","# of electrons per event @HLT",10,0.,10.);
   h_ePt_leading_HLT = ibooker_.book1D("ElePt_leading_HLT","Pt of leading electron @HLT",50,0.,100.);
   h_eEta_leading_HLT = ibooker_.book1D("EleEta_leading_HLT","Eta of leading electron @HLT",50,-5.,5.);
   h_ePhi_leading_HLT = ibooker_.book1D("ElePhi_leading_HLT","Phi of leading electron @HLT",50,-3.5,3.5);
 
   h_ePt_leading_HLT_matched = ibooker_.book1D("ElePt_leading_HLT_matched","Pt of leading electron @HLT",50,0.,100.);
   h_eEta_leading_HLT_matched = ibooker_.book1D("EleEta_leading_HLT_matched","Eta of leading electron @HLT",50,-5.,5.);
   h_ePhi_leading_HLT_matched = ibooker_.book1D("ElePhi_leading_HLT_matched","Phi of leading electron @HLT",50,-3.5,3.5);
 
   h_ePt_diff = ibooker_.book1D("ElePt_diff_matched","pT(RECO) - pT(HLT) for mathed candidates",100,-10,10.);
 
   ibooker_.cd();  
 
 }

double DQMExample_Step1::calcDeltaPhi	(	double	phi1,
		double	phi2
	)

private

Definition at line 356 of file DQMExample_Step1.cc.

References hiPixelPairStep_cff::deltaPhi.

                                                               {
   double deltaPhi = phi1 - phi2;
   if (deltaPhi < 0) deltaPhi = -deltaPhi;
   if (deltaPhi > 3.1415926) {
     deltaPhi = 2 * 3.1415926 - deltaPhi;
   }
   return deltaPhi;
 }

double DQMExample_Step1::Distance	(	const reco::Candidate &	c1,
		const reco::Candidate &	c2
	)

private

Definition at line 347 of file DQMExample_Step1.cc.

References deltaR().

Referenced by analyze().

                                                                                       {
         return  deltaR(c1,c2);
 }

double DQMExample_Step1::DistancePhi	(	const reco::Candidate &	c1,
		const reco::Candidate &	c2
	)

private

Definition at line 351 of file DQMExample_Step1.cc.

References hiPixelPairStep_cff::deltaPhi, and reco::Candidate::p4().

                                                                                          {
         return  deltaPhi(c1.p4().phi(),c2.p4().phi());
 }

void DQMExample_Step1::dqmBeginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 58 of file DQMExample_Step1.cc.

 {
   edm::LogInfo("DQMExample_Step1") <<  "DQMExample_Step1::beginRun" << std::endl;
 }

void DQMExample_Step1::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

overrideprotected

Definition at line 296 of file DQMExample_Step1.cc.

 {
   edm::LogInfo("DQMExample_Step1") <<  "DQMExample_Step1::endRun" << std::endl;
 }

bool DQMExample_Step1::MediumEle	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iESetup,
		const reco::GsfElectron &	electron
	)

private

Definition at line 368 of file DQMExample_Step1.cc.

Referenced by analyze().

 {
     
   //********* CONVERSION TOOLS
   edm::Handle<reco::ConversionCollection> conversions_h;
   iEvent.getByToken(theConversionCollection_, conversions_h);
   
   bool isMediumEle = false; 
   
   float pt = electron.pt();
   float eta = electron.eta();
     
   int isEB            = electron.isEB();
   float sigmaIetaIeta = electron.sigmaIetaIeta();
   float DetaIn        = electron.deltaEtaSuperClusterTrackAtVtx();
   float DphiIn        = electron.deltaPhiSuperClusterTrackAtVtx();
   float HOverE        = electron.hadronicOverEm();
   float ooemoop       = (1.0/electron.ecalEnergy() - electron.eSuperClusterOverP()/electron.ecalEnergy());
   
   int mishits             = electron.gsfTrack()->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
   int nAmbiguousGsfTracks = electron.ambiguousGsfTracksSize();
   
   reco::GsfTrackRef eleTrack  = electron.gsfTrack() ;
   float dxy           = eleTrack->dxy(PVPoint_);  
   float dz            = eleTrack->dz (PVPoint_);
   
   edm::Handle<reco::BeamSpot> BSHandle;
   iEvent.getByToken(theBSCollection_, BSHandle);
   const reco::BeamSpot BS = *BSHandle;
   
   bool isConverted = ConversionTools::hasMatchedConversion(electron, conversions_h, BS.position());
   
   // default
   if(  (pt > 12.) && (fabs(eta) < 2.5) &&
        ( ( (isEB == 1) && (fabs(DetaIn)  < 0.004) ) || ( (isEB == 0) && (fabs(DetaIn)  < 0.007) ) ) &&
        ( ( (isEB == 1) && (fabs(DphiIn)  < 0.060) ) || ( (isEB == 0) && (fabs(DphiIn)  < 0.030) ) ) &&
        ( ( (isEB == 1) && (sigmaIetaIeta < 0.010) ) || ( (isEB == 0) && (sigmaIetaIeta < 0.030) ) ) &&
        ( ( (isEB == 1) && (HOverE        < 0.120) ) || ( (isEB == 0) && (HOverE        < 0.100) ) ) &&
        ( ( (isEB == 1) && (fabs(ooemoop) < 0.050) ) || ( (isEB == 0) && (fabs(ooemoop) < 0.050) ) ) &&
        ( ( (isEB == 1) && (fabs(dxy)     < 0.020) ) || ( (isEB == 0) && (fabs(dxy)     < 0.020) ) ) &&
        ( ( (isEB == 1) && (fabs(dz)      < 0.100) ) || ( (isEB == 0) && (fabs(dz)      < 0.100) ) ) &&
        ( ( (isEB == 1) && (!isConverted) ) || ( (isEB == 0) && (!isConverted) ) ) &&
        ( mishits == 0 ) &&
        ( nAmbiguousGsfTracks == 0 )      
        )
     isMediumEle=true;
   
   return isMediumEle;
 }