---

EcalEleCalibLooper.h

Go to the documentation of this file.

#ifndef __CINT__
#ifndef EcalEleCalibLooper_H
#define EcalEleCalibLooper_H
#include "Calibration/EcalCalibAlgos/interface/VEcalCalibBlock.h"
#include "Calibration/Tools/interface/smartSelector.h"
#include "FWCore/Framework/interface/EDLooper.h"
#include "FWCore/Framework/interface/Event.h"

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

  private:

    typedef reco::GsfElectronCollection::const_iterator eleIterator;
    typedef edm::Handle<reco::BasicClusterShapeAssociationCollection>  HandleBasicCSAC; 

    DetId getMaxId(eleIterator EleIt,
                   HandleBasicCSAC  & barrelClShpHandle,
                   HandleBasicCSAC  & endcapClShpHandle);

    DetId findMaxHit (const std::vector<DetId> & v1,
                                  const EBRecHitCollection* EBhits,
                                  const EERecHitCollection* EEhits) ;

    void fillEBMap (EBDetId EBmax,
                    const EcalRecHitCollection * barrelHitsCollection,
                    std::map<int,double> & EBXtlMap,
                    int EBNumberOfRegion, double & pSubtract ) ;

    void fillEEMap (EEDetId EEmax,
                    const EcalRecHitCollection * endcapHitsCollection,
                    std::map<int,double> & EExtlMap,
                    int EENumberOfRegion, double & pSubtract ) ;


    int makeReport (std::string baseName="output") ;
    int evalKaliX2Num () ;

 //   int indexFinder (const int etaWorld, const int phiWorld) ;//FIXME
  //  TH1F * fillTrend (std::vector<double> const & vettore, 
    //                  const int & index) ;
   int findVoidLine (const CLHEP::HepMatrix & suspect) ;
  //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_recoWindowSide ;

    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_EBxtlNum[170][360] ;
    int m_EBxtlReg[170][360] ;
    int m_EExtlNum[100][100] ;
    int m_EExtlReg[100][100] ;

//FIXME    double m_halfXBand ;
//FIXME    double m_halfYBand ;
    int m_maxSelectedNumPerXtal ;  
//FIXME    int m_smallestFraction ;
//FIXME    int m_howManyFractions ;
//FIXME!!    smartSelector m_eventSelector ;
//    int halfSelecting ;
//    int takeOdd ;

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

    std::map<int,double> m_recalibMap ;

    unsigned int m_loops ;
    edm::InputTag m_ElectronLabel ;


  //DS numero delle regioni lungo il raggio (onion rings) (da fare divisione lungo phi)
  inline int EEregionsNum () const ;
  //DS numero delle regioni in EB
  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;
  //  std::map<int,double> m_miscalibMap;
};
#endif
#endif

---