#include <Calibration/EcalCalibAlgos/src/ElectronCalibration.cc>

Inheritance diagram for ElectronCalibration:

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

void	beginJob () override

	ElectronCalibration (const edm::ParameterSet &)

void	endJob () override

	~ElectronCalibration () override

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
EBDetId	findMaxHit (edm::Handle< EBRecHitCollection > &)

EBDetId	findMaxHit2 (const std::vector< DetId > &, const EBRecHitCollection *)

Private Attributes
std::string	calibAlgo_

CalibrationCluster	calibCluster

int	calibClusterSize

TH1F *	calibinter

TH1F *	calibinterNoCuts

TH1F *	calibs

TH1F *	calibsNoCuts

int	ClusterSize_

double	cut1_

double	cut2_

double	cut3_

double	cutCalo1_

double	cutCalo2_

double	cutEPCalo1_

double	cutEPCalo2_

double	cutEPin1_

double	cutEPin2_

double	cutESeed_

TH1F *	DiffCalibs

TH1F *	DiffCalibsNoCuts

TH1F *	e25

TH1F *	e25NoCuts

TH1F *	E25oP

TH1F *	E25oPNoCuts

TH1F *	e25OverScE

TH1F *	e25OverScENoCuts

TH1F *	e9

TH1F *	e9NoCuts

TH1F *	e9Overe25

TH1F *	e9Overe25NoCuts

int	elecclass_

edm::EDGetTokenT< reco::GsfElectronCollection >	electronToken_

double	ElePt_

std::vector< float >	EnergyVector

std::vector< float >	EnergyVectorNoCuts

TH1F *	EoP

TH1F *	EoP_all

TH1F *	EoPNoCuts

TH1F *	Error1

TH1F *	Error1NoCuts

TH1F *	Error2

TH1F *	Error2NoCuts

TH1F *	Error3

TH1F *	Error3NoCuts

TH1F *	eSeedOverPout

TH1F *	eSeedOverPout2

TH1F *	eSeedOverPout2ESeed

TH1F *	eSeedOverPout2NoCuts

TH1F *	eSeedOverPoutNoCuts

int	etaMax

int	etaMin

int	eventcrystal [171][360]

std::vector< std::vector< float > >	EventMatrix

std::vector< std::vector< float > >	EventMatrixNoCuts

TFile *	f

TH2F *	GeneralMap

TH1F *	hadOverEm

TH1F *	hadOverEmESeed

TH1F *	hadOverEmNoCuts

int	keventweight

int	keventweight_

TH2F *	Map

TH2F *	Map3Dcalib

TH2F *	Map3DcalibNoCuts

TH2F *	MapCor1

TH2F *	MapCor10

TH2F *	MapCor10ESeed

TH2F *	MapCor10NoCuts

TH2F *	MapCor11

TH2F *	MapCor11ESeed

TH2F *	MapCor11NoCuts

TH2F *	MapCor1ESeed

TH2F *	MapCor1NoCuts

TH2F *	MapCor2

TH2F *	MapCor2ESeed

TH2F *	MapCor2NoCuts

TH2F *	MapCor3

TH2F *	MapCor3ESeed

TH2F *	MapCor3NoCuts

TH2F *	MapCor4

TH2F *	MapCor4ESeed

TH2F *	MapCor4NoCuts

TH2F *	MapCor5

TH2F *	MapCor5ESeed

TH2F *	MapCor5NoCuts

TH2F *	MapCor6

TH2F *	MapCor6ESeed

TH2F *	MapCor6NoCuts

TH2F *	MapCor7

TH2F *	MapCor7ESeed

TH2F *	MapCor7NoCuts

TH2F *	MapCor8

TH2F *	MapCor8ESeed

TH2F *	MapCor8NoCuts

TH2F *	MapCor9

TH2F *	MapCor9ESeed

TH2F *	MapCor9NoCuts

TH2F *	MapCorCalib

TH2F *	MapCorCalibNoCuts

TH2F *	MapNoCuts

std::vector< int >	MaxCCeta

std::vector< int >	MaxCCetaNoCuts

std::vector< int >	MaxCCphi

std::vector< int >	MaxCCphiNoCuts

int	maxeta_

int	maxphi_

int	mineta_

int	minphi_

std::string	miscalibfile_

TH1F *	MisCalibs

HouseholderDecomposition *	MyHH

MinL3Algorithm *	MyL3Algo1

std::vector< float >	newCalibs

int	numevent_

std::vector< float >	oldCalibs

int	phiMax

int	phiMin

TH1F *	PinMinPout

TH1F *	PinMinPoutNoCuts

TH1F *	PinOverPout

TH1F *	PinOverPoutNoCuts

TH1F *	RatioCalibs

TH1F *	RatioCalibsNoCuts

int	read_events

edm::EDGetTokenT< EBRecHitCollection >	recHitToken_

CalibrationCluster::CalibMap	ReducedMap

std::string	rootfile_

TH1F *	scE

TH1F *	scENoCuts

std::vector< float >	solution

std::vector< float >	solutionNoCuts

edm::InputTag	trackLabel_

TH1F *	trP

TH1F *	trPNoCuts

std::vector< float >	WeightVector

std::vector< float >	WeightVectorNoCuts

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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: Perform single electron calibration (tested on TB data only).

Implementation: <Notes on="" implementation>="">

Definition at line 46 of file ElectronCalibration.h.

Constructor & Destructor Documentation

◆ ElectronCalibration()

ElectronCalibration::ElectronCalibration ( const edm::ParameterSet & iConfig )

explicit

Definition at line 33 of file ElectronCalibration.cc.

References calibAlgo_, ClusterSize_, cut1_, cut2_, cut3_, cutCalo1_, cutCalo2_, cutEPCalo1_, cutEPCalo2_, cutEPin1_, cutEPin2_, cutESeed_, elecclass_, electronToken_, ElePt_, edm::ParameterSet::getParameter(), keventweight_, maxeta_, maxphi_, mineta_, minphi_, miscalibfile_, numevent_, recHitToken_, rootfile_, AlCaHLTBitMon_QueryRunRegistry::string, and trackLabel_.

                                                                        {
   rootfile_ = iConfig.getParameter<std::string>("rootfile");
   recHitToken_ = consumes<EBRecHitCollection>(iConfig.getParameter<edm::InputTag>("ebRecHitsLabel"));
   electronToken_ = consumes<reco::GsfElectronCollection>(iConfig.getParameter<edm::InputTag>("electronLabel"));
   trackLabel_ = iConfig.getParameter<edm::InputTag>("trackLabel");
   calibAlgo_ = iConfig.getParameter<std::string>("CALIBRATION_ALGO");
   edm::LogVerbatim("EcalCalib") << " The used Algorithm is  " << calibAlgo_;
   keventweight_ = iConfig.getParameter<int>("keventweight");
   ClusterSize_ = iConfig.getParameter<int>("Clustersize");
   ElePt_ = iConfig.getParameter<double>("ElePt");
   maxeta_ = iConfig.getParameter<int>("maxeta");
   mineta_ = iConfig.getParameter<int>("mineta");
   maxphi_ = iConfig.getParameter<int>("maxphi");
   minphi_ = iConfig.getParameter<int>("minphi");
   cut1_ = iConfig.getParameter<double>("cut1");
   cut2_ = iConfig.getParameter<double>("cut2");
   cut3_ = iConfig.getParameter<double>("cut3");
   elecclass_ = iConfig.getParameter<int>("elecclass");
   edm::LogVerbatim("EcalCalib") << " The electronclass is " << elecclass_;
   numevent_ = iConfig.getParameter<int>("numevent");
   miscalibfile_ = iConfig.getParameter<std::string>("miscalibfile");
 
   cutEPCalo1_ = iConfig.getParameter<double>("cutEPCaloMin");
   cutEPCalo2_ = iConfig.getParameter<double>("cutEPCaloMax");
   cutEPin1_ = iConfig.getParameter<double>("cutEPinMin");
   cutEPin2_ = iConfig.getParameter<double>("cutEPinMax");
   cutCalo1_ = iConfig.getParameter<double>("cutCaloMin");
   cutCalo2_ = iConfig.getParameter<double>("cutCaloMax");
 
   cutESeed_ = iConfig.getParameter<double>("cutESeed");
 }

◆ ~ElectronCalibration()

