#include <MCElectronAnalyzer.h>

Inheritance diagram for MCElectronAnalyzer:

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

void	beginJob () override

void	endJob () override

	MCElectronAnalyzer (const edm::ParameterSet &)

	~MCElectronAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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 &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
float	etaTransformation (float a, float b)

float	phiNormalization (float &a)

Private Attributes
TFile *	fOutputFile_

std::string	fOutputFileName_

TH1F *	h_BremEnergy_

TH1F *	h_BremFrac_

TH1F *	h_MCEleE_

TH1F *	h_MCEleEta_

TH1F *	h_MCElePhi_

std::string	HepMCLabel

double	mcEta_

double	mcPhi_
	global variable for the MC photon More...

int	nEvt_

int	nMatched_

TProfile *	p_BremVsEta_

TProfile *	p_BremVsR_

std::string	SimHitLabel

std::string	SimTkLabel

std::string	SimVtxLabel

ElectronMCTruthFinder *	theElectronMCTruthFinder_

const TrackerGeometry *	trackerGeom

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 20 of file MCElectronAnalyzer.h.

Constructor & Destructor Documentation

◆ MCElectronAnalyzer()

MCElectronAnalyzer::MCElectronAnalyzer ( const edm::ParameterSet & pset )

explicit

Definition at line 41 of file MCElectronAnalyzer.cc.

42 : fOutputFileName_(pset.getUntrackedParameter<string>("HistOutFile", std::string("TestConversions.root"))),

43 fOutputFile_(nullptr) {}

◆ ~MCElectronAnalyzer()

MCElectronAnalyzer::~MCElectronAnalyzer ( )

override

Definition at line 45 of file MCElectronAnalyzer.cc.

45 { delete theElectronMCTruthFinder_; }

References theElectronMCTruthFinder_.

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 120 of file MCElectronAnalyzer.cc.

                                                                         {
   using namespace edm;
   //const float etaPhiDistance=0.01;
   // Fiducial region
   //const float TRK_BARL =0.9;
   //const float BARL = 1.4442; // DAQ TDR p.290
   //const float END_LO = 1.566;
   //const float END_HI = 2.5;
   // Electron mass
   //const Float_t mElec= 0.000511;
  
   nEvt_++;
   LogInfo("MCElectronAnalyzer") << "MCElectronAnalyzer Analyzing event number: " << e.id() << " Global Counter "
                                 << nEvt_ << "\n";
   //  LogDebug("MCElectronAnalyzer") << "MCElectronAnalyzer Analyzing event number: "  << e.id() << " Global Counter " << nEvt_ <<"\n";
   std::cout << "MCElectronAnalyzer Analyzing event number: " << e.id() << " Global Counter " << nEvt_ << "\n";
  
   std::cout << " MCElectronAnalyzer Looking for MC truth "
             << "\n";
  
   //get simtrack info
   std::vector<SimTrack> theSimTracks;
   std::vector<SimVertex> theSimVertices;
  
   edm::Handle<SimTrackContainer> SimTk;
   edm::Handle<SimVertexContainer> SimVtx;
   e.getByLabel("g4SimHits", SimTk);
   e.getByLabel("g4SimHits", SimVtx);
  
   theSimTracks.insert(theSimTracks.end(), SimTk->begin(), SimTk->end());
   theSimVertices.insert(theSimVertices.end(), SimVtx->begin(), SimVtx->end());
   std::cout << " MCElectronAnalyzer This Event has " << theSimTracks.size() << " sim tracks " << std::endl;
   std::cout << " MCElectronAnalyzer This Event has " << theSimVertices.size() << " sim vertices " << std::endl;
   if (theSimTracks.empty())
     std::cout << " Event number " << e.id() << " has NO sim tracks " << std::endl;
  
   std::vector<ElectronMCTruth> MCElectronctrons = theElectronMCTruthFinder_->find(theSimTracks, theSimVertices);
   std::cout << " MCElectronAnalyzer MCElectronctrons size " << MCElectronctrons.size() << std::endl;
  
   for (std::vector<ElectronMCTruth>::const_iterator iEl = MCElectronctrons.begin(); iEl != MCElectronctrons.end();
        ++iEl) {
     h_MCEleE_->Fill((*iEl).fourMomentum().e());
     h_MCEleEta_->Fill((*iEl).fourMomentum().pseudoRapidity());
     h_MCElePhi_->Fill((*iEl).fourMomentum().phi());
  
     float totBrem = 0;
     unsigned int iBrem;
     for (iBrem = 0; iBrem < (*iEl).bremVertices().size(); ++iBrem) {
       float rBrem = (*iEl).bremVertices()[iBrem].perp();
       float etaBrem = (*iEl).bremVertices()[iBrem].eta();
       if (rBrem < 120) {
         totBrem += (*iEl).bremMomentum()[iBrem].e();
         p_BremVsR_->Fill(rBrem, (*iEl).bremMomentum()[iBrem].e());
         p_BremVsEta_->Fill(etaBrem, (*iEl).bremMomentum()[iBrem].e());
       }
     }
  
     h_BremFrac_->Fill(totBrem / (*iEl).fourMomentum().e());
     h_BremEnergy_->Fill(totBrem);
   }
 }

