EcalEleCalibLooper.h

Go to the documentation of this file.

#ifndef __CINT__
#ifndef EcalEleCalibLooper_H
#define EcalEleCalibLooper_H
#include "Calibration/EcalCalibAlgos/interface/VEcalCalibBlock.h"
#include "FWCore/Framework/interface/EDLooper.h"
#include "FWCore/Framework/interface/Event.h"
#include "Calibration/EcalCalibAlgos/interface/VFillMap.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include <string>
#include <vector>

#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EgammaReco/interface/BasicClusterShapeAssociation.h"
#include "DataFormats/EgammaCandidates/interface/ElectronFwd.h"
#include "DataFormats/EgammaCandidates/interface/Electron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"

#include "CLHEP/Matrix/GenMatrix.h"
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/Vector.h"
class EcalEleCalibLooper : public edm::EDLooper {
public:
  explicit EcalEleCalibLooper(const edm::ParameterSet &);
  ~EcalEleCalibLooper() 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;

private:
  //DS to divide in Regions
  int EBRegionId(const int, const int) const;
  int EERegionId(const int, const int) const;
  //DS to define the regions for each cristal
  void EBRegionDefinition();
  void EERegionDefinition();
  //DS defines the limit for the tan of phi
  double giveLimit(int);

  //DS checks if the values of ics and ips are in EE or not
  int EEregionCheck(const int, const int) const;
  int EBregionCheck(const int eta, const int phi) const;

  int etaShifter(const int) const;

private:
  const edm::InputTag m_barrelAlCa;
  const edm::InputTag m_endcapAlCa;
  edm::InputTag m_ElectronLabel;

  const int m_recoWindowSidex;
  const int m_recoWindowSidey;

  const int m_etaWidth;  //PG sub matrix size and borders
  //    int m_etaBorder ; //FIXME
  const int m_phiWidthEB;
  //    int m_phiBorderEB //FIXME;

  const int m_etaStart;  //PG ECAL region to be calibrated
  const int m_etaEnd;
  const int m_phiStartEB;
  const int m_phiEndEB;
  const int m_radStart;
  const int m_radEnd;
  const int m_radWidth;
  //FIXME    int m_radBorder ;
  const int m_phiStartEE;
  const int m_phiEndEE;
  const int m_phiWidthEE;

  const int m_maxSelectedNumPerXtal;

  std::vector<VEcalCalibBlock *> m_EcalCalibBlocks;
  const double m_minEnergyPerCrystal;
  const double m_maxEnergyPerCrystal;
  const double m_minCoeff;
  const double m_maxCoeff;
  const int m_usingBlockSolver;

  EcalIntercalibConstantMap m_barrelMap;
  EcalIntercalibConstantMap m_endcapMap;

  const unsigned int m_loops;
  //The map Filler
  VFillMap *m_MapFiller;

  //DS number of regions in the Endcap
  inline int EEregionsNum() const;
  //DS number of regions in the barrel
  inline int EBregionsNum() const;

  std::vector<int> m_regions;

  std::vector<DetId> m_barrelCells;
  std::vector<DetId> m_endcapCells;

  std::map<int, int> m_xtalRegionId;
  std::map<int, int> m_xtalPositionInRegion;
  std::map<int, int> m_xtalNumOfHits;

  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