#include <TrackingTruthValid.h>

Inheritance diagram for TrackingTruthValid:

Public Types
typedef std::vector< TrackingParticle >	TrackingParticleCollection

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

void	endJob () override

	TrackingTruthValid (const edm::ParameterSet &conf)

	~TrackingTruthValid () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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 edm::one::EDProducerBase
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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

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

inlineoverride

Definition at line 21 of file TrackingTruthValid.h.

References analyze(), bookHistograms(), endJob(), and findQualityFiles::run.

21 {} ;

Member Function Documentation

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

override

Definition at line 82 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(), lumiQTWidget::t, and vec_TrackingParticle_Token_.

Referenced by ~TrackingTruthValid().

                                                                              {
 
   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::bookHistograms	(	DQMStore::IBooker &	ibooker,
		const edm::Run &	run,
		const edm::EventSetup &	es
	)

override

Definition at line 38 of file TrackingTruthValid.cc.

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

Referenced by ~TrackingTruthValid().

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

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 183 of file TrackingTruthValid.cc.

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

Referenced by ~TrackingTruthValid().

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