ElectronCalibration.h

Go to the documentation of this file.

#ifndef CALIBRATION_ECALCALIBALGOS_ELECTRONCALIBRATION
#define CALIBRATION_ECALCALIBALGOS_ELECTRONCALIBRATION

// -*- C++ -*-
//
// Package:    ElectronCalibration
// Class:      ElectronCalibration
//
//
// Original Author:  Lorenzo AGOSTINO, Radek Ofierzynski
//         Created:  Tue Jul 18 12:17:01 CEST 2006
//
//

// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "Calibration/Tools/interface/HouseholderDecomposition.h"
#include "Calibration/Tools/interface/MinL3Algorithm.h"
#include "Calibration/Tools/interface/CalibrationCluster.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "TFile.h"
#include "TH1.h"
#include "TH2.h"

// class declaration

class ElectronCalibration : public edm::one::EDAnalyzer<> {
public:
  explicit ElectronCalibration(const edm::ParameterSet&);
  ~ElectronCalibration() override;

  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void beginJob() override;
  void endJob() override;

private:
  EBDetId findMaxHit(edm::Handle<EBRecHitCollection>&);
  EBDetId findMaxHit2(const std::vector<DetId>&, const EBRecHitCollection*);

  // ----------member data ---------------------------
  std::string rootfile_;
  edm::EDGetTokenT<EBRecHitCollection> recHitToken_;
  edm::EDGetTokenT<reco::GsfElectronCollection> electronToken_;
  edm::InputTag trackLabel_;
  std::string calibAlgo_;
  std::string miscalibfile_;
  int keventweight_;
  double ElePt_;
  int maxeta_;
  int mineta_;
  int maxphi_;
  int minphi_;
  double cut1_;
  double cut2_;
  double cut3_;
  int elecclass_;
  double cutEPCalo1_;
  double cutEPCalo2_;
  double cutEPin1_;
  double cutEPin2_;
  double cutCalo1_;
  double cutCalo2_;
  double cutESeed_;
  int ClusterSize_;

  CalibrationCluster calibCluster;
  CalibrationCluster::CalibMap ReducedMap;

  int read_events;

  int calibClusterSize, keventweight;
  int etaMin, etaMax, phiMin, phiMax;
  std::vector<float> EnergyVector;
  std::vector<float> WeightVector;
  std::vector<std::vector<float> > EventMatrix;
  std::vector<int> MaxCCeta;
  std::vector<int> MaxCCphi;
  std::vector<float> EnergyVectorNoCuts;
  std::vector<float> WeightVectorNoCuts;
  std::vector<std::vector<float> > EventMatrixNoCuts;
  std::vector<int> MaxCCetaNoCuts;
  std::vector<int> MaxCCphiNoCuts;
  MinL3Algorithm* MyL3Algo1;
  std::vector<float> solution;
  std::vector<float> solutionNoCuts;
  std::vector<float> newCalibs;
  std::vector<float> oldCalibs;

  int eventcrystal[171][360];
  int numevent_;

  TFile* f;

  TH1F* e25;
  TH1F* e9;
  TH1F* scE;
  TH1F* trP;
  TH1F* EoP;
  TH1F* EoP_all;
  TH1F* calibs;
  TH1F* e9Overe25;
  TH1F* e25OverScE;
  TH2F* Map;
  TH1F* E25oP;

  TH1F* PinOverPout;
  TH1F* eSeedOverPout;
  TH1F* MisCalibs;
  TH1F* RatioCalibs;
  TH1F* DiffCalibs;
  TH1F* RatioCalibsNoCuts;
  TH1F* DiffCalibsNoCuts;

  TH1F* e25NoCuts;
  TH1F* e9NoCuts;
  TH1F* scENoCuts;
  TH1F* trPNoCuts;
  TH1F* EoPNoCuts;
  TH1F* calibsNoCuts;
  TH1F* e9Overe25NoCuts;
  TH1F* e25OverScENoCuts;
  TH2F* MapNoCuts;
  TH1F* E25oPNoCuts;
  TH2F* Map3Dcalib;
  TH2F* Map3DcalibNoCuts;
  TH1F* calibinter;
  TH1F* calibinterNoCuts;
  HouseholderDecomposition* MyHH;
  TH1F* PinOverPoutNoCuts;
  TH1F* eSeedOverPoutNoCuts;

  TH2F* GeneralMap;

  TH2F* MapCor1;
  TH2F* MapCor2;
  TH2F* MapCor3;
  TH2F* MapCor4;
  TH2F* MapCor5;
  TH2F* MapCor6;
  TH2F* MapCor7;
  TH2F* MapCor8;
  TH2F* MapCor9;
  TH2F* MapCor10;
  TH2F* MapCor11;
  TH2F* MapCorCalib;

  TH2F* MapCor1NoCuts;
  TH2F* MapCor2NoCuts;
  TH2F* MapCor3NoCuts;
  TH2F* MapCor4NoCuts;
  TH2F* MapCor5NoCuts;
  TH2F* MapCor6NoCuts;
  TH2F* MapCor7NoCuts;
  TH2F* MapCor8NoCuts;
  TH2F* MapCor9NoCuts;
  TH2F* MapCor10NoCuts;
  TH2F* MapCor11NoCuts;
  TH2F* MapCorCalibNoCuts;

  TH2F* MapCor1ESeed;
  TH2F* MapCor2ESeed;
  TH2F* MapCor3ESeed;
  TH2F* MapCor4ESeed;
  TH2F* MapCor5ESeed;
  TH2F* MapCor6ESeed;
  TH2F* MapCor7ESeed;
  TH2F* MapCor8ESeed;
  TH2F* MapCor9ESeed;
  TH2F* MapCor10ESeed;
  TH2F* MapCor11ESeed;

  TH1F* PinMinPout;
  TH1F* PinMinPoutNoCuts;

  TH1F* Error1;
  TH1F* Error2;
  TH1F* Error3;
  TH1F* Error1NoCuts;
  TH1F* Error2NoCuts;
  TH1F* Error3NoCuts;

  TH1F* eSeedOverPout2;
  TH1F* eSeedOverPout2NoCuts;
  TH1F* eSeedOverPout2ESeed;

  TH1F* hadOverEm;
  TH1F* hadOverEmNoCuts;
  TH1F* hadOverEmESeed;
};
#endif