#include <MuonTruth.h>

Public Types
typedef edm::DetSetVector< StripDigiSimLink >	DigiSimLinks

typedef edm::DetSet< StripDigiSimLink >	LayerLinks

typedef std::pair< uint32_t, EncodedEventId >	SimHitIdpr

typedef edm::DetSetVector< StripDigiSimLink >	WireDigiSimLinks

Public Member Functions
void	analyze (const CSCRecHit2D &recHit)

void	analyze (const CSCStripDigi &stripDigi, int rawDetIdCorrespondingToCSCLayer)

void	analyze (const CSCWireDigi &wireDigi, int rawDetIdCorrespondingToCSCLayer)

void	initEvent (const edm::Event &, const edm::EventSetup &)

float	muonFraction ()
	analyze() must be called before any of the following More...

std::vector< PSimHit >	muonHits ()

	MuonTruth (const edm::Event &, const edm::EventSetup &, const edm::ParameterSet &, edm::ConsumesCollector &)

	MuonTruth (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

std::vector< PSimHit >	simHits ()

Public Attributes
const CSCBadChambers *	cscBadChambers

Private Member Functions
void	addChannel (const LayerLinks &layerLinks, int channel, float weight=1.)

std::vector< PSimHit >	hitsFromSimTrack (SimHitIdpr truthId)

int	particleType (SimHitIdpr truthId)

Private Attributes
const edm::ESGetToken< CSCBadChambers, CSCBadChambersRcd >	badToken_

const bool	crossingframe

const CSCGeometry *	cscgeom

const edm::InputTag	CSCsimHitsTag

const edm::InputTag	CSCsimHitsXFTag

const edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	geomToken_

const edm::InputTag	linksTag

const edm::EDGetTokenT< DigiSimLinks >	linksToken_

edm::EDGetTokenT< edm::PSimHitContainer >	simHitsToken_

edm::EDGetTokenT< CrossingFrame< PSimHit > >	simHitsXFToken_

std::map< SimHitIdpr, float >	theChargeMap

unsigned int	theDetId

const DigiSimLinks *	theDigiSimLinks

std::map< unsigned int, edm::PSimHitContainer >	theSimHitMap

float	theTotalCharge

const DigiSimLinks *	theWireDigiSimLinks

const edm::InputTag	wireLinksTag

const edm::EDGetTokenT< DigiSimLinks >	wireLinksToken_

Detailed Description

Definition at line 25 of file MuonTruth.h.

Member Typedef Documentation

◆ DigiSimLinks

typedef edm::DetSetVector<StripDigiSimLink> MuonTruth::DigiSimLinks

Definition at line 27 of file MuonTruth.h.

◆ LayerLinks

typedef edm::DetSet<StripDigiSimLink> MuonTruth::LayerLinks

Definition at line 29 of file MuonTruth.h.

◆ SimHitIdpr

typedef std::pair<uint32_t, EncodedEventId> MuonTruth::SimHitIdpr

Definition at line 30 of file MuonTruth.h.

◆ WireDigiSimLinks

typedef edm::DetSetVector<StripDigiSimLink> MuonTruth::WireDigiSimLinks

Definition at line 28 of file MuonTruth.h.

Constructor & Destructor Documentation

◆ MuonTruth() [1/2]

MuonTruth::MuonTruth	(	const edm::Event &	event,
		const edm::EventSetup &	setup,
		const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	iC
	)

Definition at line 7 of file MuonTruth.cc.

References initEvent(), and singleTopDQM_cfi::setup.

     : theDigiSimLinks(nullptr),
       theWireDigiSimLinks(nullptr),
       linksTag(conf.getParameter<edm::InputTag>("CSClinksTag")),
       wireLinksTag(conf.getParameter<edm::InputTag>("CSCwireLinksTag")),
       // CrossingFrame used or not ?
       crossingframe(conf.getParameter<bool>("crossingframe")),
       CSCsimHitsTag(conf.getParameter<edm::InputTag>("CSCsimHitsTag")),
       CSCsimHitsXFTag(conf.getParameter<edm::InputTag>("CSCsimHitsXFTag")),
       geomToken_(iC.esConsumes<CSCGeometry, MuonGeometryRecord>()),
       badToken_(iC.esConsumes<CSCBadChambers, CSCBadChambersRcd>()),
       linksToken_(iC.consumes<DigiSimLinks>(linksTag)),
       wireLinksToken_(iC.consumes<DigiSimLinks>(wireLinksTag)) {
   initEvent(event, setup);
 }

◆ MuonTruth() [2/2]

MuonTruth::MuonTruth	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	iC
	)

