#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 &)

	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
bool	crossingframe

const CSCGeometry *	cscgeom

edm::InputTag	CSCsimHitsTag

edm::InputTag	CSCsimHitsXFTag

edm::InputTag	linksTag

std::map< SimHitIdpr, float >	theChargeMap

unsigned int	theDetId

const DigiSimLinks *	theDigiSimLinks

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

float	theTotalCharge

const DigiSimLinks *	theWireDigiSimLinks

edm::InputTag	wireLinksTag

Detailed Description

Definition at line 24 of file MuonTruth.h.

Member Typedef Documentation

◆ DigiSimLinks

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

Definition at line 26 of file MuonTruth.h.

◆ LayerLinks

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

Definition at line 28 of file MuonTruth.h.

◆ SimHitIdpr

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

Definition at line 29 of file MuonTruth.h.

◆ WireDigiSimLinks

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

Definition at line 27 of file MuonTruth.h.

Constructor & Destructor Documentation

◆ MuonTruth() [1/2]

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

Definition at line 8 of file MuonTruth.cc.

     : 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"))
  
 {
   initEvent(event, setup);
 }

References initEvent(), and singleTopDQM_cfi::setup.

◆ MuonTruth() [2/2]

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

Definition at line 22 of file MuonTruth.cc.

     : 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"))
  
 {
   iC.consumes<DigiSimLinks>(linksTag);
   iC.consumes<DigiSimLinks>(wireLinksTag);
   if (crossingframe) {
     iC.consumes<CrossingFrame<PSimHit>>(CSCsimHitsXFTag);
   } else if (!CSCsimHitsTag.label().empty()) {
     iC.consumes<edm::PSimHitContainer>(CSCsimHitsTag);
   }
 }

References crossingframe, CSCsimHitsTag, CSCsimHitsXFTag, edm::InputTag::label(), linksTag, and wireLinksTag.

Member Function Documentation

◆ addChannel()

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

private

Definition at line 206 of file MuonTruth.cc.

                                                                                   {
   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;
       }
     }
   }
 }

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

Referenced by analyze().

◆ analyze() [1/3]

void MuonTruth::analyze ( const CSCRecHit2D & recHit )

Definition at line 156 of file MuonTruth.cc.

                                                  {
   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);
     }
   }
 }

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

Referenced by CSCTruthTest::analyze().

◆ analyze() [2/3]

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

Definition at line 179 of file MuonTruth.cc.

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

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

◆ analyze() [3/3]

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

Definition at line 190 of file MuonTruth.cc.

                                                                                         {
   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.);
   }
 }

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

◆ hitsFromSimTrack()

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

private

Definition at line 130 of file MuonTruth.cc.

                                                                           {
   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;
 }

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

Referenced by particleType(), and simHits().

◆ initEvent()

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

Definition at line 42 of file MuonTruth.cc.

                                                                            {
   edm::Handle<DigiSimLinks> digiSimLinks;
   LogTrace("MuonTruth") << "getting CSC Strip DigiSimLink collection - " << linksTag;
   event.getByLabel(linksTag, digiSimLinks);
   theDigiSimLinks = digiSimLinks.product();
  
   edm::Handle<DigiSimLinks> wireDigiSimLinks;
   LogTrace("MuonTruth") << "getting CSC Wire DigiSimLink collection - " << wireLinksTag;
   event.getByLabel(wireLinksTag, wireDigiSimLinks);
   theWireDigiSimLinks = wireDigiSimLinks.product();
  
   // get CSC Geometry to use CSCLayer methods
   edm::ESHandle<CSCGeometry> mugeom;
   setup.get<MuonGeometryRecord>().get(mugeom);
   cscgeom = &*mugeom;
  
   // get CSC Bad Chambers (ME4/2)
   edm::ESHandle<CSCBadChambers> badChambers;
   setup.get<CSCBadChambersRcd>().get(badChambers);
   cscBadChambers = badChambers.product();
  
   theSimHitMap.clear();
  
   if (crossingframe) {
     edm::Handle<CrossingFrame<PSimHit>> cf;
     LogTrace("MuonTruth") << "getting CrossingFrame<PSimHit> collection - " << CSCsimHitsXFTag;
     event.getByLabel(CSCsimHitsXFTag, cf);
  
     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()) {
     edm::Handle<edm::PSimHitContainer> CSCsimhits;
     LogTrace("MuonTruth") << "getting PSimHit collection - " << CSCsimHitsTag;
     event.getByLabel(CSCsimHitsTag, CSCsimhits);
     LogTrace("MuonTruth") << "... size = " << CSCsimhits->size();
  
     for (edm::PSimHitContainer::const_iterator hitItr = CSCsimhits->begin(); hitItr != CSCsimhits->end(); ++hitItr) {
       theSimHitMap[hitItr->detUnitId()].push_back(*hitItr);
     }
   }
 }

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

Referenced by MuonTruth().

◆ muonFraction()

float MuonTruth::muonFraction ( )

analyze() must be called before any of the following

Definition at line 89 of file MuonTruth.cc.

                               {
   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;
 }

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

Referenced by CSCTruthTest::analyze().

◆ muonHits()

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

Definition at line 117 of file MuonTruth.cc.

                                        {
   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;
 }

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

◆ particleType()

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

private

Definition at line 147 of file MuonTruth.cc.

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

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

Referenced by muonFraction().

◆ simHits()

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

Definition at line 105 of file MuonTruth.cc.

                                       {
   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;
 }

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

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

Member Data Documentation

◆ crossingframe

bool MuonTruth::crossingframe

private

Definition at line 67 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ cscBadChambers

const CSCBadChambers* MuonTruth::cscBadChambers

Definition at line 47 of file MuonTruth.h.

Referenced by initEvent().

◆ cscgeom

const CSCGeometry* MuonTruth::cscgeom

private

Definition at line 73 of file MuonTruth.h.

Referenced by analyze(), and initEvent().

◆ CSCsimHitsTag

edm::InputTag MuonTruth::CSCsimHitsTag

private

Definition at line 68 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ CSCsimHitsXFTag

edm::InputTag MuonTruth::CSCsimHitsXFTag

private

Definition at line 69 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ linksTag

edm::InputTag MuonTruth::linksTag

private

Definition at line 64 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

◆ theChargeMap

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

private

Definition at line 56 of file MuonTruth.h.

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

◆ theDetId

unsigned int MuonTruth::theDetId

private

Definition at line 59 of file MuonTruth.h.

Referenced by analyze(), and hitsFromSimTrack().

◆ theDigiSimLinks

const DigiSimLinks* MuonTruth::theDigiSimLinks

private

Definition at line 61 of file MuonTruth.h.

Referenced by analyze(), and initEvent().

◆ theSimHitMap

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

private

Definition at line 71 of file MuonTruth.h.

Referenced by hitsFromSimTrack(), and initEvent().

◆ theTotalCharge

float MuonTruth::theTotalCharge

private

Definition at line 57 of file MuonTruth.h.

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

◆ theWireDigiSimLinks

const DigiSimLinks* MuonTruth::theWireDigiSimLinks

private

Definition at line 62 of file MuonTruth.h.

Referenced by analyze(), and initEvent().

◆ wireLinksTag

edm::InputTag MuonTruth::wireLinksTag

private

Definition at line 65 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

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

◆ crossingframe

◆ cscBadChambers

◆ cscgeom

◆ CSCsimHitsTag

◆ CSCsimHitsXFTag

◆ linksTag

◆ theChargeMap

◆ theDetId

◆ theDigiSimLinks

◆ theSimHitMap

◆ theTotalCharge

◆ theWireDigiSimLinks

◆ wireLinksTag