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 "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:
  edm::InputTag m_barrelAlCa;
  edm::InputTag m_endcapAlCa;

  int m_recoWindowSidex;
  int m_recoWindowSidey;

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

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

  int m_maxSelectedNumPerXtal;

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

  EcalIntercalibConstantMap m_barrelMap;
  EcalIntercalibConstantMap m_endcapMap;

  unsigned int m_loops;
  edm::InputTag m_ElectronLabel;
  //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_;
};
#endif
#endif