ShallowRechitClustersProducer Class Reference

#include <ShallowRechitClustersProducer.h>

Inheritance diagram for ShallowRechitClustersProducer:

Public Member Functions
	ShallowRechitClustersProducer (const edm::ParameterSet &)
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	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
std::vector< edm::InputTag >	inputTags
std::string	Prefix
std::string	Suffix
edm::InputTag	theClustersLabel

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)
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)

Detailed Description

Definition at line 8 of file ShallowRechitClustersProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 16 of file ShallowRechitClustersProducer.cc.

References Prefix, and Suffix.

  :  Suffix       ( iConfig.getParameter<std::string>("Suffix") ),
     Prefix       ( iConfig.getParameter<std::string>("Prefix") ),
     theClustersLabel( iConfig.getParameter<edm::InputTag>("Clusters")),
     inputTags    ( iConfig.getParameter<std::vector<edm::InputTag> >("InputTags"))
{
  produces <std::vector<float> >        ( Prefix + "strip"      + Suffix );   
  produces <std::vector<float> >        ( Prefix + "merr"       + Suffix );   
  produces <std::vector<float> >        ( Prefix + "localx"     + Suffix );   
  produces <std::vector<float> >        ( Prefix + "localy"     + Suffix );   
  produces <std::vector<float> >        ( Prefix + "localxerr"  + Suffix );   
  produces <std::vector<float> >        ( Prefix + "localyerr"  + Suffix );   
  produces <std::vector<float> >        ( Prefix + "globalx"    + Suffix );   
  produces <std::vector<float> >        ( Prefix + "globaly"    + Suffix );   
  produces <std::vector<float> >        ( Prefix + "globalz"    + Suffix );   
}

Member Function Documentation

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

privatevirtual

Implements edm::EDProducer.

Definition at line 34 of file ShallowRechitClustersProducer.cc.

References SiStripRecHit2D::cluster(), edm::hlt::Exception, TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByLabel(), i, collect_tpl::input, inputTags, BaseSiTrackerRecHit2DLocalPos::localPosition(), BaseSiTrackerRecHit2DLocalPos::localPositionError(), shallow::make_cluster_map(), Topology::measurementError(), Prefix, edm::Event::put(), DetId::rawId(), findQualityFiles::size, StripGeomDetUnit::specificTopology(), mathSSE::sqrt(), strip(), StripTopology::strip(), Suffix, theClustersLabel, GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), ExpressReco_HICollisions_FallBack::x, LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), ExpressReco_HICollisions_FallBack::y, LocalError::yy(), and z.

Referenced by python.JSONExport.JsonExport::export(), and python.HTMLExport.HTMLExport::export().

                                                       {
  shallow::CLUSTERMAP clustermap = shallow::make_cluster_map(iEvent,theClustersLabel);

  int size = clustermap.size();
  std::auto_ptr<std::vector<float> >  strip       ( new std::vector<float>(size,  -10000  ));   
  std::auto_ptr<std::vector<float> >  merr        ( new std::vector<float>(size,  -10000  ));   
  std::auto_ptr<std::vector<float> >  localx      ( new std::vector<float>(size,  -10000  ));   
  std::auto_ptr<std::vector<float> >  localy      ( new std::vector<float>(size,  -10000  ));   
  std::auto_ptr<std::vector<float> >  localxerr   ( new std::vector<float>(size,  -1  ));   
  std::auto_ptr<std::vector<float> >  localyerr   ( new std::vector<float>(size,  -1  ));     
  std::auto_ptr<std::vector<float> >  globalx     ( new std::vector<float>(size,  -10000  ));   
  std::auto_ptr<std::vector<float> >  globaly     ( new std::vector<float>(size,  -10000  ));   
  std::auto_ptr<std::vector<float> >  globalz     ( new std::vector<float>(size,  -10000  ));   

  edm::ESHandle<TrackerGeometry> theTrackerGeometry; iSetup.get<TrackerDigiGeometryRecord>().get( theTrackerGeometry );

  BOOST_FOREACH(const edm::InputTag& input, inputTags ) {  edm::Handle<SiStripRecHit2DCollection> recHits; iEvent.getByLabel(input, recHits);
    BOOST_FOREACH( const SiStripRecHit2DCollection::value_type& ds, *recHits) {
      BOOST_FOREACH( const SiStripRecHit2D& hit, ds) {
    
    shallow::CLUSTERMAP::iterator cluster = clustermap.find( std::make_pair(hit.geographicalId().rawId(), hit.cluster()->firstStrip()   ) );
    if(cluster != clustermap.end() ) {
      const StripGeomDetUnit* theStripDet = dynamic_cast<const StripGeomDetUnit*>( theTrackerGeometry->idToDet( hit.geographicalId() ) );
      unsigned int i = cluster->second;
      strip->at(i)  =    theStripDet->specificTopology().strip(hit.localPosition());
      merr->at(i)   =    sqrt(theStripDet->specificTopology().measurementError(hit.localPosition(), hit.localPositionError()).uu());
      localx->at(i) =    hit.localPosition().x();
      localy->at(i) =    hit.localPosition().y();
      localxerr->at(i) = sqrt(hit.localPositionError().xx());
      localyerr->at(i) = sqrt(hit.localPositionError().yy());
      globalx->at(i) =   theStripDet->toGlobal(hit.localPosition()).x();
      globaly->at(i) =   theStripDet->toGlobal(hit.localPosition()).y();
      globalz->at(i) =   theStripDet->toGlobal(hit.localPosition()).z();
    }
    else {throw cms::Exception("cluster not found");}
      }
    }
  }
  
  iEvent.put( strip,       Prefix + "strip"      + Suffix );   
  iEvent.put( merr,        Prefix + "merr"       + Suffix );   
  iEvent.put( localx ,     Prefix + "localx"     + Suffix );   
  iEvent.put( localy ,     Prefix + "localy"     + Suffix );   
  iEvent.put( localxerr ,  Prefix + "localxerr"  + Suffix );   
  iEvent.put( localyerr ,  Prefix + "localyerr"  + Suffix );   
  iEvent.put( globalx ,    Prefix + "globalx"    + Suffix );   
  iEvent.put( globaly ,    Prefix + "globaly"    + Suffix );   
  iEvent.put( globalz ,    Prefix + "globalz"    + Suffix );   
}

Member Data Documentation

std::vector<edm::InputTag> ShallowRechitClustersProducer::inputTags

private

Definition at line 15 of file ShallowRechitClustersProducer.h.

Referenced by python.Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputEventContent(), produce(), python.Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::readConnections(), and python.Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::uses().

std::string ShallowRechitClustersProducer::Prefix

private

Definition at line 13 of file ShallowRechitClustersProducer.h.

Referenced by produce(), and ShallowRechitClustersProducer().

std::string ShallowRechitClustersProducer::Suffix

private

Definition at line 12 of file ShallowRechitClustersProducer.h.

Referenced by produce(), and ShallowRechitClustersProducer().

edm::InputTag ShallowRechitClustersProducer::theClustersLabel

private

Definition at line 14 of file ShallowRechitClustersProducer.h.

Referenced by produce().