ElectronCalibration::~ElectronCalibration ( )

override

Definition at line 65 of file ElectronCalibration.cc.

65 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 452 of file ElectronCalibration.cc.

                                                                                      {
   //=================================================================================
   using namespace edm;
 
   // Get EBRecHits
   const Handle<EBRecHitCollection>& phits = iEvent.getHandle(recHitToken_);
   if (!phits.isValid()) {
     edm::LogError("EcalCalib") << "Error! can't get the product EBRecHitCollection: ";
   }
 
   const EBRecHitCollection* hits = phits.product();  // get a ptr to the product
 
   // Get pixelElectrons
   const Handle<reco::GsfElectronCollection>& pElectrons = iEvent.getHandle(electronToken_);
   if (!pElectrons.isValid()) {
     edm::LogError("EcalCalib") << "Error! can't get the product ElectronCollection: ";
   }
 
   const reco::GsfElectronCollection* electronCollection = pElectrons.product();
   read_events++;
   if (read_events % 1000 == 0)
     edm::LogVerbatim("EcalCalib") << "read_events = " << read_events << std::endl;
 
   if (!hits)
     return;
   if (hits->empty())
     return;
   if (!electronCollection)
     return;
   if (electronCollection->empty())
     return;
 
   //                          START HERE....
   reco::GsfElectronCollection::const_iterator eleIt = electronCollection->begin();
 
   reco::GsfElectron highPtElectron;
 
   float highestElePt = 0.;
   bool found = false;
   for (eleIt = electronCollection->begin(); eleIt != electronCollection->end(); eleIt++) {
     //Comments
     if (fabs(eleIt->eta()) > (maxeta_ + 3) * 0.0175)
       continue;
     if (eleIt->eta() < (mineta_ - 3) * 0.0175)
       continue;
 
     if (eleIt->pt() > highestElePt) {
       highestElePt = eleIt->pt();
       highPtElectron = *eleIt;
       found = true;
     }
   }
   if (highestElePt < ElePt_)
     return;
   if (!found)
     return;
   const reco::SuperCluster& sc = *(highPtElectron.superCluster());
   if (fabs(sc.eta()) > (maxeta_ + 3) * 0.0175) {
     edm::LogVerbatim("EcalCalib") << "++++ Problem with electron, electron eta is " << highPtElectron.eta()
                                   << " while SC is " << sc.eta() << std::endl;
     return;
   }
 
   std::vector<DetId> v1;
   //Loop to fill the vector of DetIds
   for (std::vector<std::pair<DetId, float> >::const_iterator idsIt = sc.hitsAndFractions().begin();
        idsIt != sc.hitsAndFractions().end();
        ++idsIt) {
     v1.push_back(idsIt->first);
   }
 
   //Change function name
   EBDetId maxHitId;
 
   maxHitId = findMaxHit2(v1, hits);
 
   if (maxHitId.null()) {
     edm::LogVerbatim("EcalCalib") << " Null " << std::endl;
     return;
   }
 
   int maxCC_Eta = maxHitId.ieta();
   int maxCC_Phi = maxHitId.iphi();
 
   if (maxCC_Eta > maxeta_)
     return;
   if (maxCC_Eta < mineta_)
     return;
   if (maxCC_Phi > maxphi_)
     return;
   if (maxCC_Phi < minphi_)
     return;
 
   // number of events per crystal is set
   if (numevent_ > 0) {
     eventcrystal[maxCC_Eta + 85][maxCC_Phi - 1] += 1;
     if (eventcrystal[maxCC_Eta + 85][maxCC_Phi - 1] > numevent_)
       return;
   }
 
   std::vector<EBDetId> Xtals5x5 = calibCluster.get5x5Id(maxHitId);
 
   if ((int)Xtals5x5.size() != ClusterSize_ * ClusterSize_)
     return;
 
   // fill cluster energy
   std::vector<float> energy;
   float energy3x3 = 0.;
   float energy5x5 = 0.;
 
   for (int icry = 0; icry < ClusterSize_ * ClusterSize_; icry++) {
     EBRecHitCollection::const_iterator itrechit;
     if (Xtals5x5[icry].subdetId() != 1)
       continue;
     itrechit = hits->find(Xtals5x5[icry]);
     if (itrechit == hits->end()) {
       edm::LogVerbatim("EcalCalib") << "DetId not is e25";
       continue;
     }
 
     if (edm::isNotFinite(itrechit->energy()))
       return;
     energy.push_back(itrechit->energy());
     energy5x5 += energy[icry];
 
     if (icry == 6 || icry == 7 || icry == 8 || icry == 11 || icry == 12 || icry == 13 || icry == 16 || icry == 17 ||
         icry == 18) {
       energy3x3 += energy[icry];
     }
   }
   if ((int)energy.size() != ClusterSize_ * ClusterSize_)
     return;
   // Once we have the matrix 5x5, we have to correct for gaps/cracks/umbrella and maincontainement
 
   GeneralMap->Fill(maxCC_Eta, maxCC_Phi);
 
   EoP_all->Fill(highPtElectron.eSuperClusterOverP());
 
   if (highPtElectron.classification() == elecclass_ || elecclass_ == -1) {
     float Ptrack_in =
         sqrt(pow(highPtElectron.trackMomentumAtVtx().X(), 2) + pow(highPtElectron.trackMomentumAtVtx().Y(), 2) +
              pow(highPtElectron.trackMomentumAtVtx().Z(), 2));
 
     float UncorrectedPatCalo =
         sqrt(pow(highPtElectron.trackMomentumAtCalo().X(), 2) + pow(highPtElectron.trackMomentumAtCalo().Y(), 2) +
              pow(highPtElectron.trackMomentumAtCalo().Z(), 2));
 
     float Ptrack_out =
         sqrt(pow(highPtElectron.trackMomentumOut().X(), 2) + pow(highPtElectron.trackMomentumOut().Y(), 2) +
              pow(highPtElectron.trackMomentumOut().Z(), 2));
 
     EventMatrixNoCuts.push_back(energy);
     EnergyVectorNoCuts.push_back(UncorrectedPatCalo);
 
     MaxCCetaNoCuts.push_back(maxCC_Eta);
     MaxCCphiNoCuts.push_back(maxCC_Phi);
 
     WeightVectorNoCuts.push_back(energy5x5 / UncorrectedPatCalo);
 
     //---------------------------------------------------No Cuts-------------------------------------------------------
     e9NoCuts->Fill(energy3x3);
     e25NoCuts->Fill(energy5x5);
     e9Overe25NoCuts->Fill(energy3x3 / energy5x5);
     scENoCuts->Fill(sc.energy());
 
     trPNoCuts->Fill(UncorrectedPatCalo);
 
     EoPNoCuts->Fill(highPtElectron.eSuperClusterOverP());
     e25OverScENoCuts->Fill(energy5x5 / sc.energy());
 
     E25oPNoCuts->Fill(energy5x5 / UncorrectedPatCalo);
 
     MapNoCuts->Fill(maxCC_Eta, maxCC_Phi);
     PinOverPoutNoCuts->Fill(
         sqrt(pow(highPtElectron.trackMomentumAtVtx().X(), 2) + pow(highPtElectron.trackMomentumAtVtx().Y(), 2) +
              pow(highPtElectron.trackMomentumAtVtx().Z(), 2)) /
         sqrt(pow(highPtElectron.trackMomentumOut().X(), 2) + pow(highPtElectron.trackMomentumOut().Y(), 2) +
              pow(highPtElectron.trackMomentumOut().Z(), 2)));
     eSeedOverPoutNoCuts->Fill(highPtElectron.eSuperClusterOverP());
 
     MapCor1NoCuts->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / Ptrack_in);
     MapCor2NoCuts->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSuperClusterOverP());
     MapCor3NoCuts->Fill(energy5x5 / UncorrectedPatCalo, Ptrack_out / Ptrack_in);
     MapCor4NoCuts->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / highPtElectron.p());
     MapCor5NoCuts->Fill(energy5x5 / UncorrectedPatCalo, UncorrectedPatCalo / Ptrack_out);
     MapCor6NoCuts->Fill(Ptrack_out / Ptrack_in, energy5x5 / Ptrack_in);
     MapCor7NoCuts->Fill(Ptrack_out / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
     MapCor8NoCuts->Fill(energy5x5 / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
     MapCor9NoCuts->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSeedClusterOverPout());
     MapCor10NoCuts->Fill(highPtElectron.eSeedClusterOverPout(), Ptrack_out / Ptrack_in);
     MapCor11NoCuts->Fill(highPtElectron.eSeedClusterOverPout(), energy5x5 / Ptrack_in);
 
     PinMinPoutNoCuts->Fill((Ptrack_in - Ptrack_out) / Ptrack_in);
 
     Error1NoCuts->Fill(highPtElectron.trackMomentumError() / Ptrack_in);
     Error2NoCuts->Fill(highPtElectron.trackMomentumError() / Ptrack_out);
     Error3NoCuts->Fill(highPtElectron.trackMomentumError() / UncorrectedPatCalo);
     eSeedOverPout2NoCuts->Fill(highPtElectron.eSeedClusterOverPout());
 
     hadOverEmNoCuts->Fill(highPtElectron.hadronicOverEm());
 
     //------------------------------------------------Cuts-----------------------------------------------------
     //Cuts!
     if ((energy3x3 / energy5x5) < cut1_)
       return;
 
     if ((Ptrack_out / Ptrack_in) < cut2_ || (Ptrack_out / Ptrack_in) > cut3_)
       return;
     if ((energy5x5 / Ptrack_in) < cutEPin1_ || (energy5x5 / Ptrack_in) > cutEPin2_)
       return;
 
     e9->Fill(energy3x3);
     e25->Fill(energy5x5);
     e9Overe25->Fill(energy3x3 / energy5x5);
     scE->Fill(sc.energy());
     trP->Fill(UncorrectedPatCalo);
 
     EoP->Fill(highPtElectron.eSuperClusterOverP());
     e25OverScE->Fill(energy5x5 / sc.energy());
 
     E25oP->Fill(energy5x5 / UncorrectedPatCalo);
 
     Map->Fill(maxCC_Eta, maxCC_Phi);
     PinOverPout->Fill(
         sqrt(pow(highPtElectron.trackMomentumAtVtx().X(), 2) + pow(highPtElectron.trackMomentumAtVtx().Y(), 2) +
              pow(highPtElectron.trackMomentumAtVtx().Z(), 2)) /
         sqrt(pow(highPtElectron.trackMomentumOut().X(), 2) + pow(highPtElectron.trackMomentumOut().Y(), 2) +
              pow(highPtElectron.trackMomentumOut().Z(), 2)));
     eSeedOverPout->Fill(highPtElectron.eSuperClusterOverP());
 
     MapCor1->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / Ptrack_in);
     MapCor2->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSuperClusterOverP());
     MapCor3->Fill(energy5x5 / UncorrectedPatCalo, Ptrack_out / Ptrack_in);
     MapCor4->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / highPtElectron.p());
     MapCor5->Fill(energy5x5 / UncorrectedPatCalo, UncorrectedPatCalo / Ptrack_out);
     MapCor6->Fill(Ptrack_out / Ptrack_in, energy5x5 / Ptrack_in);
     MapCor7->Fill(Ptrack_out / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
     MapCor8->Fill(energy5x5 / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
     MapCor9->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSeedClusterOverPout());
     MapCor10->Fill(highPtElectron.eSeedClusterOverPout(), Ptrack_out / Ptrack_in);
     MapCor11->Fill(highPtElectron.eSeedClusterOverPout(), energy5x5 / Ptrack_in);
 
     PinMinPout->Fill((Ptrack_in - Ptrack_out) / Ptrack_in);
 
     Error1->Fill(highPtElectron.trackMomentumError() / Ptrack_in);
     Error2->Fill(highPtElectron.trackMomentumError() / Ptrack_out);
     Error3->Fill(highPtElectron.trackMomentumError() / UncorrectedPatCalo);
 
     eSeedOverPout2->Fill(highPtElectron.eSeedClusterOverPout());
     hadOverEm->Fill(highPtElectron.hadronicOverEm());
 
     EventMatrix.push_back(energy);
     EnergyVector.push_back(UncorrectedPatCalo);
     MaxCCeta.push_back(maxCC_Eta);
     MaxCCphi.push_back(maxCC_Phi);
 
     WeightVector.push_back(energy5x5 / UncorrectedPatCalo);
 
     //-------------------------------------------------------Extra Cut-----------------------------------------------------
     if (highPtElectron.eSeedClusterOverPout() < cutESeed_)
       return;
 
     MapCor1ESeed->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / Ptrack_in);
     MapCor2ESeed->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSuperClusterOverP());
     MapCor3ESeed->Fill(energy5x5 / UncorrectedPatCalo, Ptrack_out / Ptrack_in);
     MapCor4ESeed->Fill(energy5x5 / UncorrectedPatCalo, energy5x5 / highPtElectron.p());
     MapCor5ESeed->Fill(energy5x5 / UncorrectedPatCalo, UncorrectedPatCalo / Ptrack_out);
     MapCor6ESeed->Fill(Ptrack_out / Ptrack_in, energy5x5 / Ptrack_in);
     MapCor7ESeed->Fill(Ptrack_out / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
     MapCor8ESeed->Fill(energy5x5 / Ptrack_in, UncorrectedPatCalo / Ptrack_out);
     MapCor9ESeed->Fill(energy5x5 / UncorrectedPatCalo, highPtElectron.eSeedClusterOverPout());
     MapCor10ESeed->Fill(highPtElectron.eSeedClusterOverPout(), Ptrack_out / Ptrack_in);
     MapCor11ESeed->Fill(highPtElectron.eSeedClusterOverPout(), energy5x5 / Ptrack_in);
 
     eSeedOverPout2ESeed->Fill(highPtElectron.eSeedClusterOverPout());
 
     hadOverEmESeed->Fill(highPtElectron.hadronicOverEm());
 
   } else {
     return;
   }
 }

