larg4::ISCalcAna Class Reference

Inheritance diagram for larg4::ISCalcAna:

Public Member Functions
	ISCalcAna (fhicl::ParameterSet const &p)
	ISCalcAna (ISCalcAna const &)=delete
	ISCalcAna (ISCalcAna &&)=delete
ISCalcAna &	operator= (ISCalcAna const &)=delete
ISCalcAna &	operator= (ISCalcAna &&)=delete
void	analyze (art::Event const &event) override
void	beginJob () override
void	endJob () override

Private Attributes
std::unique_ptr< ISCalc >	fISAlg
art::InputTag	fEDepTag
art::InputTag	calcTag
CLHEP::HepRandomEngine &	fEngine
TNtuple *	fNtuple

Detailed Description

Definition at line 37 of file ISCalcAna_module.cc.

Constructor & Destructor Documentation

larg4::ISCalcAna::ISCalcAna ( fhicl::ParameterSet const &  p )

explicit

Definition at line 57 of file ISCalcAna_module.cc.

    : EDAnalyzer(pset)
    , fEDepTag{pset.get<art::InputTag>("SimulationLabel")}
    , calcTag{pset.get<art::InputTag>("ISCalcAlg")}
    , fEngine(art::ServiceHandle<rndm::NuRandomService>()
                ->createEngine(*this, "HepJamesRandom", "NEST", pset, "SeedNEST"))
  {
    std::cout << "ISCalcAna constructor." << std::endl;
  }

larg4::ISCalcAna::ISCalcAna ( ISCalcAna const &  )

delete

larg4::ISCalcAna::ISCalcAna ( ISCalcAna &&  )

delete

Member Function Documentation

void larg4::ISCalcAna::analyze ( art::Event const &  event )

override

Definition at line 97 of file ISCalcAna_module.cc.

  {
    art::Handle<std::vector<sim::SimEnergyDeposit>> edep_handle;
    if (!event.getByLabel(fEDepTag, edep_handle)) {
      std::cout << "PDFastSimPAR Module Cannot getByLabel: " << fEDepTag << std::endl;
      return;
    }

    auto const detProp =
      art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(event);
    for (auto const& edepi : *edep_handle) {
      auto const [energyDeposit, nElectrons, nPhotons, scintYieldRatio] =
        fISAlg->CalcIonAndScint(detProp, edepi);
      fNtuple->Fill(event.run(),
                    event.event(),
                    edepi.T(),
                    edepi.X(),
                    edepi.Y(),
                    edepi.Z(),
                    edepi.StepLength(),
                    edepi.Energy(),
                    edepi.TrackID(),
                    edepi.PdgCode(),
                    energyDeposit,
                    nElectrons,
                    nPhotons,
                    scintYieldRatio);
    }

    std::cout << "ISCalcAna analyze completed." << std::endl;
  }

void larg4::ISCalcAna::beginJob ( )

override

Definition at line 68 of file ISCalcAna_module.cc.

  {
    std::cout << "ISCalcAna beginJob." << std::endl;
    std::cout << "Using " << calcTag.label() << " algorithm to calculate IS." << std::endl;

    if (calcTag.label() == "Separate")
      fISAlg = std::make_unique<larg4::ISCalcSeparate>();
    else if (calcTag.label() == "Correlated") {
      auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService>()->DataForJob();
      fISAlg = std::make_unique<larg4::ISCalcCorrelated>(detProp, fEngine);
    }
    else if (calcTag.label() == "NEST")
      fISAlg = std::make_unique<larg4::ISCalcNESTLAr>(fEngine);
    else
      mf::LogWarning("IonAndScint") << "No ISCalculation set, this can't be good.";

    art::ServiceHandle<art::TFileService const> tfs;
    fNtuple = tfs->make<TNtuple>(
      "nt_is",
      "EDep IS Calc Ntuple",
      "run:event:t:x:y:z:ds:e:trackid:pdg:e_deposit:n_electron:n_photon:scintyield");
  }

void larg4::ISCalcAna::endJob ( )

override

Definition at line 91 of file ISCalcAna_module.cc.

  {
    std::cout << "ISCalcAna endJob." << std::endl;
  }

ISCalcAna& larg4::ISCalcAna::operator= ( ISCalcAna const &  )

delete

ISCalcAna& larg4::ISCalcAna::operator= ( ISCalcAna &&  )

delete

Member Data Documentation

art::InputTag larg4::ISCalcAna::calcTag

private

Definition at line 52 of file ISCalcAna_module.cc.

art::InputTag larg4::ISCalcAna::fEDepTag

private

Definition at line 51 of file ISCalcAna_module.cc.

CLHEP::HepRandomEngine& larg4::ISCalcAna::fEngine

private

Definition at line 53 of file ISCalcAna_module.cc.

std::unique_ptr<ISCalc> larg4::ISCalcAna::fISAlg

private

Definition at line 50 of file ISCalcAna_module.cc.

TNtuple* larg4::ISCalcAna::fNtuple

private

Definition at line 54 of file ISCalcAna_module.cc.

---

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

• ISCalcAna_module.cc