#include <HITSiStripRawToClustersRoI.h>

Inheritance diagram for HITSiStripRawToClustersRoI:

Public Types
typedef edm::LazyGetter < SiStripCluster >	LazyGetter

typedef SiStripRegionCabling::Position	Position

typedef edm::RefGetter < SiStripCluster >	RefGetter

typedef SiStripRegionCabling::SubDet	SubDet

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Member Functions
virtual void	beginJob ()

virtual void	endJob ()

	HITSiStripRawToClustersRoI (const edm::ParameterSet &)

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

	~HITSiStripRawToClustersRoI ()

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
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
bool	physicalLayer (SubDet &, uint32_t &) const

void	ptracks (const trigger::TriggerFilterObjectWithRefs &, RefGetter &, edm::Handle< LazyGetter > &) const

void	random (RefGetter &, edm::Handle< LazyGetter > &) const

void	taujets (const l1extra::L1JetParticleCollection &, RefGetter &, edm::Handle< LazyGetter > &) const

Private Attributes
std::vector< uint32_t >	allregions_

edm::ESHandle < SiStripRegionCabling >	cabling_

bool	global_

int	nlayers_

bool	ptrack_

double	ptrackEta_

double	ptrackPhi_

bool	random_

double	taujetdeta_

double	taujetdphi_

bool	taujets_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	tok_ptrack_

edm::EDGetTokenT< LazyGetter >	tok_siStrip_

edm::EDGetTokenT < l1extra::L1JetParticleCollection >	tok_tauL1_

Additional Inherited Members
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

Definition at line 36 of file HITSiStripRawToClustersRoI.h.

Member Typedef Documentation

typedef edm::LazyGetter<SiStripCluster> HITSiStripRawToClustersRoI::LazyGetter

Definition at line 40 of file HITSiStripRawToClustersRoI.h.

typedef SiStripRegionCabling::Position HITSiStripRawToClustersRoI::Position

Definition at line 42 of file HITSiStripRawToClustersRoI.h.

typedef edm::RefGetter<SiStripCluster> HITSiStripRawToClustersRoI::RefGetter

Definition at line 41 of file HITSiStripRawToClustersRoI.h.

typedef SiStripRegionCabling::SubDet HITSiStripRawToClustersRoI::SubDet

Definition at line 43 of file HITSiStripRawToClustersRoI.h.

Constructor & Destructor Documentation

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

Definition at line 10 of file HITSiStripRawToClustersRoI.cc.

References edm::ParameterSet::getParameter(), tok_ptrack_, tok_siStrip_, and tok_tauL1_.

                                                                                   :
 
   cabling_(),
   allregions_(),
   nlayers_(conf.getUntrackedParameter<int>("layers",-1)),
   global_(conf.getUntrackedParameter<bool>("doGlobal",true)),
   random_(conf.getUntrackedParameter<bool>("random",false)),
   taujets_(conf.getUntrackedParameter<bool>("useTauJets",false)),
   ptrack_(conf.getUntrackedParameter<bool>("usePixelTracks",true)),
   taujetdeta_(conf.getUntrackedParameter<double>("tjetEtaWindow",0.2)),
   taujetdphi_(conf.getUntrackedParameter<double>("tjetPhiWindow",0.2)),
   ptrackEta_(conf.getUntrackedParameter<double>("ptrackEtaWindow",0.3)),
   ptrackPhi_(conf.getUntrackedParameter<double>("ptrackPhiWindow",0.3))
        
 {
   tok_siStrip_ = consumes<LazyGetter>(conf.getParameter<edm::InputTag>("siStripLazyGetter"));
   tok_tauL1_ = consumes<l1extra::L1JetParticleCollection>(conf.getParameter<edm::InputTag>("l1tauJetLabel"));
   tok_ptrack_ = consumes<trigger::TriggerFilterObjectWithRefs>(conf.getParameter<edm::InputTag>("pixelTrackLabel"));
 
   produces< RefGetter >();
 }

HITSiStripRawToClustersRoI::~HITSiStripRawToClustersRoI ( )

Definition at line 32 of file HITSiStripRawToClustersRoI.cc.

32 {}

Member Function Documentation

