HGCalTriggerGeometryImp1 Class Reference

Inheritance diagram for HGCalTriggerGeometryImp1:

Public Member Functions
	HGCalTriggerGeometryImp1 (const edm::ParameterSet &conf)
void	initialize (const edm::ESHandle< CaloGeometry > &) 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 std::string &	name () const
virtual	~HGCalTriggerGeometryBase ()

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)
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 12 of file HGCalTriggerGeometryImp1.cc.

Constructor & Destructor Documentation

HGCalTriggerGeometryImp1::HGCalTriggerGeometryImp1

(

const edm::ParameterSet &

conf

)

Definition at line 25 of file HGCalTriggerGeometryImp1.cc.

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

Member Function Documentation

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

finalvirtual

Implements HGCalTriggerGeometryBase.

Definition at line 34 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(), HGCalTriggerGeometryBase::setCaloGeometry(), HGCalTriggerGeometryGenericMapping::trigger_cells_, HGCalTriggerGeometryGenericMapping::trigger_cells_to_modules_, HGCTriggerDetId::UndefinedCell(), z, and ecaldqm::zside().

{
    // 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);
    //
    // 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)) );
    }
}

Member Data Documentation

edm::FileInPath HGCalTriggerGeometryImp1::l1tCellsMapping_

private

Definition at line 20 of file HGCalTriggerGeometryImp1.cc.

Referenced by initialize().