Public Member Functions
	RecoHists ()

	RecoHists (const std::string &subdir)

Public Attributes
TH1F *	fHstartx

TH1F *	fHstarty

TH1F *	fHstartz

TH1F *	fHstartd

TH1F *	fHendx

TH1F *	fHendy

TH1F *	fHendz

TH1F *	fHendd

TH1F *	fHtheta

TH1F *	fHphi

TH1F *	fHtheta_xz

TH1F *	fHtheta_yz

TH1F *	fHmom

TH1F *	fHmoml

TH1F *	fHlen

TH1F *	fHlens

TH1F *	fHHitChg

TH1F *	fHHitWidth

TH1F *	fHHitPdg

TH1F *	fHHitTrkId

TH1F *	fModeFrac

TH1F *	fNTrkIdTrks

TH2F *	fNTrkIdTrks2

TH2F *	fNTrkIdTrks3

Detailed Description

Definition at line 278 of file TrackAna_module.cc.

Constructor & Destructor Documentation

trkf::TrackAna::RecoHists::RecoHists ( )

Definition at line 472 of file TrackAna_module.cc.

                                :
     //
     // Purpose: Default constructor.
     //
     fHstartx(0),
     fHstarty(0),
     fHstartz(0),
     fHstartd(0),
     fHendx(0),
     fHendy(0),
     fHendz(0),
     fHendd(0),
     fHtheta(0),
     fHphi(0),
     fHtheta_xz(0),
     fHtheta_yz(0),
     fHmom(0),
     fHmoml(0),
     fHlen(0),
     fHlens(0)
     ,fHHitChg(0)
     ,fHHitWidth(0)
     ,fHHitPdg(0)
     ,fHHitTrkId(0)
     ,fModeFrac(0)
     ,fNTrkIdTrks(0)
     ,fNTrkIdTrks2(0)
     ,fNTrkIdTrks3(0)
   {}

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

Definition at line 502 of file TrackAna_module.cc.

References geo::GeometryCore::DetHalfHeight(), geo::GeometryCore::DetHalfWidth(), geo::GeometryCore::DetLength(), dir, art::TFileDirectory::make(), and art::TFileDirectory::mkdir().

   {
     // Make sure all histogram pointers are initially zero.
 
     *this = RecoHists();
 
     // Get services.
 
     art::ServiceHandle<geo::Geometry> 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");
  
   }