Inheritance diagram for HGCalTriggerGeometryHexImp2:

Public Member Functions
virtual geom_set	getCellsFromModule (const unsigned) const overridefinal

virtual geom_set	getCellsFromTriggerCell (const unsigned) const overridefinal

virtual unsigned	getModuleFromCell (const unsigned) const overridefinal

virtual unsigned	getModuleFromTriggerCell (const unsigned) const overridefinal

virtual GlobalPoint	getModulePosition (const unsigned) const overridefinal

virtual geom_ordered_set	getOrderedCellsFromModule (const unsigned) const overridefinal

virtual geom_ordered_set	getOrderedTriggerCellsFromModule (const unsigned) const overridefinal

virtual unsigned	getTriggerCellFromCell (const unsigned) const overridefinal

virtual GlobalPoint	getTriggerCellPosition (const unsigned) const overridefinal

virtual geom_set	getTriggerCellsFromModule (const unsigned) const overridefinal

	HGCalTriggerGeometryHexImp2 (const edm::ParameterSet &conf)

virtual void	initialize (const es_info &) overridefinal

Public Member Functions inherited from HGCalTriggerGeometryBase
const std::string &	bhSDName () const

const std::string &	eeSDName () const

const std::string &	fhSDName () const

	HGCalTriggerGeometryBase (const edm::ParameterSet &conf)

virtual void	initialize (const es_info &)=0

const std::string &	name () const

virtual void	reset ()

virtual	~HGCalTriggerGeometryBase ()

Private Member Functions
void	fillMaps (const es_info &)

Private Attributes
std::map< std::pair< short, short >, short >	cells_to_trigger_cells_

es_info	es_info_

edm::FileInPath	l1tCellsMapping_

edm::FileInPath	l1tModulesMapping_

std::unordered_multimap< short, short >	module_to_wafers_ee_

std::unordered_multimap< short, short >	module_to_wafers_fh_

std::unordered_map< short, short >	number_cells_in_wafers_

std::unordered_map< short, short >	number_trigger_cells_in_wafers_

std::multimap< std::pair < short, short >, short >	trigger_cells_to_cells_

std::unordered_map< short, short >	wafer_to_module_ee_

std::unordered_map< short, short >	wafer_to_module_fh_

Additional Inherited Members
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

Detailed Description

Definition at line 12 of file HGCalTriggerGeometryHexImp2.cc.

Constructor & Destructor Documentation

HGCalTriggerGeometryHexImp2::HGCalTriggerGeometryHexImp2 ( const edm::ParameterSet & conf )

Definition at line 54 of file HGCalTriggerGeometryHexImp2.cc.

                                                         :
     HGCalTriggerGeometryBase(conf),
     l1tCellsMapping_(conf.getParameter<edm::FileInPath>("L1TCellsMapping")),
     l1tModulesMapping_(conf.getParameter<edm::FileInPath>("L1TModulesMapping"))
 {
 }

Member Function Documentation

void HGCalTriggerGeometryHexImp2::fillMaps ( const es_info & esInfo )

private

Definition at line 322 of file HGCalTriggerGeometryHexImp2.cc.

References cells_to_trigger_cells_, edm::FileInPath::fullPath(), HGCEE, HGCHEF, l1tCellsMapping_, l1tModulesMapping_, python.rootplot.argparse::module, module_to_wafers_ee_, module_to_wafers_fh_, number_cells_in_wafers_, number_trigger_cells_in_wafers_, trigger_cells_to_cells_, wafer_to_module_ee_, and wafer_to_module_fh_.

