#include <TrackingTruthValid.h>

Inheritance diagram for TrackingTruthValid:

Public Types
typedef std::vector < TrackingParticle >	TrackingParticleCollection

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

void	beginJob (const edm::ParameterSet &conf)

void	beginRun (const edm::Run &, const edm::EventSetup &)

void	endJob ()

	TrackingTruthValid (const edm::ParameterSet &conf)

	~TrackingTruthValid ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	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

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

edm::EDGetTokenT < TrackingParticleCollection >	vec_TrackingParticle_Token_

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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 13 of file TrackingTruthValid.h.

Member Typedef Documentation

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

Definition at line 15 of file TrackingTruthValid.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 35 of file TrackingTruthValid.cc.

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

TrackingTruthValid::~TrackingTruthValid ( )

inline

Definition at line 19 of file TrackingTruthValid.h.

19 {} ;

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

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

                                                                              {
 
   edm::Handle<TrackingParticleCollection>  TruthTrackContainer ;
   //  edm::Handle<TrackingVertexCollection>    TruthVertexContainer;
 
   event.getByToken( vec_TrackingParticle_Token_, TruthTrackContainer );
   
   const TrackingParticleCollection *tPC   = TruthTrackContainer.product();
   //  const TrackingVertexCollection   *tVC   = TruthVertexContainer.product();
 
   /*
   // Get and print HepMC event for comparison
   edm::Handle<edm::HepMCProduct> hepMC;
   event.getByLabel("source",hepMC);
   const edm::HepMCProduct *mcp = hepMC.product();
   //  const HepMC::GenEvent *genEvent = mcp -> GetEvent();
   */
   //  cout << "Found " << tPC -> size() << " tracks and " << tVC -> size() << " vertices." <<endl;
    
 
 // 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 33 of file TrackingTruthValid.cc.

33 {}

void TrackingTruthValid::beginRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 40 of file TrackingTruthValid.cc.

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

                                                                        {
   dbe_  = edm::Service<DQMStore>().operator->();
   dbe_->setCurrentFolder("Tracking/TrackingMCTruth/TrackingParticle");
   
 
   meTPMass = dbe_->book1D("TPMass","Tracking Particle Mass",100, -1,+5.);
   meTPCharge = dbe_->book1D("TPCharge","Tracking Particle Charge",10, -5, 5);
   meTPId = dbe_->book1D("TPId","Tracking Particle Id",500, -5000, 5000);
   meTPProc = dbe_->book1D("TPProc","Tracking Particle Proc",20, -0.5, 19.5);
   meTPAllHits = dbe_->book1D("TPAllHits", "Tracking Particle All Hits", 200, -0.5, 199.5);
   meTPMatchedHits = dbe_->book1D("TPMatchedHits", "Tracking Particle Matched Hits", 100, -0.5, 99.5);
   meTPPt = dbe_->book1D("TPPt", "Tracking Particle Pt",100, 0, 100.);
   meTPEta = dbe_->book1D("TPEta", "Tracking Particle Eta",100, -7., 7.);
   meTPPhi = dbe_->book1D("TPPhi", "Tracking Particle Phi",100, -4., 4);
   meTPVtxX = dbe_->book1D("TPVtxX", "Tracking Particle VtxX",100, -100, 100.);
   meTPVtxY = dbe_->book1D("TPVtxY", "Tracking Particle VtxY",100, -100, 100.);
   meTPVtxZ = dbe_->book1D("TPVtxZ", "Tracking Particle VtxZ",100, -100, 100.);
   meTPtip = dbe_->book1D("TPtip", "Tracking Particle tip",100, 0, 1000.);
   meTPlip = dbe_->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 )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 197 of file TrackingTruthValid.cc.

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

                                { 
 
   if ( outputFile.size() != 0 && dbe_ ) dbe_->save(outputFile);
 
 } 

Member Data Documentation

DQMStore* TrackingTruthValid::dbe_

private

Definition at line 30 of file TrackingTruthValid.h.

Referenced by beginRun(), and endJob().

MonitorElement* TrackingTruthValid::meTPAllHits

private

Definition at line 35 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPCharge

private

Definition at line 32 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPEta

private

Definition at line 38 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPId

private

Definition at line 33 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPlip

private

Definition at line 44 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPMass

private

Definition at line 31 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPMatchedHits

private

Definition at line 36 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPPhi

private

Definition at line 39 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPProc

private

Definition at line 34 of file TrackingTruthValid.h.

Referenced by beginRun().

MonitorElement* TrackingTruthValid::meTPPt

private

Definition at line 37 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPtip

private

Definition at line 43 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPVtxX

private

Definition at line 40 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPVtxY

private

Definition at line 41 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* TrackingTruthValid::meTPVtxZ

private

Definition at line 42 of file TrackingTruthValid.h.

Referenced by analyze(), and beginRun().

std::string TrackingTruthValid::outputFile

private

Definition at line 28 of file TrackingTruthValid.h.

Referenced by endJob().

edm::EDGetTokenT<TrackingParticleCollection> TrackingTruthValid::vec_TrackingParticle_Token_

private

Definition at line 46 of file TrackingTruthValid.h.

Referenced by analyze().

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation