MuonTruth Class Reference

#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::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()

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

Definition at line 7 of file MuonTruth.cc.

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

    : cscBadChambers(nullptr),
      theTotalCharge(0.0),
      theDetId(0),
      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)),
      cscgeom(nullptr) {
  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 193 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 143 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 166 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 177 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 117 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 31 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_, theChargeMap, theDetId, theDigiSimLinks, theSimHitMap, theTotalCharge, theWireDigiSimLinks, wireLinksTag, and wireLinksToken_.

Referenced by CSCTruthTest::analyze().

                                                                           {
  theChargeMap.clear();
  theTotalCharge = 0.0;
  theDetId = 0;

  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 76 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 104 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 134 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 92 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 72 of file MuonTruth.h.

Referenced by initEvent().

crossingframe

const 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 81 of file MuonTruth.h.

Referenced by analyze(), and initEvent().

CSCsimHitsTag

const edm::InputTag MuonTruth::CSCsimHitsTag

private

Definition at line 68 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

CSCsimHitsXFTag

const edm::InputTag MuonTruth::CSCsimHitsXFTag

private

Definition at line 69 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

geomToken_

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

private

Definition at line 71 of file MuonTruth.h.

Referenced by initEvent().

linksTag

const edm::InputTag MuonTruth::linksTag

private

Definition at line 64 of file MuonTruth.h.

Referenced by initEvent().

linksToken_

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

private

Definition at line 74 of file MuonTruth.h.

Referenced by initEvent().

simHitsToken_

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

private

Definition at line 77 of file MuonTruth.h.

Referenced by initEvent(), and MuonTruth().

simHitsXFToken_

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

private

Definition at line 76 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(), initEvent(), muonFraction(), and simHits().

theDetId

unsigned int MuonTruth::theDetId

private

Definition at line 59 of file MuonTruth.h.

Referenced by analyze(), hitsFromSimTrack(), and initEvent().

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 79 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(), initEvent(), and muonFraction().

theWireDigiSimLinks

const DigiSimLinks* MuonTruth::theWireDigiSimLinks

private

Definition at line 62 of file MuonTruth.h.

Referenced by analyze(), and initEvent().

wireLinksTag

const edm::InputTag MuonTruth::wireLinksTag

private

Definition at line 65 of file MuonTruth.h.

Referenced by initEvent().

wireLinksToken_

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

private

Definition at line 75 of file MuonTruth.h.

Referenced by initEvent().