SeedingLayerSetsBuilder Class Reference

#include <SeedingLayerSetsBuilder.h>

Classes
struct	LayerSpec

Public Types
typedef unsigned short	LayerSetIndex
using	SeedingLayerId = std::tuple< GeomDetEnumerators::SubDetector, ctfseeding::SeedingLayer::Side, int >

Public Member Functions
bool	check (const edm::EventSetup &es)
void	hits (const edm::Event &ev, const edm::EventSetup &es, std::vector< unsigned int > &indices, ctfseeding::SeedingLayer::Hits &hits) const
const std::vector< const DetLayer * > &	layerDets () const
const std::vector< std::string > &	layerNames () const
ctfseeding::SeedingLayerSets	layers (const edm::EventSetup &es)
const std::vector< LayerSetIndex > &	layerSetIndices () const
unsigned short	numberOfLayers () const
unsigned short	numberOfLayersInSet () const
	SeedingLayerSetsBuilder ()
	SeedingLayerSetsBuilder (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)
	SeedingLayerSetsBuilder (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
void	updateEventSetup (const edm::EventSetup &es)
	~SeedingLayerSetsBuilder ()

Static Public Member Functions
static std::vector< std::vector< std::string > >	layerNamesInSets (const std::vector< std::string > &namesPSet)
static SeedingLayerId	nameToEnumId (const std::string &name)

Private Member Functions
edm::ParameterSet	layerConfig (const std::string &nameLayer, const edm::ParameterSet &cfg) const

Private Attributes
edm::ESWatcher< TrackerRecoGeometryRecord >	geometryWatcher_
std::vector< const DetLayer * >	theLayerDets
std::vector< std::string >	theLayerNames
std::vector< LayerSpec >	theLayers
std::vector< LayerSetIndex >	theLayerSetIndices
unsigned short	theNumberOfLayersInSet
std::vector< const TransientTrackingRecHitBuilder * >	theTTRHBuilders
edm::ESWatcher< TransientRecHitRecord >	trhWatcher_

Detailed Description

Definition at line 20 of file SeedingLayerSetsBuilder.h.

Member Typedef Documentation

typedef unsigned short SeedingLayerSetsBuilder::LayerSetIndex

Definition at line 34 of file SeedingLayerSetsBuilder.h.

using SeedingLayerSetsBuilder::SeedingLayerId = std::tuple<GeomDetEnumerators::SubDetector, ctfseeding::SeedingLayer::Side, int>

Definition at line 43 of file SeedingLayerSetsBuilder.h.

Constructor & Destructor Documentation

SeedingLayerSetsBuilder::SeedingLayerSetsBuilder

(

)

Definition at line 182 of file SeedingLayerSetsBuilder.cc.

{}

SeedingLayerSetsBuilder::SeedingLayerSetsBuilder	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC
	)

Definition at line 186 of file SeedingLayerSetsBuilder.cc.

References Exception, spr::find(), runEdmFileComparison::found, edm::ParameterSet::getParameter(), layerConfig(), layerNamesInSets(), hpstanc_transforms::max, harvestTrackValidationPlots::str, AlCaHLTBitMon_QueryRunRegistry::string, theLayerDets, theLayerNames, theLayers, theLayerSetIndices, theNumberOfLayersInSet, and theTTRHBuilders.

{
  std::vector<std::string> namesPset = cfg.getParameter<std::vector<std::string> >("layerList");
  std::vector<std::vector<std::string> > layerNamesInSets = this->layerNamesInSets(namesPset);

  // debug printout of layers
  typedef std::vector<std::string>::const_iterator IS;
  typedef std::vector<std::vector<std::string> >::const_iterator IT;
  std::ostringstream str;
  // The following should not be set to cout
//  for (IT it = layerNamesInSets.begin(); it != layerNamesInSets.end(); it++) {
//    str << "SET: ";
//    for (IS is = it->begin(); is != it->end(); is++)  str << *is <<" ";  
//    str << std::endl;
//  }
//  std::cout << str.str() << std::endl;
  if(layerNamesInSets.size() == 0)
    theNumberOfLayersInSet = 0;
  else
    theNumberOfLayersInSet = layerNamesInSets[0].size();


  for (IT it = layerNamesInSets.begin(); it != layerNamesInSets.end(); it++) {
    if(it->size() != theNumberOfLayersInSet)
      throw cms::Exception("Configuration") << "Assuming all SeedingLayerSets to have same number of layers. LayerSet " << (it-layerNamesInSets.begin()) << " has " << it->size() << " while 0th has " << theNumberOfLayersInSet;
    for(const std::string& layerName: *it) {
      auto found = std::find(theLayerNames.begin(), theLayerNames.end(), layerName);
      unsigned short layerIndex = 0;
      if(found != theLayerNames.end()) {
        layerIndex = found-theLayerNames.begin();
      }
      else {
        if(std::numeric_limits<unsigned short>::max() == theLayerNames.size()) {
          throw cms::Exception("Assert") << "Too many layers in " << __FILE__ << ":" << __LINE__ << ", we may have to enlarge the index type from unsigned short to unsigned int";
        }

        layerIndex = theLayers.size();
        theLayers.emplace_back(theLayerNames.size(), layerName, layerConfig(layerName, cfg), iC);
        theLayerNames.push_back(layerName);
      }
      theLayerSetIndices.push_back(layerIndex);
    }
  }
  theLayerDets.resize(theLayers.size());
  theTTRHBuilders.resize(theLayers.size());

  // debug printout
  // The following should not be set to cout
  //for(const LayerSpec& layer: theLayers) {
  //  std::cout << layer.print(theLayerNames) << std::endl;
  //}
}

SeedingLayerSetsBuilder::SeedingLayerSetsBuilder	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &&	iC
	)

Definition at line 183 of file SeedingLayerSetsBuilder.cc.

                                                                                                      :
  SeedingLayerSetsBuilder(cfg, iC)
{}

SeedingLayerSetsBuilder::~SeedingLayerSetsBuilder

(

)

Definition at line 239 of file SeedingLayerSetsBuilder.cc.

{}

Member Function Documentation

bool SeedingLayerSetsBuilder::check

(

const edm::EventSetup &

es

)

Definition at line 356 of file SeedingLayerSetsBuilder.cc.

References edm::ESWatcher< T >::check(), geometryWatcher_, and trhWatcher_.

Referenced by QuadrupletSeedMerger::mergeTriplets(), and SeedingLayersEDProducer::produce().

                                                           {
  // We want to evaluate both in the first invocation (to properly
  // initialize ESWatcher), and this way we avoid one branch compared
  // to || (should be tiny effect)
  return geometryWatcher_.check(es) | trhWatcher_.check(es);
}

void SeedingLayerSetsBuilder::hits	(	const edm::Event &	ev,
		const edm::EventSetup &	es,
		std::vector< unsigned int > &	indices,
		ctfseeding::SeedingLayer::Hits &	hits
	)		const

Definition at line 365 of file SeedingLayerSetsBuilder.cc.

References mps_fire::i, eostools::move(), theLayers, theTTRHBuilders, and tmp.

Referenced by ntupleDataFormat._TrackingParticleMatchAdaptor::bestMatchingTrackingParticle(), and SeedingLayersEDProducer::produce().

                                                                                              {
  indices.reserve(theLayers.size());
  for(unsigned int i=0; i<theLayers.size(); ++i) {
    // The index of the first hit of this layer
    indices.push_back(hits.size());

    // Obtain and copy the hits
    ctfseeding::SeedingLayer::Hits && tmp = theLayers[i].extractor->hits((const TkTransientTrackingRecHitBuilder &)(*theTTRHBuilders[i]), ev, es);
    std::move(tmp.begin(), tmp.end(), std::back_inserter(hits));
  }
}

edm::ParameterSet SeedingLayerSetsBuilder::layerConfig ( const std::string &  nameLayer, const edm::ParameterSet &  cfg  ) const

private

Definition at line 241 of file SeedingLayerSetsBuilder.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), dataset::name, and mps_fire::result.

Referenced by SeedingLayerSetsBuilder().

{
  edm::ParameterSet result;
   
  for (string::size_type iEnd=nameLayer.size(); iEnd > 0; --iEnd) {
    string name = nameLayer.substr(0,iEnd);
    if (cfg.exists(name)) return cfg.getParameter<edm::ParameterSet>(name);
  } 
  edm::LogError("SeedingLayerSetsBuilder") <<"configuration for layer: "<<nameLayer<<" not found, job will probably crash!";
  return result;
}

