HGCalTriggerGeometryImp1 Class Reference

Inheritance diagram for HGCalTriggerGeometryImp1:

Public Member Functions
	HGCalTriggerGeometryImp1 (const edm::ParameterSet &conf)
void	initialize (const edm::ESHandle< CaloGeometry > &) final
void	initialize (const edm::ESHandle< HGCalGeometry > &, const edm::ESHandle< HGCalGeometry > &, const edm::ESHandle< HGCalGeometry > &) final
Public Member Functions inherited from HGCalTriggerGeometryGenericMapping
bool	disconnectedModule (const unsigned module_id) const final
geom_set	getCellsFromModule (const unsigned cell_det_id) const final
geom_set	getCellsFromTriggerCell (const unsigned cell_det_id) const final
unsigned	getModuleFromCell (const unsigned cell_det_id) const final
unsigned	getModuleFromTriggerCell (const unsigned trigger_cell_det_id) const final
GlobalPoint	getModulePosition (const unsigned module_det_id) const final
geom_set	getNeighborsFromTriggerCell (const unsigned trigger_cell_det_id) const final
geom_ordered_set	getOrderedCellsFromModule (const unsigned cell_det_id) const final
geom_ordered_set	getOrderedTriggerCellsFromModule (const unsigned trigger_cell_det_id) const final
unsigned	getTriggerCellFromCell (const unsigned cell_det_id) const final
GlobalPoint	getTriggerCellPosition (const unsigned trigger_cell_det_id) const final
geom_set	getTriggerCellsFromModule (const unsigned trigger_cell_det_id) const final
	HGCalTriggerGeometryGenericMapping (const edm::ParameterSet &conf)
void	reset () final
unsigned	triggerLayer (const unsigned id) const final
bool	validTriggerCell (const unsigned trigger_cell_det_id) const final
	~HGCalTriggerGeometryGenericMapping () override
Public Member Functions inherited from HGCalTriggerGeometryBase
const HcalGeometry *	bhGeometry () const
const HcalTopology &	bhTopology () const
const edm::ESHandle< CaloGeometry > &	caloGeometry () const
const HGCalGeometry *	eeGeometry () const
const HGCalTopology &	eeTopology () const
const HGCalGeometry *	fhGeometry () const
const HGCalTopology &	fhTopology () const
	HGCalTriggerGeometryBase (const edm::ParameterSet &conf)
const HGCalGeometry *	hscGeometry () const
const HGCalTopology &	hscTopology () const
const HGCalGeometry *	hsiGeometry () const
const HGCalTopology &	hsiTopology () const
bool	isV9Geometry () const
const std::string &	name () const
virtual	~HGCalTriggerGeometryBase ()

Private Member Functions
void	buildMaps ()

Private Attributes
edm::FileInPath	l1tCellsMapping_

Additional Inherited Members
Public Types inherited from HGCalTriggerGeometryGenericMapping
typedef std::unordered_map< unsigned, std::unique_ptr< const HGCalTriggerGeometry::Module > >	module_map
typedef std::unordered_map< unsigned, std::unique_ptr< const HGCalTriggerGeometry::TriggerCell > >	trigger_cell_map
Public Types inherited from HGCalTriggerGeometryBase
typedef std::unordered_map< unsigned, unsigned >	geom_map
typedef std::set< unsigned >	geom_ordered_set
typedef std::unordered_set< unsigned >	geom_set
Protected Member Functions inherited from HGCalTriggerGeometryBase
void	setCaloGeometry (const edm::ESHandle< CaloGeometry > &geom)
void	setEEGeometry (const edm::ESHandle< HGCalGeometry > &geom)
void	setHScGeometry (const edm::ESHandle< HGCalGeometry > &geom)
void	setHSiGeometry (const edm::ESHandle< HGCalGeometry > &geom)
Protected Attributes inherited from HGCalTriggerGeometryGenericMapping
geom_map	cells_to_trigger_cells_
module_map	modules_
trigger_cell_map	trigger_cells_
geom_map	trigger_cells_to_modules_

Detailed Description

Definition at line 11 of file HGCalTriggerGeometryImp1.cc.

Constructor & Destructor Documentation

HGCalTriggerGeometryImp1::HGCalTriggerGeometryImp1

