PFRecHitProducerPS Class Reference

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 (edm::Run &run, const edm::EventSetup &es)
virtual void	endRun ()
	PFRecHitProducer (const edm::ParameterSet &)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup)
	~PFRecHitProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

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
typedef WorkerT< EDProducer >	WorkerType
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::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)
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.

{}

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, CaloNavigator< T >::east(), PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, geometry, PFLayer::HCAL_BARREL1, PFLayer::HCAL_ENDCAP, CaloNavigator< T >::home(), i, reco::PFRecHit::layer(), north, CaloNavigator< T >::north(), PFLayer::PS1, PFLayer::PS2, DetId::rawId(), south, CaloNavigator< T >::south(), west, and CaloNavigator< T >::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().