◆ beginJob()

void ElectronCalibration::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 68 of file ElectronCalibration.cc.

References calibAlgo_, calibCluster, calibClusterSize, calibinter, calibinterNoCuts, calibs, calibsNoCuts, ClusterSize_, DiffCalibs, DiffCalibsNoCuts, e25, e25NoCuts, E25oP, E25oPNoCuts, e25OverScE, e25OverScENoCuts, e9, e9NoCuts, e9Overe25, e9Overe25NoCuts, elecclass_, EoP, EoP_all, EoPNoCuts, Error1, Error1NoCuts, Error2, Error2NoCuts, Error3, Error3NoCuts, eSeedOverPout, eSeedOverPout2, eSeedOverPout2ESeed, eSeedOverPout2NoCuts, eSeedOverPoutNoCuts, PVValHelper::eta, etaMax, etaMin, eventcrystal, f, GeneralMap, CalibrationCluster::getMap(), hadOverEm, hadOverEmESeed, hadOverEmNoCuts, keventweight_, Map, Map3Dcalib, Map3DcalibNoCuts, MapCor1, MapCor10, MapCor10ESeed, MapCor10NoCuts, MapCor11, MapCor11ESeed, MapCor11NoCuts, MapCor1ESeed, MapCor1NoCuts, MapCor2, MapCor2ESeed, MapCor2NoCuts, MapCor3, MapCor3ESeed, MapCor3NoCuts, MapCor4, MapCor4ESeed, MapCor4NoCuts, MapCor5, MapCor5ESeed, MapCor5NoCuts, MapCor6, MapCor6ESeed, MapCor6NoCuts, MapCor7, MapCor7ESeed, MapCor7NoCuts, MapCor8, MapCor8ESeed, MapCor8NoCuts, MapCor9, MapCor9ESeed, MapCor9NoCuts, MapCorCalib, MapCorCalibNoCuts, MapNoCuts, maxeta_, maxphi_, mineta_, minphi_, MisCalibs, MyHH, MyL3Algo1, oldCalibs, phi, phiMax, phiMin, PinMinPout, PinMinPoutNoCuts, PinOverPout, PinOverPoutNoCuts, RatioCalibs, RatioCalibsNoCuts, read_events, ReducedMap, rootfile_, scE, scENoCuts, trP, and trPNoCuts.

                                    {
   //========================================================================
   f = new TFile(rootfile_.c_str(), "RECREATE");
 
   // Book histograms
   e9 = new TH1F("e9", "E9 energy", 300, 0., 150.);
   e25 = new TH1F("e25", "E25 energy", 300, 0., 150.);
   scE = new TH1F("scE", "SC energy", 300, 0., 150.);
   trP = new TH1F("trP", "Trk momentum", 300, 0., 150.);
   EoP = new TH1F("EoP", "EoP", 600, 0., 3.);
   EoP_all = new TH1F("EoP_all", "EoP_all", 600, 0., 3.);
 
   if (elecclass_ == 0 || elecclass_ == -1) {
     calibs = new TH1F("calib", "Calibration constants", 4000, 0.5, 2.);
   } else {
     calibs = new TH1F("calib", "Calibration constants", 800, 0.5, 2.);
   }
 
   e25OverScE = new TH1F("e25OverscE", "E25 / SC energy", 400, 0., 2.);
   E25oP = new TH1F("E25oP", "E25 / P", 1200, 0., 1.5);
 
   Map = new TH2F("Map", "Nb Events in Crystal", 85, 1, 85, 70, 5, 75);
   e9Overe25 = new TH1F("e9Overe25", "E9 / E25", 400, 0., 2.);
   Map3Dcalib = new TH2F("3Dcalib", "3Dcalib", 85, 1, 85, 70, 5, 75);
 
   MapCor1 = new TH2F("MapCor1", "Correlation E25/Pcalo versus E25/Pin", 100, 0., 5., 100, 0., 5.);
   MapCor2 = new TH2F("MapCor2", "Correlation E25/Pcalo versus E/P", 100, 0., 5., 100, 0., 5.);
   MapCor3 = new TH2F("MapCor3", "Correlation E25/Pcalo versus Pout/Pin", 100, 0., 5., 100, 0., 5.);
   MapCor4 = new TH2F("MapCor4", "Correlation E25/Pcalo versus E25/highestP", 100, 0., 5., 100, 0., 5.);
   MapCor5 = new TH2F("MapCor5", "Correlation E25/Pcalo versus Pcalo/Pout", 100, 0., 5., 100, 0., 5.);
   MapCor6 = new TH2F("MapCor6", "Correlation Pout/Pin versus E25/Pin", 100, 0., 5., 100, 0., 5.);
   MapCor7 = new TH2F("MapCor7", "Correlation Pout/Pin versus Pcalo/Pout", 100, 0., 5., 100, 0., 5.);
   MapCor8 = new TH2F("MapCor8", "Correlation E25/Pin versus Pcalo/Pout", 100, 0., 5., 100, 0., 5.);
   MapCor9 = new TH2F("MapCor9", "Correlation  E25/Pcalo versus Eseed/Pout", 100, 0., 5., 100, 0., 5.);
   MapCor10 = new TH2F("MapCor10", "Correlation Eseed/Pout versus Pout/Pin", 100, 0., 5., 100, 0., 5.);
   MapCor11 = new TH2F("MapCor11", "Correlation Eseed/Pout versus E25/Pin", 100, 0., 5., 100, 0., 5.);
   MapCorCalib =
       new TH2F("MapCorCalib", "Correlation Miscalibration versus Calibration constants", 100, 0.5, 1.5, 100, 0.5, 1.5);
 
   PinMinPout = new TH1F("PinMinPout", "(Pin - Pout)/Pin", 600, -2.0, 2.0);
 
   if (elecclass_ == 0 || elecclass_ == -1) {
     calibinter = new TH1F("calibinter", "internal calibration constants", 2000, 0.5, 2.);
     PinOverPout = new TH1F("PinOverPout", "pinOverpout", 600, 0., 3.);
     eSeedOverPout = new TH1F("eSeedOverPout", "eSeedOverpout ", 600, 0., 3.);
     MisCalibs = new TH1F("MisCalibs", "Miscalibration constants", 4000, 0.5, 2.);
     RatioCalibs = new TH1F("RatioCalibs", "Ratio in Calibration Constants", 4000, 0.5, 2.0);
     DiffCalibs = new TH1F("DiffCalibs", "Difference in Calibration constants", 4000, -1.0, 1.0);
   } else {
     calibinter = new TH1F("calibinter", "internal calibration constants", 400, 0.5, 2.);
     PinOverPout = new TH1F("PinOverPout", "pinOverpout", 600, 0., 3.);
     eSeedOverPout = new TH1F("eSeedOverPout", "eSeedOverpout ", 600, 0., 3.);
     MisCalibs = new TH1F("MisCalibs", "Miscalibration constants", 800, 0.5, 2.);
     RatioCalibs = new TH1F("RatioCalibs", "Ratio in Calibration Constants", 800, 0.5, 2.0);
     DiffCalibs = new TH1F("DiffCalibs", "Difference in Calibration constants", 800, -1.0, 1.0);
   }
   Error1 = new TH1F("Error1", "DeltaP/Pin", 800, -1.0, 1.0);
   Error2 = new TH1F("Error2", "DeltaP/Pout", 800, -1.0, 1.0);
   Error3 = new TH1F("Error3", "DeltaP/Pcalo", 800, -1.0, 1.0);
   eSeedOverPout2 = new TH1F("eSeedOverPout2", "eSeedOverpout (No Supercluster)", 600, 0., 4.);
   hadOverEm = new TH1F("hadOverEm", "Had/EM distribution", 600, -2., 2.);
 
   // Book histograms
   Map3DcalibNoCuts = new TH2F("3DcalibNoCuts", "3Dcalib (Before Cuts)", 85, 1, 85, 70, 5, 75);
   e9NoCuts = new TH1F("e9NoCuts", "E9 energy (Before Cuts)", 300, 0., 150.);
   e25NoCuts = new TH1F("e25NoCuts", "E25 energy (Before Cuts)", 300, 0., 150.);
   scENoCuts = new TH1F("scENoCuts", "SC energy (Before Cuts)", 300, 0., 150.);
   trPNoCuts = new TH1F("trPNoCuts", "Trk momentum (Before Cuts)", 300, 0., 150.);
   EoPNoCuts = new TH1F("EoPNoCuts", "EoP (Before Cuts)", 600, 0., 3.);
   if (elecclass_ == 0 || elecclass_ == -1) {
     calibsNoCuts = new TH1F("calibNoCuts", "Calibration constants (Before Cuts)", 4000, 0., 2.);
   } else {
     calibsNoCuts = new TH1F("calibNoCuts", "Calibration constants (Before Cuts)", 800, 0., 2.);
   }
   e25OverScENoCuts = new TH1F("e25OverscENoCuts", "E25 / SC energy (Before Cuts)", 400, 0., 2.);
   E25oPNoCuts = new TH1F("E25oPNoCuts", "E25 / P (Before Cuts)", 1200, 0., 1.5);
   MapNoCuts = new TH2F("MapNoCuts", "Nb Events in Crystal (Before Cuts)", 85, 1, 85, 70, 5, 75);
   e9Overe25NoCuts = new TH1F("e9Overe25NoCuts", "E9 / E25 (Before Cuts)", 400, 0., 2.);
   PinOverPoutNoCuts = new TH1F("PinOverPoutNoCuts", "pinOverpout (Before Cuts)", 600, 0., 3.);
   eSeedOverPoutNoCuts = new TH1F(" eSeedOverPoutNoCuts", "eSeedOverpout (Before Cuts) ", 600, 0., 4.);
   PinMinPoutNoCuts = new TH1F("PinMinPoutNoCuts", "(Pin - Pout)/Pin (Before Cuts)", 600, -2.0, 2.0);
 
   RatioCalibsNoCuts = new TH1F("RatioCalibsNoCuts", "Ratio in Calibration Constants (Before Cuts)", 4000, 0.5, 2.0);
   DiffCalibsNoCuts = new TH1F("DiffCalibsNoCuts", "Difference in Calibration constants (Before Cuts)", 4000, -1.0, 1.0);
   calibinterNoCuts = new TH1F("calibinterNoCuts", "internal calibration constants", 2000, 0.5, 2.);
 
   MapCor1NoCuts =
       new TH2F("MapCor1NoCuts", "Correlation E25/PatCalo versus E25/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor2NoCuts =
       new TH2F("MapCor2NoCuts", "Correlation E25/PatCalo versus E/P (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor3NoCuts =
       new TH2F("MapCor3NoCuts", "Correlation E25/PatCalo versus Pout/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor4NoCuts =
       new TH2F("MapCor4NoCuts", "Correlation E25/PatCalo versus E25/highestP (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor5NoCuts =
       new TH2F("MapCor5NoCuts", "Correlation E25/Pcalo versus Pcalo/Pout (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor6NoCuts =
       new TH2F("MapCor6NoCuts", "Correlation Pout/Pin versus E25/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor7NoCuts =
       new TH2F("MapCor7NoCuts", "Correlation Pout/Pin versus Pcalo/Pout (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor8NoCuts =
       new TH2F("MapCor8NoCuts", "Correlation E25/Pin versus Pcalo/Pout (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor9NoCuts =
       new TH2F("MapCor9NoCuts", "Correlation  E25/Pcalo versus Eseed/Pout (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor10NoCuts =
       new TH2F("MapCor10NoCuts", "Correlation Eseed/Pout versus Pout/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCor11NoCuts =
       new TH2F("MapCor11NoCuts", "Correlation Eseed/Pout versus E25/Pin (Before Cuts)", 100, 0., 5., 100, 0., 5.);
   MapCorCalibNoCuts = new TH2F("MapCorCalibNoCuts",
                                "Correlation Miscalibration versus Calibration constants (Before Cuts)",
                                100,
                                0.,
                                3.,
                                100,
                                0.,
                                3.);
 
   Error1NoCuts = new TH1F("Eror1NoCuts", "DeltaP/Pin (Before Cuts)", 800, -1.0, 1.0);
   Error2NoCuts = new TH1F("Error2NoCuts", "DeltaP/Pout (Before Cuts)", 800, -1.0, 1.0);
   Error3NoCuts = new TH1F("Error3NoCuts", "DeltaP/Pcalo (Before Cuts)", 800, -1.0, 1.0);
   eSeedOverPout2NoCuts = new TH1F("eSeedOverPout2NoCuts", "eSeedOverpout (No Supercluster, Before Cuts)", 600, 0., 4.);
   hadOverEmNoCuts = new TH1F("hadOverEmNoCuts", "Had/EM distribution (Before Cuts)", 600, -2., 2.);
 
   //Book histograms after ESeed cut
   MapCor1ESeed =
       new TH2F("MapCor1ESeed", "Correlation E25/Pcalo versus E25/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor2ESeed =
       new TH2F("MapCor2ESeed", "Correlation E25/Pcalo versus E/P (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor3ESeed = new TH2F(
       "MapCor3ESeed", "Correlation E25/Pcalo versus Pout/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor4ESeed = new TH2F(
       "MapCor4ESeed", "Correlation E25/Pcalo versus E25/highestP (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor5ESeed = new TH2F(
       "MapCor5ESeed", "Correlation E25/Pcalo versus Pcalo/Pout (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor6ESeed =
       new TH2F("MapCor6ESeed", "Correlation Pout/Pin versus E25/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor7ESeed = new TH2F(
       "MapCor7ESeed", "Correlation Pout/Pin versus Pcalo/Pout (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor8ESeed = new TH2F(
       "MapCor8ESeed", "Correlation E25/Pin versus Pcalo/Pout (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor9ESeed = new TH2F(
       "MapCor9ESeed", "Correlation  E25/Pcalo versus Eseed/Pout (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor10ESeed = new TH2F(
       "MapCor10ESeed", "Correlation Eseed/Pout versus Pout/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
   MapCor11ESeed = new TH2F(
       "MapCor11ESeed", "Correlation Eseed/Pout versus E25/Pin (after Eseed/Pout cut)", 100, 0., 5., 100, 0., 5.);
 
   eSeedOverPout2ESeed =
       new TH1F("eSeedOverPout2ESeed", "eSeedOverpout (No Supercluster, after Eseed/Pout cut)", 600, 0., 4.);
 
   hadOverEmESeed = new TH1F("hadOverEmESeed", "Had/EM distribution (after Eseed/Pout cut)", 600, -2., 2.);
 
   //Book histograms without any cut
   GeneralMap = new TH2F("GeneralMap", "Map without any cuts", 85, 1, 85, 70, 5, 75);
 
   calibClusterSize = ClusterSize_;
   etaMin = mineta_;
   etaMax = maxeta_;
   phiMin = minphi_;
   phiMax = maxphi_;
   if (calibAlgo_ == "L3") {
     MyL3Algo1 = new MinL3Algorithm(keventweight_, calibClusterSize, etaMin, etaMax, phiMin, phiMax);
   } else {
     if (calibAlgo_ == "HH" || calibAlgo_ == "HHReg") {
       MyHH = new HouseholderDecomposition(calibClusterSize, etaMin, etaMax, phiMin, phiMax);
     } else {
       edm::LogError("EcalCalib") << " Name of Algorithm is not recognize " << calibAlgo_
                                  << " Should be either L3, HH or HHReg. Abort! ";
     }
   }
   read_events = 0;
 
   // get Region to be calibrated
   ReducedMap = calibCluster.getMap(etaMin, etaMax, phiMin, phiMax);
 
   oldCalibs.resize(ReducedMap.size(), 0.);
 
   // table is set to zero
   for (int phi = 0; phi < 360; phi++) {
     for (int eta = 0; eta < 171; eta++) {
       eventcrystal[eta][phi] = 0;
     }
   }
 
   edm::LogVerbatim("EcalCalib") << " Begin JOB ";
 }

◆ endJob()

void ElectronCalibration::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 257 of file ElectronCalibration.cc.

References calibAlgo_, calibinter, calibinterNoCuts, calibs, calibsNoCuts, DiffCalibs, DiffCalibsNoCuts, EnergyVector, EnergyVectorNoCuts, PVValHelper::eta, EBDetId::ETAPHIMODE, EventMatrix, EventMatrixNoCuts, f, hcalRecHitTable_cff::ieta, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), cuy::ii, createfilelist::int, hcalRecHitTable_cff::iphi, MinL3Algorithm::iterate(), HouseholderDecomposition::iterate(), Map3Dcalib, Map3DcalibNoCuts, MapCorCalib, MapCorCalibNoCuts, MaxCCeta, MaxCCetaNoCuts, MaxCCphi, MaxCCphiNoCuts, maxeta_, maxphi_, mineta_, minphi_, miscalibfile_, MisCalibs, MyHH, MyL3Algo1, newCalibs, phi, RatioCalibs, RatioCalibsNoCuts, read_events, ReducedMap, HouseholderDecomposition::runRegional(), solution, solutionNoCuts, and calibXMLwriter::writeLine().

Referenced by o2olib.O2ORunMgr::executeJob().

                                  {
   //========================================================================
 
   int nIterations = 10;
   if (calibAlgo_ == "L3") {
     solution = MyL3Algo1->iterate(EventMatrix, MaxCCeta, MaxCCphi, EnergyVector, nIterations);
   } else {
     if (calibAlgo_ == "HH") {
       solution = MyHH->iterate(EventMatrix, MaxCCeta, MaxCCphi, EnergyVector, 1, false);
     } else {
       if (calibAlgo_ == "HHReg") {
         solution = MyHH->runRegional(EventMatrix, MaxCCeta, MaxCCphi, EnergyVector, 2);
       } else {
         edm::LogError("EcalCalib") << " Calibration not run due to problem in Algo Choice...";
         return;
       }
     }
   }
   for (int ii = 0; ii < (int)solution.size(); ii++) {
     edm::LogVerbatim("EcalCalib") << "solution[" << ii << "] = " << solution[ii];
     calibs->Fill(solution[ii]);
   }
 
   newCalibs.resize(ReducedMap.size(), 0.);
 
   calibXMLwriter write_calibrations;
 
   FILE* MisCalib;
   MisCalib = fopen(miscalibfile_.c_str(), "r");
   int fileStatus = 1;
   int eta = -1;
   int phi = -1;
   float coeff = -1;
 
   std::map<EBDetId, float> OldCoeff;
 
   while (fileStatus != EOF) {
     fileStatus = fscanf(MisCalib, "%d %d %f\n", &eta, &phi, &coeff);
     if (eta != -1 && phi != -1 && coeff != -1) {
       OldCoeff.insert(std::make_pair(EBDetId(eta, phi, EBDetId::ETAPHIMODE), coeff));
     }
   }
 
   fclose(MisCalib);
 
   int icry = 0;
   CalibrationCluster::CalibMap::iterator itmap;
   for (itmap = ReducedMap.begin(); itmap != ReducedMap.end(); itmap++) {
     newCalibs[icry] = solution[icry];
 
     write_calibrations.writeLine(itmap->first, newCalibs[icry]);
     float Compare = 1.;
     std::map<EBDetId, float>::iterator iter = OldCoeff.find(itmap->first);
     if (iter != OldCoeff.end())
       Compare = iter->second;
 
     if ((itmap->first).ieta() > mineta_ && (itmap->first).ieta() < maxeta_ && (itmap->first).iphi() > minphi_ &&
         (itmap->first).iphi() < maxphi_) {
       Map3Dcalib->Fill((itmap->first).ieta(), (itmap->first).iphi(), newCalibs[icry] * Compare);
       MisCalibs->Fill(Compare);
     }
     if ((itmap->first).ieta() < mineta_ + 2) {
       icry++;
       continue;
     }
     if ((itmap->first).ieta() > maxeta_ - 2) {
       icry++;
       continue;
     }
     if ((itmap->first).iphi() < minphi_ + 2) {
       icry++;
       continue;
     }
     if ((itmap->first).iphi() > maxphi_ - 2) {
       icry++;
       continue;
     }
 
     calibinter->Fill(newCalibs[icry]);
     DiffCalibs->Fill(newCalibs[icry] - 1. / Compare);
     RatioCalibs->Fill(newCalibs[icry] * Compare);
     MapCorCalib->Fill(1. / Compare, newCalibs[icry]);
     icry++;
   }
 
   if (calibAlgo_ == "L3") {
     solutionNoCuts =
         MyL3Algo1->iterate(EventMatrixNoCuts, MaxCCetaNoCuts, MaxCCphiNoCuts, EnergyVectorNoCuts, nIterations);
   } else {
     if (calibAlgo_ == "HH") {
       solutionNoCuts = MyHH->iterate(EventMatrixNoCuts, MaxCCetaNoCuts, MaxCCphiNoCuts, EnergyVectorNoCuts, 1, false);
     } else {
       if (calibAlgo_ == "HHReg") {
         solutionNoCuts = MyHH->runRegional(EventMatrixNoCuts, MaxCCetaNoCuts, MaxCCphiNoCuts, EnergyVectorNoCuts, 2);
       } else {
         edm::LogError("EcalCalib") << " Calibration not run due to problem in AlgoChoice...";
         return;
       }
     }
   }
   for (int ii = 0; ii < (int)solutionNoCuts.size(); ii++) {
     calibsNoCuts->Fill(solutionNoCuts[ii]);
   }
   int icryp = 0;
   CalibrationCluster::CalibMap::iterator itmapp;
   for (itmapp = ReducedMap.begin(); itmapp != ReducedMap.end(); itmapp++) {
     newCalibs[icryp] = solutionNoCuts[icryp];
     float Compare2 = 1.;
     std::map<EBDetId, float>::iterator iter2 = OldCoeff.find(itmapp->first);
     if (iter2 != OldCoeff.end())
       Compare2 = iter2->second;
 
     if ((itmapp->first).ieta() > mineta_ && (itmapp->first).ieta() < maxeta_ && (itmapp->first).iphi() > minphi_ &&
         (itmapp->first).iphi() < maxphi_)
       Map3DcalibNoCuts->Fill((itmapp->first).ieta(), (itmapp->first).iphi(), newCalibs[icryp] * Compare2);
     if ((itmapp->first).ieta() < mineta_ + 2) {
       icryp++;
       continue;
     }
     if ((itmapp->first).ieta() > maxeta_ - 2) {
       icryp++;
       continue;
     }
     if ((itmapp->first).iphi() < minphi_ + 2) {
       icryp++;
       continue;
     }
     if ((itmapp->first).iphi() > maxphi_ - 2) {
       icryp++;
       continue;
     }
     calibinterNoCuts->Fill(newCalibs[icryp]);
     DiffCalibsNoCuts->Fill(newCalibs[icryp] - 1. / (Compare2));
     RatioCalibsNoCuts->Fill(newCalibs[icryp] * Compare2);
     MapCorCalibNoCuts->Fill(1. / Compare2, newCalibs[icryp]);
     icryp++;
   }
 
 
   edm::LogVerbatim("EcalCalib") << "\n************* STATISTICS **************";
   edm::LogVerbatim("EcalCalib") << " Events Studied " << read_events;
 
 
   f->Write();
 
   f->Close();
 }

◆ findMaxHit()

EBDetId ElectronCalibration::findMaxHit ( edm::Handle< EBRecHitCollection > & phits )

private

Definition at line 408 of file ElectronCalibration.cc.

References edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, AlCaHLTBitMon_ParallelJobs::p, phits, and cuy::save.

                                                                             {
   //=================================================================================
 
   EcalRecHitCollection ecrh = *phits;
   EcalRecHitCollection::iterator it;
   EBDetId save;
   float en_save = 0;
   for (it = ecrh.begin(); it != ecrh.end(); it++) {
     EBDetId p = EBDetId(it->id().rawId());
     if (it->energy() > en_save) {
       en_save = it->energy();
       save = p;
     }
   }
   return save;
 }

◆ findMaxHit2()

EBDetId ElectronCalibration::findMaxHit2	(	const std::vector< DetId > &	v1,
		const EBRecHitCollection *	hits
	)

private

Definition at line 426 of file ElectronCalibration.cc.

References hfClusterShapes_cfi::hits, and MTVHistoProducerAlgoForTrackerBlock_cfi::maxHit.

Referenced by analyze().

                                                                                                    {
   //=================================================================================
 
   double currEnergy = 0.;
   EBDetId maxHit;
 
   for (std::vector<DetId>::const_iterator idsIt = v1.begin(); idsIt != v1.end(); ++idsIt) {
     if (idsIt->subdetId() != 1)
       continue;
     EBRecHitCollection::const_iterator itrechit;
     itrechit = hits->find(*idsIt);
 
     if (itrechit == hits->end()) {
       edm::LogVerbatim("EcalCalib") << "ElectronCalibration::findMaxHit2: rechit not found! ";
       continue;
     }
     if (itrechit->energy() > currEnergy) {
       currEnergy = itrechit->energy();
       maxHit = *idsIt;
     }
   }
 
   return maxHit;
 }

Member Data Documentation

◆ calibAlgo_

std::string ElectronCalibration::calibAlgo_

private

Definition at line 64 of file ElectronCalibration.h.

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

◆ calibCluster

CalibrationCluster ElectronCalibration::calibCluster

private

Definition at line 85 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ calibClusterSize

int ElectronCalibration::calibClusterSize

private

Definition at line 90 of file ElectronCalibration.h.

Referenced by beginJob().

◆ calibinter

TH1F* ElectronCalibration::calibinter

private

Definition at line 145 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ calibinterNoCuts

TH1F* ElectronCalibration::calibinterNoCuts

private

Definition at line 146 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ calibs

TH1F* ElectronCalibration::calibs

private

Definition at line 119 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ calibsNoCuts

TH1F* ElectronCalibration::calibsNoCuts

private

Definition at line 138 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ ClusterSize_

int ElectronCalibration::ClusterSize_

private

Definition at line 83 of file ElectronCalibration.h.

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

◆ cut1_

double ElectronCalibration::cut1_

private

Definition at line 72 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ cut2_

double ElectronCalibration::cut2_

private

Definition at line 73 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ cut3_

double ElectronCalibration::cut3_

private

Definition at line 74 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ cutCalo1_

double ElectronCalibration::cutCalo1_

private

Definition at line 80 of file ElectronCalibration.h.

Referenced by ElectronCalibration().

◆ cutCalo2_

double ElectronCalibration::cutCalo2_

private

Definition at line 81 of file ElectronCalibration.h.

Referenced by ElectronCalibration().

◆ cutEPCalo1_

double ElectronCalibration::cutEPCalo1_

private

Definition at line 76 of file ElectronCalibration.h.

Referenced by ElectronCalibration().

◆ cutEPCalo2_

double ElectronCalibration::cutEPCalo2_

private

Definition at line 77 of file ElectronCalibration.h.

Referenced by ElectronCalibration().

◆ cutEPin1_

double ElectronCalibration::cutEPin1_

private

Definition at line 78 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ cutEPin2_

double ElectronCalibration::cutEPin2_

private

Definition at line 79 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ cutESeed_

double ElectronCalibration::cutESeed_

private

Definition at line 82 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ DiffCalibs

TH1F* ElectronCalibration::DiffCalibs

private

Definition at line 129 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ DiffCalibsNoCuts

TH1F* ElectronCalibration::DiffCalibsNoCuts

private

Definition at line 131 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ e25

TH1F* ElectronCalibration::e25

private

Definition at line 113 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ e25NoCuts

TH1F* ElectronCalibration::e25NoCuts

private

Definition at line 133 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ E25oP

TH1F* ElectronCalibration::E25oP

private

Definition at line 123 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ E25oPNoCuts

TH1F* ElectronCalibration::E25oPNoCuts

private

Definition at line 142 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ e25OverScE

TH1F* ElectronCalibration::e25OverScE

private

Definition at line 121 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ e25OverScENoCuts

TH1F* ElectronCalibration::e25OverScENoCuts

private

Definition at line 140 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ e9

TH1F* ElectronCalibration::e9

private

Definition at line 114 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ e9NoCuts

TH1F* ElectronCalibration::e9NoCuts

private

Definition at line 134 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ e9Overe25

TH1F* ElectronCalibration::e9Overe25

private

Definition at line 120 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ e9Overe25NoCuts

TH1F* ElectronCalibration::e9Overe25NoCuts

private

Definition at line 139 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ elecclass_

int ElectronCalibration::elecclass_

private

Definition at line 75 of file ElectronCalibration.h.

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

◆ electronToken_

edm::EDGetTokenT<reco::GsfElectronCollection> ElectronCalibration::electronToken_

private

Definition at line 62 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ ElePt_

double ElectronCalibration::ElePt_

private

Definition at line 67 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ EnergyVector

std::vector<float> ElectronCalibration::EnergyVector

private

Definition at line 92 of file ElectronCalibration.h.

Referenced by analyze(), and endJob().

◆ EnergyVectorNoCuts

std::vector<float> ElectronCalibration::EnergyVectorNoCuts

private

Definition at line 97 of file ElectronCalibration.h.

Referenced by analyze(), and endJob().

◆ EoP

TH1F* ElectronCalibration::EoP

private

Definition at line 117 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ EoP_all

TH1F* ElectronCalibration::EoP_all

private

Definition at line 118 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ EoPNoCuts

TH1F* ElectronCalibration::EoPNoCuts

private

Definition at line 137 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ Error1

TH1F* ElectronCalibration::Error1

private

Definition at line 194 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ Error1NoCuts

TH1F* ElectronCalibration::Error1NoCuts

private

Definition at line 197 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ Error2

TH1F* ElectronCalibration::Error2

private

Definition at line 195 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ Error2NoCuts

TH1F* ElectronCalibration::Error2NoCuts

private

Definition at line 198 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ Error3

TH1F* ElectronCalibration::Error3

private

Definition at line 196 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ Error3NoCuts

TH1F* ElectronCalibration::Error3NoCuts

private

Definition at line 199 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ eSeedOverPout

TH1F* ElectronCalibration::eSeedOverPout

private

Definition at line 126 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ eSeedOverPout2

TH1F* ElectronCalibration::eSeedOverPout2

private

Definition at line 201 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ eSeedOverPout2ESeed

TH1F* ElectronCalibration::eSeedOverPout2ESeed

private

Definition at line 203 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ eSeedOverPout2NoCuts

TH1F* ElectronCalibration::eSeedOverPout2NoCuts

private

Definition at line 202 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ eSeedOverPoutNoCuts

TH1F* ElectronCalibration::eSeedOverPoutNoCuts

private

Definition at line 149 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ etaMax

int ElectronCalibration::etaMax

private

Definition at line 91 of file ElectronCalibration.h.

Referenced by beginJob().

◆ etaMin

int ElectronCalibration::etaMin

private

Definition at line 91 of file ElectronCalibration.h.

Referenced by beginJob().

◆ eventcrystal

int ElectronCalibration::eventcrystal[171][360]

private

Definition at line 108 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ EventMatrix

std::vector<std::vector<float> > ElectronCalibration::EventMatrix

private

Definition at line 94 of file ElectronCalibration.h.

Referenced by analyze(), and endJob().

◆ EventMatrixNoCuts

std::vector<std::vector<float> > ElectronCalibration::EventMatrixNoCuts

private

Definition at line 99 of file ElectronCalibration.h.

Referenced by analyze(), and endJob().

◆ f

TFile* ElectronCalibration::f

private

Definition at line 111 of file ElectronCalibration.h.

Referenced by beginJob(), endJob(), ztail.Decoder::follow(), ztail.Decoder::initial_synchronize(), DQMIO2histo.DQMIO::print_index(), electronCompare.flushfile::write(), electronStore.flushfile::write(), and DQMIO2histo.DQMIO::write_to_file().

◆ GeneralMap

TH2F* ElectronCalibration::GeneralMap

private

Definition at line 151 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ hadOverEm

TH1F* ElectronCalibration::hadOverEm

private

Definition at line 205 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ hadOverEmESeed

TH1F* ElectronCalibration::hadOverEmESeed

private

Definition at line 207 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ hadOverEmNoCuts

TH1F* ElectronCalibration::hadOverEmNoCuts

private

Definition at line 206 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ keventweight

int ElectronCalibration::keventweight

private

Definition at line 90 of file ElectronCalibration.h.

◆ keventweight_

int ElectronCalibration::keventweight_

private

Definition at line 66 of file ElectronCalibration.h.

Referenced by beginJob(), and ElectronCalibration().

◆ Map

TH2F* ElectronCalibration::Map

private

Definition at line 122 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ Map3Dcalib

TH2F* ElectronCalibration::Map3Dcalib

private

Definition at line 143 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ Map3DcalibNoCuts

TH2F* ElectronCalibration::Map3DcalibNoCuts

private

Definition at line 144 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ MapCor1

TH2F* ElectronCalibration::MapCor1

private

Definition at line 153 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor10

TH2F* ElectronCalibration::MapCor10

private

Definition at line 162 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor10ESeed

TH2F* ElectronCalibration::MapCor10ESeed

private

Definition at line 188 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor10NoCuts

TH2F* ElectronCalibration::MapCor10NoCuts

private

Definition at line 175 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor11

TH2F* ElectronCalibration::MapCor11

private

Definition at line 163 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor11ESeed

TH2F* ElectronCalibration::MapCor11ESeed

private

Definition at line 189 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor11NoCuts

TH2F* ElectronCalibration::MapCor11NoCuts

private

Definition at line 176 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor1ESeed

TH2F* ElectronCalibration::MapCor1ESeed

private

Definition at line 179 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor1NoCuts

TH2F* ElectronCalibration::MapCor1NoCuts

private

Definition at line 166 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor2

TH2F* ElectronCalibration::MapCor2

private

Definition at line 154 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor2ESeed

TH2F* ElectronCalibration::MapCor2ESeed

private

Definition at line 180 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor2NoCuts

TH2F* ElectronCalibration::MapCor2NoCuts

private

Definition at line 167 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor3

TH2F* ElectronCalibration::MapCor3

private

Definition at line 155 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor3ESeed

TH2F* ElectronCalibration::MapCor3ESeed

private

Definition at line 181 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor3NoCuts

TH2F* ElectronCalibration::MapCor3NoCuts

private

Definition at line 168 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor4

TH2F* ElectronCalibration::MapCor4

private

Definition at line 156 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor4ESeed

TH2F* ElectronCalibration::MapCor4ESeed

private

Definition at line 182 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor4NoCuts

TH2F* ElectronCalibration::MapCor4NoCuts

private

Definition at line 169 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor5

TH2F* ElectronCalibration::MapCor5

private

Definition at line 157 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor5ESeed

TH2F* ElectronCalibration::MapCor5ESeed

private

Definition at line 183 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor5NoCuts

TH2F* ElectronCalibration::MapCor5NoCuts

private

Definition at line 170 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor6

TH2F* ElectronCalibration::MapCor6

private

Definition at line 158 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor6ESeed

TH2F* ElectronCalibration::MapCor6ESeed

private

Definition at line 184 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor6NoCuts

TH2F* ElectronCalibration::MapCor6NoCuts

private

Definition at line 171 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor7

TH2F* ElectronCalibration::MapCor7

private

Definition at line 159 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor7ESeed

TH2F* ElectronCalibration::MapCor7ESeed

private

Definition at line 185 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor7NoCuts

TH2F* ElectronCalibration::MapCor7NoCuts

private

Definition at line 172 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor8

TH2F* ElectronCalibration::MapCor8

private

Definition at line 160 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor8ESeed

TH2F* ElectronCalibration::MapCor8ESeed

private

Definition at line 186 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor8NoCuts

TH2F* ElectronCalibration::MapCor8NoCuts

private

Definition at line 173 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor9

TH2F* ElectronCalibration::MapCor9

private

Definition at line 161 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor9ESeed

TH2F* ElectronCalibration::MapCor9ESeed

private

Definition at line 187 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCor9NoCuts

TH2F* ElectronCalibration::MapCor9NoCuts

private

Definition at line 174 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MapCorCalib

TH2F* ElectronCalibration::MapCorCalib

private

Definition at line 164 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ MapCorCalibNoCuts

TH2F* ElectronCalibration::MapCorCalibNoCuts

private

Definition at line 177 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ MapNoCuts

TH2F* ElectronCalibration::MapNoCuts

private

Definition at line 141 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ MaxCCeta

std::vector<int> ElectronCalibration::MaxCCeta

private

Definition at line 95 of file ElectronCalibration.h.

Referenced by analyze(), and endJob().

◆ MaxCCetaNoCuts

std::vector<int> ElectronCalibration::MaxCCetaNoCuts

private

Definition at line 100 of file ElectronCalibration.h.

Referenced by analyze(), and endJob().

◆ MaxCCphi

std::vector<int> ElectronCalibration::MaxCCphi

private

Definition at line 96 of file ElectronCalibration.h.

Referenced by analyze(), and endJob().

◆ MaxCCphiNoCuts

std::vector<int> ElectronCalibration::MaxCCphiNoCuts

private

Definition at line 101 of file ElectronCalibration.h.

Referenced by analyze(), and endJob().

◆ maxeta_

int ElectronCalibration::maxeta_

private

Definition at line 68 of file ElectronCalibration.h.

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

◆ maxphi_

int ElectronCalibration::maxphi_

private

Definition at line 70 of file ElectronCalibration.h.

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

◆ mineta_

int ElectronCalibration::mineta_

private

Definition at line 69 of file ElectronCalibration.h.

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

◆ minphi_

int ElectronCalibration::minphi_

private

Definition at line 71 of file ElectronCalibration.h.

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

◆ miscalibfile_

std::string ElectronCalibration::miscalibfile_

private

Definition at line 65 of file ElectronCalibration.h.

Referenced by ElectronCalibration(), and endJob().

◆ MisCalibs

TH1F* ElectronCalibration::MisCalibs

private

Definition at line 127 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ MyHH

HouseholderDecomposition* ElectronCalibration::MyHH

private

Definition at line 147 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ MyL3Algo1

MinL3Algorithm* ElectronCalibration::MyL3Algo1

private

Definition at line 102 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ newCalibs

std::vector<float> ElectronCalibration::newCalibs

private

Definition at line 105 of file ElectronCalibration.h.

Referenced by endJob().

◆ numevent_

int ElectronCalibration::numevent_

private

Definition at line 109 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ oldCalibs

std::vector<float> ElectronCalibration::oldCalibs

private

Definition at line 106 of file ElectronCalibration.h.

Referenced by beginJob().

◆ phiMax

int ElectronCalibration::phiMax

private

Definition at line 91 of file ElectronCalibration.h.

Referenced by beginJob().

◆ phiMin

int ElectronCalibration::phiMin

private

Definition at line 91 of file ElectronCalibration.h.

Referenced by beginJob().

◆ PinMinPout

TH1F* ElectronCalibration::PinMinPout

private

Definition at line 191 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ PinMinPoutNoCuts

TH1F* ElectronCalibration::PinMinPoutNoCuts

private

Definition at line 192 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ PinOverPout

TH1F* ElectronCalibration::PinOverPout

private

Definition at line 125 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ PinOverPoutNoCuts

TH1F* ElectronCalibration::PinOverPoutNoCuts

private

Definition at line 148 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ RatioCalibs

TH1F* ElectronCalibration::RatioCalibs

private

Definition at line 128 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ RatioCalibsNoCuts

TH1F* ElectronCalibration::RatioCalibsNoCuts

private

Definition at line 130 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ read_events

int ElectronCalibration::read_events

private

Definition at line 88 of file ElectronCalibration.h.

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

◆ recHitToken_

edm::EDGetTokenT<EBRecHitCollection> ElectronCalibration::recHitToken_

private

Definition at line 61 of file ElectronCalibration.h.

Referenced by analyze(), and ElectronCalibration().

◆ ReducedMap

CalibrationCluster::CalibMap ElectronCalibration::ReducedMap

private

Definition at line 86 of file ElectronCalibration.h.

Referenced by beginJob(), and endJob().

◆ rootfile_

std::string ElectronCalibration::rootfile_

private

Definition at line 60 of file ElectronCalibration.h.

Referenced by beginJob(), and ElectronCalibration().

◆ scE

TH1F* ElectronCalibration::scE

private

Definition at line 115 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ scENoCuts

TH1F* ElectronCalibration::scENoCuts

private

Definition at line 135 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ solution

std::vector<float> ElectronCalibration::solution

private

Definition at line 103 of file ElectronCalibration.h.

Referenced by endJob().

◆ solutionNoCuts

std::vector<float> ElectronCalibration::solutionNoCuts

private

Definition at line 104 of file ElectronCalibration.h.

Referenced by endJob().

◆ trackLabel_

edm::InputTag ElectronCalibration::trackLabel_

private

Definition at line 63 of file ElectronCalibration.h.

Referenced by ElectronCalibration().

◆ trP

TH1F* ElectronCalibration::trP

private

Definition at line 116 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ trPNoCuts

TH1F* ElectronCalibration::trPNoCuts

private

Definition at line 136 of file ElectronCalibration.h.

Referenced by analyze(), and beginJob().

◆ WeightVector

std::vector<float> ElectronCalibration::WeightVector

private

Definition at line 93 of file ElectronCalibration.h.

Referenced by analyze().

◆ WeightVectorNoCuts

std::vector<float> ElectronCalibration::WeightVectorNoCuts

private

Definition at line 98 of file ElectronCalibration.h.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ElectronCalibration()

◆ ~ElectronCalibration()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ findMaxHit()

◆ findMaxHit2()

Member Data Documentation

◆ calibAlgo_

◆ calibCluster

◆ calibClusterSize

◆ calibinter

◆ calibinterNoCuts

◆ calibs

◆ calibsNoCuts

◆ ClusterSize_

◆ cut1_

◆ cut2_

◆ cut3_

◆ cutCalo1_

◆ cutCalo2_

◆ cutEPCalo1_

◆ cutEPCalo2_

◆ cutEPin1_

◆ cutEPin2_

◆ cutESeed_

◆ DiffCalibs

◆ DiffCalibsNoCuts

◆ e25

◆ e25NoCuts

◆ E25oP

◆ E25oPNoCuts

◆ e25OverScE

◆ e25OverScENoCuts

◆ e9

◆ e9NoCuts

◆ e9Overe25

◆ e9Overe25NoCuts

◆ elecclass_

◆ electronToken_

◆ ElePt_

◆ EnergyVector

◆ EnergyVectorNoCuts

◆ EoP

◆ EoP_all

◆ EoPNoCuts

◆ Error1

◆ Error1NoCuts

◆ Error2

◆ Error2NoCuts

◆ Error3

◆ Error3NoCuts

◆ eSeedOverPout

◆ eSeedOverPout2

◆ eSeedOverPout2ESeed

◆ eSeedOverPout2NoCuts

◆ eSeedOverPoutNoCuts

◆ etaMax

◆ etaMin

◆ eventcrystal

◆ EventMatrix

◆ EventMatrixNoCuts

◆ f

◆ GeneralMap

◆ hadOverEm

◆ hadOverEmESeed

◆ hadOverEmNoCuts

◆ keventweight

◆ keventweight_

◆ Map

◆ Map3Dcalib

◆ Map3DcalibNoCuts

◆ MapCor1