trk::DemoShell Class Reference

#include <DemoShell.h>

Inheritance diagram for trk::DemoShell:

List of all members.

Public Member Functions
	DemoShell (const char *version)
void	Update (const cfg::Config &c)
jobc::Result	Reco (edm::EventHandle &evt)
jobc::Result	Ana (const edm::EventHandle &evt)
Private Attributes
int	fDebug
	Controls debug output 0=none, 99=scream.

---

Constructor & Destructor Documentation

DemoShell::DemoShell (  const char *  version  )

Definition at line 20 of file DemoShell.cxx. References jobc::Module::SetCfgVersion(). : jobc::Module("trk::DemoShell"), fDebug(0) { this->SetCfgVersion(version); }

---

Member Function Documentation

jobc::Result DemoShell::Ana (  const edm::EventHandle &  evt  )  [virtual]

The analysis method is intended to perform checks of the performance of the reconstruction for the purposes of optimization, testing, and validation. Parameters: evt  - Read-only access to the event data Reimplemented from jobc::Module. Definition at line 45 of file DemoShell.cxx. { return jobc::kPassed; }

jobc::Result DemoShell::Reco (  edm::EventHandle &  evt  )  [virtual]

Reconstruction method for this module. Provides read-write access to the event data. Parameters: evt  - Read-write access to the event data Reimplemented from jobc::Module. Definition at line 58 of file DemoShell.cxx. References recobase::Prong::Add(), recobase::Cluster::Add(), edm::EventHandle::Cal(), recobase::CellHit::Cell(), fDebug, util::EDMUtils::GetGeometry(), recobase::CellHit::Plane(), geo::Geometry::Plane(), edm::EventHandle::Reco(), recobase::Prong::SetDir(), recobase::Prong::SetStart(), recobase::Prong::SetStop(), and recobase::CellHit::View(). { // // "hits" is an C++ STL vector of CellHits (ok, really pointer to // CellHits...). The vector will be populated by data from the event // std::vector<const recobase::CellHit*> hits; // Try to pull the CellHits out of the event. std::string fInputCalHits = "./"; try { evt.Cal().Get(fInputCalHits.c_str(), hits); } catch (...) { // Failed to find any CalHits. Fail the event and keep processing if (fDebug>0) { std::cout << __FILE__ << ":" << __LINE__ << "Failed to find any CalHits" << std::endl; } return jobc::kFailed; } // If debugging is turned on, report that we found some hits to work // with if (fDebug>0) { std::cout << __FILE__ << ":" << __LINE__ << " Found " << hits.size() << " hits in event." << std::endl; } // // To give an example of how to work with the CellHit I'll show how // to do a least squares fit. We'll need some information about the // cell locations. That comes from the Geometry object. // const geo::Geometry& geom = util::EDMUtils::GetGeometry(evt); // We now have the basics of what we need. For example, I can fit a // straight line double xzViewSx = 0.0; double xzViewSx2 = 0.0; double xzViewSy = 0.0; double xzViewSxy = 0.0; double xzViewN = 0.0; double yzViewSx = 0.0; double yzViewSx2 = 0.0; double yzViewSy = 0.0; double yzViewSxy = 0.0; double yzViewN = 0.0; double xyz1[3] = { 0,0, 99.E99}; // Upstream point in hit cluster double xyz2[3] = { 0,0,-99.E99}; // Downstream point in hit cluster for (unsigned int i=0; i<hits.size(); ++i) { const recobase::CellHit* h = hits[i]; // Get the address of the ith hit once // Load the xyz location of the cell center into "xyz". Since we // have no information about the transverse position leave the // optional "localz" argument off double xyz[3]; // xyz[0]=x, xyz[1]=y, xyz[2]=z geom.Plane(h->Plane()).Cell(h->Cell()).GetCenter(xyz); // Store the direction fit information into a prong object // Keep track of point which is most upstream in z and most downstream in z if (xyz[2]<xyz1[2]) { xyz1[0] = xyz[0]; xyz1[1] = xyz[1]; xyz1[2] = xyz[2]; } if (xyz[2]>xyz2[2]) { xyz2[0] = xyz[0]; xyz2[1] = xyz[1]; xyz2[2] = xyz[2]; } // Accumulate sums for least squares fit in the X and Y views if (h->View()==geo::kX) { xzViewSx += xyz[2]; // in the x-z view z is the independent variable xzViewSx2 += xyz[2]*xyz[2]; xzViewSy += xyz[0]; // in the x-z view x is the dependent variable xzViewSxy += xyz[2]*xyz[0]; xzViewN += 1.0; } if (h->View()==geo::kY) { yzViewSx += xyz[2]; // in the y-z view z is the independent variable yzViewSx2 += xyz[2]*xyz[2]; yzViewSy += xyz[1]; // in the y-z view y is the dependent variable yzViewSxy += xyz[2]*xyz[1]; yzViewN += 1.0; } } double xzDenom = xzViewN*xzViewSx2 - xzViewSx*xzViewSx; double yzDenom = yzViewN*yzViewSx2 - yzViewSx*xzViewSx; // Intercepts double xzA = (xzViewSy*xzViewSx2-xzViewSx*xzViewSxy)/xzDenom; double yzA = (yzViewSy*yzViewSx2-yzViewSx*yzViewSxy)/yzDenom; // Slopes double xzB = (xzViewN*xzViewSxy - xzViewSx*xzViewSy)/xzDenom; double yzB = (yzViewN*yzViewSxy - yzViewSx*xzViewSy)/yzDenom; if (fDebug>0) { std::cout << __FILE__ << ":" << __LINE__ << " XZ(a,b)=(" << xzA << "," << xzB << ")" << std::endl; std::cout << __FILE__ << ":" << __LINE__ << " YZ(a,b)=(" << yzA << "," << yzB << ")" << std::endl; } // Now pack the data into objects other modules can use // Store the collection of hits used as a cluster. In this simple // case, I just add every hits to the cluster. A real algorithm // would have to make some more intelligent decisions. recobase::Cluster clust; for (unsigned int i=0; i<hits.size(); ++i) clust.Add(hits[i]); // Store information about the prong we've just formed recobase::Prong prong; prong.Add(&clust); // Set the start and stop points based on the line fit xyz1[0] = xzA + xzB*xyz1[2]; xyz1[1] = yzA + yzB*xyz1[2]; xyz2[0] = xzA + xzB*xyz2[2]; xyz2[1] = yzA + yzB*xyz2[2]; prong.SetStart(xyz1); prong.SetStop (xyz2); double dir[3]; double norm = sqrt(xzB*xzB + yzB*yzB + 1.0); dir[0] = xzB/norm; dir[1] = yzB/norm; dir[2] = 1.0/norm; prong.SetDir(dir); // And finally, make a track. Currently this adds little beyond what // is already in the prong, but will eventually be further // developed. recobase::Track track(&prong); // Finally, store the results in the event if (evt.Reco().GetFolder("./demoshell")==0) { evt.Reco().MakeFolder("./demoshell"); } evt.Reco().Put(clust,"./demoshell"); evt.Reco().Put(prong,"./demoshell"); evt.Reco().Put(track,"./demoshell"); return jobc::kPassed; }

void DemoShell::Update (  const cfg::Config &  c  )  [virtual]

Implements cfg::Observer. Definition at line 29 of file DemoShell.cxx. References fDebug. { c("Debug").Get(fDebug); // Unpack other configuration data following the above example // ... }

---

Member Data Documentation

int trk::DemoShell::fDebug [private]

Controls debug output 0=none, 99=scream. Definition at line 26 of file DemoShell.h. Referenced by Reco(), and Update().

---

The documentation for this class was generated from the following files:

• DemoShell.h
• DemoShell.cxx

---