Inheritance diagram for PixelClusterShapeExtractor:

Public Member Functions
void	analyze (edm::StreamID, const edm::Event &evt, const edm::EventSetup &) const override

void	endJob () override

	PixelClusterShapeExtractor (const edm::ParameterSet &pset)

Public Member Functions inherited from edm::global::EDAnalyzer<>
	EDAnalyzer ()=default

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Types
using	Lock = std::unique_lock< std::mutex >

Private Member Functions
void	analyzeRecTracks (const edm::Event &ev, const edm::EventSetup &es) const

void	analyzeSimHits (const edm::Event &ev, const edm::EventSetup &es) const

bool	checkSimHits (const TrackingRecHit &recHit, TrackerHitAssociator const &theAssociator, PSimHit &simHit, pair< unsigned int, float > &key, unsigned int &ss) const

void	init ()

bool	isSuitable (const PSimHit &simHit, const GeomDetUnit &gdu) const

void	processPixelRecHits (SiPixelRecHitCollection::DataContainer const &recHits, TrackerHitAssociator const &theAssociator, ClusterShapeHitFilter const &theFilter, SiPixelClusterShapeCache const &clusterShapeCache, const TrackerTopology &tkTpl) const

void	processRec (const SiPixelRecHit &recHit, ClusterShapeHitFilter const &theFilter, LocalVector ldir, const SiPixelClusterShapeCache &clusterShapeCache, const vector< TH2F * > &histo) const

void	processSim (const SiPixelRecHit &recHit, ClusterShapeHitFilter const &theClusterFilter, const PSimHit &simHit, const SiPixelClusterShapeCache &clusterShapeCache, const vector< TH2F * > &histo) const

Private Attributes
const edm::EDGetTokenT< SiPixelClusterShapeCache >	clusterShapeCache_token

TFile *	file

const bool	hasRecTracks

const bool	hasSimHits

std::vector< TH2F * >	hrpc

std::vector< TH2F * >	hspc

const bool	noBPIX1

const edm::EDGetTokenT< edmNew::DetSetVector< SiPixelRecHit > >	pixelRecHits_token

std::unique_ptr< std::mutex[]>	theMutex

const TrackerHitAssociator::Config	trackerHitAssociatorConfig_

const edm::EDGetTokenT< reco::TrackCollection >	tracks_token

Additional Inherited Members
Public Types inherited from edm::global::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::global::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 52 of file PixelClusterShapeExtractor.cc.

Member Typedef Documentation

using PixelClusterShapeExtractor::Lock = std::unique_lock<std::mutex>

private

Definition at line 103 of file PixelClusterShapeExtractor.cc.

Constructor & Destructor Documentation

PixelClusterShapeExtractor::PixelClusterShapeExtractor ( const edm::ParameterSet & pset )

explicit

Definition at line 143 of file PixelClusterShapeExtractor.cc.

References file, and init().

     : hasSimHits(pset.getParameter<bool>("hasSimHits")),
       hasRecTracks(pset.getParameter<bool>("hasRecTracks")),
       noBPIX1(pset.getParameter<bool>("noBPIX1")),
       tracks_token(hasRecTracks ? consumes<reco::TrackCollection>(pset.getParameter<edm::InputTag>("tracks"))
                                 : edm::EDGetTokenT<reco::TrackCollection>()),
       pixelRecHits_token(consumes<edmNew::DetSetVector<SiPixelRecHit>>(edm::InputTag("siPixelRecHits"))),
       clusterShapeCache_token(
           consumes<SiPixelClusterShapeCache>(pset.getParameter<edm::InputTag>("clusterShapeCacheSrc"))),
       trackerHitAssociatorConfig_(pset, consumesCollector()) {
   file = new TFile("clusterShape.root", "RECREATE");
   file->cd();
   init();
 }

Member Function Documentation

void PixelClusterShapeExtractor::analyze	(	edm::StreamID	,
		const edm::Event &	evt,
		const edm::EventSetup &	es
	)		const

override

Definition at line 380 of file PixelClusterShapeExtractor.cc.

References analyzeRecTracks(), analyzeSimHits(), DEFINE_FWK_MODULE, hasRecTracks, hasSimHits, and LogTrace.

                                                                                                        {
   if (hasSimHits) {
     LogTrace("MinBiasTracking") << " [ClusterShape] analyze simHits, recHits";
     analyzeSimHits(ev, es);
   }
 
   if (hasRecTracks) {
     LogTrace("MinBiasTracking") << " [ClusterShape] analyze recHits on recTracks";
     analyzeRecTracks(ev, es);
   }
 }

