Public Member Functions
	RecoHists (const std::string &subdir)

Public Attributes
TH1F *	fHstartx {nullptr}

TH1F *	fHstarty {nullptr}

TH1F *	fHstartz {nullptr}

TH1F *	fHstartd {nullptr}

TH1F *	fHendx {nullptr}

TH1F *	fHendy {nullptr}

TH1F *	fHendz {nullptr}

TH1F *	fHendd {nullptr}

TH1F *	fHtheta {nullptr}

TH1F *	fHphi {nullptr}

TH1F *	fHtheta_xz {nullptr}

TH1F *	fHtheta_yz {nullptr}

TH1F *	fHmom {nullptr}

TH1F *	fHmoml {nullptr}

TH1F *	fHlen {nullptr}

TH1F *	fHlens {nullptr}

TH1F *	fHHitChg {nullptr}

TH1F *	fHHitWidth {nullptr}

TH1F *	fHHitPdg {nullptr}

TH1F *	fHHitTrkId {nullptr}

TH1F *	fModeFrac {nullptr}

TH1F *	fNTrkIdTrks {nullptr}

TH2F *	fNTrkIdTrks2 {nullptr}

TH2F *	fNTrkIdTrks3 {nullptr}

Detailed Description

Definition at line 256 of file TrackAna_module.cc.

Constructor & Destructor Documentation

trkf::TrackAna::RecoHists::RecoHists ( const std::string & subdir )

explicit

Definition at line 439 of file TrackAna_module.cc.

References geo::GeometryCore::DetHalfHeight(), geo::GeometryCore::DetHalfWidth(), geo::GeometryCore::DetLength(), and dir.

   {
     // Get services.
 
     art::ServiceHandle<geo::Geometry const> geom;
     art::ServiceHandle<art::TFileService> tfs;
 
     // Make histogram directory.
 
     art::TFileDirectory topdir = tfs->mkdir("trkana", "TrackAna histograms");
     art::TFileDirectory dir = topdir.mkdir(subdir);
 
     // Book histograms.
 
     fHstartx = dir.make<TH1F>(
       "xstart", "X Start Position", 100, -2. * geom->DetHalfWidth(), 4. * geom->DetHalfWidth());
     fHstarty = dir.make<TH1F>(
       "ystart", "Y Start Position", 100, -geom->DetHalfHeight(), geom->DetHalfHeight());
     fHstartz = dir.make<TH1F>("zstart", "Z Start Position", 100, 0., geom->DetLength());
     fHstartd = dir.make<TH1F>(
       "dstart", "Start Position Distance to Boundary", 100, -10., geom->DetHalfWidth());
     fHendx = dir.make<TH1F>(
       "xend", "X End Position", 100, -2. * geom->DetHalfWidth(), 4. * geom->DetHalfWidth());
     fHendy =
       dir.make<TH1F>("yend", "Y End Position", 100, -geom->DetHalfHeight(), geom->DetHalfHeight());
     fHendz = dir.make<TH1F>("zend", "Z End Position", 100, 0., geom->DetLength());
     fHendd =
       dir.make<TH1F>("dend", "End Position Distance to Boundary", 100, -10., geom->DetHalfWidth());
     fHtheta = dir.make<TH1F>("theta", "Theta", 100, 0., 3.142);
     fHphi = dir.make<TH1F>("phi", "Phi", 100, -3.142, 3.142);
     fHtheta_xz = dir.make<TH1F>("theta_xz", "Theta_xz", 100, -3.142, 3.142);
     fHtheta_yz = dir.make<TH1F>("theta_yz", "Theta_yz", 100, -3.142, 3.142);
     fHmom = dir.make<TH1F>("mom", "Momentum", 100, 0., 10.);
     fHmoml = dir.make<TH1F>("moml", "Momentum", 100, 0., 1.);
     fHlen = dir.make<TH1F>("len", "Track Length", 100, 0., 1.1 * geom->DetLength());
     fHlens = dir.make<TH1F>("lens", "Track Length", 100, 0., 0.1 * geom->DetLength());
     fHHitChg = dir.make<TH1F>("hchg", "Hit Charge (ADC counts)", 100, 0., 4000.);
     fHHitWidth = dir.make<TH1F>("hwid", "Hit Width (ticks)", 40, 0., 20.);
     fHHitPdg = dir.make<TH1F>("hpdg", "Hit Pdg code", 5001, -2500.5, +2500.5);
     fHHitTrkId = dir.make<TH1F>("htrkid", "Hit Track ID", 401, -200.5, +200.5);
     fModeFrac =
       dir.make<TH1F>("hmodefrac",
                      "quasi-Purity: Fraction of component tracks with the Track mode value",
                      20,
                      0.01,
                      1.01);
     fNTrkIdTrks =
       dir.make<TH1F>("hntrkids",
                      "quasi-Efficiency: Number of stitched tracks in which TrkId appears",
                      20,
                      0.,
                      +10.0);
     fNTrkIdTrks2 = dir.make<TH2F>("hntrkids2",
                                   "Number of stitched tracks in which TrkId appears vs KE [GeV]",
                                   20,
                                   0.,
                                   +10.0,
                                   20,
                                   0.0,
                                   1.5);
     fNTrkIdTrks3 = dir.make<TH2F>("hntrkids3",
                                   "MC Track vs Reco Track, wtd by nhits on Collection Plane",
                                   10,
                                   -0.5,
                                   9.5,
                                   10,
                                   -0.5,
                                   9.5);
     fNTrkIdTrks3->GetXaxis()->SetTitle("Sorted-by-Descending-CPlane-Hits outer Track Number");
     fNTrkIdTrks3->GetYaxis()->SetTitle("Sorted-by-Descending-True-Length G4Track");
   }