453 art::ServiceHandle<geo::Geometry const> geom;
454 art::ServiceHandle<art::TFileService>
tfs;
458 art::TFileDirectory topdir = tfs->mkdir(
"trkana",
"TrackAna histograms");
459 art::TFileDirectory
dir = topdir.mkdir(subdir);
464 "xstart",
"X Start Position", 100, -2. * geom->DetHalfWidth(), 4. * geom->DetHalfWidth());
466 "ystart",
"Y Start Position", 100, -geom->DetHalfHeight(), geom->DetHalfHeight());
467 fHstartz = dir.make<TH1F>(
"zstart",
"Z Start Position", 100, 0., geom->DetLength());
469 "dstart",
"Start Position Distance to Boundary", 100, -10., geom->DetHalfWidth());
471 "xend",
"X End Position", 100, -2. * geom->DetHalfWidth(), 4. * geom->DetHalfWidth());
473 dir.make<TH1F>(
"yend",
"Y End Position", 100, -geom->DetHalfHeight(), geom->DetHalfHeight());
474 fHendz = dir.make<TH1F>(
"zend",
"Z End Position", 100, 0., geom->DetLength());
476 dir.make<TH1F>(
"dend",
"End Position Distance to Boundary", 100, -10., geom->DetHalfWidth());
477 fHtheta = dir.make<TH1F>(
"theta",
"Theta", 100, 0., 3.142);
478 fHphi = dir.make<TH1F>(
"phi",
"Phi", 100, -3.142, 3.142);
479 fHtheta_xz = dir.make<TH1F>(
"theta_xz",
"Theta_xz", 100, -3.142, 3.142);
480 fHtheta_yz = dir.make<TH1F>(
"theta_yz",
"Theta_yz", 100, -3.142, 3.142);
481 fHmom = dir.make<TH1F>(
"mom",
"Momentum", 100, 0., 10.);
482 fHmoml = dir.make<TH1F>(
"moml",
"Momentum", 100, 0., 1.);
483 fHlen = dir.make<TH1F>(
"len",
"Track Length", 100, 0., 1.1 * geom->DetLength());
484 fHlens = dir.make<TH1F>(
"lens",
"Track Length", 100, 0., 0.1 * geom->DetLength());
485 fHHitChg = dir.make<TH1F>(
"hchg",
"Hit Charge (ADC counts)", 100, 0., 4000.);
486 fHHitWidth = dir.make<TH1F>(
"hwid",
"Hit Width (ticks)", 40, 0., 20.);
487 fHHitPdg = dir.make<TH1F>(
"hpdg",
"Hit Pdg code", 5001, -2500.5, +2500.5);
488 fHHitTrkId = dir.make<TH1F>(
"htrkid",
"Hit Track ID", 401, -200.5, +200.5);
490 dir.make<TH1F>(
"hmodefrac",
491 "quasi-Purity: Fraction of component tracks with the Track mode value",
496 dir.make<TH1F>(
"hntrkids",
497 "quasi-Efficiency: Number of stitched tracks in which TrkId appears",
502 "Number of stitched tracks in which TrkId appears vs KE [GeV]",
510 "MC Track vs Reco Track, wtd by nhits on Collection Plane",
517 fNTrkIdTrks3->GetXaxis()->SetTitle(
"Sorted-by-Descending-CPlane-Hits outer Track Number");
518 fNTrkIdTrks3->GetYaxis()->SetTitle(
"Sorted-by-Descending-True-Length G4Track");
art::ServiceHandle< art::TFileService > tfs