void PixelClusterShapeExtractor::analyzeRecTracks	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)		const

private

Definition at line 329 of file PixelClusterShapeExtractor.cc.

References SiPixelPhase1OnlineDQM_cff::clusterShapeCache, clusterShapeCache_token, edm::EventSetup::get(), edm::Event::getByToken(), h, hcalSimParameters_cfi::hb, reco::TrackBase::highPurity, hrpc, fastTrackerRecHitType::isPixel(), noBPIX1, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, processRec(), edm::ESHandle< T >::product(), rpcPointValidation_cfi::recHit, HLT_2018_cff::track, PDWG_EXOHSCP_cff::tracks, and tracks_token.

Referenced by analyze().

                                                                                                    {
   edm::ESHandle<ClusterShapeHitFilter> shape;
   es.get<CkfComponentsRecord>().get("ClusterShapeHitFilter", shape);
   auto const& theClusterShape = *shape.product();
 
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   auto const& tkTpl = *tTopoHandle;
 
   // Get tracks
   edm::Handle<reco::TrackCollection> tracks;
   ev.getByToken(tracks_token, tracks);
 
   edm::Handle<SiPixelClusterShapeCache> clusterShapeCache;
   ev.getByToken(clusterShapeCache_token, clusterShapeCache);
 
   for (auto const& track : *tracks) {
     if (!track.quality(reco::Track::highPurity))
       continue;
     if (track.numberOfValidHits() < 8)
       continue;
     auto const& trajParams = track.extra()->trajParams();
     assert(trajParams.size() == track.recHitsSize());
     auto hb = track.recHitsBegin();
     for (unsigned int h = 0; h < track.recHitsSize(); h++) {
       auto recHit = *(hb + h);
       if (!recHit->isValid())
         continue;
       auto id = recHit->geographicalId();
       if (noBPIX1 && tkTpl.pxbLayer(id) == 1)
         continue;
 
       // check that we are in the pixel
       auto subdetid = (id.subdetId());
       bool isPixel = subdetid == PixelSubdetector::PixelBarrel || subdetid == PixelSubdetector::PixelEndcap;
 
       auto const& ltp = trajParams[h];
       auto ldir = ltp.momentum() / ltp.momentum().mag();
 
       if (isPixel) {
         // Pixel
         const SiPixelRecHit* pixelRecHit = dynamic_cast<const SiPixelRecHit*>(recHit);
 
         if (pixelRecHit != nullptr)
           processRec(*pixelRecHit, theClusterShape, ldir, *clusterShapeCache, hrpc);
       }
     }
   }
 }

void PixelClusterShapeExtractor::analyzeSimHits	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)		const

private

Definition at line 300 of file PixelClusterShapeExtractor.cc.

References SiPixelPhase1OnlineDQM_cff::clusterShapeCache, clusterShapeCache_token, coll, edmNew::DetSetVector< T >::data(), ev, edm::EventSetup::get(), edm::Event::getByToken(), pixelRecHits_token, processPixelRecHits(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), FastTrackerRecHitMaskProducer_cfi::recHits, and trackerHitAssociatorConfig_.

Referenced by analyze().

                                                                                                  {
   edm::ESHandle<ClusterShapeHitFilter> shape;
   es.get<CkfComponentsRecord>().get("ClusterShapeHitFilter", shape);
   auto const& theClusterShape = *shape.product();
 
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   auto const& tkTpl = *tTopoHandle;
 
   edm::Handle<SiPixelClusterShapeCache> clusterShapeCache;
   ev.getByToken(clusterShapeCache_token, clusterShapeCache);
 
   // Get associator
   auto theHitAssociator = std::make_unique<TrackerHitAssociator>(ev, trackerHitAssociatorConfig_);
 
   // Pixel hits
   {
     edm::Handle<SiPixelRecHitCollection> coll;
     ev.getByToken(pixelRecHits_token, coll);
 
     edm::Handle<SiPixelClusterShapeCache> clusterShapeCache;
     ev.getByToken(clusterShapeCache_token, clusterShapeCache);
 
     auto const& recHits = coll.product()->data();
     processPixelRecHits(recHits, *theHitAssociator, theClusterShape, *clusterShapeCache, tkTpl);
   }
 }

bool PixelClusterShapeExtractor::checkSimHits	(	const TrackingRecHit &	recHit,
		TrackerHitAssociator const &	theAssociator,
		PSimHit &	simHit,
		pair< unsigned int, float > &	key,
		unsigned int &	ss
	)		const

private

Definition at line 236 of file PixelClusterShapeExtractor.cc.

