Producer for particle flow rechits (PFRecHit) More...

#include <PFRecHitProducerPS.h>

Inheritance diagram for PFRecHitProducerPS:

Public Member Functions
	PFRecHitProducerPS (const edm::ParameterSet &)

	~PFRecHitProducerPS ()

Public Member Functions inherited from PFRecHitProducer
virtual void	beginRun (const edm::Run &run, const edm::EventSetup &es) override

	PFRecHitProducer (const edm::ParameterSet &)

void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

	~PFRecHitProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
void	createRecHits (std::vector< reco::PFRecHit > &rechits, std::vector< reco::PFRecHit > &rechitsCleaned, edm::Event &, const edm::EventSetup &)

void	findRecHitNeighbours (reco::PFRecHit &rh, const std::map< unsigned, unsigned > &sortedHits, const CaloSubdetectorTopology &barrelTopo, const CaloSubdetectorGeometry &barrelGeom, const CaloSubdetectorTopology &endcapTopo, const CaloSubdetectorGeometry &endcapGeom)

Private Attributes
edm::InputTag	inputTagEcalRecHitsES_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Types inherited from PFRecHitProducer
typedef std::map< unsigned, unsigned >::const_iterator	IDH

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from PFRecHitProducer
const HcalPFCorrs *	myPFCorr

const EcalChannelStatus *	theEcalChStatus

const HcalChannelQuality *	theHcalChStatus

const CaloTowerConstituentsMap *	theTowerConstituentsMap

double	thresh_Barrel_
	rechits with E < threshold will not give rise to a PFRecHit More...

double	thresh_Endcap_

bool	verbose_
	verbose ? More...

Detailed Description

Producer for particle flow rechits (PFRecHit)

Author: Colin Bernet

Date: february 2008

Definition at line 32 of file PFRecHitProducerPS.h.

Constructor & Destructor Documentation

PFRecHitProducerPS::PFRecHitProducerPS ( const edm::ParameterSet & iConfig )

explicit

Definition at line 32 of file PFRecHitProducerPS.cc.

References edm::ParameterSet::getParameter(), and inputTagEcalRecHitsES_.

  : PFRecHitProducer(iConfig) {
 
   // access to the collections of rechits
   
   inputTagEcalRecHitsES_ = 
     iConfig.getParameter<InputTag>("ecalRecHitsES");
 }

PFRecHitProducerPS::~PFRecHitProducerPS ( )

Definition at line 45 of file PFRecHitProducerPS.cc.

45 {}

Member Function Documentation

