#include <TrackingTruthValid.h>

Inheritance diagram for TrackingTruthValid:

Public Types
typedef std::vector < TrackingParticle >	TrackingParticleCollection

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

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

void	beginJob (const edm::ParameterSet &conf)

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

void	endJob ()

	TrackingTruthValid (const edm::ParameterSet &conf)

	~TrackingTruthValid ()

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

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

Detailed Description

Definition at line 15 of file TrackingTruthValid.h.

Member Typedef Documentation

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

Definition at line 17 of file TrackingTruthValid.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 34 of file TrackingTruthValid.cc.

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

TrackingTruthValid::~TrackingTruthValid ( )

inline

Definition at line 21 of file TrackingTruthValid.h.

21 {} ;

Member Function Documentation

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

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 84 of file TrackingTruthValid.cc.

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

                                                                              {
 
   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
   */
 
 
 }

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

Definition at line 32 of file TrackingTruthValid.cc.

32 {}

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

virtual

Implements DQMEDAnalyzer.

Definition at line 40 of file TrackingTruthValid.cc.

References DQMStore::IBooker::book1D(), dbe_, meTPAllHits, meTPCharge, meTPEta, meTPId, meTPlip, meTPMass, meTPMatchedHits, meTPPhi, meTPProc, meTPPt, meTPtip, meTPVtxX, meTPVtxY, meTPVtxZ, cppFunctionSkipper::operator, MonitorElement::setBinLabel(), and DQMStore::IBooker::setCurrentFolder().

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

void TrackingTruthValid::endJob ( void )

Definition at line 185 of file TrackingTruthValid.cc.

References dbe_, outputFile, runStandalone, and DQMStore::save().

                                { 
   //Only in standalone mode save local root file 
   if (runStandalone &&  outputFile.size() != 0 && dbe_ ){dbe_->save(outputFile);}
 
 } 