Definition at line 26 of file MuonTruth.cc.

References crossingframe, CSCsimHitsTag, CSCsimHitsXFTag, edm::InputTag::label(), simHitsToken_, and simHitsXFToken_.

     : theDigiSimLinks(nullptr),
       theWireDigiSimLinks(nullptr),
       linksTag(conf.getParameter<edm::InputTag>("CSClinksTag")),
       wireLinksTag(conf.getParameter<edm::InputTag>("CSCwireLinksTag")),
       // CrossingFrame used or not ?
       crossingframe(conf.getParameter<bool>("crossingframe")),
       CSCsimHitsTag(conf.getParameter<edm::InputTag>("CSCsimHitsTag")),
       CSCsimHitsXFTag(conf.getParameter<edm::InputTag>("CSCsimHitsXFTag")),
       linksToken_(iC.consumes<DigiSimLinks>(linksTag)),
       wireLinksToken_(iC.consumes<DigiSimLinks>(wireLinksTag)) {
   if (crossingframe) {
     simHitsXFToken_ = iC.consumes<CrossingFrame<PSimHit>>(CSCsimHitsXFTag);
   } else if (!CSCsimHitsTag.label().empty()) {
     simHitsToken_ = iC.consumes<edm::PSimHitContainer>(CSCsimHitsTag);
   }
 }

Member Function Documentation

◆ addChannel()

void MuonTruth::addChannel	(	const LayerLinks &	layerLinks,
		int	channel,
		float	weight = `1.`
	)

private

Definition at line 202 of file MuonTruth.cc.

References edm::DetSet< T >::begin(), ALCARECOTkAlJpsiMuMu_cff::charge, edm::DetSet< T >::end(), theChargeMap, theTotalCharge, and mps_merge::weight.

Referenced by analyze().

                                                                                   {
   LayerLinks::const_iterator linkItr = layerLinks.begin(), lastLayerLink = layerLinks.end();
 
   for (; linkItr != lastLayerLink; ++linkItr) {
     int linkChannel = linkItr->channel();
     if (linkChannel == channel) {
       float charge = linkItr->fraction() * weight;
       theTotalCharge += charge;
       // see if it's in the map
       SimHitIdpr truthId(linkItr->SimTrackId(), linkItr->eventId());
       std::map<SimHitIdpr, float>::const_iterator chargeMapItr = theChargeMap.find(truthId);
       if (chargeMapItr == theChargeMap.end()) {
         theChargeMap[truthId] = charge;
       } else {
         theChargeMap[truthId] += charge;
       }
     }
   }
 }

◆ analyze() [1/3]

void MuonTruth::analyze ( const CSCRecHit2D & recHit )

Definition at line 152 of file MuonTruth.cc.

References addChannel(), CSCLayerGeometry::channel(), cscgeom, edm::DetSetVector< T >::end(), edm::DetSetVector< T >::find(), relativeConstraints::geometry, CSCGeometry::layer(), rpcPointValidation_cfi::recHit, theChargeMap, theDetId, theDigiSimLinks, and theTotalCharge.

Referenced by CSCTruthTest::analyze().

                                                  {
   theChargeMap.clear();
   theTotalCharge = 0.;
   theDetId = recHit.geographicalId().rawId();
 
   int nchannels = recHit.nStrips();
   const CSCLayerGeometry *laygeom = cscgeom->layer(recHit.cscDetId())->geometry();
 
   for (int idigi = 0; idigi < nchannels; ++idigi) {
     // strip and readout channel numbers may differ in ME1/1A
     int istrip = recHit.channels(idigi);
     int channel = laygeom->channel(istrip);
     float weight = recHit.adcs(idigi, 0);  // DL: I think this is wrong before and after...seems to
                                            // assume one time binadcContainer[idigi];
 
     DigiSimLinks::const_iterator layerLinks = theDigiSimLinks->find(theDetId);
 
     if (layerLinks != theDigiSimLinks->end()) {
       addChannel(*layerLinks, channel, weight);
     }
   }
 }