Referenced by initialize().

 {
     //
     // read module mapping file
     std::ifstream l1tModulesMappingStream(l1tModulesMapping_.fullPath());
     if(!l1tModulesMappingStream.is_open()) edm::LogError("HGCalTriggerGeometry") << "Cannot open L1TModulesMapping file\n";
     short subdet  = 0;
     short wafer   = 0;
     short module  = 0;
     for(; l1tModulesMappingStream>>subdet>>wafer>>module; )
     {
         switch(subdet)
         {
             case ForwardSubdetector::HGCEE:
             {
                 // fill module <-> wafers mappings
                 wafer_to_module_ee_.emplace(wafer,module);
                 module_to_wafers_ee_.emplace(module, wafer);
                 // fill number of cells for a given wafer type
                 // translate wafer type 1/2 to 1/-1 
                 int wafer_type = esInfo.topo_ee->dddConstants().waferTypeT(wafer)==1?1:-1;
                 number_cells_in_wafers_.emplace(wafer_type, esInfo.topo_ee->dddConstants().numberCellsHexagon(wafer));
                 break;
             }
             case ForwardSubdetector::HGCHEF:
             {
                 // fill module <-> wafers mappings
                 wafer_to_module_fh_.emplace(wafer,module);
                 module_to_wafers_fh_.emplace(module, wafer);
                 // fill number of cells for a given wafer type
                 // translate wafer type 1/2 to 1/-1
                 int wafer_type = esInfo.topo_fh->dddConstants().waferTypeT(wafer)==1?1:-1;
                 number_cells_in_wafers_.emplace(wafer_type, esInfo.topo_fh->dddConstants().numberCellsHexagon(wafer));
                 break;
             }
             default:
                 edm::LogWarning("HGCalTriggerGeometry") << "Unsupported subdetector number ("<<subdet<<") in L1TModulesMapping file\n";
                 break;
         }
     }
     if(!l1tModulesMappingStream.eof()) edm::LogWarning("HGCalTriggerGeometry") << "Error reading L1TModulesMapping '"<<wafer<<" "<<module<<"' \n";
     l1tModulesMappingStream.close();
     // read trigger cell mapping file
     std::ifstream l1tCellsMappingStream(l1tCellsMapping_.fullPath());
     if(!l1tCellsMappingStream.is_open()) edm::LogError("HGCalTriggerGeometry") << "Cannot open L1TCellsMapping file\n";
     short waferType   = 0;
     short cell        = 0;
     short triggerCell = 0;
     for(; l1tCellsMappingStream>>waferType>>cell>>triggerCell; )
     {
         // fill cell <-> trigger cell mappings
         cells_to_trigger_cells_.emplace(std::make_pair((waferType?1:-1),cell), triggerCell);
         trigger_cells_to_cells_.emplace(std::make_pair((waferType?1:-1),triggerCell), cell);
         // fill number of cells for a given wafer type
         auto itr_insert = number_trigger_cells_in_wafers_.emplace((waferType?1:-1), 0);
         if(triggerCell+1 > itr_insert.first->second) itr_insert.first->second = triggerCell+1;
     }
     if(!l1tCellsMappingStream.eof()) edm::LogWarning("HGCalTriggerGeometry") << "Error reading L1TCellsMapping'"<<waferType<<" "<<cell<<" "<<triggerCell<<"' \n";
     l1tCellsMappingStream.close();
 }

HGCalTriggerGeometryBase::geom_set HGCalTriggerGeometryHexImp2::getCellsFromModule ( const unsigned module_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 150 of file HGCalTriggerGeometryHexImp2.cc.

References HGCEE, HGCHEF, HGCalDetId::layer(), module_to_wafers_ee_, module_to_wafers_fh_, number_cells_in_wafers_, DetId::subdetId(), HGCalDetId::wafer(), HGCalDetId::waferType(), and HGCalDetId::zside().

Referenced by getModulePosition().

 {
 
     HGCalDetId module_det_id(module_id);
     unsigned module = module_det_id.wafer();
     int wafer_type = module_det_id.waferType();
     unsigned subdet = module_det_id.subdetId();
     std::pair<std::unordered_multimap<short, short>::const_iterator,
         std::unordered_multimap<short, short>::const_iterator> wafer_itrs;
     switch(subdet)
     {
         case ForwardSubdetector::HGCEE:
             wafer_itrs = module_to_wafers_ee_.equal_range(module);
             break;
         case ForwardSubdetector::HGCHEF:
             wafer_itrs = module_to_wafers_fh_.equal_range(module);
             break;
         default:
             edm::LogError("HGCalTriggerGeometry") << "Unknown module->wafers mapping for subdet "<<subdet<<"\n";
             return geom_set();
     };
     geom_set cell_det_ids;
     for(auto wafer_itr=wafer_itrs.first; wafer_itr!=wafer_itrs.second; wafer_itr++)
     {
         // loop on the cells in each wafer
         for(int cell=0; cell<number_cells_in_wafers_.at(wafer_type); cell++)
         {
             cell_det_ids.emplace(HGCalDetId((ForwardSubdetector)module_det_id.subdetId(), module_det_id.zside(), module_det_id.layer(), module_det_id.waferType(), wafer_itr->second, cell).rawId());
         }
     }
     return cell_det_ids;
 }

HGCalTriggerGeometryBase::geom_set HGCalTriggerGeometryHexImp2::getCellsFromTriggerCell ( const unsigned trigger_cell_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 133 of file HGCalTriggerGeometryHexImp2.cc.

References HGCalDetId::cell(), HGCalDetId::layer(), DetId::subdetId(), trigger_cells_to_cells_, HGCalDetId::wafer(), HGCalDetId::waferType(), and HGCalDetId::zside().

Referenced by getTriggerCellPosition().

 {
     HGCalDetId trigger_cell_det_id(trigger_cell_id);
     int wafer_type = trigger_cell_det_id.waferType();
     unsigned trigger_cell = trigger_cell_det_id.cell();
     // FIXME: better way to do this TC->cell mapping?
     const auto& cell_range = trigger_cells_to_cells_.equal_range(std::make_pair(wafer_type,trigger_cell));
     geom_set cell_det_ids;
     for(auto tc_c_itr=cell_range.first; tc_c_itr!=cell_range.second; tc_c_itr++)
     {
         cell_det_ids.emplace(HGCalDetId((ForwardSubdetector)trigger_cell_det_id.subdetId(), trigger_cell_det_id.zside(), trigger_cell_det_id.layer(), trigger_cell_det_id.waferType(), trigger_cell_det_id.wafer(), tc_c_itr->second).rawId());
     }
     return cell_det_ids;
 }

unsigned HGCalTriggerGeometryHexImp2::getModuleFromCell ( const unsigned cell_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 87 of file HGCalTriggerGeometryHexImp2.cc.

References HGCEE, HGCHEF, HGCalDetId::kHGCalCellMask, HGCalDetId::layer(), python.rootplot.argparse::module, DetId::subdetId(), HGCalDetId::wafer(), wafer_to_module_ee_, wafer_to_module_fh_, HGCalDetId::waferType(), and HGCalDetId::zside().

 {
     HGCalDetId cell_det_id(cell_id);
     unsigned wafer = cell_det_id.wafer();
     unsigned subdet = cell_det_id.subdetId();
     unsigned module = 0;
     switch(subdet)
     {
         case ForwardSubdetector::HGCEE:
             module = wafer_to_module_ee_.at(wafer);
             break;
         case ForwardSubdetector::HGCHEF:
             module = wafer_to_module_fh_.at(wafer);
             break;
         default:
             edm::LogError("HGCalTriggerGeometry") << "Unknown wafer->module mapping for subdet "<<subdet<<"\n";
             return 0;
     };
     return HGCalDetId((ForwardSubdetector)cell_det_id.subdetId(), cell_det_id.zside(), cell_det_id.layer(), cell_det_id.waferType(), module, HGCalDetId::kHGCalCellMask).rawId();
 }

unsigned HGCalTriggerGeometryHexImp2::getModuleFromTriggerCell ( const unsigned trigger_cell_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 110 of file HGCalTriggerGeometryHexImp2.cc.

References HGCEE, HGCHEF, HGCalDetId::kHGCalCellMask, HGCalDetId::layer(), python.rootplot.argparse::module, DetId::subdetId(), HGCalDetId::wafer(), wafer_to_module_ee_, wafer_to_module_fh_, HGCalDetId::waferType(), and HGCalDetId::zside().

 {
     HGCalDetId trigger_cell_det_id(trigger_cell_id);
     unsigned wafer = trigger_cell_det_id.wafer();
     unsigned subdet = trigger_cell_det_id.subdetId();
     unsigned module = 0;
     switch(subdet)
     {
         case ForwardSubdetector::HGCEE:
             module = wafer_to_module_ee_.at(wafer);
             break;
         case ForwardSubdetector::HGCHEF:
             module = wafer_to_module_fh_.at(wafer);
             break;
         default:
             edm::LogError("HGCalTriggerGeometry") << "Unknown wafer->module mapping for subdet "<<subdet<<"\n";
             return 0;
     };
     return HGCalDetId((ForwardSubdetector)trigger_cell_det_id.subdetId(), trigger_cell_det_id.zside(), trigger_cell_det_id.layer(), trigger_cell_det_id.waferType(), module, HGCalDetId::kHGCalCellMask).rawId();
 }

GlobalPoint HGCalTriggerGeometryHexImp2::getModulePosition ( const unsigned module_det_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 306 of file HGCalTriggerGeometryHexImp2.cc.

References es_info_, HGCalTriggerGeometryBase::es_info::geom_ee, HGCalTriggerGeometryBase::es_info::geom_fh, getCellsFromModule(), HGCEE, and DetId::subdetId().

 {
     // Position: barycenter of the module.
     Basic3DVector<float> moduleVector(0.,0.,0.);
     const auto cell_ids = getCellsFromModule(module_det_id);
     for(const auto& cell : cell_ids)
     {
         HGCalDetId cellDetId(cell);
         moduleVector += (cellDetId.subdetId()==ForwardSubdetector::HGCEE ? es_info_.geom_ee->getPosition(cellDetId) :  es_info_.geom_fh->getPosition(cellDetId)).basicVector();
     }
     return GlobalPoint( moduleVector/cell_ids.size() );
 }

HGCalTriggerGeometryBase::geom_ordered_set HGCalTriggerGeometryHexImp2::getOrderedCellsFromModule ( const unsigned module_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 185 of file HGCalTriggerGeometryHexImp2.cc.

References HGCEE, HGCHEF, HGCalDetId::layer(), module_to_wafers_ee_, module_to_wafers_fh_, number_cells_in_wafers_, DetId::subdetId(), HGCalDetId::wafer(), HGCalDetId::waferType(), and HGCalDetId::zside().

 {
     HGCalDetId module_det_id(module_id);
     unsigned module = module_det_id.wafer();
     int wafer_type = module_det_id.waferType();
     unsigned subdet = module_det_id.subdetId();
     std::pair<std::unordered_multimap<short, short>::const_iterator,
         std::unordered_multimap<short, short>::const_iterator> wafer_itrs;
     switch(subdet)
     {
         case ForwardSubdetector::HGCEE:
             wafer_itrs = module_to_wafers_ee_.equal_range(module);
             break;
         case ForwardSubdetector::HGCHEF:
             wafer_itrs = module_to_wafers_fh_.equal_range(module);
             break;
         default:
             edm::LogError("HGCalTriggerGeometry") << "Unknown module->wafers mapping for subdet "<<subdet<<"\n";
             return geom_ordered_set();
     };
     geom_ordered_set cell_det_ids;
     for(auto wafer_itr=wafer_itrs.first; wafer_itr!=wafer_itrs.second; wafer_itr++)
     {
         // loop on the cells in each wafer
         for(int cell=0; cell<number_cells_in_wafers_.at(wafer_type); cell++)
         {
             cell_det_ids.emplace(HGCalDetId((ForwardSubdetector)module_det_id.subdetId(), module_det_id.zside(), module_det_id.layer(), module_det_id.waferType(), wafer_itr->second, cell).rawId());
         }
     }
     return cell_det_ids;
 }

HGCalTriggerGeometryBase::geom_ordered_set HGCalTriggerGeometryHexImp2::getOrderedTriggerCellsFromModule ( const unsigned module_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 254 of file HGCalTriggerGeometryHexImp2.cc.

References HGCEE, HGCHEF, HGCalDetId::layer(), module_to_wafers_ee_, module_to_wafers_fh_, number_trigger_cells_in_wafers_, DetId::subdetId(), HGCalDetId::wafer(), HGCalDetId::waferType(), and HGCalDetId::zside().

 {
     HGCalDetId module_det_id(module_id);
     unsigned module = module_det_id.wafer();
     int wafer_type = module_det_id.waferType();
     unsigned subdet = module_det_id.subdetId();
     std::pair<std::unordered_multimap<short, short>::const_iterator,
         std::unordered_multimap<short, short>::const_iterator> wafer_itrs;
     switch(subdet)
     {
         case ForwardSubdetector::HGCEE:
             wafer_itrs = module_to_wafers_ee_.equal_range(module);
             break;
         case ForwardSubdetector::HGCHEF:
             wafer_itrs = module_to_wafers_fh_.equal_range(module);
             break;
         default:
             edm::LogError("HGCalTriggerGeometry") << "Unknown module->wafers mapping for subdet "<<subdet<<"\n";
             return geom_ordered_set();
     };
     geom_ordered_set trigger_cell_det_ids;
     // loop on the wafers included in the module
     for(auto wafer_itr=wafer_itrs.first; wafer_itr!=wafer_itrs.second; wafer_itr++)
     {
         // loop on the trigger cells in each wafer
         for(int trigger_cell=0; trigger_cell<number_trigger_cells_in_wafers_.at(wafer_type); trigger_cell++)
         {
             trigger_cell_det_ids.emplace(HGCalDetId((ForwardSubdetector)module_det_id.subdetId(), module_det_id.zside(), module_det_id.layer(), module_det_id.waferType(), wafer_itr->second, trigger_cell).rawId());
         }
     }
     return trigger_cell_det_ids;
 }

unsigned HGCalTriggerGeometryHexImp2::getTriggerCellFromCell ( const unsigned cell_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 75 of file HGCalTriggerGeometryHexImp2.cc.

References HGCalDetId::cell(), cells_to_trigger_cells_, HGCalDetId::layer(), DetId::subdetId(), HGCalDetId::wafer(), HGCalDetId::waferType(), and HGCalDetId::zside().

 {
     HGCalDetId cell_det_id(cell_id);
     int wafer_type = cell_det_id.waferType();
     unsigned cell = cell_det_id.cell();
     // FIXME: better way to do this cell->TC mapping?
     unsigned trigger_cell = cells_to_trigger_cells_.at(std::make_pair(wafer_type,cell));
     return HGCalDetId((ForwardSubdetector)cell_det_id.subdetId(), cell_det_id.zside(), cell_det_id.layer(), cell_det_id.waferType(), cell_det_id.wafer(), trigger_cell).rawId();
 }

GlobalPoint HGCalTriggerGeometryHexImp2::getTriggerCellPosition ( const unsigned trigger_cell_det_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 290 of file HGCalTriggerGeometryHexImp2.cc.

References es_info_, HGCalTriggerGeometryBase::es_info::geom_ee, HGCalTriggerGeometryBase::es_info::geom_fh, getCellsFromTriggerCell(), HGCEE, and DetId::subdetId().

 {
     // Position: barycenter of the trigger cell.
     Basic3DVector<float> triggerCellVector(0.,0.,0.);
     const auto cell_ids = getCellsFromTriggerCell(trigger_cell_det_id);
     for(const auto& cell : cell_ids)
     {
         HGCalDetId cellDetId(cell);
         triggerCellVector += (cellDetId.subdetId()==ForwardSubdetector::HGCEE ? es_info_.geom_ee->getPosition(cellDetId) :  es_info_.geom_fh->getPosition(cellDetId)).basicVector();
     }
     return GlobalPoint( triggerCellVector/cell_ids.size() );
 
 }

HGCalTriggerGeometryBase::geom_set HGCalTriggerGeometryHexImp2::getTriggerCellsFromModule ( const unsigned module_id ) const

finaloverridevirtual

Implements HGCalTriggerGeometryBase.

Definition at line 219 of file HGCalTriggerGeometryHexImp2.cc.

References HGCEE, HGCHEF, HGCalDetId::layer(), module_to_wafers_ee_, module_to_wafers_fh_, number_trigger_cells_in_wafers_, DetId::subdetId(), HGCalDetId::wafer(), HGCalDetId::waferType(), and HGCalDetId::zside().

 {
     HGCalDetId module_det_id(module_id);
     unsigned module = module_det_id.wafer();
     int wafer_type = module_det_id.waferType();
     unsigned subdet = module_det_id.subdetId();
     std::pair<std::unordered_multimap<short, short>::const_iterator,
         std::unordered_multimap<short, short>::const_iterator> wafer_itrs;
     switch(subdet)
     {
         case ForwardSubdetector::HGCEE:
             wafer_itrs = module_to_wafers_ee_.equal_range(module);
             break;
         case ForwardSubdetector::HGCHEF:
             wafer_itrs = module_to_wafers_fh_.equal_range(module);
             break;
         default:
             edm::LogError("HGCalTriggerGeometry") << "Unknown module->wafers mapping for subdet "<<subdet<<"\n";
             return geom_set();
     };
     geom_set trigger_cell_det_ids;
     // loop on the wafers included in the module
     for(auto wafer_itr=wafer_itrs.first; wafer_itr!=wafer_itrs.second; wafer_itr++)
     {
         // loop on the trigger cells in each wafer
         for(int trigger_cell=0; trigger_cell<number_trigger_cells_in_wafers_.at(wafer_type); trigger_cell++)
         {
             trigger_cell_det_ids.emplace(HGCalDetId((ForwardSubdetector)module_det_id.subdetId(), module_det_id.zside(), module_det_id.layer(), module_det_id.waferType(), wafer_itr->second, trigger_cell).rawId());
         }
     }
     return trigger_cell_det_ids;
 }

void HGCalTriggerGeometryHexImp2::initialize ( const es_info & esInfo )

finaloverridevirtual

Definition at line 64 of file HGCalTriggerGeometryHexImp2.cc.

References es_info_, and fillMaps().

 {
     edm::LogWarning("HGCalTriggerGeometry") << "WARNING: This HGCal trigger geometry is incomplete.\n"\
                                             << "WARNING: There is no neighbor information.\n";
     es_info_ = esInfo;
     fillMaps(esInfo);
 
 }