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();
    void beginOfJob() ;
    void endOfJob();
    void startingNewLoop(unsigned int) ;
    Status duringLoop(const edm::Event&, const edm::EventSetup&) ;
    Status endOfLoop(const edm::EventSetup&,unsigned int iCounter) ;
    

  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