◆ analyze() [2/3]

void MuonTruth::analyze	(	const CSCStripDigi &	stripDigi,
		int	rawDetIdCorrespondingToCSCLayer
	)

Definition at line 175 of file MuonTruth.cc.

References addChannel(), edm::DetSetVector< T >::end(), edm::DetSetVector< T >::find(), CSCStripDigi::getStrip(), theChargeMap, theDetId, theDigiSimLinks, and theTotalCharge.

                                                                                           {
   theDetId = rawDetIdCorrespondingToCSCLayer;
   theChargeMap.clear();
   theTotalCharge = 0.;
 
   DigiSimLinks::const_iterator layerLinks = theDigiSimLinks->find(theDetId);
   if (layerLinks != theDigiSimLinks->end()) {
     addChannel(*layerLinks, stripDigi.getStrip(), 1.);
   }
 }

◆ analyze() [3/3]

void MuonTruth::analyze	(	const CSCWireDigi &	wireDigi,
		int	rawDetIdCorrespondingToCSCLayer
	)

Definition at line 186 of file MuonTruth.cc.

References addChannel(), edm::DetSetVector< T >::end(), edm::DetSetVector< T >::find(), CSCWireDigi::getWireGroup(), theChargeMap, theDetId, theDigiSimLinks, theTotalCharge, and theWireDigiSimLinks.

                                                                                         {
   theDetId = rawDetIdCorrespondingToCSCLayer;
   theChargeMap.clear();
   theTotalCharge = 0.;
 
   WireDigiSimLinks::const_iterator layerLinks = theWireDigiSimLinks->find(theDetId);
 
   if (layerLinks != theDigiSimLinks->end()) {
     // In the simulation digis, the channel labels for wires and strips must be
     // distinct, therefore:
     int wireDigiInSimulation = wireDigi.getWireGroup() + 100;
     //
     addChannel(*layerLinks, wireDigiInSimulation, 1.);
   }
 }

◆ hitsFromSimTrack()

std::vector< PSimHit > MuonTruth::hitsFromSimTrack ( MuonTruth::SimHitIdpr truthId )

private

Definition at line 126 of file MuonTruth.cc.

References newFWLiteAna::found, mps_fire::result, theDetId, and theSimHitMap.

Referenced by particleType(), and simHits().

                                                                           {
   std::vector<PSimHit> result;
 
   auto found = theSimHitMap.find(theDetId);
   if (found != theSimHitMap.end()) {
     for (auto const &hit : found->second) {
       unsigned int hitTrack = hit.trackId();
       EncodedEventId hitEvId = hit.eventId();
 
       if (hitTrack == truthId.first && hitEvId == truthId.second) {
         result.push_back(hit);
       }
     }
   }
   return result;
 }

◆ initEvent()

void MuonTruth::initEvent	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

Definition at line 44 of file MuonTruth.cc.

References badToken_, crossingframe, cscBadChambers, cscgeom, CSCsimHitsTag, CSCsimHitsXFTag, TTClusterAssociation_cfi::digiSimLinks, geomToken_, edm::InputTag::label(), linksTag, linksToken_, LogTrace, edm::Handle< T >::product(), edm::ESHandle< T >::product(), singleTopDQM_cfi::setup, simHitsToken_, simHitsXFToken_, theDigiSimLinks, theSimHitMap, theWireDigiSimLinks, wireLinksTag, and wireLinksToken_.

