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
Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDProducer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept
	Producer (fhicl::ParameterSet const &)
	Producer (Producer const &)=delete
	Producer (Producer &&)=delete
Producer &	operator= (Producer const &)=delete
Producer &	operator= (Producer &&)=delete
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Modifier
	~Modifier () noexcept
	Modifier ()
	Modifier (Modifier const &)=delete
	Modifier (Modifier &&)=delete
Modifier &	operator= (Modifier const &)=delete
Modifier &	operator= (Modifier &&)=delete
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

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

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >
Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >
Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 33 of file ShiftEdepSCE_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 59 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 ( )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 68 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 )

overridevirtual

Implements art::EDProducer.

Definition at line 80 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 51 of file ShiftEdepSCE_module.cc.

larg4::ISCalcSeparate spacecharge::ShiftEdepSCE::fISAlg

private

Definition at line 56 of file ShiftEdepSCE_module.cc.

bool spacecharge::ShiftEdepSCE::fMakeAnaTree

private

Definition at line 52 of file ShiftEdepSCE_module.cc.

TNtuple* spacecharge::ShiftEdepSCE::fNtEdepAna

private

Definition at line 53 of file ShiftEdepSCE_module.cc.

---

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

• larsim/larsim/ElectronDrift/ShiftEdepSCE_module.cc