SimG4HGCalValidation Class Reference

Inheritance diagram for SimG4HGCalValidation:

Public Member Functions
void	beginRun (edm::EventSetup const &) override
void	produce (edm::Event &, const edm::EventSetup &) override
void	registerConsumes (edm::ConsumesCollector) override
	SimG4HGCalValidation (const edm::ParameterSet &p)
	~SimG4HGCalValidation () override
Public Member Functions inherited from SimProducer
const SimProducer &	operator= (const SimProducer &)=delete
void	registerProducts (edm::ProducesCollector producesCollector)
	SimProducer ()
	SimProducer (const SimProducer &)=delete
Public Member Functions inherited from SimWatcher
bool	isMT () const
const SimWatcher &	operator= (const SimWatcher &)=delete
	SimWatcher ()
	SimWatcher (const SimWatcher &)=delete
virtual	~SimWatcher ()
Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const G4Step *>
	Observer ()
void	slotForUpdate (const G4Step * iT)
virtual	~Observer ()

Private Member Functions
void	clear ()
void	init ()
void	layerAnalysis (PHGCalValidInfo &)
const SimG4HGCalValidation &	operator= (const SimG4HGCalValidation &)
	SimG4HGCalValidation (const SimG4HGCalValidation &)
void	update (const BeginOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const G4Step *step) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
unsigned int	count_
std::vector< edm::ESGetToken< HGCalDDDConstants, IdealGeometryRecord > >	ddconsToken_
std::vector< unsigned int >	dets_
std::vector< int >	detTypes_
double	edepEE_
double	edepHEB_
double	edepHEF_
std::vector< HGCalNumberingScheme * >	hgcalNumbering_
std::vector< double >	hgcEEedep_
std::vector< double >	hgcHEBedep_
std::vector< double >	hgcHEFedep_
std::vector< unsigned int >	hgchitDets_
std::vector< unsigned int >	hgchitIndex_
std::vector< double >	hgchitX_
std::vector< double >	hgchitY_
std::vector< double >	hgchitZ_
std::vector< HGCNumberingScheme * >	hgcNumbering_
std::string	labelLayer_
int	levelT1_
int	levelT2_
std::vector< std::string >	names_
std::vector< std::string >	nameXs_
std::vector< int >	subdet_
std::vector< int >	types_
int	verbosity_

Additional Inherited Members
Protected Member Functions inherited from SimProducer
template<class T >
void	produces ()
template<class T >
void	produces (const std::string &instanceName)
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 51 of file SimG4HGCalValidation.cc.

Constructor & Destructor Documentation

SimG4HGCalValidation() [1/2]

SimG4HGCalValidation::SimG4HGCalValidation

(

const edm::ParameterSet &

p

)

Definition at line 98 of file SimG4HGCalValidation.cc.

References detTypes_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), dqmdumpme::k, labelLayer_, names_, AlCaHLTBitMon_ParallelJobs::p, AlCaHLTBitMon_QueryRunRegistry::string, types_, and verbosity_.

                                                                   : levelT1_(999), levelT2_(999), count_(0) {
  edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("SimG4HGCalValidation");
  names_ = m_Anal.getParameter<std::vector<std::string> >("Names");
  types_ = m_Anal.getParameter<std::vector<int> >("Types");
  detTypes_ = m_Anal.getParameter<std::vector<int> >("DetTypes");
  labelLayer_ = m_Anal.getParameter<std::string>("LabelLayerInfo");
  verbosity_ = m_Anal.getUntrackedParameter<int>("Verbosity", 0);

  produces<PHGCalValidInfo>(labelLayer_);

  if (verbosity_ > 0) {
    edm::LogVerbatim("ValidHGCal") << "HGCalTestAnalysis:: Initialised as "
                                   << "observer of begin events and of G4step "
                                   << "with Parameter values: \n\tLabel : " << labelLayer_ << " and with "
                                   << names_.size() << " detectors";
    for (unsigned int k = 0; k < names_.size(); ++k)
      edm::LogVerbatim("ValidHGCal") << " [" << k << "] " << names_[k] << " Type " << types_[k] << " DetType "
                                     << detTypes_[k];
  }
}

~SimG4HGCalValidation()

SimG4HGCalValidation::~SimG4HGCalValidation ( )

override

Definition at line 119 of file SimG4HGCalValidation.cc.

References hgcalNumbering_, hgcNumbering_, and contentValuesFiles::number.

                                            {
  for (auto number : hgcNumbering_)
    delete number;
  for (auto number : hgcalNumbering_)
    delete number;
}