References TrackerHitAssociator::associateHit(), TrackingRecHit::detUnit(), isSuitable(), FastTrackerRecHitCombiner_cfi::simHits, PSimHit::timeOfFlight(), and PSimHit::trackId().

Referenced by processPixelRecHits().

                                                                       {
   auto const& simHits = theHitAssociator.associateHit(recHit);
 
   //std::cout << "simHits.size() = " << simHits.size() << std::endl;
   for (auto const& sh : simHits) {
     if (isSuitable(sh, *recHit.detUnit())) {
       simHit = sh;
       key = {simHit.trackId(), simHit.timeOfFlight()};
       ss = simHits.size();
       return true;
     }
   }
 
   return false;
 }

void PixelClusterShapeExtractor::endJob ( void )

overridevirtual

Reimplemented from edm::global::EDAnalyzerBase.

Definition at line 159 of file PixelClusterShapeExtractor.cc.

References file, h, hrpc, and hspc.

                                         {
   file->cd();
 
   // simulated
   for (auto h = hspc.begin(); h != hspc.end(); h++)
     (*h)->Write();
 
   // reconstructed
   for (auto h = hrpc.begin(); h != hrpc.end(); h++)
     (*h)->Write();
 
   file->Close();
 }

void PixelClusterShapeExtractor::init ( void )

private

Definition at line 110 of file PixelClusterShapeExtractor.cc.

References exMax, eyMax, and mutex.

Referenced by PixelClusterShapeExtractor().

                                       {
   // Declare histograms
   char histName[256];
 
   // pixel
   for (int subdet = 0; subdet <= 1; subdet++) {
     for (int ex = 0; ex <= exMax; ex++)
       for (int ey = 0; ey <= eyMax; ey++) {
         sprintf(histName, "hspc_%d_%d_%d", subdet, ex, ey);
         hspc.push_back(new TH2F(histName,
                                 histName,
                                 10 * 2 * (exMax + 2),
                                 -(exMax + 2),
                                 (exMax + 2),
                                 10 * 2 * (eyMax + 2),
                                 -(eyMax + 2),
                                 (eyMax + 2)));
 
         sprintf(histName, "hrpc_%d_%d_%d", subdet, ex, ey);
         hrpc.push_back(new TH2F(histName,
                                 histName,
                                 10 * 2 * (exMax + 2),
                                 -(exMax + 2),
                                 (exMax + 2),
                                 10 * 2 * (eyMax + 2),
                                 -(eyMax + 2),
                                 (eyMax + 2)));
       }
   }
   theMutex.reset(new std::mutex[hspc.size()]);
 }

bool PixelClusterShapeExtractor::isSuitable	(	const PSimHit &	simHit,
		const GeomDetUnit &	gdu
	)		const

private

?????

Definition at line 174 of file PixelClusterShapeExtractor.cc.

References funct::abs(), constexpr, PSimHit::entryPoint(), PSimHit::exitPoint(), f, trackerHitRTTI::isFast(), PSimHit::momentumAtEntry(), PV3DBase< T, PVType, FrameType >::perp2(), GeomDet::position(), PSimHit::processType(), GeomDet::toLocal(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by checkSimHits().

                                                                                                {
   // Outgoing?
   // very expensive....
   GlobalVector gvec = gdu.position() - GlobalPoint(0, 0, 0);
   LocalVector lvec = gdu.toLocal(gvec);
   LocalVector ldir = simHit.exitPoint() - simHit.entryPoint();
 
   // cut on size as well (pixel is 285um thick...
   bool isOutgoing = std::abs(ldir.z()) > 0.01f && (lvec.z() * ldir.z() > 0);
 
   const bool isRelevant = RelevantProcesses.count(simHit.processType());
   // From a relevant process? primary or decay
   //bool isRelevant = (simHit.processType() == 2 ||
   //                   simHit.processType() == 4);
 
   constexpr float ptCut2 = 0.2 * 0.2;  //  0.050*0.050;
   // Fast enough? pt > 50 MeV/c   FIXME (at least 200MeV....
   bool isFast = (simHit.momentumAtEntry().perp2() > ptCut2);
 
   //std::cout << "isOutgoing = " << isOutgoing << ", isRelevant = " << simHit.processType() << ", isFast = " << isFast << std::endl;
   return (isOutgoing && isRelevant && isFast);
 }

void PixelClusterShapeExtractor::processPixelRecHits	(	SiPixelRecHitCollection::DataContainer const &	recHits,
		TrackerHitAssociator const &	theAssociator,
		ClusterShapeHitFilter const &	theFilter,
		SiPixelClusterShapeCache const &	clusterShapeCache,
		const TrackerTopology &	tkTpl
	)		const

private

Definition at line 257 of file PixelClusterShapeExtractor.cc.

References checkSimHits(), HTMLExport::elem(), hspc, crabWrapper::key, SiStripPI::max, noBPIX1, processSim(), TrackerTopology::pxbLayer(), rpcPointValidation_cfi::recHit, rpcPointValidation_cfi::simHit, findQualityFiles::size, and contentValuesCheck::ss.

Referenced by analyzeSimHits().

                                                                                          {
   struct Elem {
     const SiPixelRecHit* rhit;
     PSimHit shit;
     unsigned int size;
   };
   std::map<pair<unsigned int, float>, Elem> simHitMap;
 
   PSimHit simHit;
   pair<unsigned int, float> key;
   unsigned int ss;
 
   for (auto const& recHit : recHits) {
     if (noBPIX1 && tkTpl.pxbLayer(recHit.geographicalId()) == 1)
       continue;
     if (!checkSimHits(recHit, theHitAssociator, simHit, key, ss))
       continue;
     // Fill map
     if (simHitMap.count(key) == 0) {
       simHitMap[key] = {&recHit, simHit, ss};
     } else if (recHit.cluster()->size() > simHitMap[key].rhit->cluster()->size())
       simHitMap[key] = {&recHit, simHit, std::max(ss, simHitMap[key].size)};
   }
   for (auto const& elem : simHitMap) {
     /* irrelevant
    auto const rh = *elem.second.rhit;
    auto const& topol = reinterpret_cast<const PixelGeomDetUnit*>(rh.detUnit())->specificTopology();
    auto const & cl = *rh.cluster();
    if (cl.minPixelCol()==0) continue;
    if (cl.maxPixelCol()+1==topol.ncolumns()) continue;
    if (cl.minPixelRow()==0) continue; 
    if (cl.maxPixelRow()+1==topol.nrows()) continue;
    */
     if (elem.second.size == 1)
       processSim(*elem.second.rhit, theFilter, elem.second.shit, clusterShapeCache, hspc);
   }
 }

void PixelClusterShapeExtractor::processRec	(	const SiPixelRecHit &	recHit,
		ClusterShapeHitFilter const &	theFilter,
		LocalVector	ldir,
		const SiPixelClusterShapeCache &	clusterShapeCache,
		const vector< TH2F * > &	histo
	)		const

private

Definition at line 199 of file PixelClusterShapeExtractor.cc.

References funct::abs(), gather_cfg::cout, exMax, eyMax, edm::VecArray< T, N >::front(), TrackingRecHit::geographicalId(), ClusterShapeHitFilter::getSizes(), mps_fire::i, CommonMethods::lock(), PV3DBase< T, PVType, FrameType >::mag(), edm::VecArray< T, N >::size(), and theMutex.

Referenced by analyzeRecTracks(), and processSim().

                                                                               {
   int part;
   ClusterData::ArrayType meas;
   pair<float, float> pred;
 
   if (theClusterShape.getSizes(recHit, ldir, clusterShapeCache, part, meas, pred))
     if (meas.size() == 1)
       if (meas.front().first <= exMax && meas.front().second <= eyMax) {
         int i = (part * (exMax + 1) + meas.front().first) * (eyMax + 1) + meas.front().second;
 #ifdef DO_DEBUG
         if (meas.front().second == 0 && std::abs(pred.second) > 3) {
           Lock lock(theMutex[0]);
           int id = recHit.geographicalId();
           std::cout << id << " bigpred " << meas.front().first << '/' << meas.front().second << ' ' << pred.first << '/'
                     << pred.second << ' ' << ldir << ' ' << ldir.mag() << std::endl;
         }
 #endif
         Lock lock(theMutex[i]);
         histo[i]->Fill(pred.first, pred.second);
       }
 }

void PixelClusterShapeExtractor::processSim	(	const SiPixelRecHit &	recHit,
		ClusterShapeHitFilter const &	theClusterFilter,
		const PSimHit &	simHit,
		const SiPixelClusterShapeCache &	clusterShapeCache,
		const vector< TH2F * > &	histo
	)		const

private

Definition at line 226 of file PixelClusterShapeExtractor.cc.

References PSimHit::entryPoint(), PSimHit::exitPoint(), and processRec().

Referenced by processPixelRecHits().

                                                                               {
   LocalVector ldir = simHit.exitPoint() - simHit.entryPoint();
   processRec(recHit, theClusterFilter, ldir, clusterShapeCache, histo);
 }