#include <SiClusterTranslator.h>

Inheritance diagram for SiClusterTranslator:

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

virtual void	produce (edm::Event &e, const edm::EventSetup &c)

	SiClusterTranslator (const edm::ParameterSet &conf)

virtual	~SiClusterTranslator ()

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	loadPixelClusters (std::map< DetId, std::vector< SiPixelCluster > > &theClusters, SiPixelClusterCollectionNew &theClusterCollection) const

void	loadStripClusters (std::map< DetId, std::vector< SiStripCluster > > &theClusters, edmNew::DetSetVector< SiStripCluster > &theClusterCollection) const

Private Attributes
const TrackerGeometry *	geometry

std::vector< std::pair < PSimHit, int > >	theNewSimHitList

std::vector< edm::DetSet < PixelDigiSimLink > >	thePixelDigiLinkVector

std::vector< edm::DetSet < StripDigiSimLink > >	theStripDigiLinkVector

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

Author: : Douglas Orbaker : Dougl.nosp@m.as.O.nosp@m.rbake.nosp@m.r@ce.nosp@m.rn.ch Clusters Produced by the FastSim and converts them to the standard cluster type.

Definition at line 35 of file SiClusterTranslator.h.

Constructor & Destructor Documentation

SiClusterTranslator::SiClusterTranslator ( const edm::ParameterSet & conf )

explicit

SiClusterTranslator.cc

Description: see SiClusterTranslator.h Authors: R. Ranieri (CERN)

History: Sep 27, 2006 - initial version

Definition at line 56 of file SiClusterTranslator.cc.

 {
   produces<edmNew::DetSetVector<SiStripCluster> >();
   produces<edmNew::DetSetVector<SiPixelCluster> >();
   produces<edm::DetSetVector<StripDigiSimLink> >();
   produces<edm::DetSetVector<PixelDigiSimLink> >();
 }

SiClusterTranslator::~SiClusterTranslator ( )

virtual

Definition at line 65 of file SiClusterTranslator.cc.

65 {}

Member Function Documentation

void SiClusterTranslator::beginRun	(	edm::Run &	run,
		const edm::EventSetup &	es
	)

virtual

Reimplemented from edm::EDProducer.

Definition at line 68 of file SiClusterTranslator.cc.

References geometry, and edm::EventSetup::get().

                                                                     {
 
   // Initialize the Tracker Geometry
   edm::ESHandle<TrackerGeometry> theGeometry;
   es.get<TrackerDigiGeometryRecord> ().get (theGeometry);
   geometry = &(*theGeometry);
 }

void SiClusterTranslator::loadPixelClusters	(	std::map< DetId, std::vector< SiPixelCluster > > &	theClusters,
		SiPixelClusterCollectionNew &	theClusterCollection
	)		const

private

Definition at line 290 of file SiClusterTranslator.cc.

References edmNew::DetSetVector< T >::FastFiller::push_back().

Referenced by produce().

 {
   
   std::map<DetId,std::vector<SiPixelCluster> >::const_iterator 
     it = theClusters.begin();
   std::map<DetId,std::vector<SiPixelCluster> >::const_iterator 
     lastCluster = theClusters.end();
   for( ; it != lastCluster ; ++it ) { 
     edmNew::DetSetVector<SiPixelCluster>::FastFiller spc(theClusterCollection, it->first);
     
     std::vector<SiPixelCluster>::const_iterator clust_it = it->second.begin();
     std::vector<SiPixelCluster>::const_iterator clust_end = it->second.end();
     for( ; clust_it != clust_end ; ++clust_it)  {
       spc.push_back(*clust_it);
     }
   }
 }

void SiClusterTranslator::loadStripClusters	(	std::map< DetId, std::vector< SiStripCluster > > &	theClusters,
		edmNew::DetSetVector< SiStripCluster > &	theClusterCollection
	)		const

private

Definition at line 270 of file SiClusterTranslator.cc.

References edmNew::DetSetVector< T >::FastFiller::push_back().

Referenced by produce().

 {
   std::map<DetId,std::vector<SiStripCluster> >::const_iterator 
     it = theClusters.begin();
   std::map<DetId,std::vector<SiStripCluster> >::const_iterator 
     lastDet = theClusters.end();
   for( ; it != lastDet ; ++it ) { 
     edmNew::DetSetVector<SiStripCluster>::FastFiller cluster_col(theClusterCollection, it->first);
     
     std::vector<SiStripCluster>::const_iterator clust_it = it->second.begin();
     std::vector<SiStripCluster>::const_iterator clust_end = it->second.end();
     for( ; clust_it != clust_end ; ++clust_it)  {
       cluster_col.push_back(*clust_it);
     }
   }
 }

void SiClusterTranslator::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

virtual

Implements edm::EDProducer.

Definition at line 77 of file SiClusterTranslator.cc.

References MixCollection< T >::begin(), DeDxDiscriminatorTools::charge(), ClusterParameterEstimator< T >::clearParameters(), edm::DetSet< T >::data, ClusterParameterEstimator< T >::enterLocalParameters(), error, edm::hlt::Exception, edm::EventSetup::get(), edm::Event::getByLabel(), edm::Event::getByType(), edm::DetSet< T >::id, TrackerGeometry::idToDetUnit(), loadPixelClusters(), loadStripClusters(), PixelTopology::pixel(), PixelChannelIdentifier::pixelToChannel(), position, edm::Handle< T >::product(), edm::Event::put(), DetId::rawId(), RectangularStripTopology::strip(), RadialStripTopology::strip(), DetId::subdetId(), theNewSimHitList, thePixelDigiLinkVector, theStripDigiLinkVector, GeomDetType::topology(), PixelGeomDetUnit::topology(), patCandidatesForDimuonsSequences_cff::tracker, and StripGeomDetUnit::type().

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

 {
   // Step A: Get Inputs (FastGSRecHit's)
   edm::Handle<FastTrackerClusterCollection> theFastClusters; 
   e.getByType(theFastClusters);
   
   edm::ESHandle<TrackerGeometry> tkgeom;
   es.get<TrackerDigiGeometryRecord>().get( tkgeom ); 
   const TrackerGeometry &tracker(*tkgeom);
   
   edm::ESHandle<PixelClusterParameterEstimator> pixelCPE;
   es.get<TkPixelCPERecord>().get("FastPixelCPE",pixelCPE);
   const PixelClusterParameterEstimator &pixelcpe(*pixelCPE);
   
   edm::ESHandle<StripClusterParameterEstimator> stripCPE;
   es.get<TkStripCPERecord>().get("FastStripCPE", stripCPE); 
   const StripClusterParameterEstimator &stripcpe(*stripCPE);
     
   edm::Handle<CrossingFrame<PSimHit> > cf_simhit;
   std::vector<const CrossingFrame<PSimHit> *> cf_simhitvec;
   e.getByLabel("mix","famosSimHitsTrackerHits", cf_simhit);
   cf_simhitvec.push_back(cf_simhit.product());
   
   std::auto_ptr<MixCollection<PSimHit> > allTrackerHits(new MixCollection<PSimHit>(cf_simhitvec));
   int counter =0;
   
   //Clearing vector to hopefully make it run faster.
   theNewSimHitList.clear();
   thePixelDigiLinkVector.clear();
   theStripDigiLinkVector.clear();
 
   for(MixCollection<PSimHit>::iterator it = allTrackerHits->begin(); it!= allTrackerHits->end();it++){
     counter++;
     theNewSimHitList.push_back(std::make_pair((*it), counter));
   }
 
   // Step B: fill a temporary full Cluster collection from the fast Cluster collection
   FastTrackerClusterCollection::const_iterator aCluster = theFastClusters->begin();
   FastTrackerClusterCollection::const_iterator theLastHit = theFastClusters->end();
   std::map< DetId, std::vector<SiPixelCluster> > temporaryPixelClusters;
   std::map< DetId, std::vector<SiStripCluster> > temporaryStripClusters;
   
   //Clearing CPE maps from previous event.
   stripcpe.clearParameters();
   pixelcpe.clearParameters();
   
   int ClusterNum = 0;
   
   // loop on Fast GS Hits
   for ( ; aCluster != theLastHit; ++aCluster ) {
     ClusterNum++;
     
     //Finding SubDet Id of cluster: 1 & 2 = Pixel, 3 & 4 & 5 & 6 = Strip, >6 = ?
     DetId det = aCluster->id();
     unsigned int subdet   = det.subdetId();
     int sim_counter = 0;
     for (std::vector<std::pair<PSimHit,int> >::const_iterator 
        simcount = theNewSimHitList.begin() ; simcount != theNewSimHitList.end(); simcount ++){
       if((aCluster->simtrackId() == (int)(*simcount).first.trackId())&&(det.rawId() == (*simcount).first.detUnitId())&&(aCluster->eeId() == (*simcount).first.eventId().rawId()))
     sim_counter = (*simcount).second;
     }
     if (sim_counter == 0)  throw cms::Exception("SiClusterTranslator") << "No Matching SimHit found.";
     
     //Separating into Pixels and Strips
     
     //Pixel
     if (subdet < 3) {
       //Here is the hard part. From the position of the FastSim Cluster I need to figure out the Pixel location for the Cluster.
       LocalPoint position = aCluster->localPosition();
       LocalError error = aCluster->localPositionError();
       //std::cout << "The pixel charge is " << aCluster->charge() << std::endl;
       int charge = (int)(aCluster->charge() + 0.5);
       //std::cout << "The pixel charge after integer conversion is " << charge << std::endl;
 
       //std::vector<int> digi_vec;
       //while (charge > 255) { 
       //digi_vec.push_back(charge);
       //charge -= 256;
       //}
       //digi_vec.push_back(charge);
       
       const GeomDetUnit *  geoDet = tracker.idToDetUnit(det);
       const PixelGeomDetUnit * pixelDet=(const PixelGeomDetUnit*)(geoDet);
       const PixelTopology& topol=(PixelTopology&)pixelDet->topology();
       
       //Out of pixel is float, but input of pixel is int. Hopeful it works...
       std::pair<float,float> pixelPos_out = topol.pixel(position);
       SiPixelCluster::PixelPos pixelPos((int)pixelPos_out.first, (int)pixelPos_out.second);
       
       //Filling Pixel CPE with information.
       std::pair<int,int> row_col((int)pixelPos_out.first,(int)pixelPos_out.second);
       pixelcpe.enterLocalParameters((unsigned int) det.rawId() , row_col, std::make_pair(position,error));
       
       unsigned int ch = PixelChannelIdentifier::pixelToChannel((int)pixelPos_out.first, (int)pixelPos_out.second);
       
       //Creating a new pixel cluster.
       SiPixelCluster temporaryPixelCluster(pixelPos, charge);
       temporaryPixelClusters[det].push_back(temporaryPixelCluster);
       
       //Making a PixelDigiSimLink.
       edm::DetSet<PixelDigiSimLink> pixelDetSet;
       pixelDetSet.id = det.rawId();
       pixelDetSet.data.push_back(PixelDigiSimLink(ch,
                           aCluster->simtrackId(),
                           EncodedEventId(aCluster->eeId()),
                           1.0));
       thePixelDigiLinkVector.push_back(pixelDetSet);
     } 
     
     //Strips
     else if ((subdet > 2) && (subdet < 7)) {
       //Getting pos/err info from the cluster
       LocalPoint position = aCluster->localPosition();
       LocalError error = aCluster->localPositionError();
       
       //Will have to make charge into ADC eventually...
       uint16_t charge = (uint16_t)(aCluster->charge() + 0.5);
       
       //std::cout << "The charge is " << charge << std::endl;
       
       uint16_t strip_num = 0;
       std::vector<uint16_t> digi_vec;
       while (charge > 255) { 
     digi_vec.push_back(255);
     charge -= 255;
       }
       if (charge > 0) digi_vec.push_back(charge);
       //std::cout << "The digi_vec size is " << digi_vec.size() << std::endl;
       //int totcharge = 0;
       //for(int i = 0; i < digi_vec.size(); ++i) {
       //totcharge += digi_vec[i];
       //} 
       const GeomDetUnit *  geoDet = tracker.idToDetUnit(det);
       const StripGeomDetUnit * stripDet = (const StripGeomDetUnit*)(geoDet);
       
       //3 = TIB, 4 = TID, 5 = TOB, 6 = TEC
       if((subdet == 3) || (subdet == 5)) {
     const RectangularStripTopology& topol=(RectangularStripTopology&)stripDet->type().topology();
     strip_num = (uint16_t)topol.strip(position);
       } else if ((subdet == 4) || (subdet == 6)) {
     const RadialStripTopology& topol=(RadialStripTopology&)stripDet->type().topology();
     strip_num = (uint16_t)topol.strip(position);
       }
       
       //Filling Strip CPE with info.
       stripcpe.enterLocalParameters(det.rawId(), strip_num, std::make_pair(position,error));
       
       //Creating a new strip cluster
       SiStripCluster temporaryStripCluster(det.rawId(), strip_num, digi_vec.begin(), digi_vec.end());
       temporaryStripClusters[det].push_back(temporaryStripCluster);
       
       //Making a StripDigiSimLink
       edm::DetSet<StripDigiSimLink> stripDetSet;
       stripDetSet.id = det.rawId();
       stripDetSet.data.push_back(StripDigiSimLink(strip_num,
                           aCluster->simtrackId(),
                           sim_counter,
                           EncodedEventId(aCluster->eeId()),
                           1.0));
       theStripDigiLinkVector.push_back(stripDetSet);
     } 
     
     //?????
     else {
       throw cms::Exception("SiClusterTranslator") <<
     "Trying to build a cluster that is not in the SiStripTracker or Pixels.\n";
     }
     
   }//Cluster loop
 
   // Step C: from the temporary Cluster collections, create the real ones.
   
   //Pixels
   std::auto_ptr<edmNew::DetSetVector<SiPixelCluster> >
     siPixelClusterCollection(new edmNew::DetSetVector<SiPixelCluster>);
   loadPixelClusters(temporaryPixelClusters, *siPixelClusterCollection);
   std::auto_ptr<edm::DetSetVector<StripDigiSimLink> > stripoutputlink(new edm::DetSetVector<StripDigiSimLink>(theStripDigiLinkVector) );
   
   
   //Strips
   std::auto_ptr<edmNew::DetSetVector<SiStripCluster> > 
     siStripClusterCollection(new edmNew::DetSetVector<SiStripCluster>);
   loadStripClusters(temporaryStripClusters, *siStripClusterCollection);
   std::auto_ptr<edm::DetSetVector<PixelDigiSimLink> > pixeloutputlink(new edm::DetSetVector<PixelDigiSimLink>(thePixelDigiLinkVector) );
   
   // Step D: write output to file
   e.put(siPixelClusterCollection);
   e.put(siStripClusterCollection);
   e.put(stripoutputlink);
   e.put(pixeloutputlink);
 }

Member Data Documentation

const TrackerGeometry* SiClusterTranslator::geometry

private

Definition at line 58 of file SiClusterTranslator.h.

Referenced by beginRun(), and Vispa.Gui.ConnectableWidget.ConnectableWidget::leaveEvent().

std::vector<std::pair<PSimHit, int > > SiClusterTranslator::theNewSimHitList

private

Definition at line 51 of file SiClusterTranslator.h.

Referenced by produce().

std::vector<edm::DetSet<PixelDigiSimLink> > SiClusterTranslator::thePixelDigiLinkVector

private

Definition at line 50 of file SiClusterTranslator.h.

Referenced by produce().

std::vector<edm::DetSet<StripDigiSimLink> > SiClusterTranslator::theStripDigiLinkVector

private

Definition at line 49 of file SiClusterTranslator.h.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

SiClusterTranslator.cc

History: Sep 27, 2006 - initial version

Member Function Documentation

Member Data Documentation