void HITSiStripRawToClustersRoI::beginJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 34 of file HITSiStripRawToClustersRoI.cc.

34 {}

void HITSiStripRawToClustersRoI::endJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 36 of file HITSiStripRawToClustersRoI.cc.

36 {}

bool HITSiStripRawToClustersRoI::physicalLayer	(	SubDet &	subdet,
		uint32_t &	layer
	)		const

private

Defines allowed physical layer numbers

Definition at line 91 of file HITSiStripRawToClustersRoI.cc.

References SiStripRegionCabling::ALLLAYERS, nlayers_, SiStripRegionCabling::TEC, SiStripRegionCabling::TIB, SiStripRegionCabling::TIBLAYERS, SiStripRegionCabling::TID, SiStripRegionCabling::TIDLAYERS, and SiStripRegionCabling::TOB.

Referenced by ptracks(), random(), and taujets().

                                                                                     {
   uint32_t signedlayer = 0;
   if (subdet == SiStripRegionCabling::TIB || subdet == SiStripRegionCabling::TID) signedlayer = layer;
   else if (subdet == SiStripRegionCabling::TOB) signedlayer = SiStripRegionCabling::TIBLAYERS + layer;
   else if (subdet == SiStripRegionCabling::TEC) signedlayer = SiStripRegionCabling::TIDLAYERS + layer;
   else signedlayer = SiStripRegionCabling::ALLLAYERS;
 
   return (nlayers_ == -1 || int(signedlayer) < nlayers_) ? true : false;
 }

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

virtual

All regions

Random region number. Starts from 0.

Seeded by L1 taujets.

seeded by isol tracks

Implements edm::EDProducer.

Definition at line 38 of file HITSiStripRawToClustersRoI.cc.

References allregions_, cabling_, runEdmFileComparison::collection, SiStripRegionCabling::elementIndex(), edm::EventSetup::get(), global_, cmsHarvester::index, ptrack_, ptracks(), random(), random_, taujets(), taujets_, tok_ptrack_, tok_siStrip_, and tok_tauL1_.

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

                                                                                     {
 
   setup.get<SiStripRegionCablingRcd>().get(cabling_);
   allregions_.reserve(cabling_->getRegionCabling().size());
   for (uint32_t iregion=0;iregion<cabling_->getRegionCabling().size();iregion++) {
     for (uint32_t isubdet=0;isubdet<cabling_->getRegionCabling()[iregion].size();isubdet++) {
       for (uint32_t ilayer=0;ilayer<cabling_->getRegionCabling()[iregion][isubdet].size();ilayer++) {
         uint32_t index = SiStripRegionCabling::elementIndex(iregion,static_cast<SubDet>(isubdet),ilayer);
         allregions_.push_back(index);
       }
     }
   }
 
   
   edm::Handle< LazyGetter > lazygetter;
   event.getByToken(tok_siStrip_,lazygetter);
 
   
   if (global_) {
     std::auto_ptr<RefGetter> globalrefgetter(new RefGetter(lazygetter,allregions_));
     /*//unpack
       for (RefGetter::const_iterator iReg = globalrefgetter->begin(); iReg != globalrefgetter->end(); iReg++) *iReg;*/
     event.put(globalrefgetter);
     return;
   }
   
   std::auto_ptr<RefGetter> refgetter(new RefGetter(allregions_.size()));
   refgetter->reserve(10000);
   
   
   if (random_) {random(*refgetter,lazygetter);}
   
 
   if (taujets_) {
     edm::Handle<l1extra::L1JetParticleCollection> collection;
     event.getByToken(tok_tauL1_,collection);
     taujets(*collection,*refgetter,lazygetter);
     }
 
 
   if (ptrack_) {
       edm::Handle<trigger::TriggerFilterObjectWithRefs> collection;
       event.getByToken(tok_ptrack_,collection);
       ptracks(*collection,*refgetter,lazygetter);
     }
   
   event.put(refgetter);
 }

void HITSiStripRawToClustersRoI::ptracks	(	const trigger::TriggerFilterObjectWithRefs &	collection,
		RefGetter &	refgetter,
		edm::Handle< LazyGetter > &	lazygetter
	)		const

private

seeding by isolated pixel tracks