References gather_cfg::cout, MillePedeFileConverter_cfg::e, ElectronMCTruthFinder::find(), h_BremEnergy_, h_BremFrac_, h_MCEleE_, h_MCEleEta_, h_MCElePhi_, nEvt_, p_BremVsEta_, p_BremVsR_, and theElectronMCTruthFinder_.

◆ beginJob()

void MCElectronAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 47 of file MCElectronAnalyzer.cc.

                                   {
   nEvt_ = 0;
  
   theElectronMCTruthFinder_ = new ElectronMCTruthFinder();
  
   fOutputFile_ = new TFile(fOutputFileName_.c_str(), "RECREATE");
  
   h_MCEleE_ = new TH1F("MCEleE", "MC ele energy", 100, 0., 200.);
   h_MCElePhi_ = new TH1F("MCElePhi", "MC ele phi", 40, -3.14, 3.14);
   h_MCEleEta_ = new TH1F("MCEleEta", "MC ele eta", 40, -3., 3.);
   h_BremFrac_ = new TH1F("bremFrac", "brem frac ", 100, 0., 1.);
   h_BremEnergy_ = new TH1F("BremE", "Brem energy", 100, 0., 200.);
  
   p_BremVsR_ = new TProfile("BremVsR", " Mean Brem energy vs R ", 48, 0., 120.);
   p_BremVsEta_ = new TProfile("BremVsEta", " Mean Brem energy vs Eta ", 50, -2.5, 2.5);
  
   return;
 }

References fOutputFile_, fOutputFileName_, h_BremEnergy_, h_BremFrac_, h_MCEleE_, h_MCEleEta_, h_MCElePhi_, nEvt_, p_BremVsEta_, p_BremVsR_, and theElectronMCTruthFinder_.

◆ endJob()

void MCElectronAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 183 of file MCElectronAnalyzer.cc.

                                 {
   fOutputFile_->Write();
   fOutputFile_->Close();
  
   edm::LogInfo("MCElectronAnalyzer") << "Analyzed " << nEvt_ << "\n";
   std::cout << "MCElectronAnalyzer::endJob Analyzed " << nEvt_ << " events "
             << "\n";
  
   return;
 }

References gather_cfg::cout, fOutputFile_, and nEvt_.

◆ etaTransformation()

float MCElectronAnalyzer::etaTransformation	(	float	a,
		float	b
	)

private

Definition at line 67 of file MCElectronAnalyzer.cc.

                                                                             {
   //---Definitions
   const float PI = 3.1415927;
   //const float TWOPI = 2.0*PI;
  
   //---Definitions for ECAL
   const float R_ECAL = 136.5;
   const float Z_Endcap = 328.0;
   const float etaBarrelEndcap = 1.479;
  
   //---ETA correction
  
   float Theta = 0.0;
   float ZEcal = R_ECAL * sinh(EtaParticle) + Zvertex;
  
   if (ZEcal != 0.0)
     Theta = atan(R_ECAL / ZEcal);
   if (Theta < 0.0)
     Theta = Theta + PI;
   float ETA = -log(tan(0.5 * Theta));
  
   if (std::abs(ETA) > etaBarrelEndcap) {
     float Zend = Z_Endcap;
     if (EtaParticle < 0.0)
       Zend = -Zend;
     float Zlen = Zend - Zvertex;
     float RR = Zlen / sinh(EtaParticle);
     Theta = atan(RR / Zend);
     if (Theta < 0.0)
       Theta = Theta + PI;
     ETA = -log(tan(0.5 * Theta));
   }
   //---Return the result
   return ETA;
   //---end
 }

References funct::abs(), ETA, etaBarrelEndcap, dqm-mbProfile::log, PI, R_ECAL, funct::tan(), Z_Endcap, and ZEcal.

◆ phiNormalization()

float MCElectronAnalyzer::phiNormalization ( float & a )

private

Definition at line 104 of file MCElectronAnalyzer.cc.

                                                      {
   //---Definitions
   const float PI = 3.1415927;
   const float TWOPI = 2.0 * PI;
  
   if (phi > PI) {
     phi = phi - TWOPI;
   }
   if (phi < -PI) {
     phi = phi + TWOPI;
   }
  
   //  cout << " Float_t PHInormalization out " << PHI << endl;
   return phi;
 }

References phi, PI, and TWOPI.

Member Data Documentation

◆ fOutputFile_

TFile* MCElectronAnalyzer::fOutputFile_

private

Definition at line 40 of file MCElectronAnalyzer.h.

Referenced by beginJob(), and endJob().

◆ fOutputFileName_

std::string MCElectronAnalyzer::fOutputFileName_

private

Definition at line 39 of file MCElectronAnalyzer.h.

Referenced by beginJob().

◆ h_BremEnergy_

TH1F* MCElectronAnalyzer::h_BremEnergy_

private

Definition at line 58 of file MCElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h_BremFrac_

TH1F* MCElectronAnalyzer::h_BremFrac_

private

Definition at line 57 of file MCElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h_MCEleE_

TH1F* MCElectronAnalyzer::h_MCEleE_

private

Definition at line 54 of file MCElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h_MCEleEta_

TH1F* MCElectronAnalyzer::h_MCEleEta_

private

Definition at line 55 of file MCElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ h_MCElePhi_

TH1F* MCElectronAnalyzer::h_MCElePhi_

private

Definition at line 56 of file MCElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ HepMCLabel

std::string MCElectronAnalyzer::HepMCLabel

private

Definition at line 49 of file MCElectronAnalyzer.h.

◆ mcEta_

double MCElectronAnalyzer::mcEta_

private

Definition at line 47 of file MCElectronAnalyzer.h.

◆ mcPhi_

double MCElectronAnalyzer::mcPhi_

private

global variable for the MC photon

Definition at line 46 of file MCElectronAnalyzer.h.

◆ nEvt_

int MCElectronAnalyzer::nEvt_

private

Definition at line 42 of file MCElectronAnalyzer.h.

Referenced by analyze(), beginJob(), and endJob().

◆ nMatched_

int MCElectronAnalyzer::nMatched_

private

Definition at line 43 of file MCElectronAnalyzer.h.

◆ p_BremVsEta_

TProfile* MCElectronAnalyzer::p_BremVsEta_

private

Definition at line 61 of file MCElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ p_BremVsR_

TProfile* MCElectronAnalyzer::p_BremVsR_

private

Definition at line 60 of file MCElectronAnalyzer.h.

Referenced by analyze(), and beginJob().

◆ SimHitLabel

std::string MCElectronAnalyzer::SimHitLabel

private

Definition at line 52 of file MCElectronAnalyzer.h.

◆ SimTkLabel

std::string MCElectronAnalyzer::SimTkLabel

private

Definition at line 50 of file MCElectronAnalyzer.h.

◆ SimVtxLabel

std::string MCElectronAnalyzer::SimVtxLabel

private

Definition at line 51 of file MCElectronAnalyzer.h.

◆ theElectronMCTruthFinder_

ElectronMCTruthFinder* MCElectronAnalyzer::theElectronMCTruthFinder_

private

Definition at line 35 of file MCElectronAnalyzer.h.

Referenced by analyze(), beginJob(), and ~MCElectronAnalyzer().

◆ trackerGeom

const TrackerGeometry* MCElectronAnalyzer::trackerGeom

private

Definition at line 37 of file MCElectronAnalyzer.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MCElectronAnalyzer()

◆ ~MCElectronAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ etaTransformation()

◆ phiNormalization()

Member Data Documentation

◆ fOutputFile_

◆ fOutputFileName_

◆ h_BremEnergy_

◆ h_BremFrac_

◆ h_MCEleE_

◆ h_MCEleEta_

◆ h_MCElePhi_

◆ HepMCLabel

◆ mcEta_

◆ mcPhi_

◆ nEvt_

◆ nMatched_

◆ p_BremVsEta_

◆ p_BremVsR_

◆ SimHitLabel

◆ SimTkLabel

◆ SimVtxLabel

◆ theElectronMCTruthFinder_

◆ trackerGeom