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/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 "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 decleration
//

class ElectronCalibration : public edm::EDAnalyzer {
   public:
      explicit ElectronCalibration(const edm::ParameterSet&);
      ~ElectronCalibration();

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

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

      // ----------member data ---------------------------
      std::string rootfile_;
      edm::InputTag recHitLabel_;
      edm::InputTag electronLabel_;
      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