const std::vector<const DetLayer *>& SeedingLayerSetsBuilder::layerDets ( ) const

inline

Definition at line 40 of file SeedingLayerSetsBuilder.h.

References ev, and hfClusterShapes_cfi::hits.

Referenced by SeedingLayersEDProducer::produce().

{ return theLayerDets; }

const std::vector<std::string>& SeedingLayerSetsBuilder::layerNames ( ) const

inline

Definition at line 39 of file SeedingLayerSetsBuilder.h.

Referenced by SeedingLayersEDProducer::produce().

{ return theLayerNames; }

vector< vector< string > > SeedingLayerSetsBuilder::layerNamesInSets ( const std::vector< std::string > &  namesPSet )

static

Definition at line 253 of file SeedingLayerSetsBuilder.cc.

References geometryCSVtoXML::line, and mps_fire::result.

Referenced by MTVHistoProducerAlgoForTracker::MTVHistoProducerAlgoForTracker(), and SeedingLayerSetsBuilder().

{
  std::vector<std::vector<std::string> > result; 
  for (std::vector<std::string>::const_iterator is=namesPSet.begin(); is < namesPSet.end(); ++is) {
    vector<std::string> layersInSet;
    string line = *is;
    string::size_type pos=0;
    while (pos != string::npos ) {
      pos=line.find("+");
      string layer = line.substr(0,pos);
      layersInSet.push_back(layer);
      line=line.substr(pos+1,string::npos); 
    }
    result.push_back(layersInSet);
  }
  return result;
}

SeedingLayerSets SeedingLayerSetsBuilder::layers

(

const edm::EventSetup &

es

)

Definition at line 335 of file SeedingLayerSetsBuilder.cc.

References SeedingLayerSetsBuilder::LayerSpec::extractor, mps_fire::i, gen::n, mps_fire::result, theLayerDets, theLayerNames, theLayers, theLayerSetIndices, theNumberOfLayersInSet, theTTRHBuilders, and updateEventSetup().

Referenced by QuadrupletSeedMerger::mergeTriplets().

{
  updateEventSetup(es);

  typedef std::vector<SeedingLayer> Set;
  SeedingLayerSets  result;

  for(size_t i=0, n=theLayerSetIndices.size(); i<n; i += theNumberOfLayersInSet) {
    Set set;
    for(size_t j=0; j<theNumberOfLayersInSet; ++j) {
      const unsigned short layerIndex = theLayerSetIndices[i+j];
      const LayerSpec& layer = theLayers[layerIndex];
      const DetLayer *detLayer = theLayerDets[layerIndex];

      set.push_back( SeedingLayer( theLayerNames[layerIndex], layerIndex, detLayer, theTTRHBuilders[layerIndex], layer.extractor.get()));
    }
    result.push_back(set);
  }
  return result;
}

const std::vector<LayerSetIndex>& SeedingLayerSetsBuilder::layerSetIndices ( ) const

inline

Definition at line 36 of file SeedingLayerSetsBuilder.h.

Referenced by SeedingLayersEDProducer::produce().

{ return theLayerSetIndices; }

SeedingLayerSetsBuilder::SeedingLayerId SeedingLayerSetsBuilder::nameToEnumId ( const std::string &  name )

static

Definition at line 32 of file SeedingLayerSetsBuilder.cc.

References ctfseeding::SeedingLayer::Barrel, diffTreeTool::index, GeomDetEnumerators::invalidDet, ctfseeding::SeedingLayer::NegEndcap, GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, ctfseeding::SeedingLayer::PosEndcap, GeomDetEnumerators::TEC, GeomDetEnumerators::TIB, GeomDetEnumerators::TID, and GeomDetEnumerators::TOB.

Referenced by SeedingLayerSetsBuilder::LayerSpec::LayerSpec(), and MTVHistoProducerAlgoForTracker::MTVHistoProducerAlgoForTracker().

                                                                                                 {
    GeomDetEnumerators::SubDetector subdet = GeomDetEnumerators::invalidDet;
    SeedingLayer::Side side = SeedingLayer::Barrel;
    int idLayer = 0;

    size_t index;
    //
    // BPIX
    //
    if ((index = name.find("BPix")) != string::npos) {
      subdet = GeomDetEnumerators::PixelBarrel;
      side = SeedingLayer::Barrel;
      idLayer = atoi(name.substr(index+4,1).c_str());
    }
    //
    // FPIX
    //
    else if ((index = name.find("FPix")) != string::npos) {
      subdet = GeomDetEnumerators::PixelEndcap;
      idLayer = atoi(name.substr(index+4).c_str());
      if ( name.find("pos") != string::npos ) {
        side = SeedingLayer::PosEndcap;
      } else {
        side = SeedingLayer::NegEndcap;
      }
    }
    //
    // TIB
    //
    else if ((index = name.find("TIB")) != string::npos) {
      subdet = GeomDetEnumerators::TIB;
      side = SeedingLayer::Barrel;
      idLayer = atoi(name.substr(index+3,1).c_str());
    }
    //
    // TID
    //
    else if ((index = name.find("TID")) != string::npos) {
      subdet = GeomDetEnumerators::TID;
      idLayer = atoi(name.substr(index+3,1).c_str());
      if ( name.find("pos") != string::npos ) {
        side = SeedingLayer::PosEndcap;
      } else {
        side = SeedingLayer::NegEndcap;
      }
    }
    //
    // TOB
    //
    else if ((index = name.find("TOB")) != string::npos) {
      subdet = GeomDetEnumerators::TOB;
      side = SeedingLayer::Barrel;
      idLayer = atoi(name.substr(index+3,1).c_str());
    }
    //
    // TEC
    //
    else if ((index = name.find("TEC")) != string::npos) {
      subdet = GeomDetEnumerators::TEC;
      idLayer = atoi(name.substr(index+3,1).c_str());
      if ( name.find("pos") != string::npos ) {
        side = SeedingLayer::PosEndcap;
      } else {
        side = SeedingLayer::NegEndcap;
      }
    }
    return std::make_tuple(subdet, side, idLayer);
}

unsigned short SeedingLayerSetsBuilder::numberOfLayers ( ) const

inline

Definition at line 38 of file SeedingLayerSetsBuilder.h.

{ return theLayers.size(); }

unsigned short SeedingLayerSetsBuilder::numberOfLayersInSet ( ) const

inline

Definition at line 35 of file SeedingLayerSetsBuilder.h.

Referenced by SeedingLayersEDProducer::produce().

{ return theNumberOfLayersInSet; }

void SeedingLayerSetsBuilder::updateEventSetup

(

const edm::EventSetup &

es

)

Definition at line 271 of file SeedingLayerSetsBuilder.cc.

