Inheritance diagram for EcalZmassTask:

Public Member Functions
	EcalZmassTask (const edm::ParameterSet &)

	~EcalZmassTask () override

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

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

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

	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)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

Private Attributes
const edm::EDGetTokenT< reco::GsfElectronCollection >	electronCollectionToken_

MonitorElement *	h_95_ee_invMass_BB

MonitorElement *	h_95_ee_invMass_EB

MonitorElement *	h_95_ee_invMass_EE

MonitorElement *	h_e1_et

MonitorElement *	h_e1_eta

MonitorElement *	h_e1_phi

MonitorElement *	h_e2_et

MonitorElement *	h_e2_eta

MonitorElement *	h_e2_phi

MonitorElement *	h_ee_invMass

MonitorElement *	h_ee_invMass_BB

MonitorElement *	h_ee_invMass_EB

MonitorElement *	h_ee_invMass_EE

const std::string	prefixME_

const edm::EDGetTokenT< reco::GsfTrackCollection >	trackCollectionToken_

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

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 45 of file EcalZmassTask.cc.

Constructor & Destructor Documentation

EcalZmassTask::EcalZmassTask ( const edm::ParameterSet & parameters )

explicit

Definition at line 76 of file EcalZmassTask.cc.

                                                                 :
   electronCollectionToken_(consumes<reco::GsfElectronCollection>(parameters.getParameter < edm::InputTag > ("electronCollection"))),
   trackCollectionToken_(consumes<reco::GsfTrackCollection>(parameters.getParameter<edm::InputTag>("trackCollection"))),
   prefixME_(parameters.getUntrackedParameter < std::string > ("prefixME", ""))
 {
 }

EcalZmassTask::~EcalZmassTask ( )

inlineoverride

Definition at line 49 of file EcalZmassTask.cc.

References analyze(), and bookHistograms().

49 {}

Member Function Documentation

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

overrideprivate

Definition at line 85 of file EcalZmassTask.cc.

References ElectronMcFakeValidator_cfi::electronCollection, electronCollectionToken_, MonitorElement::Fill(), edm::Event::getByToken(), h_ee_invMass_BB, h_ee_invMass_EB, h_ee_invMass_EE, HE, edm::HandleBase::isValid(), reco::HitPattern::MISSING_INNER_HITS, and trackCollectionToken_.

Referenced by ~EcalZmassTask().

 {
   using namespace edm;
   Handle < reco::GsfElectronCollection > electronCollection;
   iEvent.getByToken (electronCollectionToken_, electronCollection);
   if (!electronCollection.isValid ()) return;
 
   //get GSF Tracks
   Handle < reco::GsfTrackCollection > gsftracks_h;
   iEvent.getByToken (trackCollectionToken_, gsftracks_h);
 
   bool isIsolatedBarrel;
   bool isIDBarrel;
   bool isConvertedBarrel;
   bool isIsolatedEndcap;
   bool isIDEndcap;
   bool isConvertedEndcap;
 
   int elIsAccepted=0;
   int elIsAcceptedEB=0;
   int elIsAcceptedEE=0;
 
   std::vector<TLorentzVector> LV;
 
   for (reco::GsfElectronCollection::const_iterator recoElectron =
      electronCollection->begin ();
        recoElectron != electronCollection->end (); recoElectron++)
     {
 
       if (recoElectron->et () <= 25)
     continue;
 
       // Define Isolation variables
       double IsoTrk = (recoElectron->dr03TkSumPt () / recoElectron->et ());
       double IsoEcal =
     (recoElectron->dr03EcalRecHitSumEt () / recoElectron->et ());
       double IsoHcal =
     (recoElectron->dr03HcalTowerSumEt () / recoElectron->et ());
       double HE = (recoElectron->hcalOverEcal ());
 
       //Define ID variables
 
       float DeltaPhiTkClu = recoElectron->deltaPhiSuperClusterTrackAtVtx ();
       float DeltaEtaTkClu = recoElectron->deltaEtaSuperClusterTrackAtVtx ();
       float sigmaIeIe = recoElectron->sigmaIetaIeta ();
 
       //Define Conversion Rejection Variables
 
       float Dcot = recoElectron->convDcot ();
       float Dist = recoElectron->convDist ();
       int NumberOfExpectedInnerHits = recoElectron->gsfTrack()->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
 
       //quality flags
 
       isIsolatedBarrel = false;
       isIDBarrel = false;
       isConvertedBarrel = false;
       isIsolatedEndcap = false;
       isIDEndcap = false;
       isConvertedEndcap = false;
 
 
       /***** Barrel WP80 Cuts *****/
 
       if (fabs (recoElectron->eta ()) <= 1.4442)
     {
 
       /* Isolation */
       if (IsoTrk < 0.09 && IsoEcal < 0.07 && IsoHcal < 0.10)
         {
           isIsolatedBarrel = true;
         }
 
       /* Identification */
       if (fabs (DeltaEtaTkClu) < 0.004 && fabs (DeltaPhiTkClu) < 0.06
           && sigmaIeIe < 0.01 && HE < 0.04)
         {
           isIDBarrel = true;
         }
 
       /* Conversion Rejection */
       if ((fabs (Dist) >= 0.02 || fabs (Dcot) >= 0.02)
           && NumberOfExpectedInnerHits <= 1.0)
         {
           isConvertedBarrel = true;
         }
     }
 
       if (isIsolatedBarrel && isIDBarrel && isConvertedBarrel) {
     elIsAccepted++;
     elIsAcceptedEB++;
     TLorentzVector b_e2(recoElectron->momentum ().x (),recoElectron->momentum ().y (),recoElectron->momentum ().z (), recoElectron->p ());
     LV.push_back(b_e2);
       }
 
       /***** Endcap WP80 Cuts *****/
 
       if (fabs (recoElectron->eta ()) >= 1.5660
       && fabs (recoElectron->eta ()) <= 2.5000)
     {
 
       /* Isolation */
       if (IsoTrk < 0.04 && IsoEcal < 0.05 && IsoHcal < 0.025)
         {
           isIsolatedEndcap = true;
         }
 
       /* Identification */
       if (fabs (DeltaEtaTkClu) < 0.007 && fabs (DeltaPhiTkClu) < 0.03
           && sigmaIeIe < 0.031 && HE < 0.15)
         {
           isIDEndcap = true;
         }
 
       /* Conversion Rejection */
       if ((fabs (Dcot) > 0.02 || fabs (Dist) > 0.02)
           && NumberOfExpectedInnerHits <= 1.0)
         {
           isConvertedEndcap = true;
         }
     }
       if (isIsolatedEndcap && isIDEndcap && isConvertedEndcap) {
     elIsAccepted++;
     elIsAcceptedEE++;
     TLorentzVector e_e2(recoElectron->momentum ().x (),recoElectron->momentum ().y (),recoElectron->momentum ().z (), recoElectron->p ());
     LV.push_back(e_e2);
       }
 
     }
 
   // Calculate the Z invariant masses
 
   if (elIsAccepted>1){
     double e_ee_invMass=0; 
     if (LV.size()==2){
       TLorentzVector e_pair = LV[0] + LV[1];
       e_ee_invMass = e_pair.M ();
     }  
               
     if (elIsAcceptedEB==2){
       h_ee_invMass_BB->Fill(e_ee_invMass);
     }
     if (elIsAcceptedEE==2){
       h_ee_invMass_EE->Fill(e_ee_invMass);
     }
     if (elIsAcceptedEB==1 && elIsAcceptedEE==1){
       h_ee_invMass_EB->Fill(e_ee_invMass);
     }
               
     LV.clear();
                   
   }
 }

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 241 of file EcalZmassTask.cc.

References DQMStore::IBooker::book1D(), DEFINE_FWK_MODULE, h_ee_invMass_BB, h_ee_invMass_EB, h_ee_invMass_EE, LogTrace, prefixME_, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~EcalZmassTask().

 {
   std::string logTraceName("EcalZmassTask");
 
   h_ee_invMass_EB = nullptr;
   h_ee_invMass_EE = nullptr;
   h_ee_invMass_BB = nullptr;
 
   LogTrace (logTraceName) << "Parameters initialization";
 
   iBooker.setCurrentFolder (prefixME_ + "/Zmass");  // Use folder with name of PAG
 
   h_ee_invMass_EB =
     iBooker.book1D ("Z peak - WP80 EB-EE",
                     "Z peak - WP80 EB-EE;InvMass (GeV)", 60, 60.0, 120.0);
   h_ee_invMass_EE =
     iBooker.book1D ("Z peak - WP80 EE-EE",
                     "Z peak - WP80 EE-EE;InvMass (Gev)", 60, 60.0, 120.0);
   h_ee_invMass_BB =
     iBooker.book1D ("Z peak - WP80 EB-EB",
                     "Z peak - WP80 EB-EB;InvMass (Gev)", 60, 60.0, 120.0);
 }

Member Data Documentation

const edm::EDGetTokenT<reco::GsfElectronCollection> EcalZmassTask::electronCollectionToken_

private

Definition at line 55 of file EcalZmassTask.cc.

Referenced by analyze().

MonitorElement* EcalZmassTask::h_95_ee_invMass_BB

private

Definition at line 72 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_95_ee_invMass_EB

private

Definition at line 70 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_95_ee_invMass_EE

private

Definition at line 71 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_e1_et

private

Definition at line 64 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_e1_eta

private

Definition at line 66 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_e1_phi

private

Definition at line 68 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_e2_et

private

Definition at line 65 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_e2_eta

private

Definition at line 67 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_e2_phi

private

Definition at line 69 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_ee_invMass

private

Definition at line 63 of file EcalZmassTask.cc.

MonitorElement* EcalZmassTask::h_ee_invMass_BB

private

Definition at line 62 of file EcalZmassTask.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* EcalZmassTask::h_ee_invMass_EB

private

Definition at line 60 of file EcalZmassTask.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* EcalZmassTask::h_ee_invMass_EE

private

Definition at line 61 of file EcalZmassTask.cc.

Referenced by analyze(), and bookHistograms().

const std::string EcalZmassTask::prefixME_

private

Definition at line 58 of file EcalZmassTask.cc.

Referenced by bookHistograms().

const edm::EDGetTokenT<reco::GsfTrackCollection> EcalZmassTask::trackCollectionToken_

private

Definition at line 56 of file EcalZmassTask.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation