ECAL TB 2006 calibration with matrix inversion technique. More...

#include <InvRingCalib.h>

Inheritance diagram for InvRingCalib:

Public Member Functions
void	beginOfJob () override
	BeginOfJob. More...

Status	duringLoop (const edm::Event &, const edm::EventSetup &) override
	duringLoop More...

void	endOfJob () override

Status	endOfLoop (const edm::EventSetup &, unsigned int iCounter) override

	InvRingCalib (const edm::ParameterSet &)
	ctor More...

void	startingNewLoop (unsigned int) override
	startingNewLoop More...

	~InvRingCalib () override
	dtor More...

Public Member Functions inherited from edm::EDLooper
	EDLooper ()

	EDLooper (EDLooper const &)=delete

EDLooper &	operator= (EDLooper const &)=delete

	~EDLooper () override

Public Member Functions inherited from edm::EDLooperBase
virtual void	attachTo (ActivityRegistry &)
	Override this method if you need to monitor the state of the processing. More...

void	beginOfJob (EventSetupImpl const &)

virtual void	beginOfJob (EventSetup const &)

void	copyInfo (ScheduleInfo const &)

void	doBeginLuminosityBlock (LuminosityBlockPrincipal &, EventSetupImpl const &, ProcessContext *)

void	doBeginRun (RunPrincipal &, EventSetupImpl const &, ProcessContext *)

Status	doDuringLoop (EventPrincipal &eventPrincipal, EventSetupImpl const &es, ProcessingController &, StreamContext *)

void	doEndLuminosityBlock (LuminosityBlockPrincipal &, EventSetupImpl const &, ProcessContext *)

Status	doEndOfLoop (EventSetupImpl const &es)

void	doEndRun (RunPrincipal &, EventSetupImpl const &, ProcessContext *)

void	doStartingNewLoop ()

	EDLooperBase ()

	EDLooperBase (EDLooperBase const &)=delete

void	esPrefetchAsync (WaitingTaskHolder iTask, EventSetupImpl const &iImpl, Transition iTrans, ServiceToken const &iToken) const

virtual std::set< eventsetup::EventSetupRecordKey >	modifyingRecords () const

EDLooperBase &	operator= (EDLooperBase const &)=delete

void	prefetchAsync (WaitingTaskHolder iTask, ServiceToken const &token, Transition iTrans, Principal const &iPrincipal, EventSetupImpl const &iImpl) const

void	prepareForNextLoop (eventsetup::EventSetupProvider *esp)

void	setActionTable (ExceptionToActionTable const *actionTable)

void	setModuleChanger (ModuleChanger *)

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

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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
int	EBRegId (const int)
	Gives back in which region you are: More...

void	EBRegionDef ()
	Defines the regions in the barrel. More...

int	EBRegionNum () const
	Number of regions in EB. More...

int	EERegId (int)
	gives back in which region of the endcap you are. More...

void	EERegionDef ()
	Defines the regions in the endcap. More...

int	EERegionNum () const
	The number of regions in EE. More...

void	EERingDef (const edm::EventSetup &)
	Defines the rins in the endcap. More...

void	RegPrepare ()
	Prepares the EB regions;. More...

Private Attributes
bool	isfirstcall_

const edm::InputTag	m_barrelAlCa
	EcalBarrel Input Collection name. More...

std::vector< DetId >	m_barrelCells
	geometry things used all over the file More...

EcalIntercalibConstantMap	m_barrelMap

std::map< int, int >	m_cellPhi

std::map< int, GlobalPoint >	m_cellPos
	position of the cell, borders, coords etc... More...

const std::string	m_EBcoeffFile
	coeffs filenames More...

const edm::EDGetTokenT< EBRecHitCollection >	m_ebRecHitToken
	ED token. More...

const std::string	m_EEcoeffFile

const edm::EDGetTokenT< EERecHitCollection >	m_eeRecHitToken

const int	m_EEZone
	endcap zone to be calibrated More...

const edm::InputTag	m_ElectronLabel
	To take the electrons. More...

const edm::InputTag	m_endcapAlCa
	EcalEndcap Input Collection name. More...

std::vector< DetId >	m_endcapCells

EcalIntercalibConstantMap	m_endcapMap

const int	m_endRing

const int	m_etaEnd
	eta end of the zone of interest More...

const int	m_etaStart
	eta start of the zone of interest More...

const int	m_etaWidth
	eta size of the regions More...

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	m_geometryToken
	ES token. More...

const edm::EDGetTokenT< reco::GsfElectronCollection >	m_gsfElectronToken

std::vector< VEcalCalibBlock * >	m_IMACalibBlocks
	single blocks calibrators More...

unsigned int	m_loops
	LP sets the number of loops to do. More...

VFillMap *	m_MapFiller
	The class that fills the map! More...

std::string	m_mapFillerType

const double	m_maxCoeff
	maximum coefficient accepted (RAW) More...

const double	m_maxEnergyPerCrystal
	maximum energy per crystal cut More...

const int	m_maxSelectedNumPerRing
	maximum number of events per Ring More...

const double	m_minCoeff
	minimum coefficient accepted (RAW) More...

const double	m_minEnergyPerCrystal
	minimum energy per crystal cut More...

const int	m_recoWindowSidex
	reconstruction window size More...

const int	m_recoWindowSidey

std::map< int, int >	m_Reg
	EB regions vs. eta index. More...

std::map< int, int >	m_RinginRegion
	association map between raw detIds and the number of the ring inside the region More...

std::map< int, int >	m_RingNumOfHits
	number of events already read per Ring More...

const int	m_startRing
	LP define the EE region to calibrate. More...

const int	m_usingBlockSolver
	to exclude the blocksolver More...

std::map< int, int >	m_xtalRegionId
	association map between raw detIds and Region More...

std::map< int, int >	m_xtalRing
	association map between Raw detIds and Rings More...

Additional Inherited Members
Public Types inherited from edm::EDLooperBase
enum	Status { kContinue, kStop }

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDLooperBase
ModuleChanger *	moduleChanger ()
	This only returns a non-zero value during the call to endOfLoop. More...

ScheduleInfo const *	scheduleInfo () const
	This returns a non-zero value after the constructor has been called. More...

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

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

ECAL TB 2006 calibration with matrix inversion technique.

Author

Definition at line 30 of file InvRingCalib.h.

Constructor & Destructor Documentation

◆ InvRingCalib()

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

explicit

ctor

Definition at line 26 of file InvRingCalib.cc.

References funct::abs(), qcdUeDQM_cfi::algorithm, cms::cuda::assert(), EBRegionNum(), EERegionNum(), beamvalidation::exit(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, createfilelist::int, m_EEZone, m_endRing, m_etaEnd, m_etaStart, m_etaWidth, m_IMACalibBlocks, m_loops, m_mapFillerType, m_startRing, and AlCaHLTBitMon_QueryRunRegistry::string.

     : m_barrelAlCa(iConfig.getParameter<edm::InputTag>("barrelAlca")),
       m_endcapAlCa(iConfig.getParameter<edm::InputTag>("endcapAlca")),
       m_ElectronLabel(iConfig.getParameter<edm::InputTag>("ElectronLabel")),
       m_recoWindowSidex(iConfig.getParameter<int>("recoWindowSidex")),
       m_recoWindowSidey(iConfig.getParameter<int>("recoWindowSidey")),
       m_minEnergyPerCrystal(iConfig.getParameter<double>("minEnergyPerCrystal")),
       m_maxEnergyPerCrystal(iConfig.getParameter<double>("maxEnergyPerCrystal")),
       m_etaStart(iConfig.getParameter<int>("etaStart")),
       m_etaEnd(iConfig.getParameter<int>("etaEnd")),
       m_etaWidth(iConfig.getParameter<int>("etaWidth")),
       m_maxSelectedNumPerRing(iConfig.getParameter<int>("maxNumPerRing")),
       m_minCoeff(iConfig.getParameter<double>("minCoeff")),
       m_maxCoeff(iConfig.getParameter<double>("maxCoeff")),
       m_usingBlockSolver(iConfig.getParameter<int>("usingBlockSolver")),
       m_startRing(iConfig.getParameter<int>("startRing")),
       m_endRing(iConfig.getParameter<int>("endRing")),
       m_EBcoeffFile(iConfig.getParameter<std::string>("EBcoeffs")),
       m_EEcoeffFile(iConfig.getParameter<std::string>("EEcoeffs")),
       m_EEZone(iConfig.getParameter<int>("EEZone")),
       m_ebRecHitToken(consumes<EBRecHitCollection>(m_barrelAlCa)),
       m_eeRecHitToken(consumes<EERecHitCollection>(m_endcapAlCa)),
       m_gsfElectronToken(consumes<reco::GsfElectronCollection>(m_ElectronLabel)),
       m_geometryToken(esConsumes()) {
   //controls if the parameters inputed are correct
   if ((m_etaEnd * m_etaStart) > 0)
     assert(!((m_etaEnd - m_etaStart) % m_etaWidth));
   if ((m_etaEnd * m_etaStart) < 0)
     assert(!((m_etaEnd - m_etaStart - 1) % m_etaWidth));
 
   assert(m_etaStart >= -85 && m_etaStart <= 86);
   assert(m_etaEnd >= m_etaStart && m_etaEnd <= 86);
   assert(m_startRing > -1 && m_startRing <= 40);
   assert(m_endRing >= m_startRing && m_endRing <= 40);
 
   assert(!((m_endRing - m_startRing) % m_etaWidth));
   assert((abs(m_EEZone) <= 1));
 
   m_loops = (unsigned int)iConfig.getParameter<int>("loops") - 1;
   //LP CalibBlock vector instantiation
   edm::LogInfo("IML") << "[InvRingCalib][ctor] Calib Block";
   std::string algorithm = iConfig.getParameter<std::string>("algorithm");
   m_mapFillerType = iConfig.getParameter<std::string>("FillType");
   int eventWeight = iConfig.getUntrackedParameter<int>("L3EventWeight", 1);
 
   for (int i = 0; i < EBRegionNum(); ++i) {
     if (algorithm == "IMA")
       m_IMACalibBlocks.push_back(new IMACalibBlock(m_etaWidth));
     else if (algorithm == "L3")
       m_IMACalibBlocks.push_back(new L3CalibBlock(m_etaWidth, eventWeight));
     else {
       edm::LogError("building") << algorithm << " is not a valid calibration algorithm";
       exit(1);
     }
   }
   int EEBlocks = 0;
   if (m_EEZone == 0)
     EEBlocks = 2 * EERegionNum();
   if (m_EEZone == 1 || m_EEZone == -1)
     EEBlocks = EERegionNum();
 
   for (int i = 0; i < EEBlocks; ++i) {
     if (algorithm == "IMA")
       m_IMACalibBlocks.push_back(new IMACalibBlock(m_etaWidth));
     else if (algorithm == "L3")
       m_IMACalibBlocks.push_back(new L3CalibBlock(m_etaWidth, eventWeight));
     else {
       edm::LogError("building") << algorithm << " is not a valid calibration algorithm";
       exit(1);
     }
   }
   edm::LogInfo("IML") << " [InvRingCalib][ctor] end of creator";
 }

◆ ~InvRingCalib()

InvRingCalib::~InvRingCalib ( )

override

dtor

destructor

Definition at line 103 of file InvRingCalib.cc.

103 {}

Member Function Documentation

◆ beginOfJob()

void InvRingCalib::beginOfJob ( )

overridevirtual

BeginOfJob.

Reimplemented from edm::EDLooperBase.

Definition at line 108 of file InvRingCalib.cc.

References isfirstcall_.

108 { isfirstcall_ = true; }

InvRingCalib::isfirstcall_

bool isfirstcall_

Definition: InvRingCalib.h:128

◆ duringLoop()

edm::EDLooper::Status InvRingCalib::duringLoop	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

duringLoop

Implements edm::EDLooper.

Definition at line 129 of file InvRingCalib.cc.

References funct::abs(), EBRegionDef(), DetId::Ecal, EcalBarrel, EcalEndcap, EERegionDef(), EERingDef(), HcalObjRepresent::Fill(), VFillMap::fillMap(), dqmdumpme::first, CaloMiscalibMapEcal::get(), edm::EventSetup::getData(), EcalCondObjectContainer< T >::getMap(), SiStripPI::getMaximum(), EBDetId::ieta(), iEvent, EBDetId::iphi(), isfirstcall_, edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), edm::EDLooperBase::kContinue, m_barrelCells, m_barrelMap, m_EBcoeffFile, m_ebRecHitToken, m_EEcoeffFile, m_eeRecHitToken, m_endcapCells, m_endcapMap, m_geometryToken, m_gsfElectronToken, m_IMACalibBlocks, m_MapFiller, m_mapFillerType, m_maxEnergyPerCrystal, m_minEnergyPerCrystal, m_recoWindowSidex, m_recoWindowSidey, m_RinginRegion, m_RingNumOfHits, m_xtalRegionId, m_xtalRing, METSkim_cff::Max, MillePedeFileConverter_cfg::out, MiscalibReaderFromXML::parseXMLMiscalibFile(), CaloMiscalibMapEcal::prefillMap(), edm::Handle< T >::product(), DetId::rawId(), and EEDetId::zside().

                                                                                               {
   if (isfirstcall_) {
     edm::LogInfo("IML") << "[InvRingCalib][beginOfJob]";
     //gets the geometry from the event setup
     const auto& geometry = iSetup.getData(m_geometryToken);
     edm::LogInfo("IML") << "[InvRingCalib] Event Setup read";
     //fills a vector with all the cells
     m_barrelCells = geometry.getValidDetIds(DetId::Ecal, EcalBarrel);
     m_endcapCells = geometry.getValidDetIds(DetId::Ecal, EcalEndcap);
     //Defines the EB regions
     edm::LogInfo("IML") << "[InvRingCalib] Defining Barrel Regions";
     EBRegionDef();
     //Defines what is a ring in the EE
     edm::LogInfo("IML") << "[InvRingCalib] Defining endcap Rings";
     EERingDef(iSetup);
     //Defines the regions in the EE
     edm::LogInfo("IML") << "[InvRingCalib] Defining endcap Regions";
     EERegionDef();
     if (m_mapFillerType == "Cluster")
       m_MapFiller = new ClusterFillMap(m_recoWindowSidex,
                                        m_recoWindowSidey,
                                        m_xtalRegionId,
                                        m_minEnergyPerCrystal,
                                        m_maxEnergyPerCrystal,
                                        m_RinginRegion,
                                        &m_barrelMap,
                                        &m_endcapMap);
     if (m_mapFillerType == "Matrix")
       m_MapFiller = new MatrixFillMap(m_recoWindowSidex,
                                       m_recoWindowSidey,
                                       m_xtalRegionId,
                                       m_minEnergyPerCrystal,
                                       m_maxEnergyPerCrystal,
                                       m_RinginRegion,
                                       &m_barrelMap,
                                       &m_endcapMap);
     edm::LogInfo("IML") << "[InvRingCalib] Initializing the coeffs";
     //Sets the initial coefficients to 1.
     //Graphs to check ring, regions and so on, not needed in the final version
     TH2F EBRegion("EBRegion", "EBRegion", 171, -85, 86, 360, 1, 361);
     TH2F EBRing("EBRing", "EBRing", 171, -85, 86, 360, 1, 361);
     for (std::vector<DetId>::const_iterator it = m_barrelCells.begin(); it != m_barrelCells.end(); ++it) {
       EBDetId eb(*it);
       EBRing.Fill(eb.ieta(), eb.iphi(), m_RinginRegion[it->rawId()]);
       EBRegion.Fill(eb.ieta(), eb.iphi(), m_xtalRegionId[it->rawId()]);
     }
 
     TH2F EEPRegion("EEPRegion", "EEPRegion", 100, 1, 101, 100, 1, 101);
     TH2F EEPRing("EEPRing", "EEPRing", 100, 1, 101, 100, 1, 101);
     TH2F EEPRingReg("EEPRingReg", "EEPRingReg", 100, 1, 101, 100, 1, 101);
     TH2F EEMRegion("EEMRegion", "EEMRegion", 100, 1, 101, 100, 1, 101);
     TH2F EEMRing("EEMRing", "EEMRing", 100, 1, 101, 100, 1, 101);
     TH2F EEMRingReg("EEMRingReg", "EEMRingReg", 100, 1, 101, 100, 1, 101);
     //  TH1F eta ("eta","eta",250,-85,165);
     for (std::vector<DetId>::const_iterator it = m_endcapCells.begin(); it != m_endcapCells.end(); ++it) {
       EEDetId ee(*it);
       if (ee.zside() > 0) {
         EEPRegion.Fill(ee.ix(), ee.iy(), m_xtalRegionId[ee.rawId()]);
         EEPRing.Fill(ee.ix(), ee.iy(), m_xtalRing[ee.rawId()]);
         EEPRingReg.Fill(ee.ix(), ee.iy(), m_RinginRegion[ee.rawId()]);
       }
       if (ee.zside() < 0) {
         EEMRegion.Fill(ee.ix(), ee.iy(), m_xtalRegionId[ee.rawId()]);
         EEMRing.Fill(ee.ix(), ee.iy(), m_xtalRing[ee.rawId()]);
         EEMRingReg.Fill(ee.ix(), ee.iy(), m_RinginRegion[ee.rawId()]);
       }
     }
 
     //  for (std::map<int,float>::iterator it=m_eta.begin();
     //        it!=m_eta.end();++it)
     //         eta.Fill(it->first,it->second);
     TFile out("EBZone.root", "recreate");
     EBRegion.Write();
     EBRing.Write();
     EEPRegion.Write();
     EEPRing.Write();
     EEPRingReg.Write();
     EEMRegion.Write();
     EEMRing.Write();
     //  eta.Write();
     EEMRingReg.Write();
     out.Close();
     edm::LogInfo("IML") << "[InvRingCalib] Start to acquire the coeffs";
     CaloMiscalibMapEcal EBmap;
     EBmap.prefillMap();
     MiscalibReaderFromXMLEcalBarrel barrelreader(EBmap);
     if (!m_EBcoeffFile.empty())
       barrelreader.parseXMLMiscalibFile(m_EBcoeffFile);
     EcalIntercalibConstants costants(EBmap.get());
     m_barrelMap = costants.getMap();
     CaloMiscalibMapEcal EEmap;
     EEmap.prefillMap();
     MiscalibReaderFromXMLEcalEndcap endcapreader(EEmap);
     if (!m_EEcoeffFile.empty())
       endcapreader.parseXMLMiscalibFile(m_EEcoeffFile);
     EcalIntercalibConstants EEcostants(EEmap.get());
     m_endcapMap = EEcostants.getMap();
 
     isfirstcall_ = false;
   }  // if isfirstcall
 
   //gets the barrel recHits
   double pSubtract = 0.;
   double pTk = 0.;
   const EcalRecHitCollection* barrelHitsCollection = nullptr;
   edm::Handle<EBRecHitCollection> barrelRecHitsHandle;
   iEvent.getByToken(m_ebRecHitToken, barrelRecHitsHandle);
   barrelHitsCollection = barrelRecHitsHandle.product();
 
   if (!barrelRecHitsHandle.isValid()) {
     edm::LogError("IML") << "[EcalEleCalibLooper] barrel rec hits not found";
     return kContinue;
   }
 
   //gets the endcap recHits
   const EcalRecHitCollection* endcapHitsCollection = nullptr;
   edm::Handle<EERecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(m_eeRecHitToken, endcapRecHitsHandle);
   endcapHitsCollection = endcapRecHitsHandle.product();
 
   if (!endcapRecHitsHandle.isValid()) {
     edm::LogError("IML") << "[EcalEleCalibLooper] endcap rec hits not found";
     return kContinue;
   }
 
   //gets the electrons
   edm::Handle<reco::GsfElectronCollection> pElectrons;
   iEvent.getByToken(m_gsfElectronToken, pElectrons);
 
   if (!pElectrons.isValid()) {
     edm::LogError("IML") << "[EcalEleCalibLooper] electrons not found";
     return kContinue;
   }
 
   //loops over the electrons in the event
   for (reco::GsfElectronCollection::const_iterator eleIt = pElectrons->begin(); eleIt != pElectrons->end(); ++eleIt) {
     pSubtract = 0;
     pTk = eleIt->trackMomentumAtVtx().R();
     std::map<int, double> xtlMap;
     DetId Max = 0;
     if (std::abs(eleIt->eta()) < 1.49)
       Max = EcalClusterTools::getMaximum(eleIt->superCluster()->hitsAndFractions(), barrelHitsCollection).first;
     else
       Max = EcalClusterTools::getMaximum(eleIt->superCluster()->hitsAndFractions(), endcapHitsCollection).first;
     if (Max.det() == 0)
       continue;
     m_MapFiller->fillMap(
         eleIt->superCluster()->hitsAndFractions(), Max, barrelHitsCollection, endcapHitsCollection, xtlMap, pSubtract);
     if (m_xtalRegionId[Max.rawId()] == -1)
       continue;
     pSubtract += eleIt->superCluster()->preshowerEnergy();
     ++m_RingNumOfHits[m_xtalRing[Max.rawId()]];
     //fills the calibBlock
     m_IMACalibBlocks.at(m_xtalRegionId[Max.rawId()])->Fill(xtlMap.begin(), xtlMap.end(), pTk, pSubtract);
   }
   return kContinue;
 }  //end of duringLoop