(

const edm::ParameterSet &

conf

)

Definition at line 27 of file HGCalTriggerGeometryImp1.cc.

    : HGCalTriggerGeometryGenericMapping(conf),
      l1tCellsMapping_(conf.getParameter<edm::FileInPath>("L1TCellsMapping"))
/*****************************************************************/
{}

Member Function Documentation

void HGCalTriggerGeometryImp1::buildMaps ( )

private

Definition at line 55 of file HGCalTriggerGeometryImp1.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), HGCalTriggerGeometryGenericMapping::cells_to_trigger_cells_, DEFINE_EDM_PLUGIN, HGCalTriggerGeometryBase::eeGeometry(), edm::FileInPath::fullPath(), HGCalGeometry::getPosition(), HGCEE, HGCTrigger, l1tCellsMapping_, HGCalTriggerGeometryGenericMapping::modules_, eostools::move(), HGCalTriggerGeometryGenericMapping::trigger_cells_, HGCalTriggerGeometryGenericMapping::trigger_cells_to_modules_, HGCTriggerDetId::UndefinedCell(), z, and ecaldqm::zside().

{
  //
  // read trigger cell mapping file
  std::ifstream l1tCellsMappingStream(l1tCellsMapping_.fullPath());
  if (!l1tCellsMappingStream.is_open())
    edm::LogError("HGCalTriggerGeometry") << "Cannot open L1TCellsMapping file\n";
  short layer = 0;
  short subsector = 0;
  short cell = 0;
  short module = 0;
  short triggercell = 0;
  for (; l1tCellsMappingStream >> layer >> subsector >> cell >> module >> triggercell;) {
    if (layer > 30 || layer <= 0) {
      edm::LogWarning("HGCalTriggerGeometry") << "Bad layer index in L1TCellsMapping\n";
      continue;
    }
    // Loop on all sectors
    // FIXME:  Number of sectors in each zside should not be hardcoded
    for (unsigned z = 0; z <= 1; z++) {
      int zside = (z == 0 ? -1 : 1);
      for (unsigned sector = 1; sector <= 18; sector++) {
        HGCEEDetId detid(HGCEE, zside, layer, sector, subsector, cell);
        //
        // Fill cell -> trigger cell mapping
        HGCTriggerDetId triggerDetid(HGCTrigger, zside, layer, sector, module, triggercell);
        const auto& ret = cells_to_trigger_cells_.insert(std::make_pair(detid, triggerDetid));
        if (!ret.second)
          edm::LogWarning("HGCalTriggerGeometry") << "Duplicate cell in L1TCellsMapping\n";
        // Fill trigger cell -> module mapping
        HGCTriggerDetId moduleDetid(HGCTrigger, zside, layer, sector, module, HGCTriggerDetId::UndefinedCell());
        trigger_cells_to_modules_.insert(
            std::make_pair(triggerDetid, moduleDetid));  // do nothing if trigger cell has already been inserted
      }
    }
  }
  if (!l1tCellsMappingStream.eof())
    edm::LogWarning("HGCalTriggerGeometry")
        << "Error reading L1TCellsMapping'" << layer << " " << cell << " " << triggercell << " " << subsector << "' \n";
  l1tCellsMappingStream.close();
  //
  // Build trigger cells and fill map
  typedef HGCalTriggerGeometry::TriggerCell::list_type list_cells;
  // make list of cells in trigger cells
  std::map<unsigned, list_cells> trigger_cells_to_cells;
  for (const auto& cell_triggercell : cells_to_trigger_cells_) {
    unsigned cell = cell_triggercell.first;
    unsigned triggercell = cell_triggercell.second;
    trigger_cells_to_cells.insert(std::make_pair(triggercell, list_cells()));
    trigger_cells_to_cells.at(triggercell).insert(cell);
  }
  for (const auto& triggercell_cells : trigger_cells_to_cells) {
    unsigned triggercellId = triggercell_cells.first;
    list_cells cellIds = triggercell_cells.second;
    // Position: for the moment, barycenter of the trigger cell.
    Basic3DVector<float> triggercellVector(0., 0., 0.);
    for (const auto& cell : cellIds) {
      HGCTriggerDetId cellId(cell);
      triggercellVector += eeGeometry()->getPosition(cellId).basicVector();
    }
    GlobalPoint triggercellPoint(triggercellVector / cellIds.size());
    const auto& tc2mItr = trigger_cells_to_modules_.find(triggercellId);
    unsigned moduleId =
        (tc2mItr != trigger_cells_to_modules_.end() ? tc2mItr->second
                                                    : 0);  // 0 if the trigger cell doesn't belong to a module
    //unsigned moduleId = trigger_cells_to_modules_.at(triggercellId);
    // FIXME: empty neighbours
    std::unique_ptr<const HGCalTriggerGeometry::TriggerCell> triggercellPtr(
        new HGCalTriggerGeometry::TriggerCell(triggercellId, moduleId, triggercellPoint, list_cells(), cellIds));
    trigger_cells_.insert(std::make_pair(triggercellId, std::move(triggercellPtr)));
  }
  //
  // Build modules and fill map
  typedef HGCalTriggerGeometry::Module::list_type list_triggercells;
  typedef HGCalTriggerGeometry::Module::tc_map_type tc_map_to_cells;
  // make list of trigger cells in modules
  std::map<unsigned, list_triggercells> modules_to_trigger_cells;
  for (const auto& triggercell_module : trigger_cells_to_modules_) {
    unsigned triggercell = triggercell_module.first;
    unsigned module = triggercell_module.second;
    modules_to_trigger_cells.insert(std::make_pair(module, list_triggercells()));
    modules_to_trigger_cells.at(module).insert(triggercell);
  }
  for (const auto& module_triggercell : modules_to_trigger_cells) {
    unsigned moduleId = module_triggercell.first;
    list_triggercells triggercellIds = module_triggercell.second;
    tc_map_to_cells cellsInTriggerCells;
    // Position: for the moment, barycenter of the module, from trigger cell positions
    Basic3DVector<float> moduleVector(0., 0., 0.);
    for (const auto& triggercell : triggercellIds) {
      const auto& cells_in_tc = trigger_cells_to_cells[triggercell];
      for (const unsigned cell : cells_in_tc) {
        cellsInTriggerCells.emplace(triggercell, cell);
      }
      moduleVector += trigger_cells_.at(triggercell)->position().basicVector();
    }
    GlobalPoint modulePoint(moduleVector / triggercellIds.size());
    // FIXME: empty neighbours
    std::unique_ptr<const HGCalTriggerGeometry::Module> modulePtr(new HGCalTriggerGeometry::Module(
        moduleId, modulePoint, list_triggercells(), triggercellIds, cellsInTriggerCells));
    modules_.insert(std::make_pair(moduleId, std::move(modulePtr)));
  }
}

