forked from celeritas-project/celeritas
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Add optical generator data and update pre-generators (celeritas-proje…
…ct#1182) * Add optical generator data and update pre-generators * Address @stognini's feedback * Address @sethrj feedback
- Loading branch information
1 parent
fda9820
commit bb6fc2a
Showing
8 changed files
with
394 additions
and
254 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,157 @@ | ||
| //----------------------------------*-C++-*----------------------------------// | ||
| // Copyright 2024 UT-Battelle, LLC, and other Celeritas developers. | ||
| // See the top-level COPYRIGHT file for details. | ||
| // SPDX-License-Identifier: (Apache-2.0 OR MIT) | ||
| //---------------------------------------------------------------------------// | ||
| //! \file celeritas/optical/OpticalGenData.hh | ||
| //---------------------------------------------------------------------------// | ||
| #pragma once | ||
|
|
||
| #include "corecel/Macros.hh" | ||
| #include "corecel/Types.hh" | ||
| #include "corecel/data/Collection.hh" | ||
| #include "celeritas/Quantities.hh" | ||
| #include "celeritas/Types.hh" | ||
| #include "celeritas/optical/CerenkovData.hh" | ||
| #include "celeritas/optical/OpticalDistributionData.hh" | ||
| #include "celeritas/optical/OpticalPropertyData.hh" | ||
| #include "celeritas/optical/ScintillationData.hh" | ||
|
|
||
| namespace celeritas | ||
| { | ||
| //---------------------------------------------------------------------------// | ||
| /*! | ||
| * Current sizes of the buffers of distribution data. | ||
| * | ||
| * These sizes are updated by value on the host at each step. | ||
| */ | ||
| struct OpticalBufferSize | ||
| { | ||
| size_type cerenkov{0}; | ||
| size_type scintillation{0}; | ||
| }; | ||
|
|
||
| //---------------------------------------------------------------------------// | ||
| /*! | ||
| * Immutable problem data for generating optical photon distributions. | ||
| */ | ||
| template<Ownership W, MemSpace M> | ||
| struct OpticalGenParamsData | ||
| { | ||
| //// DATA //// | ||
|
|
||
| bool cerenkov{false}; //!< Whether Cerenkov is enabled | ||
| bool scintillation{false}; //!< Whether scintillation is enabled | ||
| real_type capacity{0}; //!< Distribution data buffer capacity | ||
|
|
||
| //// METHODS //// | ||
|
|
||
| //! True if all params are assigned | ||
| explicit CELER_FUNCTION operator bool() const | ||
| { | ||
| return (cerenkov || scintillation) && capacity > 0; | ||
| } | ||
|
|
||
| //! Assign from another set of data | ||
| template<Ownership W2, MemSpace M2> | ||
| OpticalGenParamsData& operator=(OpticalGenParamsData<W2, M2> const& other) | ||
| { | ||
| CELER_EXPECT(other); | ||
| cerenkov = other.cerenkov; | ||
| scintillation = other.scintillation; | ||
| capacity = other.capacity; | ||
| return *this; | ||
| } | ||
| }; | ||
|
|
||
| //---------------------------------------------------------------------------// | ||
| /*! | ||
| * Pre-step data needed to generate optical photon distributions. | ||
| */ | ||
| struct OpticalPreStepData | ||
| { | ||
| units::LightSpeed speed; | ||
| Real3 pos{}; | ||
| real_type time{}; | ||
|
|
||
| //! Check whether the data are assigned | ||
| explicit CELER_FUNCTION operator bool() const | ||
| { | ||
| return speed > zero_quantity(); | ||
| } | ||
| }; | ||
|
|
||
| //---------------------------------------------------------------------------// | ||
| /*! | ||
| * Optical photon distribution data. | ||
| */ | ||
| template<Ownership W, MemSpace M> | ||
| struct OpticalGenStateData | ||
| { | ||
| //// TYPES //// | ||
|
|
||
| template<class T> | ||
| using StateItems = StateCollection<T, W, M>; | ||
| template<class T> | ||
| using Items = Collection<T, W, M>; | ||
|
|
||
| //// DATA //// | ||
|
|
||
| // Pre-step data for generating optical photon distributions | ||
| StateItems<OpticalPreStepData> step; | ||
|
|
||
| // Buffers of distribution data for generating optical primaries | ||
| Items<OpticalDistributionData> cerenkov; | ||
| Items<OpticalDistributionData> scintillation; | ||
|
|
||
| //// METHODS //// | ||
|
|
||
| //! Number of states | ||
| CELER_FUNCTION size_type size() const { return step.size(); } | ||
|
|
||
| //! Whether all data are assigned and valid | ||
| explicit CELER_FUNCTION operator bool() const | ||
| { | ||
| return !step.empty() && !(cerenkov.empty() && scintillation.empty()); | ||
| } | ||
|
|
||
| //! Assign from another set of data | ||
| template<Ownership W2, MemSpace M2> | ||
| OpticalGenStateData& operator=(OpticalGenStateData<W2, M2>& other) | ||
| { | ||
| CELER_EXPECT(other); | ||
| step = other.step; | ||
| cerenkov = other.cerenkov; | ||
| scintillation = other.scintillation; | ||
| return *this; | ||
| } | ||
| }; | ||
|
|
||
| //---------------------------------------------------------------------------// | ||
| /*! | ||
| * Resize optical states. | ||
| */ | ||
| template<MemSpace M> | ||
| void resize(OpticalGenStateData<Ownership::value, M>* state, | ||
| HostCRef<OpticalGenParamsData> const& params, | ||
| StreamId, | ||
| size_type size) | ||
| { | ||
| CELER_EXPECT(params); | ||
| CELER_EXPECT(size > 0); | ||
|
|
||
| resize(&state->step, size); | ||
| if (params.cerenkov) | ||
| { | ||
| resize(&state->cerenkov, params.capacity); | ||
| } | ||
| if (params.scintillation) | ||
| { | ||
| resize(&state->scintillation, params.capacity); | ||
| } | ||
|
|
||
| CELER_ENSURE(*state); | ||
| } | ||
|
|
||
| //---------------------------------------------------------------------------// | ||
| } // namespace celeritas |
Oops, something went wrong.