SamplingFractionInLayers.h

Go to the documentation of this file.

#ifndef DETSTUDIES_SAMPLINGFRACTIONINLAYERS_H
#define DETSTUDIES_SAMPLINGFRACTIONINLAYERS_H
 
// GAUDI
#include "Gaudi/Algorithm.h"
#include "GaudiKernel/ServiceHandle.h"
 
// FCCSW
#include "k4FWCore/DataHandle.h"
class IGeoSvc;
 
// datamodel
namespace edm4hep {
class SimCalorimeterHitCollection;
}
 
class TH1F;
class ITHistSvc;
class SamplingFractionInLayers : public Gaudi::Algorithm {
public:
  explicit SamplingFractionInLayers(const std::string&, ISvcLocator*);
  virtual ~SamplingFractionInLayers();
  virtual StatusCode initialize() final;
  virtual StatusCode execute(const EventContext&) const final;
  virtual StatusCode finalize() final;
 
private:
  ServiceHandle<ITHistSvc> m_histSvc;
  ServiceHandle<IGeoSvc> m_geoSvc;
  mutable DataHandle<edm4hep::SimCalorimeterHitCollection> m_deposits{"rec/caloHits", Gaudi::DataHandle::Reader, this};
  Gaudi::Property<std::string> m_activeFieldName{this, "activeFieldName", "", "Identifier of active material"};
  Gaudi::Property<int> m_activeFieldValue{this, "activeFieldValue", 0, "Value of identifier for active material"};
  Gaudi::Property<std::string> m_layerFieldName{this, "layerFieldName", "", "Identifier of layers"};
  Gaudi::Property<uint> m_numLayers{this, "numLayers", 8, "Number of layers"};
  Gaudi::Property<uint> m_firstLayerId{this, "firstLayerId", 0, "ID of first layer"};
  Gaudi::Property<std::string> m_readoutName{this, "readoutName", "", "Name of the detector readout"};
  // Maximum energy for the axis range
  Gaudi::Property<double> m_energy{this, "energyAxis", 500, "Maximum energy for axis range"};
  // Histograms of total deposited energy within layer
  // Layers are numbered starting at 1. Layer 0 includes total energy deposited in cryostat and bath (in front and
  // behind calo)
  std::vector<TH1F*> m_totalEnLayers;
  // Histogram of total deposited energy in the calorimeter (in active and passive material, excluding cryostat and
  // bath)
  TH1F* m_totalEnergy;
  // Histograms of energy deposited in the active material within layer
  std::vector<TH1F*> m_activeEnLayers;
  // Histogram of energy deposited in the active material of the calorimeter
  TH1F* m_totalActiveEnergy;
  // Histograms of sampling fraction (active/total energy) calculated within layer
  std::vector<TH1F*> m_sfLayers;
  // Histogram of sampling fraction (active/total energy) calculated for the calorimeter (excluding cryostat and bath)
  TH1F* m_sf;
};
#endif /* DETSTUDIES_SAMPLINGFRACTIONINLAYERS_H */