void HGCalTriggerGeometryImp1::initialize ( const edm::ESHandle< CaloGeometry > &  calo_geometry )

finalvirtual

Implements HGCalTriggerGeometryBase.

Definition at line 34 of file HGCalTriggerGeometryImp1.cc.

References HGCalTriggerGeometryBase::setCaloGeometry().

{
  // FIXME: !!!Only for HGCEE for the moment!!!
  edm::LogWarning("HGCalTriggerGeometry") << "WARNING: This HGCal trigger geometry is incomplete.\n"
                                          << "WARNING: Only the EE part is covered.\n"
                                          << "WARNING: There is no neighbor information.\n";
  setCaloGeometry(calo_geometry);
}

void HGCalTriggerGeometryImp1::initialize ( const edm::ESHandle< HGCalGeometry > &  hgc_ee_geometry, const edm::ESHandle< HGCalGeometry > &  hgc_hsi_geometry, const edm::ESHandle< HGCalGeometry > &  hgc_hsc_geometry  )

finalvirtual

Implements HGCalTriggerGeometryBase.

Definition at line 45 of file HGCalTriggerGeometryImp1.cc.

References Exception.

{
  throw cms::Exception("BadGeometry")
      << "HGCalTriggerGeometryImp1 geometry cannot be initialized with the V9 HGCAL geometry";
}

Member Data Documentation

edm::FileInPath HGCalTriggerGeometryImp1::l1tCellsMapping_

private

Definition at line 21 of file HGCalTriggerGeometryImp1.cc.

Referenced by buildMaps().