SimG4HGCalValidation() [2/2]

SimG4HGCalValidation::SimG4HGCalValidation ( const SimG4HGCalValidation &  )

private

Member Function Documentation

beginRun()

void SimG4HGCalValidation::beginRun ( edm::EventSetup const &  es )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 172 of file SimG4HGCalValidation.cc.

References ddconsToken_, dets_, detTypes_, Exception, edm::EventSetup::getHandle(), hgcalNumbering_, hgcEEedep_, hgcHEBedep_, hgcHEFedep_, hgcNumbering_, mps_fire::i, edm::ESHandleBase::isValid(), hgcalTopologyTester_cfi::layers, HGCalDDDConstants::layers(), levelT1_, levelT2_, HGCalDDDConstants::levelTop(), nameXs_, subdet_, HGCalDDDConstants::tileTrapezoid(), types_, verbosity_, and HGCalDDDConstants::waferHexagon6().

                                                           {
  for (unsigned int type = 0; type < types_.size(); ++type) {
    int layers(0);
    int detType = detTypes_[type];
    edm::ESHandle<HGCalDDDConstants> hdc = es.getHandle(ddconsToken_[type]);
    if (hdc.isValid()) {
      levelT1_ = hdc->levelTop(0);
      levelT2_ = hdc->levelTop(1);
      if (hdc->tileTrapezoid()) {
        types_[type] = -1;
        hgcalNumbering_.emplace_back(new HGCalNumberingScheme(*hdc, (DetId::Detector)(dets_[type]), nameXs_[type], ""));
      } else if (hdc->waferHexagon6()) {
        types_[type] = 1;
        hgcNumbering_.emplace_back(new HGCNumberingScheme(*hdc, nameXs_[type]));
      } else {
        types_[type] = 0;
        hgcalNumbering_.emplace_back(new HGCalNumberingScheme(*hdc, (DetId::Detector)(dets_[type]), nameXs_[type], ""));
      }
      layers = hdc->layers(false);
    } else {
      edm::LogError("ValidHGCal") << "Cannot find HGCalDDDConstants for " << nameXs_[type];
      throw cms::Exception("Unknown", "ValidHGCal") << "Cannot find HGCalDDDConstants for " << nameXs_[type] << "\n";
    }
    if (detType == 0) {
      for (int i = 0; i < layers; ++i)
        hgcEEedep_.emplace_back(0);
    } else if (detType == 1) {
      for (int i = 0; i < layers; ++i)
        hgcHEFedep_.emplace_back(0);
    } else {
      for (int i = 0; i < layers; ++i)
        hgcHEBedep_.emplace_back(0);
    }
    if (verbosity_ > 0)
      edm::LogVerbatim("ValidHGCal") << "[" << type << "]: " << nameXs_[type] << " det " << dets_[type] << " subdet "
                                     << subdet_[type] << " with " << layers << " layers";
  }
}

clear()

void SimG4HGCalValidation::clear ( void  )

private

Definition at line 344 of file SimG4HGCalValidation.cc.

References hgchitDets_, hgchitIndex_, hgchitX_, hgchitY_, and hgchitZ_.

Referenced by update().

                                 {
  hgchitDets_.erase(hgchitDets_.begin(), hgchitDets_.end());
  hgchitIndex_.erase(hgchitIndex_.begin(), hgchitIndex_.end());
  hgchitX_.erase(hgchitX_.begin(), hgchitX_.end());
  hgchitY_.erase(hgchitY_.begin(), hgchitY_.end());
  hgchitZ_.erase(hgchitZ_.begin(), hgchitZ_.end());
}

init()

void SimG4HGCalValidation::init ( )

private

layerAnalysis()

void SimG4HGCalValidation::layerAnalysis ( PHGCalValidInfo &  product )

private

Definition at line 331 of file SimG4HGCalValidation.cc.

References edepEE_, edepHEB_, edepHEF_, PHGCalValidInfo::fillhgcHits(), PHGCalValidInfo::fillhgcLayers(), hgcEEedep_, hgcHEBedep_, hgcHEFedep_, hgchitDets_, hgchitIndex_, hgchitX_, hgchitY_, hgchitZ_, and verbosity_.