References GeometricSearchTracker::barrelLayers(), Exception, edm::EventSetup::get(), SeedingLayerSetsBuilder::LayerSpec::hitBuilder, mps_fire::i, SeedingLayerSetsBuilder::LayerSpec::idLayer, diffTreeTool::index, gen::n, SeedingLayerSetsBuilder::LayerSpec::nameIndex, GeometricSearchTracker::negForwardLayers(), GeometricSearchTracker::negTecLayers(), GeometricSearchTracker::negTidLayers(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, GeometricSearchTracker::posForwardLayers(), GeometricSearchTracker::posTecLayers(), GeometricSearchTracker::posTidLayers(), edm::ESHandle< T >::product(), SeedingLayerSetsBuilder::LayerSpec::side, SeedingLayerSetsBuilder::LayerSpec::subdet, GeomDetEnumerators::TEC, theLayerDets, theLayerNames, theLayers, theTTRHBuilders, GeomDetEnumerators::TIB, GeometricSearchTracker::tibLayers(), GeomDetEnumerators::TID, GeomDetEnumerators::TOB, GeometricSearchTracker::tobLayers(), and mixOne_simraw_on_sim_cfi::tracker.

Referenced by layers(), and SeedingLayersEDProducer::produce().

                                                                      {
  edm::ESHandle<GeometricSearchTracker> htracker;
  es.get<TrackerRecoGeometryRecord>().get( htracker );
  const GeometricSearchTracker& tracker = *htracker;

  const std::vector<BarrelDetLayer const*>&  bpx  = tracker.barrelLayers();
  const std::vector<BarrelDetLayer const*>&  tib  = tracker.tibLayers();
  const std::vector<BarrelDetLayer const*>&  tob  = tracker.tobLayers();

  const std::vector<ForwardDetLayer const*>& fpx_pos = tracker.posForwardLayers();
  const std::vector<ForwardDetLayer const*>& tid_pos = tracker.posTidLayers();
  const std::vector<ForwardDetLayer const*>& tec_pos = tracker.posTecLayers();

  const std::vector<ForwardDetLayer const*>& fpx_neg = tracker.negForwardLayers();
  const std::vector<ForwardDetLayer const*>& tid_neg = tracker.negTidLayers();
  const std::vector<ForwardDetLayer const*>& tec_neg = tracker.negTecLayers();

  for(size_t i=0, n=theLayers.size(); i<n; ++i) {
    const LayerSpec& layer = theLayers[i];
    const DetLayer * detLayer = nullptr;
    int index = layer.idLayer-1;

    if (layer.subdet == GeomDetEnumerators::PixelBarrel) {
      detLayer = bpx[index];
    }
    else if (layer.subdet == GeomDetEnumerators::PixelEndcap) {
      if (layer.side == SeedingLayer::PosEndcap) {
        detLayer = fpx_pos[index];
      } else {
        detLayer = fpx_neg[index];
      }
    }
    else if (layer.subdet == GeomDetEnumerators::TIB) {
      detLayer = tib[index];
    }
    else if (layer.subdet == GeomDetEnumerators::TID) {
      if (layer.side == SeedingLayer::PosEndcap) {
        detLayer = tid_pos[index];
      } else {
        detLayer = tid_neg[index];
      }
    }
    else if (layer.subdet == GeomDetEnumerators::TOB) {
      detLayer = tob[index];
    }
    else if (layer.subdet == GeomDetEnumerators::TEC) {
      if (layer.side == SeedingLayer::PosEndcap) {
        detLayer = tec_pos[index];
      } else {
        detLayer = tec_neg[index];
      }
    }
    else {
      throw cms::Exception("Configuration") << "Did not find DetLayer for layer " << theLayerNames[layer.nameIndex];
    }

    edm::ESHandle<TransientTrackingRecHitBuilder> builder;
    es.get<TransientRecHitRecord>().get(layer.hitBuilder, builder);

    theLayerDets[i] = detLayer;
    theTTRHBuilders[i] = builder.product();
  }
}

Member Data Documentation

edm::ESWatcher<TrackerRecoGeometryRecord> SeedingLayerSetsBuilder::geometryWatcher_

private

Definition at line 50 of file SeedingLayerSetsBuilder.h.

Referenced by check().

std::vector<const DetLayer *> SeedingLayerSetsBuilder::theLayerDets

private

Definition at line 71 of file SeedingLayerSetsBuilder.h.

Referenced by layers(), SeedingLayerSetsBuilder(), and updateEventSetup().

std::vector<std::string> SeedingLayerSetsBuilder::theLayerNames

private

Definition at line 70 of file SeedingLayerSetsBuilder.h.

Referenced by layers(), SeedingLayerSetsBuilder(), and updateEventSetup().

std::vector<LayerSpec> SeedingLayerSetsBuilder::theLayers

private

Definition at line 73 of file SeedingLayerSetsBuilder.h.

Referenced by hits(), layers(), SeedingLayerSetsBuilder(), and updateEventSetup().

std::vector<LayerSetIndex> SeedingLayerSetsBuilder::theLayerSetIndices

private

Definition at line 69 of file SeedingLayerSetsBuilder.h.

Referenced by layers(), and SeedingLayerSetsBuilder().

unsigned short SeedingLayerSetsBuilder::theNumberOfLayersInSet

private

Definition at line 68 of file SeedingLayerSetsBuilder.h.

Referenced by layers(), and SeedingLayerSetsBuilder().

std::vector<const TransientTrackingRecHitBuilder *> SeedingLayerSetsBuilder::theTTRHBuilders

private

Definition at line 72 of file SeedingLayerSetsBuilder.h.

Referenced by hits(), layers(), SeedingLayerSetsBuilder(), and updateEventSetup().

edm::ESWatcher<TransientRecHitRecord> SeedingLayerSetsBuilder::trhWatcher_

private

Definition at line 51 of file SeedingLayerSetsBuilder.h.

Referenced by check().