/data/refman/pasoursint/CMSSW_4_1_8_patch13/src/RecoLocalTracker/SubCollectionProducers/src/PixelClusterSelectorTopBottom.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/MakerMacros.h"
#include "RecoLocalTracker/SubCollectionProducers/interface/PixelClusterSelectorTopBottom.h"

void PixelClusterSelectorTopBottom::produce( edm::Event& event, const edm::EventSetup& setup) {

  edm::Handle< SiPixelClusterCollectionNew > input;
  event.getByLabel(label_, input);

  edm::ESHandle<TrackerGeometry> geom;
  setup.get<TrackerDigiGeometryRecord>().get( geom );
  const TrackerGeometry& theTracker( *geom );
  
  std::auto_ptr<SiPixelClusterCollectionNew> output( new SiPixelClusterCollectionNew() );

  for (SiPixelClusterCollectionNew::const_iterator clustSet = input->begin(); clustSet!=input->end(); ++clustSet) {
    edmNew::DetSet<SiPixelCluster>::const_iterator clustIt = clustSet->begin();
    edmNew::DetSet<SiPixelCluster>::const_iterator end     = clustSet->end();
    
    DetId detIdObject( clustSet->detId() );
    edmNew::DetSetVector<SiPixelCluster>::FastFiller spc(*output, detIdObject);
    const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*> (theTracker.idToDet(detIdObject) );
    const PixelTopology * topol = (&(theGeomDet->specificTopology()));
    
    for(; clustIt!=end;++clustIt) {
      LocalPoint  lpclust = topol->localPosition(MeasurementPoint((*clustIt).x(),(*clustIt).y()));
      GlobalPoint GPclust = theGeomDet->surface().toGlobal(Local3DPoint(lpclust.x(),lpclust.y(),lpclust.z()));
      double clustY = GPclust.y();
      if ((clustY * y_) > 0) {
        spc.push_back(*clustIt);
      }
    }
  }
  event.put( output );  
}

DEFINE_FWK_MODULE( PixelClusterSelectorTopBottom );