Pi0FixedMassWindowCalibration.h

Go to the documentation of this file.

#ifndef Calibration_EcalCalibAlgos_Pi0FixedMassWindowCalibration_h
#define Calibration_EcalCalibAlgos_Pi0FixedMassWindowCalibration_h

#include <memory>
#include <string>
#include <iostream>

#include <vector>
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"

// Framework
#include "FWCore/Framework/interface/LooperFactory.h"
#include "FWCore/Framework/interface/ESProducerLooper.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "DataFormats/DetId/interface/DetId.h"

#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/HcalRecHit/interface/HBHERecHit.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EcalDetId/interface/ESDetId.h"
#include "DataFormats/DetId/interface/DetId.h"

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"


#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "RecoEcal/EgammaClusterAlgos/interface/IslandClusterAlgo.h"
#include "RecoEcal/EgammaCoreTools/interface/PositionCalc.h"
#include "RecoEcal/EgammaCoreTools/interface/ClusterShapeAlgo.h"

#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"


#include "TFile.h"
#include "TTree.h"
#include "TH1F.h"
#include "TF1.h"
#include "TGraph.h"
#include "TCanvas.h"

class Pi0FixedMassWindowCalibration :  public edm::ESProducerLooper
{

 public:

  Pi0FixedMassWindowCalibration( const edm::ParameterSet& iConfig );
  
  ~Pi0FixedMassWindowCalibration() override;

  virtual void produce(edm::Event&, const edm::EventSetup&) {};

  void beginOfJob() override;

  void endOfJob() override;

  void startingNewLoop( unsigned int iLoop ) override;

  Status endOfLoop( const edm::EventSetup&, unsigned int iLoop ) override;

  Status duringLoop( const edm::Event&, const edm::EventSetup& ) override;

 private:

  //  static const double PDGPi0Mass;

  int nevent;

  unsigned int theMaxLoops;   
  std::string ecalHitsProducer_;
  std::string barrelHits_;


  IslandClusterAlgo::VerbosityLevel verbosity;

  std::string barrelHitProducer_;
  std::string barrelHitCollection_;
  std::string barrelClusterCollection_;
  std::string clustershapecollectionEB_;
  std::string barrelClusterShapeAssociation_;

  PositionCalc posCalculator_; // position calculation algorithm
  ClusterShapeAlgo shapeAlgo_; // cluster shape algorithm
  IslandClusterAlgo * island_p;

  // Selection algorithm parameters
  double selePi0PtGammaOneMin_;
  double selePi0PtGammaTwoMin_;

  double selePi0DRBelt_;
  double selePi0DetaBelt_;

  double selePi0PtPi0Min_;

  double selePi0S4S9GammaOneMin_;
  double selePi0S4S9GammaTwoMin_;
  double selePi0S9S25GammaOneMin_;
  double selePi0S9S25GammaTwoMin_;

  double selePi0EtBeltIsoRatioMax_;
 
  double selePi0MinvMeanFixed_;
  double selePi0MinvSigmaFixed_;



  std::vector<DetId> barrelCells;

  // input calibration constants
  double oldCalibs_barl[85][360][2];
  double newCalibs_barl[85][360][2];
  double wxtals[85][360][2];
  double mwxtals[85][360][2];

  // steering parameters

  edm::ParameterSet theParameterSet;


  

  // map for all RecHits from ECal:
  std::map<DetId, EcalRecHit> *recHitsEB_map;

  const EcalRecHitCollection* ecalRecHitBarrelCollection;
  const EcalRecHitCollection* recalibEcalRecHitCollection;


  // root tree
  TFile* theFile;

  bool isfirstcall_;

};

#endif