InvRingCalib.h

Go to the documentation of this file.

#ifndef __CINT__
#ifndef InvRingCalib_H
#define InvRingCalib_H
#include "Calibration/EcalCalibAlgos/interface/VEcalCalibBlock.h"
#include "FWCore/Framework/interface/EDLooper.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/CaloRecHit/interface/CaloRecHit.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"

#include <string>
#include <vector>

#include "Calibration/EcalCalibAlgos/interface/VFillMap.h"

class InvRingCalib : public edm::EDLooper {
public:
  explicit InvRingCalib(const edm::ParameterSet &);
  ~InvRingCalib() override;
  void beginOfJob() override;
  void endOfJob() override;
  void startingNewLoop(unsigned int) override;
  Status duringLoop(const edm::Event &, const edm::EventSetup &) override;
  Status endOfLoop(const edm::EventSetup &, unsigned int iCounter) override;

  //end

private:
  inline int EERegionNum() const;
  int EBRegionNum() const;
  void EBRegionDef();
  void EERingDef(const edm::EventSetup &);
  void EERegionDef();
  void RegPrepare();
  int EBRegId(const int);
  int EERegId(int);

  VFillMap *m_MapFiller;

private:
  const edm::InputTag m_barrelAlCa;
  const edm::InputTag m_endcapAlCa;
  const edm::InputTag m_ElectronLabel;
  const int m_recoWindowSidex;
  const int m_recoWindowSidey;
  const double m_minEnergyPerCrystal;
  const double m_maxEnergyPerCrystal;
  const int m_etaStart;
  const int m_etaEnd;
  const int m_etaWidth;
  //    std::map<int,float> m_eta;
  const int m_maxSelectedNumPerRing;
  std::map<int, int> m_RingNumOfHits;
  std::vector<VEcalCalibBlock *> m_IMACalibBlocks;
  const double m_minCoeff;
  const double m_maxCoeff;
  const int m_usingBlockSolver;
  std::map<int, GlobalPoint> m_cellPos;
  std::map<int, int> m_cellPhi;
  EcalIntercalibConstantMap m_barrelMap;
  EcalIntercalibConstantMap m_endcapMap;
  unsigned int m_loops;
  const int m_startRing;
  const int m_endRing;
  std::map<int, int> m_xtalRing;
  std::map<int, int> m_xtalRegionId;
  std::map<int, int> m_RinginRegion;

  std::vector<DetId> m_barrelCells;
  std::vector<DetId> m_endcapCells;
  const std::string m_EBcoeffFile;
  const std::string m_EEcoeffFile;
  const int m_EEZone;
  std::map<int, int> m_Reg;
  std::string m_mapFillerType;
  bool isfirstcall_;

  const edm::EDGetTokenT<EBRecHitCollection> m_ebRecHitToken;
  const edm::EDGetTokenT<EERecHitCollection> m_eeRecHitToken;
  const edm::EDGetTokenT<reco::GsfElectronCollection> m_gsfElectronToken;
  const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> m_geometryToken;
};
#endif
#endif