#include <SiClusterTranslator.h>

Inheritance diagram for SiClusterTranslator:

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

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

	SiClusterTranslator (const edm::ParameterSet &conf)

virtual	~SiClusterTranslator ()

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

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) 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	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
edm::InputTag	fastTrackerClusterCollectionTag_

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

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

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

                                                                     :
   fastTrackerClusterCollectionTag_(conf.getParameter<edm::InputTag>("fastTrackerClusterCollectionTag"))
 {
   produces<edmNew::DetSetVector<SiStripCluster> >();
   produces<edmNew::DetSetVector<SiPixelCluster> >();
   produces<edm::DetSetVector<StripDigiSimLink> >();
   produces<edm::DetSetVector<PixelDigiSimLink> >();
 }

SiClusterTranslator::~SiClusterTranslator ( )

virtual

Definition at line 66 of file SiClusterTranslator.cc.

66 {}

Member Function Documentation

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

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 69 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 296 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 276 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
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 78 of file SiClusterTranslator.cc.

References counter, edm::DetSet< T >::data, edm::hlt::Exception, fastTrackerClusterCollectionTag_, edm::EventSetup::get(), edm::Event::getByLabel(), 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.getByLabel(fastTrackerClusterCollectionTag_, theFastClusters);
   
   edm::ESHandle<TrackerGeometry> tkgeom;
   es.get<TrackerDigiGeometryRecord>().get( tkgeom ); 
   const TrackerGeometry &tracker(*tkgeom);
   
   edm::ESHandle<PixelClusterParameterEstimator> pixelCPE;
   es.get<TkPixelCPERecord>().get("FastPixelCPE",pixelCPE);
   //  auto pixelcpe = pixelCPE->clone();
   
   edm::ESHandle<StripClusterParameterEstimator> stripCPE;
   es.get<TkStripCPERecord>().get("FastStripCPE", stripCPE); 
   // auto stripcpe = stripCPE->clone();
     
   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;
   
   
   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=(const 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.
       // FIXME this cannot be done. CPE is global and cannot be modified.
       // this information has to be encoded in the Pixel cluster itself w/o changing it!
       // CPE has just to retrieve it and return.
       //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;
       
       // FIXME VI there is something wrong here
       // each "digi" is a new strip that this algorithm will add at the right of "strip_num"
       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=(const RectangularStripTopology&)stripDet->type().topology();
     strip_num = (uint16_t)topol.strip(position);
       } else if ((subdet == 4) || (subdet == 6)) {
     const RadialStripTopology& topol=(const RadialStripTopology&)stripDet->type().topology();
     strip_num = (uint16_t)topol.strip(position);
       }
       
       //Filling Strip CPE with info.
       // FIXME this cannot be done. CPE is global and cannot be modified.
       // this information has to be encoded in the Strip cluster itself w/o modifing it!
       // CPE has just to retrieve it and return.
       // stripcpe->enterLocalParameters(det.rawId(), strip_num, std::make_pair(position,error));
       
       //Creating a new strip cluster
       SiStripCluster temporaryStripCluster(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

edm::InputTag SiClusterTranslator::fastTrackerClusterCollectionTag_

private

Definition at line 60 of file SiClusterTranslator.h.

Referenced by produce().

const TrackerGeometry* SiClusterTranslator::geometry

private

Definition at line 59 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 52 of file SiClusterTranslator.h.

Referenced by produce().

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

private

Definition at line 51 of file SiClusterTranslator.h.

Referenced by produce().

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

private

Definition at line 50 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