#include <SiPixelPhase1TrackClusters.h>

Inheritance diagram for SiPixelPhase1TrackClusters:

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

	SiPixelPhase1TrackClusters (const edm::ParameterSet &conf)

Public Member Functions inherited from SiPixelPhase1Base
void	bookHistograms (DQMStore::IBooker &iBooker, edm::Run const &run, edm::EventSetup const &iSetup)

	SiPixelPhase1Base (const edm::ParameterSet &iConfig)

virtual	~SiPixelPhase1Base ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 &	itemsToGetFromEvent () 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

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

virtual	~EDConsumerBase () noexcept(false)

Public Member Functions inherited from HistogramManagerHolder
	HistogramManagerHolder (const edm::ParameterSet &iConfig)

Private Types
enum	{ ONTRACK_CHARGE, ONTRACK_SIZE, ONTRACK_NCLUSTERS, ONTRACK_POSITION_B, ONTRACK_POSITION_F, OFFTRACK_CHARGE, OFFTRACK_SIZE, OFFTRACK_NCLUSTERS, OFFTRACK_POSITION_B, OFFTRACK_POSITION_F, NTRACKS, NTRACKS_VOLUME }

Private Attributes
edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > >	clustersToken_

edm::EDGetTokenT< reco::TrackCollection >	tracksToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::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 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)

Protected Attributes inherited from HistogramManagerHolder
GeometryInterface	geometryInterface

std::vector< HistogramManager >	histo

Detailed Description

Definition at line 15 of file SiPixelPhase1TrackClusters.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
ONTRACK_CHARGE
ONTRACK_SIZE
ONTRACK_NCLUSTERS
ONTRACK_POSITION_B
ONTRACK_POSITION_F
OFFTRACK_CHARGE
OFFTRACK_SIZE
OFFTRACK_NCLUSTERS
OFFTRACK_POSITION_B
OFFTRACK_POSITION_F
NTRACKS
NTRACKS_VOLUME

Definition at line 16 of file SiPixelPhase1TrackClusters.h.

        {
     ONTRACK_CHARGE,
     ONTRACK_SIZE,
     ONTRACK_NCLUSTERS,
     ONTRACK_POSITION_B,
     ONTRACK_POSITION_F,
 
     OFFTRACK_CHARGE,
     OFFTRACK_SIZE,
     OFFTRACK_NCLUSTERS,
     OFFTRACK_POSITION_B,
     OFFTRACK_POSITION_F,
 
     NTRACKS,
     NTRACKS_VOLUME
   };

Constructor & Destructor Documentation

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

explicit

Definition at line 25 of file SiPixelPhase1TrackClusters.cc.

References clustersToken_, edm::ParameterSet::getParameter(), and tracksToken_.

                                                                                      :
   SiPixelPhase1Base(iConfig) 
 {
   clustersToken_ = consumes<edmNew::DetSetVector<SiPixelCluster>>(iConfig.getParameter<edm::InputTag>("clusters"));
   tracksToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("tracks"));
 }

Member Function Documentation

void SiPixelPhase1TrackClusters::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

Definition at line 32 of file SiPixelPhase1TrackClusters.cc.

