larg::LArSimChannelAna Class Reference

Base class for creation of raw signals on wires. More...

Inheritance diagram for larg::LArSimChannelAna:

Public Member Functions
	LArSimChannelAna (fhicl::ParameterSet const &pset)

Private Member Functions
void	analyze (const art::Event &evt) override
void	ensureHists (unsigned int const nTimeSamples)

Private Attributes
std::string	fLArG4ModuleLabel
bool	initDone
TH1D *	fChargeXpos
	position of the MC Truth charge deposition More...
TH1D *	fChargeYpos
	position of the MC Truth charge deposition More...
TH1D *	fChargeZpos
	position of the MC Truth charge deposition More...
TH1D *	fTDC
	Which TDCs have activity. More...
TH1D *	fTDCsPerChannel
	Number of TDCs with activity. More...
TH1D *	fIDEsPerChannel
TH1D *	fElectrons
	Electrons in the whole channel entry. More...
TH1D *	fEnergy
	Energy in the whole channel entry. More...
TH1D *	fElectronsPerTDC
TH1D *	fEnergyPerTDC
TH1D *	fElectronsPerIDE
TH1D *	fEnergyPerIDE

Detailed Description

Base class for creation of raw signals on wires.

Definition at line 35 of file LArSimChannelAna_module.cc.

Constructor & Destructor Documentation

larg::LArSimChannelAna::LArSimChannelAna ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 77 of file LArSimChannelAna_module.cc.

    : EDAnalyzer(pset)
    , fLArG4ModuleLabel{pset.get<std::string>("LArGeantModuleLabel")}
    , initDone(false)
    , fChargeXpos()
    , fChargeYpos()
    , fChargeZpos()
    , fTDC()
    , fTDCsPerChannel()
    , fIDEsPerChannel()
    , fElectrons()
    , fEnergy()
    , fElectronsPerTDC()
    , fEnergyPerTDC()
    , fElectronsPerIDE()
    , fEnergyPerIDE()
  {}

Member Function Documentation

void larg::LArSimChannelAna::analyze ( const art::Event &  evt )

overrideprivate

Definition at line 137 of file LArSimChannelAna_module.cc.

  {
    if (evt.isRealData()) {
      throw cet::exception("LArSimChannelAna") << "Not for use on Data yet...\n";
    }

    auto const detProp =
      art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt);
    ensureHists(detProp.NumberTimeSamples());

    art::ServiceHandle<geo::Geometry const> geom;

    auto const& scVec = *evt.getValidHandle<std::vector<sim::SimChannel>>(fLArG4ModuleLabel);

    //++++++++++
    // Loop over the Chnnels and fill histograms
    //++++++++++
    unsigned int totalIDEs = 0;
    double totalElectrons = 0;
    double totalEnergy = 0;
    for (const auto& sc : scVec) {
      const auto& tdcidemap = sc.TDCIDEMap();
      fTDCsPerChannel->Fill(tdcidemap.size());

      for (auto const& [tdc, ideVec] : tdcidemap) {
        totalIDEs += ideVec.size();
        double tdcElectrons = 0;
        double tdcEnergy = 0;

        fTDC->Fill(tdc);

        for (const auto& ide : ideVec) {
          totalElectrons += ide.numElectrons;
          totalEnergy += ide.energy;
          tdcElectrons += ide.numElectrons;
          tdcEnergy += ide.energy;

          fChargeXpos->Fill(ide.x);
          fChargeYpos->Fill(ide.y);
          fChargeZpos->Fill(ide.z);
          fElectronsPerIDE->Fill(ide.numElectrons);
          fEnergyPerIDE->Fill(ide.energy);
        }
        fElectronsPerTDC->Fill(tdcElectrons);
        fEnergyPerTDC->Fill(tdcEnergy);
      }
    }
    fIDEsPerChannel->Fill(totalIDEs);
    fElectrons->Fill(totalElectrons);
    fEnergy->Fill(totalEnergy);
  } //end analyze method

