InvRingCalib Class Reference

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 ( const edm::ParameterSet &  iConfig )

explicit

ctor

Definition at line 26 of file InvRingCalib.cc.

References funct::abs(), HLT_FULL_cff::algorithm, cms::cuda::assert(), EBRegionNum(), EERegionNum(), beamvalidation::exit(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, 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 ( )

override

dtor

destructor

Definition at line 103 of file InvRingCalib.cc.

{}

Member Function Documentation

void InvRingCalib::beginOfJob ( )

overridevirtual

BeginOfJob.

Reimplemented from edm::EDLooperBase.

Definition at line 108 of file InvRingCalib.cc.

References isfirstcall_.

{ isfirstcall_ = true; }

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(), first, relativeConstraints::geometry, CaloMiscalibMapEcal::get(), edm::Event::getByToken(), edm::EventSetup::getData(), EcalCondObjectContainer< T >::getMap(), EBDetId::ieta(), 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, siStripFEDMonitor_P5_cff::Max, submitPVResolutionJobs::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

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 m_etaEnd, and m_Reg.

Referenced by EBRegionDef().

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

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;
  }
}

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;
}

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 gpuClustering::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;
}

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;
    }
  }
}

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

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

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, relativeConstraints::geometry, 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;
}

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]);
  }
}

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 lut2db_cfg::filename, 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 submitPVResolutionJobs::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;
}

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, isotrackApplyRegressor::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;
  }
}

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

bool InvRingCalib::isfirstcall_

private

Definition at line 128 of file InvRingCalib.h.

Referenced by beginOfJob(), and duringLoop().

const edm::InputTag InvRingCalib::m_barrelAlCa

private

EcalBarrel Input Collection name.

Definition at line 67 of file InvRingCalib.h.

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

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

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

private

Definition at line 100 of file InvRingCalib.h.

Referenced by EERingDef().

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

const std::string InvRingCalib::m_EBcoeffFile

private

coeffs filenames

Definition at line 121 of file InvRingCalib.h.

Referenced by duringLoop().

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

private

ED token.

Definition at line 131 of file InvRingCalib.h.

Referenced by duringLoop().

const std::string InvRingCalib::m_EEcoeffFile

private

Definition at line 122 of file InvRingCalib.h.

Referenced by duringLoop().

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

private

Definition at line 132 of file InvRingCalib.h.

Referenced by duringLoop().

const int InvRingCalib::m_EEZone

private

endcap zone to be calibrated

Definition at line 124 of file InvRingCalib.h.

Referenced by EERegId(), and InvRingCalib().

const edm::InputTag InvRingCalib::m_ElectronLabel

private

To take the electrons.

Definition at line 71 of file InvRingCalib.h.

const edm::InputTag InvRingCalib::m_endcapAlCa

private

EcalEndcap Input Collection name.

Definition at line 69 of file InvRingCalib.h.

std::vector<DetId> InvRingCalib::m_endcapCells

private

Definition at line 119 of file InvRingCalib.h.

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

EcalIntercalibConstantMap InvRingCalib::m_endcapMap

private

Definition at line 104 of file InvRingCalib.h.

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

const int InvRingCalib::m_endRing

private

Definition at line 109 of file InvRingCalib.h.

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

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

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

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

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

private

ES token.

Definition at line 135 of file InvRingCalib.h.

Referenced by duringLoop(), and EERingDef().

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

private

Definition at line 133 of file InvRingCalib.h.

Referenced by duringLoop().

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

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

VFillMap* InvRingCalib::m_MapFiller

private

The class that fills the map!

Definition at line 63 of file InvRingCalib.h.

Referenced by duringLoop().

std::string InvRingCalib::m_mapFillerType

private

Definition at line 127 of file InvRingCalib.h.

Referenced by duringLoop(), and InvRingCalib().

const double InvRingCalib::m_maxCoeff

private

maximum coefficient accepted (RAW)

Definition at line 95 of file InvRingCalib.h.

Referenced by endOfLoop().

const double InvRingCalib::m_maxEnergyPerCrystal

private

maximum energy per crystal cut

Definition at line 78 of file InvRingCalib.h.

Referenced by duringLoop().

const int InvRingCalib::m_maxSelectedNumPerRing

private

maximum number of events per Ring

Definition at line 87 of file InvRingCalib.h.

const double InvRingCalib::m_minCoeff

private

minimum coefficient accepted (RAW)

Definition at line 93 of file InvRingCalib.h.

Referenced by endOfLoop().

const double InvRingCalib::m_minEnergyPerCrystal

private

minimum energy per crystal cut

Definition at line 76 of file InvRingCalib.h.

Referenced by duringLoop().

const int InvRingCalib::m_recoWindowSidex

private

reconstruction window size

Definition at line 73 of file InvRingCalib.h.

Referenced by duringLoop().

const int InvRingCalib::m_recoWindowSidey

private

Definition at line 74 of file InvRingCalib.h.

Referenced by duringLoop().

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

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

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

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

const int InvRingCalib::m_usingBlockSolver

private

to exclude the blocksolver

Definition at line 97 of file InvRingCalib.h.

Referenced by endOfLoop().

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

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