References funct::abs(), edmNew::DetSetVector< T >::begin(), SiPixelCluster::charge(), SiPixelRecHit::cluster(), clustersToken_, allConversions_cfi::d0, DEFINE_FWK_MODULE, allConversions_cfi::dz, edmNew::DetSetVector< T >::end(), objects.autophobj::float, edm::EventSetup::get(), edm::Event::getByToken(), h, hcalSimParameters_cfi::hb, HistogramManagerHolder::histo, hcalTTPDigis_cfi::id, TrackerGeometry::idToDet(), iEvent, edm::ESHandleBase::isValid(), crabWrapper::key, Topology::localPosition(), PV3DBase< T, PVType, FrameType >::mag(), edmNew::makeRefTo(), NTRACKS, NTRACKS_VOLUME, OFFTRACK_CHARGE, OFFTRACK_NCLUSTERS, OFFTRACK_POSITION_B, OFFTRACK_POSITION_F, OFFTRACK_SIZE, ONTRACK_CHARGE, ONTRACK_NCLUSTERS, ONTRACK_POSITION_B, ONTRACK_POSITION_F, ONTRACK_SIZE, PV3DBase< T, PVType, FrameType >::phi(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, funct::pow(), SiPixelCluster::size(), mathSSE::sqrt(), funct::tan(), HiIsolationCommonParameters_cff::track, mixOne_simraw_on_sim_cfi::tracker, l1t::tracks, tracksToken_, PV3DBase< T, PVType, FrameType >::x(), SiPixelCluster::x(), PV3DBase< T, PVType, FrameType >::y(), SiPixelCluster::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                             {
 
   // get geometry
   edm::ESHandle<TrackerGeometry> tracker;
   iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
   assert(tracker.isValid());
   
   //get the map
   edm::Handle<reco::TrackCollection> tracks;
   iEvent.getByToken( tracksToken_, tracks);
   
   // get clusters
   edm::Handle< edmNew::DetSetVector<SiPixelCluster> >  clusterColl;
   iEvent.getByToken( clustersToken_, clusterColl );
   
   // we need to store some per-cluster data. Instead of a map, we use a vector,
   // exploiting the fact that all custers live in the DetSetVector and we can 
   // use the same indices to refer to them.
   // corr_charge is not strictly needed but cleaner to have it.
   std::vector<bool>  ontrack    (clusterColl->data().size(), false);
   std::vector<float> corr_charge(clusterColl->data().size(), -1.0f);
 
   for (auto const & track : *tracks) {
 
     bool isBpixtrack = false, isFpixtrack = false, crossesPixVol=false;
 
     // find out whether track crosses pixel fiducial volume (for cosmic tracks)
     double d0 = track.d0(), dz = track.dz(); 
     if(std::abs(d0)<15 && std::abs(dz)<50) crossesPixVol = true;
 
     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 hit = *(hb+h);
       if (!hit->isValid()) continue;
       DetId id = hit->geographicalId();
 
       // check that we are in the pixel
       uint32_t subdetid = (id.subdetId());
       if (subdetid == PixelSubdetector::PixelBarrel) isBpixtrack = true;
       if (subdetid == PixelSubdetector::PixelEndcap) isFpixtrack = true;
       if (subdetid != PixelSubdetector::PixelBarrel && subdetid != PixelSubdetector::PixelEndcap) continue;
       auto pixhit = dynamic_cast<const SiPixelRecHit*>(hit->hit());
       if (!pixhit) continue;
         
       // get the cluster
       auto clust = pixhit->cluster();
       if (clust.isNull()) continue; 
       ontrack[clust.key()] = true; // mark cluster as ontrack
 
 
       // correct charge for track impact angle
       auto const & ltp = trajParams[h];
       LocalVector localDir = ltp.momentum()/ltp.momentum().mag();
      
       float clust_alpha = atan2(localDir.z(), localDir.x());
       float clust_beta  = atan2(localDir.z(), localDir.y());
       double corrCharge = clust->charge() * sqrt( 1.0 / ( 1.0/pow( tan(clust_alpha), 2 ) + 
                                                           1.0/pow( tan(clust_beta ), 2 ) + 
                                                           1.0 ));
       corr_charge[clust.key()] = (float) corrCharge;
     }
 
     // statistics on tracks
     histo[NTRACKS].fill(1, DetId(0), &iEvent);
     if (isBpixtrack || isFpixtrack) 
       histo[NTRACKS].fill(2, DetId(0), &iEvent);
     if (isBpixtrack) 
       histo[NTRACKS].fill(3, DetId(0), &iEvent);
     if (isFpixtrack) 
       histo[NTRACKS].fill(4, DetId(0), &iEvent);
 
     if (crossesPixVol) {
       if (isBpixtrack || isFpixtrack)
         histo[NTRACKS_VOLUME].fill(1, DetId(0), &iEvent);
       else 
         histo[NTRACKS_VOLUME].fill(0, DetId(0), &iEvent);
     }
   }
 
   edmNew::DetSetVector<SiPixelCluster>::const_iterator it;
   for (it = clusterColl->begin(); it != clusterColl->end(); ++it) {
     auto id = DetId(it->detId());
 
     const PixelGeomDetUnit* geomdetunit = dynamic_cast<const PixelGeomDetUnit*> ( tracker->idToDet(id) );
     const PixelTopology& topol = geomdetunit->specificTopology();
 
     for(auto subit = it->begin(); subit != it->end(); ++subit) {
       // we could do subit-...->data().front() as well, but this seems cleaner.
       auto key = edmNew::makeRefTo(clusterColl, subit).key(); 
       bool is_ontrack = ontrack[key];
       float corrected_charge = corr_charge[key];
       SiPixelCluster const& cluster = *subit;
 
       LocalPoint clustlp = topol.localPosition(MeasurementPoint(cluster.x(), cluster.y()));
       GlobalPoint clustgp = geomdetunit->surface().toGlobal(clustlp);
 
       if (is_ontrack) {
         histo[ONTRACK_NCLUSTERS ].fill(id, &iEvent);
         histo[ONTRACK_CHARGE    ].fill(double(corrected_charge), id, &iEvent);
         histo[ONTRACK_SIZE      ].fill(double(cluster.size()  ), id, &iEvent);
         histo[ONTRACK_POSITION_B].fill(clustgp.z(),   clustgp.phi(),   id, &iEvent);
         histo[ONTRACK_POSITION_F].fill(clustgp.x(),   clustgp.y(),     id, &iEvent);
       } else {
         histo[OFFTRACK_NCLUSTERS ].fill(id, &iEvent);
         histo[OFFTRACK_CHARGE    ].fill(double(cluster.charge()), id, &iEvent);
         histo[OFFTRACK_SIZE      ].fill(double(cluster.size()  ), id, &iEvent);
         histo[OFFTRACK_POSITION_B].fill(clustgp.z(),   clustgp.phi(),   id, &iEvent);
         histo[OFFTRACK_POSITION_F].fill(clustgp.x(),   clustgp.y(),     id, &iEvent);
       }
     }
   }
 
   histo[ONTRACK_NCLUSTERS].executePerEventHarvesting(&iEvent);
   histo[OFFTRACK_NCLUSTERS].executePerEventHarvesting(&iEvent);
 }

Member Data Documentation

edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > SiPixelPhase1TrackClusters::clustersToken_

private

Definition at line 38 of file SiPixelPhase1TrackClusters.h.

Referenced by analyze(), and SiPixelPhase1TrackClusters().

edm::EDGetTokenT<reco::TrackCollection> SiPixelPhase1TrackClusters::tracksToken_

private

Definition at line 39 of file SiPixelPhase1TrackClusters.h.

Referenced by analyze(), and SiPixelPhase1TrackClusters().

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation