|
|
|
#include <TrackingTruthValid.h>
Public Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
| void | beginJob (const edm::ParameterSet &conf) |
| void | endJob () |
| TrackingTruthValid (const edm::ParameterSet &conf) | |
| ~TrackingTruthValid () | |
Private Attributes | |
| edm::ParameterSet | conf_ |
| 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::InputTag | src_ |
Definition at line 20 of file TrackingTruthValid.h.
| TrackingTruthValid::TrackingTruthValid | ( | const edm::ParameterSet & | conf | ) | [explicit] |
Definition at line 36 of file TrackingTruthValid.cc.
References dbe_, edm::ParameterSet::getParameter(), cppFunctionSkipper::operator, and download_sqlite_cfg::outputFile.
{
outputFile = conf.getParameter<std::string>("outputFile");
src_ = conf.getParameter<edm::InputTag>( "src" );
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,"");
}
| TrackingTruthValid::~TrackingTruthValid | ( | ) | [inline] |
Definition at line 25 of file TrackingTruthValid.h.
{} ;
| 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, edm::Handle< T >::product(), mathSSE::sqrt(), lumiQTWidget::t, and DetId::Tracker.
{
using namespace std;
edm::Handle<TrackingParticleCollection> TruthTrackContainer ;
// edm::Handle<TrackingVertexCollection> TruthVertexContainer;
event.getByLabel(src_,TruthTrackContainer );
// event.getByLabel(src_,TruthVertexContainer);
// std::cout << "Using Collection " << src_ << std::endl;
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());
std::vector<PSimHit> trackerPSimHit( t->trackPSimHit(DetId::Tracker) );
meTPAllHits->Fill(trackerPSimHit.size());
//get the process of the first hit
if(trackerPSimHit.size() !=0) meTPProc->Fill( trackerPSimHit.front().processType());
meTPMatchedHits->Fill(t->matchedHit());
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 34 of file TrackingTruthValid.cc.
{}
| void TrackingTruthValid::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 201 of file TrackingTruthValid.cc.
References dbe_, download_sqlite_cfg::outputFile, and DQMStore::save().
{
if ( outputFile.size() != 0 && dbe_ ) dbe_->save(outputFile);
}
edm::ParameterSet TrackingTruthValid::conf_ [private] |
Definition at line 34 of file TrackingTruthValid.h.
DQMStore* TrackingTruthValid::dbe_ [private] |
Definition at line 33 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPAllHits [private] |
Definition at line 42 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPCharge [private] |
Definition at line 39 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPEta [private] |
Definition at line 45 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPId [private] |
Definition at line 40 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPlip [private] |
Definition at line 51 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPMass [private] |
Definition at line 38 of file TrackingTruthValid.h.
Definition at line 43 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPPhi [private] |
Definition at line 46 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPProc [private] |
Definition at line 41 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPPt [private] |
Definition at line 44 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPtip [private] |
Definition at line 50 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPVtxX [private] |
Definition at line 47 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPVtxY [private] |
Definition at line 48 of file TrackingTruthValid.h.
MonitorElement* TrackingTruthValid::meTPVtxZ [private] |
Definition at line 49 of file TrackingTruthValid.h.
std::string TrackingTruthValid::outputFile [private] |
Definition at line 35 of file TrackingTruthValid.h.
edm::InputTag TrackingTruthValid::src_ [private] |
Definition at line 36 of file TrackingTruthValid.h.
1.7.3