Definition at line 133 of file HITSiStripRawToClustersRoI.cc.

References SiStripRegionCabling::ALLLAYERS, SiStripRegionCabling::ALLSUBDETS, cabling_, eta(), trigger::TriggerRefsCollections::getObjects(), AlCaHLTBitMon_ParallelJobs::p, phi, physicalLayer(), position, ptrackEta_, ptrackPhi_, and trigger::TriggerTrack.

Referenced by produce().

                                                                                                                                                             {
 
   std::vector< edm::Ref<reco::IsolatedPixelTrackCandidateCollection> > isoPixTrackRefs;
   collection.getObjects(trigger::TriggerTrack, isoPixTrackRefs);
   for (uint32_t p=0 ; p<isoPixTrackRefs.size(); p++) {
     Position position(isoPixTrackRefs[p]->track()->eta(),isoPixTrackRefs[p]->track()->phi());
     for (uint32_t isubdet = 0; isubdet < SiStripRegionCabling::ALLSUBDETS; isubdet++) {
       SubDet subdet = static_cast<SubDet>(isubdet);
       for (uint32_t ilayer = 0; ilayer < SiStripRegionCabling::ALLLAYERS; ilayer++) {
     if (!physicalLayer(subdet,ilayer)) break;
     cabling_->updateSiStripRefGetter<SiStripCluster>(refgetter,lazygetter,position,ptrackEta_,ptrackPhi_,subdet,ilayer);
       }
     }
     
   }
 }

void HITSiStripRawToClustersRoI::random	(	RefGetter &	refgetter,
		edm::Handle< LazyGetter > &	lazygetter
	)		const

private

Defines regions of interest randomly

Definition at line 101 of file HITSiStripRawToClustersRoI.cc.

References SiStripRegionCabling::ALLLAYERS, SiStripRegionCabling::ALLSUBDETS, cabling_, SiStripRegionCabling::elementIndex(), physicalLayer(), dumpFill::required, and pileupDistInMC::total.

Referenced by produce().

                                                                                                      {
   
   uint32_t total = cabling_->getRegionCabling().size();
   uint32_t required = static_cast<uint32_t>(CLHEP::RandFlat::shoot()*(total+1));
   for (uint32_t iregion = 0; iregion < required; iregion++) {
     for (uint32_t isubdet = 0; isubdet < SiStripRegionCabling::ALLSUBDETS; isubdet++) {
       SubDet subdet = static_cast<SubDet>(isubdet);
       for (uint32_t ilayer = 0; ilayer < SiStripRegionCabling::ALLLAYERS; ilayer++) {
     if (!physicalLayer(subdet,ilayer)) break;
     cabling_->updateSiStripRefGetter<SiStripCluster>(refgetter,lazygetter,SiStripRegionCabling::elementIndex(iregion,subdet,ilayer));
       }
     }
   }
 }

void HITSiStripRawToClustersRoI::taujets	(	const l1extra::L1JetParticleCollection &	collection,
		RefGetter &	refgetter,
		edm::Handle< LazyGetter > &	lazygetter
	)		const

private

Defines regions of interest by taujets

Definition at line 117 of file HITSiStripRawToClustersRoI.cc.

References SiStripRegionCabling::ALLLAYERS, SiStripRegionCabling::ALLSUBDETS, cabling_, physicalLayer(), position, taujetdeta_, and taujetdphi_.

Referenced by produce().

                                                                                                                                                         {
 
   l1extra::L1JetParticleCollection::const_iterator icollection = collection.begin();
   for (; icollection!=collection.end(); icollection++) {
     Position position(icollection->eta(),icollection->phi());
     for (uint32_t isubdet = 0; isubdet < SiStripRegionCabling::ALLSUBDETS; isubdet++) {
       SubDet subdet = static_cast<SubDet>(isubdet);
       for (uint32_t ilayer = 0; ilayer < SiStripRegionCabling::ALLLAYERS; ilayer++) {
     if (!physicalLayer(subdet,ilayer)) break;
     cabling_->updateSiStripRefGetter<SiStripCluster>(refgetter,lazygetter,position,taujetdeta_,taujetdphi_,subdet,ilayer);
       }
     }    
   }
 }