596 art::TFileDirectory topdir = tfs->mkdir(
"trkana",
"TrackAnaCT histograms");
597 art::TFileDirectory
dir = topdir.mkdir(subdir);
601 fHduvcosth = dir.make<TH2F>(
"duvcosth",
"Delta(uv) vs. Colinearity",
602 100, 0.95, 1., 100, 0., 1.);
603 fHcosth = dir.make<TH1F>(
"colin",
"Colinearity", 100, 0.95, 1.);
604 fHmcu =
dir.make<TH1F>(
"mcu",
"MC Truth U", 100, -5., 5.);
605 fHmcv =
dir.make<TH1F>(
"mcv",
"MC Truth V", 100, -5., 5.);
606 fHmcw =
dir.make<TH1F>(
"mcw",
"MC Truth W", 100, -20., 20.);
607 fHupull =
dir.make<TH1F>(
"dupull",
"U Pull", 100, -20., 20.);
608 fHvpull =
dir.make<TH1F>(
"dvpull",
"V Pull", 100, -20., 20.);
609 fHmcdudw =
dir.make<TH1F>(
"mcdudw",
"MC Truth U Slope", 100, -0.2, 0.2);
610 fHmcdvdw =
dir.make<TH1F>(
"mcdvdw",
"MV Truth V Slope", 100, -0.2, 0.2);
611 fHdudwpull =
dir.make<TH1F>(
"dudwpull",
"U Slope Pull", 100, -10., 10.);
612 fHdvdwpull =
dir.make<TH1F>(
"dvdwpull",
"V Slope Pull", 100, -10., 10.);
613 fHstartdx =
dir.make<TH1F>(
"startdx",
"Start Delta x", 100, -10., 10.);
614 fHstartdy =
dir.make<TH1F>(
"startdy",
"Start Delta y", 100, -10., 10.);
615 fHstartdz =
dir.make<TH1F>(
"startdz",
"Start Delta z", 100, -10., 10.);
616 fHenddx =
dir.make<TH1F>(
"enddx",
"End Delta x", 100, -10., 10.);
617 fHenddy =
dir.make<TH1F>(
"enddy",
"End Delta y", 100, -10., 10.);
618 fHenddz =
dir.make<TH1F>(
"enddz",
"End Delta z", 100, -10., 10.);
619 fHlvsl =
dir.make<TH2F>(
"lvsl",
"Reco Length vs. MC Truth Length",
620 100, 0., 1.1 * geom->Cryostat(0).Length(), 100, 0., 1.1 * geom->Cryostat(0).Length());
621 fHdl =
dir.make<TH1F>(
"dl",
"Track Length Minus MC Particle Length", 100, -50., 50.);
622 fHpvsp =
dir.make<TH2F>(
"pvsp",
"Reco Momentum vs. MC Truth Momentum",
623 100, 0., 5., 100, 0., 5.);
624 fHpvspc =
dir.make<TH2F>(
"pvspc",
"Reco Momentum vs. MC Truth Momentum (Contained Tracks)",
625 100, 0., 5., 100, 0., 5.);
626 fHdp =
dir.make<TH1F>(
"dp",
"Reco-MC Momentum Difference", 100, -5., 5.);
627 fHdpc =
dir.make<TH1F>(
"dpc",
"Reco-MC Momentum Difference (Contained Tracks)",
629 fHppull =
dir.make<TH1F>(
"ppull",
"Momentum Pull", 100, -10., 10.);
630 fHppullc =
dir.make<TH1F>(
"ppullc",
"Momentum Pull (Contained Tracks)", 100, -10., 10.);
632 fHmcstartx =
dir.make<TH1F>(
"mcxstart",
"MC X Start Position",
633 10, -2.*geom->Cryostat(0).HalfWidth(), 4.*geom->Cryostat(0).HalfWidth());
634 fHmcstarty =
dir.make<TH1F>(
"mcystart",
"MC Y Start Position",
635 10, -geom->Cryostat(0).HalfHeight(), geom->Cryostat(0).HalfHeight());
636 fHmcstartz =
dir.make<TH1F>(
"mczstart",
"MC Z Start Position",
637 10, 0., geom->Cryostat(0).Length());
638 fHmcendx =
dir.make<TH1F>(
"mcxend",
"MC X End Position",
639 10, -2.*geom->Cryostat(0).HalfWidth(), 4.*geom->Cryostat(0).HalfWidth());
640 fHmcendy =
dir.make<TH1F>(
"mcyend",
"MC Y End Position",
641 10, -geom->Cryostat(0).HalfHeight(), geom->Cryostat(0).HalfHeight());
642 fHmcendz =
dir.make<TH1F>(
"mczend",
"MC Z End Position",
643 10, 0., geom->Cryostat(0).Length());
644 fHmctheta =
dir.make<TH1F>(
"mctheta",
"MC Theta", 20, 0., 3.142);
645 fHmcphi =
dir.make<TH1F>(
"mcphi",
"MC Phi", 10, -3.142, 3.142);
646 fHmctheta_xz =
dir.make<TH1F>(
"mctheta_xz",
"MC Theta_xz", 40, -3.142, 3.142);
647 fHmctheta_yz =
dir.make<TH1F>(
"mctheta_yz",
"MC Theta_yz", 40, -3.142, 3.142);
648 fHmcmom =
dir.make<TH1F>(
"mcmom",
"MC Momentum", 10, 0., 10.);
649 fHmclen =
dir.make<TH1F>(
"mclen",
"MC Particle Length", 10, 0., 1.1 * geom->Cryostat(0).Length());
651 fHgstartx =
dir.make<TH1F>(
"gxstart",
"Good X Start Position",
652 10, -2.*geom->Cryostat(0).HalfWidth(), 4.*geom->Cryostat(0).HalfWidth());
653 fHgstarty =
dir.make<TH1F>(
"gystart",
"Good Y Start Position",
654 10, -geom->Cryostat(0).HalfHeight(), geom->Cryostat(0).HalfHeight());
655 fHgstartz =
dir.make<TH1F>(
"gzstart",
"Good Z Start Position",
656 10, 0., geom->Cryostat(0).Length());
657 fHgendx =
dir.make<TH1F>(
"gxend",
"Good X End Position",
658 10, -2.*geom->Cryostat(0).HalfWidth(), 4.*geom->Cryostat(0).HalfWidth());
659 fHgendy =
dir.make<TH1F>(
"gyend",
"Good Y End Position",
660 10, -geom->Cryostat(0).HalfHeight(), geom->Cryostat(0).HalfHeight());
661 fHgendz =
dir.make<TH1F>(
"gzend",
"Good Z End Position",
662 10, 0., geom->Cryostat(0).Length());
663 fHgtheta =
dir.make<TH1F>(
"gtheta",
"Good Theta", 20, 0., 3.142);
664 fHgphi =
dir.make<TH1F>(
"gphi",
"Good Phi", 10, -3.142, 3.142);
665 fHgtheta_xz =
dir.make<TH1F>(
"gtheta_xz",
"Good Theta_xz", 40, -3.142, 3.142);
666 fHgtheta_yz =
dir.make<TH1F>(
"gtheta_yz",
"Good Theta_yz", 40, -3.142, 3.142);
667 fHgmom =
dir.make<TH1F>(
"gmom",
"Good Momentum", 10, 0., 10.);
668 fHglen =
dir.make<TH1F>(
"glen",
"Good Particle Length", 10, 0., 1.1 * geom->Cryostat(0).Length());
670 fHestartx =
dir.make<TH1F>(
"exstart",
"Efficiency vs. X Start Position",
671 10, -2.*geom->Cryostat(0).HalfWidth(), 4.*geom->Cryostat(0).HalfWidth());
672 fHestarty =
dir.make<TH1F>(
"eystart",
"Efficiency vs. Y Start Position",
673 10, -geom->Cryostat(0).HalfHeight(), geom->Cryostat(0).HalfHeight());
674 fHestartz =
dir.make<TH1F>(
"ezstart",
"Efficiency vs. Z Start Position",
675 10, 0., geom->Cryostat(0).Length());
676 fHeendx =
dir.make<TH1F>(
"exend",
"Efficiency vs. X End Position",
677 10, -2.*geom->Cryostat(0).HalfWidth(), 4.*geom->Cryostat(0).HalfWidth());
678 fHeendy =
dir.make<TH1F>(
"eyend",
"Efficiency vs. Y End Position",
679 10, -geom->Cryostat(0).HalfHeight(), geom->Cryostat(0).HalfHeight());
680 fHeendz =
dir.make<TH1F>(
"ezend",
"Efficiency vs. Z End Position",
681 10, 0., geom->Cryostat(0).Length());
682 fHetheta =
dir.make<TH1F>(
"etheta",
"Efficiency vs. Theta", 20, 0., 3.142);
683 fHephi =
dir.make<TH1F>(
"ephi",
"Efficiency vs. Phi", 10, -3.142, 3.142);
684 fHetheta_xz =
dir.make<TH1F>(
"etheta_xz",
"Efficiency vs. Theta_xz", 40, -3.142, 3.142);
685 fHetheta_yz =
dir.make<TH1F>(
"etheta_yz",
"Efficiency vs. Theta_yz", 40, -3.142, 3.142);
686 fHemom =
dir.make<TH1F>(
"emom",
"Efficiency vs. Momentum", 10, 0., 10.);
687 fHelen =
dir.make<TH1F>(
"elen",
"Efficiency vs. Particle Length",
688 10, 0., 1.1 * geom->Cryostat(0).Length());