spacecharge::ShiftEdepSCE Class Reference

Inheritance diagram for spacecharge::ShiftEdepSCE:

Public Member Functions
	ShiftEdepSCE (fhicl::ParameterSet const &p)
	ShiftEdepSCE (ShiftEdepSCE const &)=delete
	ShiftEdepSCE (ShiftEdepSCE &&)=delete
ShiftEdepSCE &	operator= (ShiftEdepSCE const &)=delete
ShiftEdepSCE &	operator= (ShiftEdepSCE &&)=delete
void	produce (art::Event &e) override
void	beginJob () override

Private Attributes
art::InputTag	fEDepTag
bool	fMakeAnaTree
TNtuple *	fNtEdepAna
larg4::ISCalcSeparate	fISAlg

Detailed Description

Definition at line 35 of file ShiftEdepSCE_module.cc.

Constructor & Destructor Documentation

spacecharge::ShiftEdepSCE::ShiftEdepSCE ( fhicl::ParameterSet const &  p )

explicit

Definition at line 61 of file ShiftEdepSCE_module.cc.

  : EDProducer{p}
  , fEDepTag(p.get<art::InputTag>("EDepTag"))
  , fMakeAnaTree(p.get<bool>("MakeAnaTree", true))
{
  produces<std::vector<sim::SimEnergyDeposit>>();
}

spacecharge::ShiftEdepSCE::ShiftEdepSCE ( ShiftEdepSCE const &  )

delete

spacecharge::ShiftEdepSCE::ShiftEdepSCE ( ShiftEdepSCE &&  )

delete

Member Function Documentation

void spacecharge::ShiftEdepSCE::beginJob ( )

override

Definition at line 70 of file ShiftEdepSCE_module.cc.

{
  if (fMakeAnaTree) {
    art::ServiceHandle<art::TFileService const> tfs;
    fNtEdepAna = tfs->make<TNtuple>(
      "nt_edep_ana",
      "Edep PosDiff Ana Ntuple",
      "energy:orig_x:orig_y:orig_z:orig_el:orig_ph:shift_x:shift_y:shift_z:shift_el:shift_ph");
  }
}

ShiftEdepSCE& spacecharge::ShiftEdepSCE::operator= ( ShiftEdepSCE const &  )

delete

ShiftEdepSCE& spacecharge::ShiftEdepSCE::operator= ( ShiftEdepSCE &&  )

delete

void spacecharge::ShiftEdepSCE::produce ( art::Event &  e )

override

Definition at line 82 of file ShiftEdepSCE_module.cc.

{
  auto sce = lar::providerFrom<spacecharge::SpaceChargeService>();
  auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(e);

  auto const& inEdepVec = *e.getValidHandle<std::vector<sim::SimEnergyDeposit>>(fEDepTag);

  auto outEdepVecPtr = std::make_unique<std::vector<sim::SimEnergyDeposit>>();
  auto& outEdepVec = *outEdepVecPtr;
  outEdepVec.reserve(inEdepVec.size());

  geo::Vector_t posOffsetsStart{0.0, 0.0, 0.0};
  geo::Vector_t posOffsetsEnd{0.0, 0.0, 0.0};
  for (auto const& edep : inEdepVec) {
    if (sce->EnableSimSpatialSCE()) {
      posOffsetsStart = sce->GetPosOffsets({edep.StartX(), edep.StartY(), edep.StartZ()});
      posOffsetsEnd = sce->GetPosOffsets({edep.EndX(), edep.EndY(), edep.EndZ()});
      if (larsim::Utils::SCE::out_of_bounds(posOffsetsStart) ||
          larsim::Utils::SCE::out_of_bounds(posOffsetsEnd) ) { 
          continue; 
        }
    }
    auto const isData = fISAlg.CalcIonAndScint(detProp, edep);
    outEdepVec.emplace_back(
      isData.numPhotons,
      isData.numElectrons,
      0.0,
      edep.Energy(),
      geo::Point_t{(float)(edep.StartX() - posOffsetsStart.X()), //x should be subtracted
                   (float)(edep.StartY() + posOffsetsStart.Y()),
                   (float)(edep.StartZ() + posOffsetsStart.Z())},
      geo::Point_t{(float)(edep.EndX() - posOffsetsEnd.X()), //x should be subtracted
                   (float)(edep.EndY() + posOffsetsEnd.Y()),
                   (float)(edep.EndZ() + posOffsetsEnd.Z())},
      edep.StartT(),
      edep.EndT(),
      edep.TrackID(),
      edep.PdgCode());
    if (fMakeAnaTree)
      fNtEdepAna->Fill(edep.Energy(),
                       edep.X(),
                       edep.Y(),
                       edep.Z(),
                       edep.NumElectrons(),
                       edep.NumPhotons(),
                       outEdepVec.back().X(),
                       outEdepVec.back().Y(),
                       outEdepVec.back().Z(),
                       outEdepVec.back().NumElectrons(),
                       outEdepVec.back().NumPhotons());
  }

  e.put(std::move(outEdepVecPtr));
}

Member Data Documentation

art::InputTag spacecharge::ShiftEdepSCE::fEDepTag

private

Definition at line 53 of file ShiftEdepSCE_module.cc.

larg4::ISCalcSeparate spacecharge::ShiftEdepSCE::fISAlg

private

Definition at line 58 of file ShiftEdepSCE_module.cc.

bool spacecharge::ShiftEdepSCE::fMakeAnaTree

private

Definition at line 54 of file ShiftEdepSCE_module.cc.

TNtuple* spacecharge::ShiftEdepSCE::fNtEdepAna

private

Definition at line 55 of file ShiftEdepSCE_module.cc.

---

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

• ShiftEdepSCE_module.cc