Referenced by MuonTruth().

                                                                            {
   LogTrace("MuonTruth") << "getting CSC Strip DigiSimLink collection - " << linksTag;
   const edm::Handle<DigiSimLinks> &digiSimLinks = event.getHandle(linksToken_);
   theDigiSimLinks = digiSimLinks.product();
 
   LogTrace("MuonTruth") << "getting CSC Wire DigiSimLink collection - " << wireLinksTag;
   const edm::Handle<DigiSimLinks> &wireDigiSimLinks = event.getHandle(wireLinksToken_);
   theWireDigiSimLinks = wireDigiSimLinks.product();
 
   // get CSC Geometry to use CSCLayer methods
   const edm::ESHandle<CSCGeometry> &mugeom = setup.getHandle(geomToken_);
   cscgeom = mugeom.product();
 
   // get CSC Bad Chambers (ME4/2)
   const edm::ESHandle<CSCBadChambers> &badChambers = setup.getHandle(badToken_);
   cscBadChambers = badChambers.product();
 
   theSimHitMap.clear();
 
   if (crossingframe) {
     LogTrace("MuonTruth") << "getting CrossingFrame<PSimHit> collection - " << CSCsimHitsXFTag;
     const edm::Handle<CrossingFrame<PSimHit>> &cf = event.getHandle(simHitsXFToken_);
 
     std::unique_ptr<MixCollection<PSimHit>> CSCsimhits(new MixCollection<PSimHit>(cf.product()));
     LogTrace("MuonTruth") << "... size = " << CSCsimhits->size();
 
     for (MixCollection<PSimHit>::MixItr hitItr = CSCsimhits->begin(); hitItr != CSCsimhits->end(); ++hitItr) {
       theSimHitMap[hitItr->detUnitId()].push_back(*hitItr);
     }
 
   } else if (!CSCsimHitsTag.label().empty()) {
     LogTrace("MuonTruth") << "getting PSimHit collection - " << CSCsimHitsTag;
     const edm::Handle<edm::PSimHitContainer> &CSCsimhits = event.getHandle(simHitsToken_);
     LogTrace("MuonTruth") << "... size = " << CSCsimhits->size();
 
     for (edm::PSimHitContainer::const_iterator hitItr = CSCsimhits->begin(); hitItr != CSCsimhits->end(); ++hitItr) {
       theSimHitMap[hitItr->detUnitId()].push_back(*hitItr);
     }
   }
 }

◆ muonFraction()

float MuonTruth::muonFraction ( )

analyze() must be called before any of the following

Definition at line 85 of file MuonTruth.cc.

References funct::abs(), particleType(), theChargeMap, and theTotalCharge.

Referenced by CSCTruthTest::analyze().

                               {
   if (theChargeMap.empty())
     return 0.;
 
   float muonCharge = 0.;
   for (std::map<SimHitIdpr, float>::const_iterator chargeMapItr = theChargeMap.begin();
        chargeMapItr != theChargeMap.end();
        ++chargeMapItr) {
     if (abs(particleType(chargeMapItr->first)) == 13) {
       muonCharge += chargeMapItr->second;
     }
   }
 
   return muonCharge / theTotalCharge;
 }

◆ muonHits()

std::vector< PSimHit > MuonTruth::muonHits ( )

Definition at line 113 of file MuonTruth.cc.

References funct::abs(), mps_fire::result, and simHits().

                                        {
   std::vector<PSimHit> result;
   std::vector<PSimHit> allHits = simHits();
   std::vector<PSimHit>::const_iterator hitItr = allHits.begin(), lastHit = allHits.end();
 
   for (; hitItr != lastHit; ++hitItr) {
     if (abs((*hitItr).particleType()) == 13) {
       result.push_back(*hitItr);
     }
   }
   return result;
 }

◆ particleType()

int MuonTruth::particleType ( MuonTruth::SimHitIdpr truthId )

private

Definition at line 143 of file MuonTruth.cc.

References hfClusterShapes_cfi::hits, hitsFromSimTrack(), and mps_fire::result.

Referenced by muonFraction().

                                                        {
   int result = 0;
   const std::vector<PSimHit> &hits = hitsFromSimTrack(truthId);
   if (!hits.empty()) {
     result = hits[0].particleType();
   }
   return result;
 }

◆ simHits()

std::vector< PSimHit > MuonTruth::simHits ( )

Definition at line 101 of file MuonTruth.cc.

References hitsFromSimTrack(), mps_fire::result, and theChargeMap.

Referenced by ntupleDataFormat.TrackingParticle::bestMatchingTrack(), and muonHits().

                                       {
   std::vector<PSimHit> result;
   for (std::map<SimHitIdpr, float>::const_iterator chargeMapItr = theChargeMap.begin();
        chargeMapItr != theChargeMap.end();
        ++chargeMapItr) {
     std::vector<PSimHit> trackHits = hitsFromSimTrack(chargeMapItr->first);
     result.insert(result.end(), trackHits.begin(), trackHits.end());
   }
 
   return result;
 }

Member Data Documentation

◆ badToken_

const edm::ESGetToken<CSCBadChambers, CSCBadChambersRcd> MuonTruth::badToken_

private

Definition at line 73 of file MuonTruth.h.

