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/ParameterSet/interface/ConfigurationDescriptions.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;
 
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
 
  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