void larg::LArSimChannelAna::ensureHists ( unsigned int const  nTimeSamples )

private

Definition at line 97 of file LArSimChannelAna_module.cc.

  {
    if (initDone) return; // Bail if we've already done this.
    initDone = true;      // Insure that we bail later on

    // get access to the TFile service
    art::ServiceHandle<art::TFileService const> tfs;
    // geometry data.
    art::ServiceHandle<geo::Geometry const> geom;

    // assumes all TPCs are the same
    double width = 2 * geom->TPC(0).HalfWidth();
    double halfHeight = geom->TPC(0).HalfHeight();
    double length = geom->TPC(0).Length();

    // Assumes microboone dimensions. Ideally we'd fix this later...
    fChargeXpos =
      tfs->make<TH1D>("hChargeXpos", "X charge depositions;X (cm);Events", 101, 0.0, width);
    fChargeYpos = tfs->make<TH1D>(
      "hChargeYpos", "Y charge depositions;Y (cm);Events", 101, -halfHeight, halfHeight);
    fChargeZpos =
      tfs->make<TH1D>("hChargeZpos", "Z charge depositions;Z (cm);Events", 101, 0.0, length);
    fTDC = tfs->make<TH1D>("hTDC", "Active TDC;TDCs;Events;", nTimeSamples, 0, nTimeSamples);
    fTDCsPerChannel = tfs->make<TH1D>(
      "hTDCsPerChannel", "TDCs per channel entry;# TDCs;Events", 128, 0, nTimeSamples);
    fIDEsPerChannel =
      tfs->make<TH1D>("hIDEsPerChannel", "IDE per channel entry;# IDEs;Events", 100, 0, 20000);
    fElectrons =
      tfs->make<TH1D>("hElectrons", "Electrons per channel;Electrons;Events", 100, 0, 2e7);
    fEnergy = tfs->make<TH1D>("hEnergy", "Energy per channel;energy;Events", 100, 0, 2500);
    fElectronsPerIDE =
      tfs->make<TH1D>("hElectronsPerIDE", "Electrons per IDE;Electrons;Events", 100, 0, 10000);
    fEnergyPerIDE = tfs->make<TH1D>("hEnergyPerIDE", "Energy per IDE;energy;Events", 100, 0, 50);
    fElectronsPerTDC =
      tfs->make<TH1D>("hElectronsPerTDC", "Electrons per TDC;Electrons;Events", 100, 0, 10000);
    fEnergyPerTDC = tfs->make<TH1D>("hEnergyPerTDC", "Energy per YDC;energy;Events", 100, 0, 50);
  }

Member Data Documentation

TH1D* larg::LArSimChannelAna::fChargeXpos

private

position of the MC Truth charge deposition

Definition at line 56 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fChargeYpos

private

position of the MC Truth charge deposition

Definition at line 57 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fChargeZpos

private

position of the MC Truth charge deposition

Definition at line 58 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fElectrons

private

Electrons in the whole channel entry.

Definition at line 65 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fElectronsPerIDE

private

Definition at line 71 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fElectronsPerTDC

private

Definition at line 68 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fEnergy

private

Energy in the whole channel entry.

Definition at line 66 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fEnergyPerIDE

private

Definition at line 72 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fEnergyPerTDC

private

Definition at line 69 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fIDEsPerChannel

private

Definition at line 63 of file LArSimChannelAna_module.cc.

std::string larg::LArSimChannelAna::fLArG4ModuleLabel

private

Definition at line 49 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fTDC

private

Which TDCs have activity.

Definition at line 60 of file LArSimChannelAna_module.cc.

TH1D* larg::LArSimChannelAna::fTDCsPerChannel

private

Number of TDCs with activity.

Definition at line 62 of file LArSimChannelAna_module.cc.

bool larg::LArSimChannelAna::initDone

private

Definition at line 54 of file LArSimChannelAna_module.cc.

---

The documentation for this class was generated from the following file:

• LArSimChannelAna_module.cc