Referenced by initEvent().

◆ crossingframe

const bool MuonTruth::crossingframe

private

Definition at line 68 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ cscBadChambers

const CSCBadChambers* MuonTruth::cscBadChambers

Definition at line 48 of file MuonTruth.h.

Referenced by initEvent().

◆ cscgeom

const CSCGeometry* MuonTruth::cscgeom

private

Definition at line 82 of file MuonTruth.h.

Referenced by analyze(), and initEvent().

◆ CSCsimHitsTag

const edm::InputTag MuonTruth::CSCsimHitsTag

private

Definition at line 69 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ CSCsimHitsXFTag

const edm::InputTag MuonTruth::CSCsimHitsXFTag

private

Definition at line 70 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ geomToken_

const edm::ESGetToken<CSCGeometry, MuonGeometryRecord> MuonTruth::geomToken_

private

Definition at line 72 of file MuonTruth.h.

Referenced by initEvent().

◆ linksTag

const edm::InputTag MuonTruth::linksTag

private

Definition at line 65 of file MuonTruth.h.

Referenced by initEvent().

◆ linksToken_

const edm::EDGetTokenT<DigiSimLinks> MuonTruth::linksToken_

private

Definition at line 75 of file MuonTruth.h.

Referenced by initEvent().

◆ simHitsToken_

edm::EDGetTokenT<edm::PSimHitContainer> MuonTruth::simHitsToken_

private

Definition at line 78 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ simHitsXFToken_

edm::EDGetTokenT<CrossingFrame<PSimHit> > MuonTruth::simHitsXFToken_

private

Definition at line 77 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ theChargeMap

std::map<SimHitIdpr, float> MuonTruth::theChargeMap

private

Definition at line 57 of file MuonTruth.h.

Referenced by addChannel(), analyze(), muonFraction(), and simHits().

◆ theDetId

unsigned int MuonTruth::theDetId

private

Definition at line 60 of file MuonTruth.h.

Referenced by analyze(), and hitsFromSimTrack().

◆ theDigiSimLinks

const DigiSimLinks* MuonTruth::theDigiSimLinks

private

Definition at line 62 of file MuonTruth.h.

Referenced by analyze(), and initEvent().

◆ theSimHitMap

std::map<unsigned int, edm::PSimHitContainer> MuonTruth::theSimHitMap

private

Definition at line 80 of file MuonTruth.h.

Referenced by hitsFromSimTrack(), and initEvent().

◆ theTotalCharge

float MuonTruth::theTotalCharge

private

Definition at line 58 of file MuonTruth.h.

Referenced by addChannel(), analyze(), and muonFraction().

◆ theWireDigiSimLinks

const DigiSimLinks* MuonTruth::theWireDigiSimLinks

private

Definition at line 63 of file MuonTruth.h.

Referenced by analyze(), and initEvent().

◆ wireLinksTag

const edm::InputTag MuonTruth::wireLinksTag

private

Definition at line 66 of file MuonTruth.h.

Referenced by initEvent().

◆ wireLinksToken_

const edm::EDGetTokenT<DigiSimLinks> MuonTruth::wireLinksToken_

private

Definition at line 76 of file MuonTruth.h.

Referenced by initEvent().

Public Types

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ DigiSimLinks

◆ LayerLinks

◆ SimHitIdpr

◆ WireDigiSimLinks

Constructor & Destructor Documentation

◆ MuonTruth() [1/2]

◆ MuonTruth() [2/2]

Member Function Documentation

◆ addChannel()

◆ analyze() [1/3]

◆ analyze() [2/3]

◆ analyze() [3/3]

◆ hitsFromSimTrack()

◆ initEvent()

◆ muonFraction()

◆ muonHits()

◆ particleType()

◆ simHits()

Member Data Documentation

◆ badToken_

◆ crossingframe

◆ cscBadChambers

◆ cscgeom

◆ CSCsimHitsTag

◆ CSCsimHitsXFTag

◆ geomToken_

◆ linksTag

◆ linksToken_

◆ simHitsToken_

◆ simHitsXFToken_

◆ theChargeMap

◆ theDetId

◆ theDigiSimLinks

◆ theSimHitMap

◆ theTotalCharge

◆ theWireDigiSimLinks

◆ wireLinksTag

◆ wireLinksToken_