Referenced by produce().

                                                                 {
  if (verbosity_ > 0)
    edm::LogVerbatim("ValidHGCal") << "\n ===>>> SimG4HGCalValidation: Energy "
                                   << "deposit\n at EE : " << std::setw(6) << edepEE_ / CLHEP::MeV
                                   << "\n at HEF: " << std::setw(6) << edepHEF_ / CLHEP::MeV
                                   << "\n at HEB: " << std::setw(6) << edepHEB_ / CLHEP::MeV;

  //Fill HGC Variables
  product.fillhgcHits(hgchitDets_, hgchitIndex_, hgchitX_, hgchitY_, hgchitZ_);
  product.fillhgcLayers(edepEE_, edepHEF_, edepHEB_, hgcEEedep_, hgcHEFedep_, hgcHEBedep_);
}

operator=()

const SimG4HGCalValidation& SimG4HGCalValidation::operator= ( const SimG4HGCalValidation &  )

private

produce()

void SimG4HGCalValidation::produce ( edm::Event &  e, const edm::EventSetup &    )

overridevirtual

Implements SimProducer.

Definition at line 166 of file SimG4HGCalValidation.cc.

References MillePedeFileConverter_cfg::e, labelLayer_, layerAnalysis(), and eostools::move().

                                                                    {
  std::unique_ptr<PHGCalValidInfo> productLayer(new PHGCalValidInfo);
  layerAnalysis(*productLayer);
  e.put(std::move(productLayer), labelLayer_);
}

registerConsumes()

void SimG4HGCalValidation::registerConsumes ( edm::ConsumesCollector  cc )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 126 of file SimG4HGCalValidation.cc.

References edm::BeginRun, ddconsToken_, dets_, detTypes_, edm::ConsumesCollector::esConsumes(), Exception, DetId::Forward, ForwardEmpty, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCEE, HGCHEB, HGCHEF, hfnoseParametersInitialization_cfi::nameX, nameXs_, subdet_, and types_.

                                                                   {
  for (unsigned int type = 0; type < types_.size(); ++type) {
    int detType = detTypes_[type];
    G4String nameX = "HGCal";
    if (types_[type] <= 1) {
      if (types_[type] == 0) {
        subdet_.emplace_back(ForwardSubdetector::ForwardEmpty);
        if (detType == 0)
          dets_.emplace_back((int)(DetId::HGCalEE));
        else if (detType == 1)
          dets_.emplace_back((int)(DetId::HGCalHSi));
        else
          dets_.emplace_back((int)(DetId::HGCalHSc));
      } else {
        dets_.emplace_back((unsigned int)(DetId::Forward));
        if (detType == 0)
          subdet_.emplace_back((int)(ForwardSubdetector::HGCEE));
        else if (detType == 1)
          subdet_.emplace_back((int)(ForwardSubdetector::HGCHEF));
        else
          subdet_.emplace_back((int)(ForwardSubdetector::HGCHEB));
      }
      if (detType == 0)
        nameX = "HGCalEESensitive";
      else if (detType == 1)
        nameX = "HGCalHESiliconSensitive";
      else
        nameX = "HGCalHEScintillatorSensitive";
      nameXs_.emplace_back(nameX);
      ddconsToken_.emplace_back(
          cc.esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(edm::ESInputTag{"", nameX}));
      edm::LogVerbatim("ValidHGCal") << "Defines ESGetToken for type " << type << ":" << dets_.back() << ":"
                                     << subdet_.back() << ":" << nameX;
    } else {
      edm::LogError("ValidHGCal") << "Wrong Type " << types_[type];
      throw cms::Exception("Unknown", "ValidHGCal") << "Wrong Type " << types_[type] << "\n";
    }
  }
}

update() [1/2]

void SimG4HGCalValidation::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 212 of file SimG4HGCalValidation.cc.

References clear(), count_, edepEE_, edepHEB_, edepHEF_, hgcEEedep_, hgcHEBedep_, hgcHEFedep_, mps_fire::i, iev, and verbosity_.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                         {
  int iev = (*evt)()->GetEventID();
  if (verbosity_ > 0)
    edm::LogVerbatim("ValidHGCal") << "SimG4HGCalValidation: =====> Begin "
                                   << "event = " << iev;

  ++count_;
  edepEE_ = edepHEF_ = edepHEB_ = 0.;

  //HGCal variables
  for (unsigned int i = 0; i < hgcEEedep_.size(); i++)
    hgcEEedep_[i] = 0.;
  for (unsigned int i = 0; i < hgcHEFedep_.size(); i++)
    hgcHEFedep_[i] = 0.;
  for (unsigned int i = 0; i < hgcHEBedep_.size(); i++)
    hgcHEBedep_[i] = 0.;

  //Cache reset
  clear();
}

update() [2/2]

void SimG4HGCalValidation::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 234 of file SimG4HGCalValidation.cc.

References TauDecayModes::dec, dets_, detTypes_, edepEE_, edepHEB_, edepHEF_, spr::find(), hgcalNumbering_, hgcEEedep_, hgcHEBedep_, hgcHEFedep_, hgchitDets_, hgchitIndex_, hgchitX_, hgchitY_, hgchitZ_, hgcNumbering_, dqmdumpme::k, phase1PixelTopology::layer, levelT1_, levelT2_, callgraph::module, Skims_PA_cff::name, names_, subdet_, types_, verbosity_, and mps_merge::weight.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                     {
  if (aStep != nullptr) {
    G4VPhysicalVolume* curPV = aStep->GetPreStepPoint()->GetPhysicalVolume();
    const G4String& name = curPV->GetName();
    G4VSensitiveDetector* curSD = aStep->GetPreStepPoint()->GetSensitiveDetector();
    if (verbosity_ > 1)
      edm::LogVerbatim("ValidHGCal") << "ValidHGCal::Step in " << name << " at "
                                     << aStep->GetPreStepPoint()->GetPosition();

    // Only for Sensitive detector
    if (curSD != nullptr) {
      int type(-1);
      for (unsigned int k = 0; k < names_.size(); ++k) {
        if (name.find(names_[k].c_str()) != std::string::npos) {
          type = k;
          break;
        }
      }
      int detType = (type >= 0) ? detTypes_[type] : -1;
      if (verbosity_ > 0)
        edm::LogVerbatim("ValidHGCal") << "ValidHGCal::In SD " << name << " type " << type << ":" << detType;

      // Right type of SD
      if (type >= 0) {
        //Get the 32-bit index of the hit cell
        const G4TouchableHistory* touchable =
            static_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
        unsigned int index(0);
        int layer(0);
        G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
        // HGCal
        G4ThreeVector localpos = touchable->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
        float globalZ = touchable->GetTranslation(0).z();
        int iz(globalZ > 0 ? 1 : -1);
        int module(-1), cell(-1);
        if (types_[type] == 1) {
          if (touchable->GetHistoryDepth() == levelT1_) {
            layer = touchable->GetReplicaNumber(0);
          } else {
            layer = touchable->GetReplicaNumber(2);
            module = touchable->GetReplicaNumber(1);
            cell = touchable->GetReplicaNumber(0);
          }
          index =
              hgcNumbering_[type]->getUnitID((ForwardSubdetector)(subdet_[type]), layer, module, cell, iz, localpos);
        } else {
          if ((touchable->GetHistoryDepth() == levelT1_) || (touchable->GetHistoryDepth() == levelT2_)) {
            layer = touchable->GetReplicaNumber(0);
          } else {
            layer = touchable->GetReplicaNumber(3);
            module = touchable->GetReplicaNumber(2);
            cell = touchable->GetReplicaNumber(1);
          }
          double weight(0);
          index = hgcalNumbering_[type]->getUnitID(layer, module, cell, iz, hitPoint, weight);
        }
        if (verbosity_ > 1)
          edm::LogVerbatim("ValidHGCal") << "HGCal: " << name << " Layer " << layer << " Module " << module << " Cell "
                                         << cell;

        double edeposit = aStep->GetTotalEnergyDeposit();
        if (verbosity_ > 0)
          edm::LogVerbatim("ValidHGCal") << "Layer " << layer << " Index " << std::hex << index << std::dec << " Edep "
                                         << edeposit << " hit at " << hitPoint;
        if (detType == 0) {
          edepEE_ += edeposit;
          if (layer < (int)(hgcEEedep_.size()))
            hgcEEedep_[layer] += edeposit;
        } else if (detType == 1) {
          edepHEF_ += edeposit;
          if (layer < (int)(hgcHEFedep_.size()))
            hgcHEFedep_[layer] += edeposit;
        } else {
          edepHEB_ += edeposit;
          if (layer < (int)(hgcHEBedep_.size()))
            hgcHEBedep_[layer] += edeposit;
        }
        G4String nextVolume("XXX");
        if (aStep->GetTrack()->GetNextVolume() != nullptr)
          nextVolume = aStep->GetTrack()->GetNextVolume()->GetName();

        if (nextVolume.c_str() != name.c_str()) {  //save hit when it exits cell
          if (std::find(hgchitIndex_.begin(), hgchitIndex_.end(), index) == hgchitIndex_.end()) {
            hgchitDets_.emplace_back(dets_[type]);
            hgchitIndex_.emplace_back(index);
            hgchitX_.emplace_back(hitPoint.x());
            hgchitY_.emplace_back(hitPoint.y());
            hgchitZ_.emplace_back(hitPoint.z());
          }
        }
      }  // it is right type of SD
    }    // it is in a SD
  }      //if aStep!=NULL
}  //end update aStep

Member Data Documentation

count_

unsigned int SimG4HGCalValidation::count_

private

Definition at line 90 of file SimG4HGCalValidation.cc.

Referenced by update().

ddconsToken_

std::vector<edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord> > SimG4HGCalValidation::ddconsToken_

private

Definition at line 76 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), and registerConsumes().

dets_

std::vector<unsigned int> SimG4HGCalValidation::dets_

private

Definition at line 94 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), registerConsumes(), and update().

detTypes_

std::vector<int> SimG4HGCalValidation::detTypes_

private

Definition at line 82 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), registerConsumes(), SimG4HGCalValidation(), and update().

edepEE_

double SimG4HGCalValidation::edepEE_

private

Definition at line 92 of file SimG4HGCalValidation.cc.

Referenced by layerAnalysis(), and update().

edepHEB_

double SimG4HGCalValidation::edepHEB_

private

Definition at line 92 of file SimG4HGCalValidation.cc.

Referenced by layerAnalysis(), and update().

edepHEF_

double SimG4HGCalValidation::edepHEF_

private

Definition at line 92 of file SimG4HGCalValidation.cc.

Referenced by layerAnalysis(), and update().

hgcalNumbering_

std::vector<HGCalNumberingScheme*> SimG4HGCalValidation::hgcalNumbering_

private

Definition at line 78 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), update(), and ~SimG4HGCalValidation().

hgcEEedep_

std::vector<double> SimG4HGCalValidation::hgcEEedep_

private

Definition at line 93 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), layerAnalysis(), and update().

hgcHEBedep_

std::vector<double> SimG4HGCalValidation::hgcHEBedep_

private

Definition at line 93 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), layerAnalysis(), and update().

hgcHEFedep_

std::vector<double> SimG4HGCalValidation::hgcHEFedep_

private

Definition at line 93 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), layerAnalysis(), and update().

hgchitDets_

std::vector<unsigned int> SimG4HGCalValidation::hgchitDets_

private

Definition at line 94 of file SimG4HGCalValidation.cc.

Referenced by clear(), layerAnalysis(), and update().

hgchitIndex_

std::vector<unsigned int> SimG4HGCalValidation::hgchitIndex_

private

Definition at line 94 of file SimG4HGCalValidation.cc.

Referenced by clear(), layerAnalysis(), and update().

hgchitX_

std::vector<double> SimG4HGCalValidation::hgchitX_

private

Definition at line 95 of file SimG4HGCalValidation.cc.

Referenced by clear(), layerAnalysis(), and update().

hgchitY_

std::vector<double> SimG4HGCalValidation::hgchitY_

private

Definition at line 95 of file SimG4HGCalValidation.cc.

Referenced by clear(), layerAnalysis(), and update().

hgchitZ_

std::vector<double> SimG4HGCalValidation::hgchitZ_

private

Definition at line 95 of file SimG4HGCalValidation.cc.

Referenced by clear(), layerAnalysis(), and update().

hgcNumbering_

std::vector<HGCNumberingScheme*> SimG4HGCalValidation::hgcNumbering_

private

Definition at line 77 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), update(), and ~SimG4HGCalValidation().

labelLayer_

std::string SimG4HGCalValidation::labelLayer_

private

Definition at line 83 of file SimG4HGCalValidation.cc.

Referenced by produce(), and SimG4HGCalValidation().

levelT1_

int SimG4HGCalValidation::levelT1_

private

Definition at line 87 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), and update().

levelT2_

int SimG4HGCalValidation::levelT2_

private

Definition at line 87 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), and update().

names_

std::vector<std::string> SimG4HGCalValidation::names_

private

Definition at line 81 of file SimG4HGCalValidation.cc.

Referenced by SimG4HGCalValidation(), and update().

nameXs_

std::vector<std::string> SimG4HGCalValidation::nameXs_

private

Definition at line 84 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), and registerConsumes().

subdet_

std::vector<int> SimG4HGCalValidation::subdet_

private

Definition at line 82 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), registerConsumes(), and update().

types_

std::vector<int> SimG4HGCalValidation::types_

private

Definition at line 82 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), registerConsumes(), SimG4HGCalValidation(), and update().

verbosity_

int SimG4HGCalValidation::verbosity_

private

Definition at line 91 of file SimG4HGCalValidation.cc.

Referenced by beginRun(), layerAnalysis(), SimG4HGCalValidation(), and update().