◆ EBRegId()

int InvRingCalib::EBRegId ( const int ieta )

private

Gives back in which region you are:

gives the region Id given ieta

Definition at line 533 of file InvRingCalib.cc.

References LEDCalibrationChannels::ieta, m_etaEnd, and m_Reg.

Referenced by EBRegionDef().

                                         {
   if (ieta < m_etaStart || ieta >= m_etaEnd)
     return -1;
   else
     return (m_Reg[ieta]);
 }

◆ EBRegionDef()

void InvRingCalib::EBRegionDef ( )

private

Defines the regions in the barrel.

Definition at line 543 of file InvRingCalib.cc.

References EBRegId(), m_barrelCells, m_etaStart, m_etaWidth, m_RinginRegion, m_xtalRegionId, m_xtalRing, and RegPrepare().

Referenced by duringLoop().

                                {
   RegPrepare();
   for (std::vector<DetId>::const_iterator it = m_barrelCells.begin(); it != m_barrelCells.end(); ++it) {
     EBDetId eb(it->rawId());
     m_xtalRing[eb.rawId()] = eb.ieta();
     m_xtalRegionId[eb.rawId()] = EBRegId(eb.ieta());
     if (m_xtalRegionId[eb.rawId()] == -1)
       continue;
     m_RinginRegion[eb.rawId()] = (eb.ieta() - m_etaStart) % m_etaWidth;
   }
 }

◆ EBRegionNum()

int InvRingCalib::EBRegionNum ( ) const

private

Number of regions in EB.

number of Ring in EB

Definition at line 512 of file InvRingCalib.cc.

References m_etaEnd, m_etaStart, and m_etaWidth.

Referenced by EERegionDef(), and InvRingCalib().

                                     {
   if ((m_etaEnd * m_etaStart) > 0)
     return ((m_etaEnd - m_etaStart) / m_etaWidth);
 
   if ((m_etaEnd * m_etaStart) < 0)
     return ((m_etaEnd - m_etaStart - 1) / m_etaWidth);
 
   return 0;
 }

◆ EERegId()

int InvRingCalib::EERegId ( int id )

private

gives back in which region of the endcap you are.

Gives the Id of the region given the id of the xtal.

Definition at line 439 of file InvRingCalib.cc.

References l1ctLayer2EG_cff::id, m_EEZone, m_endRing, m_etaWidth, m_RinginRegion, m_startRing, m_xtalRing, relativeConstraints::ring, and EEDetId::zside().

Referenced by EERegionDef().

                                 {
   int reg;
   int ring;
   EEDetId ee(id);
   //sets the reg to -1 if the ring doesn't exist or is outside the region of interest
   if (m_xtalRing[id] == -1)
     return -1;
   //avoid the calibration in the wrong zside
   if (m_EEZone == 1) {
     if (ee.zside() < 0)
       return -1;
     ring = m_xtalRing[id] - 86;
     if (ring >= m_endRing)
       return -1;
     if (ring < m_startRing)
       return -1;
     reg = (ring - m_startRing) / m_etaWidth;
     m_RinginRegion[id] = (ring - m_startRing) % m_etaWidth;
     return reg;
   }
   if (m_EEZone == -1) {
     if (ee.zside() > 0)
       return -1;
     ring = m_xtalRing[id] - 125;
     if (ring >= m_endRing)
       return -1;
     if (ring < m_startRing)
       return -1;
     reg = (ring - m_startRing) / m_etaWidth;
     m_RinginRegion[id] = (ring - m_startRing) % m_etaWidth;
     return reg;
   }
   if (ee.zside() > 0)
     ring = m_xtalRing[id] - 86;
   else
     ring = m_xtalRing[id] - 125;
   if (ring >= m_endRing)
     return -1;
   if (ring < m_startRing)
     return -1;
   reg = (ring - m_startRing) / m_etaWidth;
   m_RinginRegion[id] = (ring - m_startRing) % m_etaWidth;
   return reg;
 }

◆ EERegionDef()

void InvRingCalib::EERegionDef ( )

private

Defines the regions in the endcap.

Loops over all the endcap xtals and sets for each xtal the value of the region the xtal is in, and the ringNumber inside the region

Definition at line 486 of file InvRingCalib.cc.

References EBRegionNum(), EERegId(), EERegionNum(), m_endcapCells, m_xtalRegionId, and EEDetId::zside().

Referenced by duringLoop().

                                {
   int reg;
   for (std::vector<DetId>::const_iterator endcapIt = m_endcapCells.begin(); endcapIt != m_endcapCells.end();
        ++endcapIt) {
     EEDetId ee(*endcapIt);
     reg = EERegId(endcapIt->rawId());
     //If the ring is not of interest saves only the region Id(-1)
     if (reg == -1)
       m_xtalRegionId[endcapIt->rawId()] = reg;
     //sums the number of region in EB or EB+EE to have different regionsId in different regions
     else {
       if (ee.zside() > 0)
         reg += EBRegionNum();
       else
         reg += EBRegionNum() + EERegionNum();
       m_xtalRegionId[endcapIt->rawId()] = reg;
     }
   }
 }

◆ EERegionNum()

int InvRingCalib::EERegionNum ( ) const

inlineprivate

The number of regions in EE.

Number of Regions in EE.

Definition at line 509 of file InvRingCalib.cc.

References m_endRing, m_etaWidth, and m_startRing.

Referenced by EERegionDef(), and InvRingCalib().

509 { return ((m_endRing - m_startRing) / m_etaWidth); }

InvRingCalib::m_etaWidth

const int m_etaWidth

eta size of the regions

Definition: InvRingCalib.h:84

InvRingCalib::m_endRing

const int m_endRing

Definition: InvRingCalib.h:109

InvRingCalib::m_startRing

const int m_startRing

LP define the EE region to calibrate.

Definition: InvRingCalib.h:108

◆ EERingDef()

void InvRingCalib::EERingDef ( const edm::EventSetup & iSetup )

private

Defines the rins in the endcap.

EE ring definition.

Definition at line 383 of file InvRingCalib.cc.

References DetId::Ecal, EcalEndcap, edm::EventSetup::getData(), CaloSubdetectorGeometry::getGeometry(), m_cellPhi, m_cellPos, m_endcapCells, m_geometryToken, m_xtalRing, DetId::rawId(), relativeConstraints::ring, EEDetId::validDetId(), EEDetId::XYMODE, and EEDetId::zside().

Referenced by duringLoop().

                                                         {
   //Gets the geometry of the endcap
   const auto& geometry = iSetup.getData(m_geometryToken);
   const CaloSubdetectorGeometry* endcapGeometry = geometry.getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
   //for every xtal gets the position Vector and the phi position
 
   // for (std::vector<DetId>::const_iterator barrelIt = m_barrelCells.begin();
   //    barrelIt!=m_barrelCells.end();
   //    ++barrelIt) {
   //     const CaloCellGeometry *cellGeometry = barrelGeometry->getGeometry(*barrelIt);
   //     GlobalPoint point;
   //     EBDetId eb (*barrelIt);
   //     point=cellGeometry->getPosition();
   //     m_eta[eb.ieta()]=point.eta() ;    //cellGeometry->getPosition().eta();
   //    }
   for (std::vector<DetId>::const_iterator endcapIt = m_endcapCells.begin(); endcapIt != m_endcapCells.end();
        ++endcapIt) {
     auto cellGeometry = endcapGeometry->getGeometry(*endcapIt);
     m_cellPos[endcapIt->rawId()] = cellGeometry->getPosition();
     m_cellPhi[endcapIt->rawId()] = cellGeometry->getPosition().phi();
   }
   //takes the first 39 xtals at a fixed y varying the x coordinate and saves their eta coordinate
   float eta_ring[39];
   for (int ring = 0; ring < 39; ring++)
     if (EEDetId::validDetId(ring, 50, 1)) {
       EEDetId det = EEDetId(ring, 50, 1, EEDetId::XYMODE);
       eta_ring[ring] = m_cellPos[det.rawId()].eta();
     }
   //defines the bonduary of the rings as the average eta of a xtal
   double etaBonduary[40];
   etaBonduary[0] = 1.49;
   etaBonduary[39] = 4.0;
   for (int ring = 1; ring < 39; ++ring)
     etaBonduary[ring] = (eta_ring[ring] + eta_ring[ring - 1]) / 2.;
   //assign to each xtal a ring number
   int CRing;
   for (int ring = 0; ring < 39; ring++)
     for (std::vector<DetId>::const_iterator endcapIt = m_endcapCells.begin(); endcapIt != m_endcapCells.end();
          ++endcapIt) {
       if (fabs(m_cellPos[endcapIt->rawId()].eta()) > etaBonduary[ring] &&
           fabs(m_cellPos[endcapIt->rawId()].eta()) < etaBonduary[ring + 1]) {
         EEDetId ee(*endcapIt);
         if (ee.zside() > 0)
           CRing = ring + 86;
         else
           CRing = ring + 125;
         m_xtalRing[endcapIt->rawId()] = CRing;
         //              m_eta[CRing]=m_cellPos[endcapIt->rawId()].eta();
       }
     }
   return;
 }

