#include <DQMExample_Step1.h>

Inheritance diagram for DQMExample_Step1:

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

virtual	~DQMExample_Step1 ()

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
	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	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 (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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

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

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

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &eSetup)

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

Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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::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

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)

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 ( )

virtual

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
	)

protectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 85 of file DQMExample_Step1.cc.

References funct::abs(), deltaR(), Distance(), reco::Particle::eta(), reco::LeafCandidate::eta(), edm::HandleBase::failedToGet(), MonitorElement::Fill(), HcalObjRepresent::Fill(), edm::Event::getByToken(), 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(), j, relativeConstraints::keys, MediumEle(), NULL, 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(), theCaloJetCollection_, theElectronCollection_, thePfMETCollection_, thePVCollection_, triggerEvent_, triggerFilter_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), triggerPath_, triggerResults_, trigNames, funct::true, and findQualityFiles::v.

 {
   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 = NULL;
   const reco::GsfElectron* ele2 = NULL;
   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 = NULL;
   const reco::CaloJet* jet2 = NULL;
   
   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 = NULL;
   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.size() >= 1 ) 
     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::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumi,
		edm::EventSetup const &	eSetup
	)

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 75 of file DQMExample_Step1.cc.

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

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 323 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 374 of file DQMExample_Step1.cc.

References SiPixelRawToDigiRegional_cfi::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 365 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 369 of file DQMExample_Step1.cc.

References SiPixelRawToDigiRegional_cfi::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::endLuminosityBlock	(	edm::LuminosityBlock const &	lumi,
		edm::EventSetup const &	eSetup
	)

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 305 of file DQMExample_Step1.cc.

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

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

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 314 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 386 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()->trackerExpectedHitsInner().numberOfHits();
   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;
 }