411 art::ServiceHandle<geo::Geometry const> geom;
412 art::ServiceHandle<art::TFileService>
tfs;
416 art::TFileDirectory topdir = tfs->mkdir(
"seedana",
"SeedAna histograms");
417 art::TFileDirectory
dir = topdir.mkdir(subdir);
422 dir.make<TH2F>(
"duvcosth",
"Delta(uv) vs. Colinearity", 100, 0.95, 1., 100, 0., 1.);
423 fHcosth = dir.make<TH1F>(
"colin",
"Colinearity", 100, 0.95, 1.);
424 fHmcu =
dir.make<TH1F>(
"mcu",
"MC Truth U", 100, -5., 5.);
425 fHmcv =
dir.make<TH1F>(
"mcv",
"MC Truth V", 100, -5., 5.);
426 fHmcw =
dir.make<TH1F>(
"mcw",
"MC Truth W", 100, -20., 20.);
427 fHmcdudw =
dir.make<TH1F>(
"mcdudw",
"MC Truth U Slope", 100, -0.2, 0.2);
428 fHmcdvdw =
dir.make<TH1F>(
"mcdvdw",
"MV Truth V Slope", 100, -0.2, 0.2);
431 "mcxstart",
"MC X Start Position", 10, -2. * geom->DetHalfWidth(), 4. * geom->DetHalfWidth());
433 "mcystart",
"MC Y Start Position", 10, -geom->DetHalfHeight(), geom->DetHalfHeight());
434 fHmcstartz =
dir.make<TH1F>(
"mczstart",
"MC Z Start Position", 10, 0., geom->DetLength());
436 "mcxend",
"MC X End Position", 10, -2. * geom->DetHalfWidth(), 4. * geom->DetHalfWidth());
438 "mcyend",
"MC Y End Position", 10, -geom->DetHalfHeight(), geom->DetHalfHeight());
439 fHmcendz =
dir.make<TH1F>(
"mczend",
"MC Z End Position", 10, 0., geom->DetLength());
440 fHmctheta =
dir.make<TH1F>(
"mctheta",
"MC Theta", 20, 0., 3.142);
441 fHmcphi =
dir.make<TH1F>(
"mcphi",
"MC Phi", 10, -3.142, 3.142);
442 fHmctheta_xz =
dir.make<TH1F>(
"mctheta_xz",
"MC Theta_xz", 40, -3.142, 3.142);
443 fHmctheta_yz =
dir.make<TH1F>(
"mctheta_yz",
"MC Theta_yz", 40, -3.142, 3.142);
444 fHmcmom =
dir.make<TH1F>(
"mcmom",
"MC Momentum", 10, 0., 10.);
445 fHmclen =
dir.make<TH1F>(
"mclen",
"MC Particle Length", 10, 0., 1.1 * geom->DetLength());
448 "Matched X Start Position",
450 -2. * geom->DetHalfWidth(),
451 4. * geom->DetHalfWidth());
453 "mystart",
"Matched Y Start Position", 10, -geom->DetHalfHeight(), geom->DetHalfHeight());
454 fHmstartz =
dir.make<TH1F>(
"mzstart",
"Matched Z Start Position", 10, 0., geom->DetLength());
456 "mxend",
"Matched X End Position", 10, -2. * geom->DetHalfWidth(), 4. * geom->DetHalfWidth());
458 "myend",
"Matched Y End Position", 10, -geom->DetHalfHeight(), geom->DetHalfHeight());
459 fHmendz =
dir.make<TH1F>(
"mzend",
"Matched Z End Position", 10, 0., geom->DetLength());
460 fHmtheta =
dir.make<TH1F>(
"mtheta",
"Matched Theta", 20, 0., 3.142);
461 fHmphi =
dir.make<TH1F>(
"mphi",
"Matched Phi", 10, -3.142, 3.142);
462 fHmtheta_xz =
dir.make<TH1F>(
"mtheta_xz",
"Matched Theta_xz", 40, -3.142, 3.142);
463 fHmtheta_yz =
dir.make<TH1F>(
"mtheta_yz",
"Matched Theta_yz", 40, -3.142, 3.142);
464 fHmmom =
dir.make<TH1F>(
"mmom",
"Matched Momentum", 10, 0., 10.);
465 fHmlen =
dir.make<TH1F>(
"mlen",
"Matched Particle Length", 10, 0., 1.1 * geom->DetLength());
468 "Good X Start Position",
470 -2. * geom->DetHalfWidth(),
471 4. * geom->DetHalfWidth());
473 "gystart",
"Good Y Start Position", 10, -geom->DetHalfHeight(), geom->DetHalfHeight());
474 fHgstartz =
dir.make<TH1F>(
"gzstart",
"Good Z Start Position", 10, 0., geom->DetLength());
476 "gxend",
"Good X End Position", 10, -2. * geom->DetHalfWidth(), 4. * geom->DetHalfWidth());
478 "gyend",
"Good Y End Position", 10, -geom->DetHalfHeight(), geom->DetHalfHeight());
479 fHgendz =
dir.make<TH1F>(
"gzend",
"Good Z End Position", 10, 0., geom->DetLength());
480 fHgtheta =
dir.make<TH1F>(
"gtheta",
"Good Theta", 20, 0., 3.142);
481 fHgphi =
dir.make<TH1F>(
"gphi",
"Good Phi", 10, -3.142, 3.142);
482 fHgtheta_xz =
dir.make<TH1F>(
"gtheta_xz",
"Good Theta_xz", 40, -3.142, 3.142);
483 fHgtheta_yz =
dir.make<TH1F>(
"gtheta_yz",
"Good Theta_yz", 40, -3.142, 3.142);
484 fHgmom =
dir.make<TH1F>(
"gmom",
"Good Momentum", 10, 0., 10.);
485 fHglen =
dir.make<TH1F>(
"glen",
"Good Particle Length", 10, 0., 1.1 * geom->DetLength());
488 "Multiplicity vs. X Start Position",
490 -2. * geom->DetHalfWidth(),
491 4. * geom->DetHalfWidth());
493 "Multiplicity vs. Y Start Position",
495 -geom->DetHalfHeight(),
496 geom->DetHalfHeight());
498 dir.make<TH1F>(
"mulzstart",
"Multiplicity vs. Z Start Position", 10, 0., geom->DetLength());
500 "Multiplicity vs. X End Position",
502 -2. * geom->DetHalfWidth(),
503 4. * geom->DetHalfWidth());
505 "Multiplicity vs. Y End Position",
507 -geom->DetHalfHeight(),
508 geom->DetHalfHeight());
510 dir.make<TH1F>(
"mulzend",
"Multiplicity vs. Z End Position", 10, 0., geom->DetLength());
511 fHmultheta =
dir.make<TH1F>(
"multheta",
"Multiplicity vs. Theta", 20, 0., 3.142);
512 fHmulphi =
dir.make<TH1F>(
"mulphi",
"Multiplicity vs. Phi", 10, -3.142, 3.142);
513 fHmultheta_xz =
dir.make<TH1F>(
"multheta_xz",
"Multiplicity vs. Theta_xz", 40, -3.142, 3.142);
514 fHmultheta_yz =
dir.make<TH1F>(
"multheta_yz",
"Multiplicity vs. Theta_yz", 40, -3.142, 3.142);
515 fHmulmom =
dir.make<TH1F>(
"mulmom",
"Multiplicity vs. Momentum", 10, 0., 10.);
517 dir.make<TH1F>(
"mullen",
"Multiplicity vs. Particle Length", 10, 0., 1.1 * geom->DetLength());
520 "Efficiency vs. X Start Position",
522 -2. * geom->DetHalfWidth(),
523 4. * geom->DetHalfWidth());
525 "Efficiency vs. Y Start Position",
527 -geom->DetHalfHeight(),
528 geom->DetHalfHeight());
530 dir.make<TH1F>(
"ezstart",
"Efficiency vs. Z Start Position", 10, 0., geom->DetLength());
532 "Efficiency vs. X End Position",
534 -2. * geom->DetHalfWidth(),
535 4. * geom->DetHalfWidth());
537 "eyend",
"Efficiency vs. Y End Position", 10, -geom->DetHalfHeight(), geom->DetHalfHeight());
538 fHeendz =
dir.make<TH1F>(
"ezend",
"Efficiency vs. Z End Position", 10, 0., geom->DetLength());
539 fHetheta =
dir.make<TH1F>(
"etheta",
"Efficiency vs. Theta", 20, 0., 3.142);
540 fHephi =
dir.make<TH1F>(
"ephi",
"Efficiency vs. Phi", 10, -3.142, 3.142);
541 fHetheta_xz =
dir.make<TH1F>(
"etheta_xz",
"Efficiency vs. Theta_xz", 40, -3.142, 3.142);
542 fHetheta_yz =
dir.make<TH1F>(
"etheta_yz",
"Efficiency vs. Theta_yz", 40, -3.142, 3.142);
543 fHemom =
dir.make<TH1F>(
"emom",
"Efficiency vs. Momentum", 10, 0., 10.);
545 dir.make<TH1F>(
"elen",
"Efficiency vs. Particle Length", 10, 0., 1.1 * geom->DetLength());
art::ServiceHandle< art::TFileService > tfs