void PFRecHitProducerPS::createRecHits	(	std::vector< reco::PFRecHit > &	rechits,
		std::vector< reco::PFRecHit > &	rechitsCleaned,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

gets PS rechits, translate them to PFRecHits, which are stored in the rechits vector

Implements PFRecHitProducer.

Definition at line 49 of file PFRecHitProducerPS.cc.

References CaloRecHit::detid(), cond::rpcobgas::detid, DetId::Ecal, EcalPreshower, CaloRecHit::energy(), relval_parameters_module::energy, edm::hlt::Exception, findRecHitNeighbours(), newFWLiteAna::found, edm::EventSetup::get(), edm::Event::getByLabel(), CaloCellGeometry::getCorners(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), i, inputTagEcalRecHitsES_, edm::HandleBase::isValid(), LogDebug, PFLayer::NONE, ESDetId::plane(), position, edm::Handle< T >::product(), PFLayer::PS1, PFLayer::PS2, DetId::rawId(), reco::PFRecHit::setNECorner(), reco::PFRecHit::setNWCorner(), reco::PFRecHit::setSECorner(), reco::PFRecHit::setSWCorner(), PFRecHitProducer::thresh_Endcap_, x, PV3DBase< T, PVType, FrameType >::x(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

                                                     {
 
   map<unsigned, unsigned > idSortedRecHits;
   typedef map<unsigned, unsigned >::iterator IDH;
 
 
   // get the ps geometry
   edm::ESHandle<CaloGeometry> geoHandle;
   iSetup.get<CaloGeometryRecord>().get(geoHandle);
     
   const CaloSubdetectorGeometry *psGeometry = 
     geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
     
 
   // ShR 28 Jul 2008: check if geometry is NULL. If so we are using
   // partial CMS gemoetry in Pilot1/2 scenarios which do not include the preshower
   if(!psGeometry) {
     LogDebug("PFRecHitProducerPS") << "No EcalPreshower geometry available. putting empty PS rechits collection in event";
     return;
   }
 
 
   // get the ps topology
   EcalPreshowerTopology psTopology(geoHandle);
 
   // process rechits
   Handle< EcalRecHitCollection >   pRecHits;
 
 
 
   bool found = iEvent.getByLabel(inputTagEcalRecHitsES_,
                  pRecHits);
 
   if(!found) {
     ostringstream err;
     err<<"could not find rechits "<<inputTagEcalRecHitsES_;
     LogError("PFRecHitProducerPS")<<err.str()<<endl;
     
     throw cms::Exception( "MissingProduct", err.str());
   }
   else {
     assert( pRecHits.isValid() );
 
     const EcalRecHitCollection& psrechits = *( pRecHits.product() );
     typedef EcalRecHitCollection::const_iterator IT;
  
     for(IT i=psrechits.begin(); i!=psrechits.end(); i++) {
       const EcalRecHit& hit = *i;
       
       double energy = hit.energy();
       if( energy < thresh_Endcap_ ) continue; 
             
       const ESDetId& detid = hit.detid();
       const CaloCellGeometry *thisCell = psGeometry->getGeometry(detid);
      
       if(!thisCell) {
     LogError("PFRecHitProducerPS")<<"warning detid "<<detid.rawId()
                      <<" not found in preshower geometry"
                      <<endl;
     return;
       }
       
       const GlobalPoint& position = thisCell->getPosition();
      
       PFLayer::Layer layer = PFLayer::NONE;
             
       switch( detid.plane() ) {
       case 1:
     layer = PFLayer::PS1;
     break;
       case 2:
     layer = PFLayer::PS2;
     break;
       default:
     LogError("PFRecHitProducerPS")
       <<"incorrect preshower plane !! plane number "
       <<detid.plane()<<endl;
     assert(0);
       }
  
       reco::PFRecHit *pfrh
     = new reco::PFRecHit( detid.rawId(), layer, energy, 
                   position.x(), position.y(), position.z(), 
                   0,0,0 );
       
       const CaloCellGeometry::CornersVec& corners = thisCell->getCorners();
       assert( corners.size() == 8 );
       
       pfrh->setNECorner( corners[0].x(), corners[0].y(),  corners[0].z() );
       pfrh->setSECorner( corners[1].x(), corners[1].y(),  corners[1].z() );
       pfrh->setSWCorner( corners[2].x(), corners[2].y(),  corners[2].z() );
       pfrh->setNWCorner( corners[3].x(), corners[3].y(),  corners[3].z() );
 
       
       // if( !pfrh ) continue; // problem with this rechit. skip it
 
       rechits.push_back( *pfrh );
       delete pfrh;
       idSortedRecHits.insert( make_pair(detid.rawId(), rechits.size()-1 ) );   
     }
   }
 
   // do navigation
   for(unsigned i=0; i<rechits.size(); i++ ) {
     
     findRecHitNeighbours( rechits[i], idSortedRecHits, 
               psTopology, 
               *psGeometry, 
               psTopology,
               *psGeometry);
   }
 }

void PFRecHitProducerPS::findRecHitNeighbours	(	reco::PFRecHit &	rh,
		const std::map< unsigned, unsigned > &	sortedHits,
		const CaloSubdetectorTopology &	barrelTopo,
		const CaloSubdetectorGeometry &	barrelGeom,
		const CaloSubdetectorTopology &	endcapTopo,
		const CaloSubdetectorGeometry &	endcapGeom
	)

private

find and set the neighbours to a given rechit this works for ecal, hcal, ps COLIN remonter cette fonction dans la classe de base

Definition at line 169 of file PFRecHitProducerPS.cc.

References reco::PFRecHit::add4Neighbour(), reco::PFRecHit::add8Neighbour(), cond::rpcobgas::detid, reco::PFRecHit::detId(), east, PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, geometry, PFLayer::HCAL_BARREL1, PFLayer::HCAL_ENDCAP, i, reco::PFRecHit::layer(), north, PFLayer::PS1, PFLayer::PS2, DetId::rawId(), south, and west.

Referenced by createRecHits().

                                                   {
   
   DetId detid( rh.detId() );
 
   const CaloSubdetectorTopology* topology = 0;
   const CaloSubdetectorGeometry* geometry = 0;
   
   switch( rh.layer() ) {
   case PFLayer::ECAL_ENDCAP: 
     topology = &endcapTopology;
     geometry = &endcapGeometry;
     break;
   case PFLayer::ECAL_BARREL: 
     topology = &barrelTopology;
     geometry = &barrelGeometry;
     break;
   case PFLayer::HCAL_ENDCAP:
     topology = &endcapTopology;
     geometry = &endcapGeometry;
     break;
   case PFLayer::HCAL_BARREL1:
     topology = &barrelTopology;
     geometry = &barrelGeometry;
     break;
   case PFLayer::PS1:
   case PFLayer::PS2:
     topology = &barrelTopology;
     geometry = &barrelGeometry;
     break;
   default:
     assert(0);
   }
   
   assert( topology && geometry );
 
   CaloNavigator<DetId> navigator(detid, topology);
 
   DetId north = navigator.north();  
   
   DetId northeast(0);
   if( north != DetId(0) ) {
     northeast = navigator.east();  
   }
   navigator.home();
 
 
   DetId south = navigator.south();
 
   
 
   DetId southwest(0); 
   if( south != DetId(0) ) {
     southwest = navigator.west();
   }
   navigator.home();
 
 
   DetId east = navigator.east();
   DetId southeast;
   if( east != DetId(0) ) {
     southeast = navigator.south(); 
   }
   navigator.home();
   DetId west = navigator.west();
   DetId northwest;
   if( west != DetId(0) ) {   
     northwest = navigator.north();  
   }
   navigator.home();
     
   IDH i = sortedHits.find( north.rawId() );
   if(i != sortedHits.end() ) 
     rh.add4Neighbour( i->second );
   
   i = sortedHits.find( northeast.rawId() );
   if(i != sortedHits.end() ) 
     rh.add8Neighbour( i->second );
   
   i = sortedHits.find( south.rawId() );
   if(i != sortedHits.end() ) 
     rh.add4Neighbour( i->second );
     
   i = sortedHits.find( southwest.rawId() );
   if(i != sortedHits.end() ) 
     rh.add8Neighbour( i->second );
     
   i = sortedHits.find( east.rawId() );
   if(i != sortedHits.end() ) 
     rh.add4Neighbour( i->second );
     
   i = sortedHits.find( southeast.rawId() );
   if(i != sortedHits.end() ) 
     rh.add8Neighbour( i->second );
     
   i = sortedHits.find( west.rawId() );
   if(i != sortedHits.end() ) 
      rh.add4Neighbour( i->second );
    
   i = sortedHits.find( northwest.rawId() );
   if(i != sortedHits.end() ) 
     rh.add8Neighbour( i->second );
     
 
 }

Member Data Documentation

edm::InputTag PFRecHitProducerPS::inputTagEcalRecHitsES_

private

Definition at line 63 of file PFRecHitProducerPS.h.

Referenced by createRecHits(), and PFRecHitProducerPS().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation