TrackingTruthValid Class Reference

#include <TrackingTruthValid.h>

Inheritance diagram for TrackingTruthValid:

Public Types
typedef std::vector< TrackingParticle >	TrackingParticleCollection
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
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 Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &es) override
	TrackingTruthValid (const edm::ParameterSet &conf)
	~TrackingTruthValid () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
DQMStore *	dbe_
MonitorElement *	meTPAllHits
MonitorElement *	meTPCharge
MonitorElement *	meTPEta
MonitorElement *	meTPId
MonitorElement *	meTPlip
MonitorElement *	meTPMass
MonitorElement *	meTPMatchedHits
MonitorElement *	meTPPhi
MonitorElement *	meTPProc
MonitorElement *	meTPPt
MonitorElement *	meTPtip
MonitorElement *	meTPVtxX
MonitorElement *	meTPVtxY
MonitorElement *	meTPVtxZ
std::string	outputFile
bool	runStandalone
edm::EDGetTokenT< TrackingParticleCollection >	vec_TrackingParticle_Token_

Additional Inherited Members
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 14 of file TrackingTruthValid.h.

Member Typedef Documentation

TrackingParticleCollection

typedef std::vector<TrackingParticle> TrackingTruthValid::TrackingParticleCollection

Definition at line 16 of file TrackingTruthValid.h.

Constructor & Destructor Documentation

TrackingTruthValid()

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

explicit

Definition at line 29 of file TrackingTruthValid.cc.

    : runStandalone(conf.getParameter<bool>("runStandalone")),
      outputFile(conf.getParameter<std::string>("outputFile")),
      dbe_(nullptr),
      vec_TrackingParticle_Token_(consumes<TrackingParticleCollection>(conf.getParameter<edm::InputTag>("src"))) {}

~TrackingTruthValid()

TrackingTruthValid::~TrackingTruthValid ( )

inlineoverride

Definition at line 20 of file TrackingTruthValid.h.

{};

Member Function Documentation

analyze()

void TrackingTruthValid::analyze ( const edm::Event &  event, const edm::EventSetup &  c  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 77 of file TrackingTruthValid.cc.

                                                                              {
  edm::Handle<TrackingParticleCollection> TruthTrackContainer;
  //  edm::Handle<TrackingVertexCollection>    TruthVertexContainer;
 
  event.getByToken(vec_TrackingParticle_Token_, TruthTrackContainer);
 
  const TrackingParticleCollection *tPC = TruthTrackContainer.product();
 
  // Loop over TrackingParticle's
  for (TrackingParticleCollection::const_iterator t = tPC->begin(); t != tPC->end(); ++t) {
    // if(t -> trackerPSimHit().size() ==0) cout << " Track with 0 SimHit " <<
    // endl;
 
    meTPMass->Fill(t->mass());
    meTPCharge->Fill(t->charge());
    meTPId->Fill(t->pdgId());
    meTPPt->Fill(sqrt(t->momentum().perp2()));
    meTPEta->Fill(t->momentum().eta());
    meTPPhi->Fill(t->momentum().Phi());
    //#warning "This file has been modified just to get it to compile without
    // any regard as to whether it still functions as intended" #ifdef
    // REMOVED_JUST_TO_GET_IT_TO_COMPILE__THIS_CODE_NEEDS_TO_BE_CHECKED
    //    std::vector<PSimHit> trackerPSimHit( t->trackPSimHit(DetId::Tracker)
    //    );
    //#endif
    meTPAllHits->Fill(t->numberOfTrackerHits());
    // get the process of the first hit
    //#warning "This file has been modified just to get it to compile without
    // any regard as to whether it still functions as intended" #ifdef
    // REMOVED_JUST_TO_GET_IT_TO_COMPILE__THIS_CODE_NEEDS_TO_BE_CHECKED
    //    if(trackerPSimHit.size() !=0) meTPProc->Fill(
    //    trackerPSimHit.front().processType());
    //#endif
 
    // there is no more the PSimHits collection !!! how to deal w/ the
    // processType ?
    //    if(t->numberOfTrackerHits() !=0) meTPProc->Fill(
    //    trackerPSimHit.front().processType());
 
    meTPMatchedHits->Fill(t->numberOfTrackerLayers());
    meTPVtxX->Fill(t->vx());
    meTPVtxY->Fill(t->vy());
    meTPVtxZ->Fill(t->vz());
    meTPtip->Fill(sqrt(t->vertex().perp2()));
    meTPlip->Fill(t->vz());
 
    /*
 // Compare momenta from sources
  cout << "T.P.   Track mass, Momentum, q , ID, & Event # "
        << t -> mass()  << " "
        << t -> p4()    << " " << t -> charge() << " "
        << t -> pdgId() << " "
        << t -> eventId().bunchCrossing() << "." << t -> eventId().event() <<
 endl;
 
  if(t->mass() < 0) cout << "======= WARNING, this particle has negative mass: "
 << t->mass()
                         << " and pdgId: " << t->pdgId() << endl;
  if(t->pdgId() == 0) cout << "======= WARNING, this particle has pdgId = 0: "
 << t->pdgId() << endl; cout << " Hits for this track: " << t ->
 trackerPSimHit().size() << endl;
  */
 
    /*
      std::cout << std::endl << "### Tracking Particle ###" << std::endl;
      std::cout << (*t) << std::endl;
      std::cout << "\t Tracker: " << t->trackerPSimHit().size() << std::endl;
      std::cout << "\t Muon: "    << t->muonPSimHit().size()    << std::endl;
      std::cout << (*t) << std::endl;
    */
  }  // End loop over TrackingParticle
 
  // Loop over TrackingVertex's
  /*
  cout << "Dumping sample vertex info" << endl;
  for (TrackingVertexCollection::const_iterator v = tVC -> begin(); v != tVC ->
  end(); ++v) { cout << " Vertex Position & Event #" << v -> position() << " "
  << v -> eventId().bunchCrossing() << "." << v -> eventId().event() << endl;
    cout << "  Associated with " << v -> daughterTracks().size() << " tracks" <<
  endl;
    // Get Geant and HepMC positions
    for (genv_iterator genV = v -> genVertices_begin(); genV != v ->
  genVertices_end(); ++genV) { cout << "  HepMC vertex position " <<
  (*(*genV)).position() << endl;
    }
    for (g4v_iterator g4V = v -> g4Vertices_begin(); g4V != v ->
  g4Vertices_end(); ++g4V) { cout << "  Geant vertex position " <<
  (*g4V).position() << endl;
      // Probably empty all the time, currently
    }
 
    // Loop over daughter track(s)
    for (tp_iterator iTP = v -> daughterTracks_begin(); iTP != v ->
  daughterTracks_end(); ++iTP) { cout << "  Daughter starts:      " <<
  (*(*iTP)).vertex(); for (g4t_iterator g4T  = (*(*iTP)).g4Track_begin(); g4T !=
  (*(*iTP)).g4Track_end(); ++g4T) { cout << " p " << g4T->momentum();
      }
      cout << endl;
    }
 
    // Loop over source track(s) (can be multiple since vertices are collapsed)
    for (tp_iterator iTP = v -> sourceTracks_begin(); iTP != v ->
  sourceTracks_end(); ++iTP) { cout << "  Source   starts: " <<
  (*(*iTP)).vertex(); for (g4t_iterator g4T  = (*iTP)->g4Track_begin(); g4T !=
  (*iTP)->g4Track_end(); ++g4T) { cout << ", p " <<  g4T ->momentum();
      }
      cout << endl;
    }
  }  // End loop over TrackingVertex
  */
}

References dqm::impl::MonitorElement::Fill(), meTPAllHits, meTPCharge, meTPEta, meTPId, meTPlip, meTPMass, meTPMatchedHits, meTPPhi, meTPPt, meTPtip, meTPVtxX, meTPVtxY, meTPVtxZ, edm::Handle< T >::product(), mathSSE::sqrt(), OrderedSet::t, and vec_TrackingParticle_Token_.

bookHistograms()

void TrackingTruthValid::bookHistograms ( DQMStore::IBooker &  ibooker, const edm::Run &  run, const edm::EventSetup &  es  )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 35 of file TrackingTruthValid.cc.

                                                                                                            {
  dbe_ = edm::Service<DQMStore>().operator->();
  ibooker.setCurrentFolder("Tracking/TrackingMCTruth/TrackingParticle");
 
  meTPMass = ibooker.book1D("TPMass", "Tracking Particle Mass", 100, -1, +5.);
  meTPCharge = ibooker.book1D("TPCharge", "Tracking Particle Charge", 10, -5, 5);
  meTPId = ibooker.book1D("TPId", "Tracking Particle Id", 500, -5000, 5000);
  meTPProc = ibooker.book1D("TPProc", "Tracking Particle Proc", 20, -0.5, 19.5);
  meTPAllHits = ibooker.book1D("TPAllHits", "Tracking Particle All Hits", 200, -0.5, 199.5);
  meTPMatchedHits = ibooker.book1D("TPMatchedHits", "Tracking Particle Matched Hits", 100, -0.5, 99.5);
  meTPPt = ibooker.book1D("TPPt", "Tracking Particle Pt", 100, 0, 100.);
  meTPEta = ibooker.book1D("TPEta", "Tracking Particle Eta", 100, -7., 7.);
  meTPPhi = ibooker.book1D("TPPhi", "Tracking Particle Phi", 100, -4., 4);
  meTPVtxX = ibooker.book1D("TPVtxX", "Tracking Particle VtxX", 100, -100, 100.);
  meTPVtxY = ibooker.book1D("TPVtxY", "Tracking Particle VtxY", 100, -100, 100.);
  meTPVtxZ = ibooker.book1D("TPVtxZ", "Tracking Particle VtxZ", 100, -100, 100.);
  meTPtip = ibooker.book1D("TPtip", "Tracking Particle tip", 100, 0, 1000.);
  meTPlip = ibooker.book1D("TPlip", "Tracking Particle lip", 100, 0, 100.);
 
  // Prepare Axes Labels for Processes
  meTPProc->setBinLabel(1, "Undefined");              // value =  0
  meTPProc->setBinLabel(2, "Unknown");                // value =  1
  meTPProc->setBinLabel(3, "Primary");                // value =  2
  meTPProc->setBinLabel(4, "Hadronic");               // value =  3
  meTPProc->setBinLabel(5, "Decay");                  // value =  4
  meTPProc->setBinLabel(6, "Compton");                // value =  5
  meTPProc->setBinLabel(7, "Annihilation");           // value =  6
  meTPProc->setBinLabel(8, "EIoni");                  // value =  7
  meTPProc->setBinLabel(9, "HIoni");                  // value =  8
  meTPProc->setBinLabel(10, "MuIoni");                // value =  9
  meTPProc->setBinLabel(11, "Photon");                // value = 10
  meTPProc->setBinLabel(12, "MuPairProd");            // value = 11
  meTPProc->setBinLabel(13, "Conversions");           // value = 12
  meTPProc->setBinLabel(14, "EBrem");                 // value = 13
  meTPProc->setBinLabel(15, "SynchrotronRadiation");  // value = 14
  meTPProc->setBinLabel(16, "MuBrem");                // value = 15
  meTPProc->setBinLabel(17, "MuNucl");                // value = 16
  meTPProc->setBinLabel(18, "");
  meTPProc->setBinLabel(19, "");
  meTPProc->setBinLabel(20, "");
}

References dqm::implementation::IBooker::book1D(), dbe_, meTPAllHits, meTPCharge, meTPEta, meTPId, meTPlip, meTPMass, meTPMatchedHits, meTPPhi, meTPProc, meTPPt, meTPtip, meTPVtxX, meTPVtxY, meTPVtxZ, dqm::impl::MonitorElement::setBinLabel(), and dqm::implementation::NavigatorBase::setCurrentFolder().

Member Data Documentation

dbe_

DQMStore* TrackingTruthValid::dbe_

private

Definition at line 30 of file TrackingTruthValid.h.

Referenced by bookHistograms().

meTPAllHits

MonitorElement* TrackingTruthValid::meTPAllHits

private

Definition at line 35 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPCharge

MonitorElement* TrackingTruthValid::meTPCharge

private

Definition at line 32 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPEta

MonitorElement* TrackingTruthValid::meTPEta

private

Definition at line 38 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPId

MonitorElement* TrackingTruthValid::meTPId

private

Definition at line 33 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPlip

MonitorElement* TrackingTruthValid::meTPlip

private

Definition at line 44 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPMass

MonitorElement* TrackingTruthValid::meTPMass

private

Definition at line 31 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPMatchedHits

MonitorElement* TrackingTruthValid::meTPMatchedHits

private

Definition at line 36 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPPhi

MonitorElement* TrackingTruthValid::meTPPhi

private

Definition at line 39 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPProc

MonitorElement* TrackingTruthValid::meTPProc

private

Definition at line 34 of file TrackingTruthValid.h.

Referenced by bookHistograms().

meTPPt

MonitorElement* TrackingTruthValid::meTPPt

private

Definition at line 37 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPtip

MonitorElement* TrackingTruthValid::meTPtip

private

Definition at line 43 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPVtxX

MonitorElement* TrackingTruthValid::meTPVtxX

private

Definition at line 40 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPVtxY

MonitorElement* TrackingTruthValid::meTPVtxY

private

Definition at line 41 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

meTPVtxZ

MonitorElement* TrackingTruthValid::meTPVtxZ

private

Definition at line 42 of file TrackingTruthValid.h.

Referenced by analyze(), and bookHistograms().

outputFile

std::string TrackingTruthValid::outputFile

private

Definition at line 28 of file TrackingTruthValid.h.

runStandalone

bool TrackingTruthValid::runStandalone

private

Definition at line 27 of file TrackingTruthValid.h.

vec_TrackingParticle_Token_

edm::EDGetTokenT<TrackingParticleCollection> TrackingTruthValid::vec_TrackingParticle_Token_

private

Definition at line 46 of file TrackingTruthValid.h.

Referenced by analyze().