Calorimeter.cc

Go to the documentation of this file.

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/EcalAlgo/interface/EcalBarrelGeometry.h"
#include "Geometry/EcalAlgo/interface/EcalEndcapGeometry.h"
#include "Geometry/EcalAlgo/interface/EcalPreshowerGeometry.h"

#include "Geometry/CaloTopology/interface/CaloTopology.h"
#include "FastSimulation/CalorimeterProperties/interface/Calorimeter.h"
#include "FastSimulation/CalorimeterProperties/interface/PreshowerLayer1Properties.h"
#include "FastSimulation/CalorimeterProperties/interface/PreshowerLayer2Properties.h"
#include "FastSimulation/CalorimeterProperties/interface/ECALBarrelProperties.h"
#include "FastSimulation/CalorimeterProperties/interface/ECALEndcapProperties.h"
#include "FastSimulation/CalorimeterProperties/interface/HCALBarrelProperties.h"
#include "FastSimulation/CalorimeterProperties/interface/HCALEndcapProperties.h"
#include "FastSimulation/CalorimeterProperties/interface/HCALForwardProperties.h"
#include "Geometry/HcalTowerAlgo/interface/HcalHardcodeGeometryLoader.h"

Calorimeter::Calorimeter():
  myPreshowerLayer1Properties_(nullptr),
  myPreshowerLayer2Properties_(nullptr),  
  myECALBarrelProperties_     (nullptr),  
  myECALEndcapProperties_     (nullptr),  
  myHCALBarrelProperties_     (nullptr),  
  myHCALEndcapProperties_     (nullptr),  
  myHCALForwardProperties_    (nullptr),
  EcalBarrelGeometry_         (nullptr),
  EcalEndcapGeometry_         (nullptr),
  HcalGeometry_               (nullptr),
  PreshowerGeometry_          (nullptr)
{
;
}

Calorimeter::Calorimeter(const edm::ParameterSet& fastCalo):
  myPreshowerLayer1Properties_(nullptr),
  myPreshowerLayer2Properties_(nullptr),  
  myECALBarrelProperties_     (nullptr),  
  myECALEndcapProperties_     (nullptr),  
  myHCALBarrelProperties_     (nullptr),  
  myHCALEndcapProperties_     (nullptr),  
  myHCALForwardProperties_    (nullptr),
  EcalBarrelGeometry_        (nullptr),
  EcalEndcapGeometry_          (nullptr),
  HcalGeometry_               (nullptr),
  PreshowerGeometry_          (nullptr)  
{
  edm::ParameterSet fastDet = fastCalo.getParameter<edm::ParameterSet>("CalorimeterProperties");
  edm::ParameterSet fastDetHF = fastCalo.getParameter<edm::ParameterSet>("ForwardCalorimeterProperties");

  myPreshowerLayer1Properties_  = new PreshowerLayer1Properties(fastDet); 
  myPreshowerLayer2Properties_  = new PreshowerLayer2Properties(fastDet);
  myECALBarrelProperties_       = new ECALBarrelProperties     (fastDet);
  myECALEndcapProperties_       = new ECALEndcapProperties     (fastDet);
  myHCALBarrelProperties_       = new HCALBarrelProperties     (fastDet);
  myHCALEndcapProperties_       = new HCALEndcapProperties     (fastDet);
  myHCALForwardProperties_      = new HCALForwardProperties    (fastDetHF);

}

Calorimeter::~Calorimeter()
{
  if(myPreshowerLayer1Properties_        )  delete myPreshowerLayer1Properties_      ;
  if(myPreshowerLayer2Properties_        )  delete myPreshowerLayer2Properties_      ;
  if(myECALBarrelProperties_             )  delete myECALBarrelProperties_           ;
  if(myECALEndcapProperties_             )  delete myECALEndcapProperties_           ;
  if(myHCALBarrelProperties_             )  delete myHCALBarrelProperties_           ;
  if(myHCALEndcapProperties_             )  delete myHCALEndcapProperties_           ;
  if(myHCALForwardProperties_            )  delete myHCALForwardProperties_          ;
}

const ECALProperties*
Calorimeter::ecalProperties(int onEcal) const {
  if ( onEcal ) {
    if ( onEcal == 1 ) 
      return myECALBarrelProperties_;
    else
      return myECALEndcapProperties_;
  } else
    return nullptr;
}

const HCALProperties*
Calorimeter::hcalProperties(int onHcal) const {
  if ( onHcal ) {
    if ( onHcal == 1 ) 
      return myHCALBarrelProperties_;
    else 
      if ( onHcal == 2 ) 
    return myHCALEndcapProperties_;
      else {
    return myHCALForwardProperties_;
    edm::LogInfo("CalorimeterProperties") << " Calorimeter::hcalProperties : set myHCALForwardProperties" << std::endl;
      }
  } else
    return nullptr;
}

const PreshowerLayer1Properties*
Calorimeter::layer1Properties(int onLayer1) const {
  if ( onLayer1 ) 
    return myPreshowerLayer1Properties_;
  else
    return nullptr;
}

const PreshowerLayer2Properties*
Calorimeter::layer2Properties(int onLayer2) const {
  if ( onLayer2 ) 
    return myPreshowerLayer2Properties_;
  else
    return nullptr;
}

void Calorimeter::setupGeometry(const CaloGeometry& pG)
{
  edm::LogInfo("CalorimeterProperties") << " setupGeometry " << std::endl;
  EcalBarrelGeometry_ = dynamic_cast<const EcalBarrelGeometry*>(pG.getSubdetectorGeometry(DetId::Ecal,EcalBarrel));
  EcalEndcapGeometry_ = dynamic_cast<const EcalEndcapGeometry*>(pG.getSubdetectorGeometry(DetId::Ecal,EcalEndcap));
  HcalGeometry_ = pG.getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
  // Takes a lot of time
  PreshowerGeometry_  = dynamic_cast<const EcalPreshowerGeometry*>(pG.getSubdetectorGeometry(DetId::Ecal,EcalPreshower));
}

void Calorimeter::setupTopology(const CaloTopology& theTopology)
{
  EcalBarrelTopology_ = theTopology.getSubdetectorTopology(DetId::Ecal,EcalBarrel);
  EcalEndcapTopology_ = theTopology.getSubdetectorTopology(DetId::Ecal,EcalEndcap);
}



const CaloSubdetectorGeometry * Calorimeter::getEcalGeometry(int subdetn) const
{
  if(subdetn==1) return EcalBarrelGeometry_;
  if(subdetn==2) return EcalEndcapGeometry_;
  if(subdetn==3) return PreshowerGeometry_;
  edm::LogWarning("Calorimeter") << "Requested an invalid ECAL subdetector geometry: " << subdetn << std::endl;
  return nullptr;
}

const CaloSubdetectorTopology * Calorimeter::getEcalTopology(int subdetn) const
{
  if(subdetn==1) return EcalBarrelTopology_;
  if(subdetn==2) return EcalEndcapTopology_;
  edm::LogWarning("Calorimeter") << "Requested an invalid ECAL subdetector topology: " << subdetn << std::endl;
  return nullptr;
}