◆ endOfJob()

void InvRingCalib::endOfJob ( )

overridevirtual

Reimplemented from edm::EDLooperBase.

Definition at line 360 of file InvRingCalib.cc.

References EcalBarrel, EcalEndcap, m_barrelCells, m_barrelMap, m_endcapCells, m_endcapMap, and calibXMLwriter::writeLine().

                             {
   edm::LogInfo("IML") << "[InvMatrixCalibLooper][endOfJob] saving calib coeffs";
   calibXMLwriter barrelWriter(EcalBarrel);
   calibXMLwriter endcapWriter(EcalEndcap);
   for (std::vector<DetId>::const_iterator barrelIt = m_barrelCells.begin(); barrelIt != m_barrelCells.end();
        ++barrelIt) {
     EBDetId eb(*barrelIt);
     barrelWriter.writeLine(eb, m_barrelMap[eb]);
   }
   for (std::vector<DetId>::const_iterator endcapIt = m_endcapCells.begin(); endcapIt != m_endcapCells.end();
        ++endcapIt) {
     EEDetId ee(*endcapIt);
     endcapWriter.writeLine(ee, m_endcapMap[ee]);
   }
 }

◆ endOfLoop()

edm::EDLooper::Status InvRingCalib::endOfLoop	(	const edm::EventSetup &	,
		unsigned int	iCounter
	)

overridevirtual

Called after the system has finished one loop over the events. Thar argument is a count of how many loops have been processed before this loo. For the first time through the events the argument will be 0.

Implements edm::EDLooperBase.

Definition at line 290 of file InvRingCalib.cc.

References corrVsCorr::filename, RemoveAddSevLevel::flag, cms::cuda::for(), edm::EDLooperBase::kContinue, edm::EDLooperBase::kStop, m_barrelCells, m_barrelMap, m_endcapCells, m_endcapMap, m_IMACalibBlocks, m_loops, m_maxCoeff, m_minCoeff, m_RinginRegion, m_RingNumOfHits, m_usingBlockSolver, m_xtalRegionId, m_xtalRing, and MillePedeFileConverter_cfg::out.

                                                                                           {
   std::map<int, double> InterRings;
   edm::LogInfo("IML") << "[InvMatrixCalibLooper][endOfLoop] Start to invert the matrixes";
   //call the autoexplaining "solve" method for every calibBlock
   for (std::vector<VEcalCalibBlock*>::iterator calibBlock = m_IMACalibBlocks.begin();
        calibBlock != m_IMACalibBlocks.end();
        ++calibBlock)
     (*calibBlock)->solve(m_usingBlockSolver, m_minCoeff, m_maxCoeff);
 
   edm::LogInfo("IML") << "[InvRingLooper][endOfLoop] Starting to write the coeffs";
   TH1F* coeffDistr = new TH1F("coeffdistr", "coeffdistr", 100, 0.7, 1.4);
   TH1F* coeffMap = new TH1F("coeffRingMap", "coeffRingMap", 250, -85, 165);
   TH1F* ringDistr = new TH1F("ringDistr", "ringDistr", 250, -85, 165);
   TH1F* RingFill = new TH1F("RingFill", "RingFill", 250, -85, 165);
   for (std::map<int, int>::const_iterator it = m_xtalRing.begin(); it != m_xtalRing.end(); ++it)
     ringDistr->Fill(it->second + 0.1);
 
   int ID;
   std::map<int, int> flag;
   for (std::map<int, int>::const_iterator it = m_xtalRing.begin(); it != m_xtalRing.end(); ++it)
     flag[it->second] = 0;
 
   for (std::vector<DetId>::const_iterator it = m_barrelCells.begin(); it != m_barrelCells.end(); ++it) {
     ID = it->rawId();
     if (m_xtalRegionId[ID] == -1)
       continue;
     if (flag[m_xtalRing[ID]])
       continue;
     flag[m_xtalRing[ID]] = 1;
     RingFill->Fill(m_xtalRing[ID], m_RingNumOfHits[m_xtalRing[ID]]);
     InterRings[m_xtalRing[ID]] = m_IMACalibBlocks.at(m_xtalRegionId[ID])->at(m_RinginRegion[ID]);
     coeffMap->Fill(m_xtalRing[ID] + 0.1, InterRings[m_xtalRing[ID]]);
     coeffDistr->Fill(InterRings[m_xtalRing[ID]]);
   }
 
   for (std::vector<DetId>::const_iterator it = m_endcapCells.begin(); it != m_endcapCells.end(); ++it) {
     ID = it->rawId();
     if (m_xtalRegionId[ID] == -1)
       continue;
     if (flag[m_xtalRing[ID]])
       continue;
     flag[m_xtalRing[ID]] = 1;
     InterRings[m_xtalRing[ID]] = m_IMACalibBlocks.at(m_xtalRegionId[ID])->at(m_RinginRegion[ID]);
     RingFill->Fill(m_xtalRing[ID], m_RingNumOfHits[m_xtalRing[ID]]);
     coeffMap->Fill(m_xtalRing[ID], InterRings[m_xtalRing[ID]]);
     coeffDistr->Fill(InterRings[m_xtalRing[ID]]);
   }
 
   char filename[80];
   sprintf(filename, "coeff%d.root", iCounter);
   TFile out(filename, "recreate");
   coeffDistr->Write();
   coeffMap->Write();
   ringDistr->Write();
   RingFill->Write();
   out.Close();
   for (std::vector<DetId>::const_iterator it = m_barrelCells.begin(); it != m_barrelCells.end(); ++it) {
     m_barrelMap[*it] *= InterRings[m_xtalRing[it->rawId()]];
   }
   for (std::vector<DetId>::const_iterator it = m_endcapCells.begin(); it != m_endcapCells.end(); ++it)
     m_endcapMap[*it] *= InterRings[m_xtalRing[it->rawId()]];
   if (iCounter < m_loops - 1)
     return kContinue;
   else
     return kStop;
 }

◆ RegPrepare()

void InvRingCalib::RegPrepare ( )

private

Prepares the EB regions;.

Divides the barrel in region, necessary to take into account the missing 0 xtal

Definition at line 523 of file InvRingCalib.cc.

References mps_fire::i, dqmdumpme::k, m_etaEnd, m_etaStart, m_etaWidth, and m_Reg.

Referenced by EBRegionDef().

                               {
   int k = 0;
   for (int i = m_etaStart; i < m_etaEnd; ++i) {
     if (i == 0)
       continue;
     m_Reg[i] = k / m_etaWidth;
     ++k;
   }
 }

◆ startingNewLoop()

void InvRingCalib::startingNewLoop ( unsigned int ciclo )

overridevirtual

startingNewLoop

Implements edm::EDLooperBase.

Definition at line 113 of file InvRingCalib.cc.

References m_IMACalibBlocks, m_RingNumOfHits, m_xtalRing, and relativeConstraints::ring.

                                                      {
   edm::LogInfo("IML") << "[InvMatrixCalibLooper][Start] entering loop " << ciclo;
   for (std::vector<VEcalCalibBlock*>::iterator calibBlock = m_IMACalibBlocks.begin();
        calibBlock != m_IMACalibBlocks.end();
        ++calibBlock) {
     //LP empties the energies vector, to fill DuringLoop.
     (*calibBlock)->reset();
   }
   for (std::map<int, int>::const_iterator ring = m_xtalRing.begin(); ring != m_xtalRing.end(); ++ring)
     m_RingNumOfHits[ring->second] = 0;
   return;
 }

Member Data Documentation

◆ isfirstcall_

bool InvRingCalib::isfirstcall_

private

Definition at line 128 of file InvRingCalib.h.

Referenced by beginOfJob(), and duringLoop().

◆ m_barrelAlCa

const edm::InputTag InvRingCalib::m_barrelAlCa

private

EcalBarrel Input Collection name.

Definition at line 67 of file InvRingCalib.h.

◆ m_barrelCells

std::vector<DetId> InvRingCalib::m_barrelCells

private

geometry things used all over the file

Definition at line 118 of file InvRingCalib.h.

Referenced by duringLoop(), EBRegionDef(), endOfJob(), and endOfLoop().

◆ m_barrelMap

EcalIntercalibConstantMap InvRingCalib::m_barrelMap

private

association map between coeff and ring coeffs for the single xtals

Definition at line 103 of file InvRingCalib.h.

Referenced by duringLoop(), endOfJob(), and endOfLoop().

◆ m_cellPhi

std::map<int, int> InvRingCalib::m_cellPhi

private

Definition at line 100 of file InvRingCalib.h.

Referenced by EERingDef().

◆ m_cellPos

std::map<int, GlobalPoint> InvRingCalib::m_cellPos

private

position of the cell, borders, coords etc...

Definition at line 99 of file InvRingCalib.h.

Referenced by EERingDef().

◆ m_EBcoeffFile

const std::string InvRingCalib::m_EBcoeffFile

private

coeffs filenames

Definition at line 121 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_ebRecHitToken

const edm::EDGetTokenT<EBRecHitCollection> InvRingCalib::m_ebRecHitToken

private

ED token.

Definition at line 131 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_EEcoeffFile

const std::string InvRingCalib::m_EEcoeffFile

private

Definition at line 122 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_eeRecHitToken

const edm::EDGetTokenT<EERecHitCollection> InvRingCalib::m_eeRecHitToken

private

Definition at line 132 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_EEZone

const int InvRingCalib::m_EEZone

private

endcap zone to be calibrated

Definition at line 124 of file InvRingCalib.h.

Referenced by EERegId(), and InvRingCalib().

◆ m_ElectronLabel

const edm::InputTag InvRingCalib::m_ElectronLabel

private

To take the electrons.

Definition at line 71 of file InvRingCalib.h.

◆ m_endcapAlCa

const edm::InputTag InvRingCalib::m_endcapAlCa

private

EcalEndcap Input Collection name.

Definition at line 69 of file InvRingCalib.h.

◆ m_endcapCells

std::vector<DetId> InvRingCalib::m_endcapCells

private

Definition at line 119 of file InvRingCalib.h.

Referenced by duringLoop(), EERegionDef(), EERingDef(), endOfJob(), and endOfLoop().

◆ m_endcapMap

EcalIntercalibConstantMap InvRingCalib::m_endcapMap

private

Definition at line 104 of file InvRingCalib.h.

Referenced by duringLoop(), endOfJob(), and endOfLoop().

◆ m_endRing

const int InvRingCalib::m_endRing

private

Definition at line 109 of file InvRingCalib.h.

Referenced by EERegId(), EERegionNum(), and InvRingCalib().

◆ m_etaEnd

const int InvRingCalib::m_etaEnd

private

eta end of the zone of interest

Definition at line 82 of file InvRingCalib.h.

Referenced by EBRegId(), EBRegionNum(), InvRingCalib(), and RegPrepare().

◆ m_etaStart

const int InvRingCalib::m_etaStart

private

eta start of the zone of interest

Definition at line 80 of file InvRingCalib.h.

Referenced by EBRegionDef(), EBRegionNum(), InvRingCalib(), and RegPrepare().

◆ m_etaWidth

const int InvRingCalib::m_etaWidth

private

eta size of the regions

Definition at line 84 of file InvRingCalib.h.

Referenced by EBRegionDef(), EBRegionNum(), EERegId(), EERegionNum(), InvRingCalib(), and RegPrepare().

◆ m_geometryToken

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> InvRingCalib::m_geometryToken

private

ES token.

Definition at line 135 of file InvRingCalib.h.

Referenced by duringLoop(), and EERingDef().

◆ m_gsfElectronToken

const edm::EDGetTokenT<reco::GsfElectronCollection> InvRingCalib::m_gsfElectronToken

private

Definition at line 133 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_IMACalibBlocks

std::vector<VEcalCalibBlock *> InvRingCalib::m_IMACalibBlocks

private

single blocks calibrators

Definition at line 91 of file InvRingCalib.h.

Referenced by duringLoop(), endOfLoop(), InvRingCalib(), and startingNewLoop().

◆ m_loops

unsigned int InvRingCalib::m_loops

private

LP sets the number of loops to do.

Definition at line 106 of file InvRingCalib.h.

Referenced by endOfLoop(), and InvRingCalib().

◆ m_MapFiller

VFillMap* InvRingCalib::m_MapFiller

private

The class that fills the map!

Definition at line 63 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_mapFillerType

std::string InvRingCalib::m_mapFillerType

private

Definition at line 127 of file InvRingCalib.h.

Referenced by duringLoop(), and InvRingCalib().

◆ m_maxCoeff

const double InvRingCalib::m_maxCoeff

private

maximum coefficient accepted (RAW)

Definition at line 95 of file InvRingCalib.h.

Referenced by endOfLoop().

◆ m_maxEnergyPerCrystal

const double InvRingCalib::m_maxEnergyPerCrystal

private

maximum energy per crystal cut

Definition at line 78 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_maxSelectedNumPerRing

const int InvRingCalib::m_maxSelectedNumPerRing

private

maximum number of events per Ring

Definition at line 87 of file InvRingCalib.h.

◆ m_minCoeff

const double InvRingCalib::m_minCoeff

private

minimum coefficient accepted (RAW)

Definition at line 93 of file InvRingCalib.h.

Referenced by endOfLoop().

◆ m_minEnergyPerCrystal

const double InvRingCalib::m_minEnergyPerCrystal

private

minimum energy per crystal cut

Definition at line 76 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_recoWindowSidex

const int InvRingCalib::m_recoWindowSidex

private

reconstruction window size

Definition at line 73 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_recoWindowSidey

const int InvRingCalib::m_recoWindowSidey

private

Definition at line 74 of file InvRingCalib.h.

Referenced by duringLoop().

◆ m_Reg

std::map<int, int> InvRingCalib::m_Reg

private

EB regions vs. eta index.

Definition at line 126 of file InvRingCalib.h.

Referenced by EBRegId(), and RegPrepare().

◆ m_RinginRegion

std::map<int, int> InvRingCalib::m_RinginRegion

private

association map between raw detIds and the number of the ring inside the region

Definition at line 115 of file InvRingCalib.h.

Referenced by duringLoop(), EBRegionDef(), EERegId(), and endOfLoop().

◆ m_RingNumOfHits

std::map<int, int> InvRingCalib::m_RingNumOfHits

private

number of events already read per Ring

Definition at line 89 of file InvRingCalib.h.

Referenced by duringLoop(), endOfLoop(), and startingNewLoop().

◆ m_startRing

const int InvRingCalib::m_startRing

private

LP define the EE region to calibrate.

Definition at line 108 of file InvRingCalib.h.

Referenced by EERegId(), EERegionNum(), and InvRingCalib().

◆ m_usingBlockSolver

const int InvRingCalib::m_usingBlockSolver

private

to exclude the blocksolver

Definition at line 97 of file InvRingCalib.h.

Referenced by endOfLoop().

◆ m_xtalRegionId

std::map<int, int> InvRingCalib::m_xtalRegionId

private

association map between raw detIds and Region

Definition at line 113 of file InvRingCalib.h.

Referenced by duringLoop(), EBRegionDef(), EERegionDef(), and endOfLoop().

◆ m_xtalRing

std::map<int, int> InvRingCalib::m_xtalRing

private

association map between Raw detIds and Rings

Definition at line 111 of file InvRingCalib.h.

Referenced by duringLoop(), EBRegionDef(), EERegId(), EERingDef(), endOfLoop(), and startingNewLoop().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ InvRingCalib()

◆ ~InvRingCalib()

Member Function Documentation

◆ beginOfJob()

◆ duringLoop()

◆ EBRegId()

◆ EBRegionDef()

◆ EBRegionNum()

◆ EERegId()

◆ EERegionDef()

◆ EERegionNum()

◆ EERingDef()

◆ endOfJob()

◆ endOfLoop()

◆ RegPrepare()

◆ startingNewLoop()

Member Data Documentation

◆ isfirstcall_

◆ m_barrelAlCa

◆ m_barrelCells

◆ m_barrelMap

◆ m_cellPhi

◆ m_cellPos

◆ m_EBcoeffFile

◆ m_ebRecHitToken

◆ m_EEcoeffFile

◆ m_eeRecHitToken

◆ m_EEZone

◆ m_ElectronLabel

◆ m_endcapAlCa

◆ m_endcapCells

◆ m_endcapMap

◆ m_endRing

◆ m_etaEnd

◆ m_etaStart

◆ m_etaWidth

◆ m_geometryToken

◆ m_gsfElectronToken

◆ m_IMACalibBlocks

◆ m_loops

◆ m_MapFiller

◆ m_mapFillerType

◆ m_maxCoeff

◆ m_maxEnergyPerCrystal

◆ m_maxSelectedNumPerRing

◆ m_minCoeff

◆ m_minEnergyPerCrystal

◆ m_recoWindowSidex

◆ m_recoWindowSidey

◆ m_Reg

◆ m_RinginRegion

◆ m_RingNumOfHits

◆ m_startRing

◆ m_usingBlockSolver

◆ m_xtalRegionId

◆ m_xtalRing