EIC issueshttps://eicweb.phy.anl.gov/groups/EIC/-/issues2022-08-09T15:13:54Zhttps://eicweb.phy.anl.gov/EIC/juggler/-/issues/102Job Failed #727010: acts@:19.3.0 fails to compile TrackProjector2022-08-09T15:13:54ZWouter DeconinckJob Failed #727010: acts@:19.3.0 fails to compile TrackProjectorJob [#727010](https://eicweb.phy.anl.gov/EIC/juggler/-/jobs/727010) failed for af1ac54e8f170dc3d4074631089a58ebc5e99af8:
Upstream: https://github.com/acts-project/acts/pull/1289Job [#727010](https://eicweb.phy.anl.gov/EIC/juggler/-/jobs/727010) failed for af1ac54e8f170dc3d4074631089a58ebc5e99af8:
Upstream: https://github.com/acts-project/acts/pull/1289Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/juggler/-/issues/101`m_BField->getField(pos,bCache).value()` compile failure in eigen with gcc@122022-06-25T00:10:05ZWouter Deconinck`m_BField->getField(pos,bCache).value()` compile failure in eigen with gcc@12```
291 inlined from 'void Eigen::PlainObjectBase<Derived>::_init1(const Eigen::DenseBase<ElseDerived>&) [with T = Eigen::Transpose<Eigen::CwiseBinaryOp<Eigen::internal::scalar_quotient_op<double, double>, const Eigen::Matrix...```
291 inlined from 'void Eigen::PlainObjectBase<Derived>::_init1(const Eigen::DenseBase<ElseDerived>&) [with T = Eigen::Transpose<Eigen::CwiseBinaryOp<Eigen::internal::scalar_quotient_op<double, double>, const Eigen::Matrix<d
ouble, 3, 1>, const Eigen::CwiseNullaryOp<Eigen::internal::scalar_constant_op<double>, const Eigen::Matrix<double, 3, 1> > > >; OtherDerived = Eigen::Transpose<Eigen::CwiseBinaryOp<Eigen::internal::scalar_quotient_op<double
, double>, const Eigen::Matrix<double, 3, 1>, const Eigen::CwiseNullaryOp<Eigen::internal::scalar_constant_op<double>, const Eigen::Matrix<double, 3, 1> > > >; Derived = Eigen::Matrix<double, 1, 3>]' at /opt/software/linux-
ubuntu22.04-skylake/gcc-12.1.0/eigen-3.4.0-xmmbrla23rhtxamudlr2l7362yihflsn/include/eigen3/Eigen/src/Core/PlainObjectBase.h:883:25,
292 inlined from 'Eigen::Matrix<_Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols>::Matrix(const T&) [with T = Eigen::Transpose<Eigen::CwiseBinaryOp<Eigen::internal::scalar_quotient_op<double, double>, const Eigen::Matrix
<double, 3, 1>, const Eigen::CwiseNullaryOp<Eigen::internal::scalar_constant_op<double>, const Eigen::Matrix<double, 3, 1> > > >; _Scalar = double; int _Rows = 1; int _Cols = 3; int _Options = 1; int _MaxRows = 1; int _MaxC
ols = 3]' at /opt/software/linux-ubuntu22.04-skylake/gcc-12.1.0/eigen-3.4.0-xmmbrla23rhtxamudlr2l7362yihflsn/include/eigen3/Eigen/src/Core/Matrix.h:332:31,
293 inlined from 'Eigen::DenseBase<Derived>::EvalReturnType Eigen::DenseBase<Derived>::eval() const [with Derived = Eigen::Transpose<Eigen::CwiseBinaryOp<Eigen::internal::scalar_quotient_op<double, double>, const Eigen::Mat
rix<double, 3, 1>, const Eigen::CwiseNullaryOp<Eigen::internal::scalar_constant_op<double>, const Eigen::Matrix<double, 3, 1> > > >]' at /opt/software/linux-ubuntu22.04-skylake/gcc-12.1.0/eigen-3.4.0-xmmbrla23rhtxamudlr2l73
62yihflsn/include/eigen3/Eigen/src/Core/DenseBase.h:412:62,
294 inlined from 'std::ostream& Eigen::operator<<(std::ostream&, const DenseBase<Derived>&) [with Derived = Transpose<CwiseBinaryOp<internal::scalar_quotient_op<double, double>, const Matrix<double, 3, 1>, const CwiseNullar
yOp<internal::scalar_constant_op<double>, const Matrix<double, 3, 1> > > >]' at /opt/software/linux-ubuntu22.04-skylake/gcc-12.1.0/eigen-3.4.0-xmmbrla23rhtxamudlr2l7362yihflsn/include/eigen3/Eigen/src/Core/IO.h:253:32,
295 inlined from 'MsgStream& operator<<(MsgStream&, const T&) [with T = Eigen::Transpose<Eigen::CwiseBinaryOp<Eigen::internal::scalar_quotient_op<double, double>, const Eigen::Matrix<double, 3, 1>, const Eigen::CwiseNullary
Op<Eigen::internal::scalar_constant_op<double>, const Eigen::Matrix<double, 3, 1> > > >]' at /opt/software/linux-ubuntu22.04-skylake/gcc-12.1.0/gaudi-36.5-scgkwadreca7hlk34q57xs7ccyotugf7/include/GaudiKernel/MsgStream.h:304
:20,
296 inlined from 'virtual StatusCode Jug::Reco::CKFTracking::initialize()' at /home/wdconinc/.spack/stage/spack-stage-juggler-master-szupytavynub76mjf5bggc7fe77kczk5/spack-src/JugTrack/src/components/CKFTracking.cpp:92:64:
>> 297 /usr/lib/gcc/x86_64-linux-gnu/12/include/emmintrin.h:131:24: error: '*(__m128d_u*)((char*)&<unnamed> + offsetof(Eigen::Matrix<double, 3, 1, 0, 3, 1>,Eigen::Matrix<double, 3, 1, 0, 3, 1>::<unnamed>.Eigen::PlainObjectBase<Eig
en::Matrix<double, 3, 1, 0, 3, 1> >::<unnamed>.Eigen::MatrixBase<Eigen::Matrix<double, 3, 1, 0, 3, 1> >::<unnamed>.Eigen::DenseBase<Eigen::Matrix<double, 3, 1, 0, 3, 1> >::<unnamed>.Eigen::DenseCoeffsBase<Eigen::Matrix<doub
le, 3, 1, 0, 3, 1>, 3>::<unnamed>.Eigen::DenseCoeffsBase<Eigen::Matrix<double, 3, 1, 0, 3, 1>, 1>::<unnamed>.Eigen::DenseCoeffsBase<Eigen::Matrix<double, 3, 1, 0, 3, 1>, 0>::<unnamed>))' may be used uninitialized [-Werror=m
aybe-uninitialized]
298 131 | return *(__m128d_u *)__P;
299 | ^~~
300 /home/wdconinc/.spack/stage/spack-stage-juggler-master-szupytavynub76mjf5bggc7fe77kczk5/spack-src/JugTrack/src/components/CKFTracking.cpp: In member function 'virtual StatusCode Jug::Reco::CKFTracking::initialize()':
301 /home/wdconinc/.spack/stage/spack-stage-juggler-master-szupytavynub76mjf5bggc7fe77kczk5/spack-src/JugTrack/src/components/CKFTracking.cpp:91:68: note: '<anonymous>' declared here
302 91 | auto b = m_BField->getField({0.0,0.0,double(z)},bCache).value()/(Acts::UnitConstants::T);
303 | ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~
```Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/juggler/-/issues/100Build all algorithms to avoid stale code2022-06-25T04:12:34ZWouter DeconinckBuild all algorithms to avoid stale codeCurrently we use explicit lists of files in [gaudi_add_module](https://eicweb.phy.anl.gov/EIC/juggler/-/blob/master/JugTrack/CMakeLists.txt#L9) calls in the CMakeLists.txt. While GLOB is not recommended (since added files do not cause a ...Currently we use explicit lists of files in [gaudi_add_module](https://eicweb.phy.anl.gov/EIC/juggler/-/blob/master/JugTrack/CMakeLists.txt#L9) calls in the CMakeLists.txt. While GLOB is not recommended (since added files do not cause a change to CMakeLists.txt), we could enable this with `CONFIGURE_DEPENDS` or make it a two step procedure.Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/detectors/ecce/-/issues/18dRICH: increase z-length, and add space for services2022-07-20T09:18:10ZChristopher DilksdRICH: increase z-length, and add space for servicesnew dimensions and z-position from <https://indico.bnl.gov/event/16314/>
![2022-07-19-203317_783x439_scrot](/uploads/8cb4dcd2541e7f072548d3d6aa3717cf/2022-07-19-203317_783x439_scrot.png)new dimensions and z-position from <https://indico.bnl.gov/event/16314/>
![2022-07-19-203317_783x439_scrot](/uploads/8cb4dcd2541e7f072548d3d6aa3717cf/2022-07-19-203317_783x439_scrot.png)Christopher DilksChristopher Dilkshttps://eicweb.phy.anl.gov/EIC/benchmarks/detector_benchmarks/-/issues/64Check ZDC particle gun acceptance2022-07-05T19:45:16ZAlex JentschCheck ZDC particle gun acceptanceAlex JentschAlex Jentschhttps://eicweb.phy.anl.gov/EIC/juggler/-/issues/98Enable -Wsign-conversion to avoid implicit conversions2022-06-22T01:24:09ZWouter DeconinckEnable -Wsign-conversion to avoid implicit conversionsYet another signed to unsigned implicit conversion caused under threshold reconstructed hit energies to overflow... !423 !415 #90Yet another signed to unsigned implicit conversion caused under threshold reconstructed hit energies to overflow... !423 !415 #90Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/irt/-/issues/7refactor `EVALUATION` for the `edm4hep` update2022-08-04T17:04:10ZChristopher Dilksrefactor `EVALUATION` for the `edm4hep` updateSince the data model has changed, we need to refactor some of the `evaluation` code. Things to do include:
- update our usage of `eic::Index`, since we now use PODIO relations instead
- should we use `edm4hep::MCParticle` to replace `dd4...Since the data model has changed, we need to refactor some of the `evaluation` code. Things to do include:
- update our usage of `eic::Index`, since we now use PODIO relations instead
- should we use `edm4hep::MCParticle` to replace `dd4pod::Geant4Particle`?
- how to match `ReconstructedParticle` to the MC generated particle? we no longer have `mcID`, do we need to add a relation in `eicd`?https://eicweb.phy.anl.gov/EIC/detectors/ecce/-/issues/17dRICH: produce IRT geometry objects2022-07-21T04:18:58ZChristopher DilksdRICH: produce IRT geometry objectsWe would like to connect the IRT geometry objects to the geometry description. We have tried doing this at the `Juggler` level, using `GeoSvc`, but find it more practical to keep the IRT geometry code integrated directly *with* the detec...We would like to connect the IRT geometry objects to the geometry description. We have tried doing this at the `Juggler` level, using `GeoSvc`, but find it more practical to keep the IRT geometry code integrated directly *with* the detector geometry code. We could automate the production of a "configuration" file containing IRT geometry objects, generated by `src/DRICH_geo.cpp` as needed.
This is how we ran for ATHENA, using https://eicweb.phy.anl.gov/EIC/detectors/athena/-/merge_requests/331 but we will need to think about how to best handle the extra configuration file in CI pipelines.Christopher DilksChristopher Dilkshttps://eicweb.phy.anl.gov/EIC/benchmarks/common_bench/-/issues/6Add JUGGLER_DETECTOR_CONFIG2022-06-15T00:40:25ZWouter DeconinckAdd JUGGLER_DETECTOR_CONFIGWe need a way to specify the geometry entry point for benchmarks.We need a way to specify the geometry entry point for benchmarks.Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/benchmarks/common_bench/-/issues/5DETECTOR_PATH should point to CMAKE_INSTALL_PREFIX subtree, not source dir2023-10-02T00:31:22ZWouter DeconinckDETECTOR_PATH should point to CMAKE_INSTALL_PREFIX subtree, not source dirRight now DETECTOR_PATH is set (in `env.sh`) to ${LOCAL_PREFIX}/detector which is where the detector geometry repository is cloned into.
```console
## detector prefix: prefix for the detector definitions
export DETECTOR_PREFIX="${LOCAL_P...Right now DETECTOR_PATH is set (in `env.sh`) to ${LOCAL_PREFIX}/detector which is where the detector geometry repository is cloned into.
```console
## detector prefix: prefix for the detector definitions
export DETECTOR_PREFIX="${LOCAL_PREFIX}/detector"
mkdir -p ${DETECTOR_PREFIX}
## detector path: actual detector definition path
export DETECTOR_PATH="${DETECTOR_PREFIX}/${JUGGLER_DETECTOR}"
```
That means any `make install` doesn't actually get used.Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/benchmarks/physics_benchmarks/-/issues/47JugFast::MatchClusters needs updating for EDM4hep migration2022-06-10T18:30:38ZSylvester JoostenJugFast::MatchClusters needs updating for EDM4hep migrationhttps://eicweb.phy.anl.gov/EIC/detectors/ecce/-/issues/14Add options for envelopes where we use assemblies2022-06-02T23:14:07ZWouter DeconinckAdd options for envelopes where we use assembliesBetter performance, better graphics rendering.
E.g.
- [x] HomogeneousCalorimeter_geo: add_disk
- [x] HomogeneousCalorimeter_geo: add_array
- [ ] ~~HomogeneousCalorimeter_geo: add_lines~~Better performance, better graphics rendering.
E.g.
- [x] HomogeneousCalorimeter_geo: add_disk
- [x] HomogeneousCalorimeter_geo: add_array
- [ ] ~~HomogeneousCalorimeter_geo: add_lines~~Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/detectors/ecce/-/issues/13Certificate issue when downloading field map: CAfile:/etc/pki/tls/certs/ca-bu...2022-07-21T15:44:47ZWouter DeconinckCertificate issue when downloading field map: CAfile:/etc/pki/tls/certs/ca-bundle.crt```
Compact INFO ++ Converted subdetector:SolenoidCoilBarrel of type athena_SolenoidCoil
Compact INFO ++ Converted subdetector:SolenoidCoilEndcapP of type ref_SolenoidEndcap
Compact INFO ++ Converted subdete...```
Compact INFO ++ Converted subdetector:SolenoidCoilBarrel of type athena_SolenoidCoil
Compact INFO ++ Converted subdetector:SolenoidCoilEndcapP of type ref_SolenoidEndcap
Compact INFO ++ Converted subdetector:SolenoidCoilEndcapN of type ref_SolenoidEndcap
FileLoader INFO Cache
FileLoader INFO Downloading fieldmaps/EIC_v.2.0.3_Magnetic_Field_Map_2021_09_28_radial_coords_cm_T.Bmap.txt as hash 6567da5a4887ad9b with curl --retry 5 -f https://eicweb.phy.anl.gov/EIC/detectors/athena/uploads/fd64db33e1c99547ab1840d8847fd1c1/EIC_v.2.0.3_Magnetic_Field_Map_2021_09_28_radial_coords_cm_T.Bmap.txt -o fieldmaps/6567da5a4887ad9b
% Total % Received % Xferd Average Speed Time Time Time Current
Dload Upload Total Spent Left Speed
0 0 0 0 0 0 0 0 --:--:-- --:--:-- 0 0 0 0 0 0 0 0 --:--:-- --:--:-- 0 0 0 0 0 0 0 0 --:--:-- --:--:----:--:-- 0
curl: (77) error setting certificate verify locations: CAfile:/etc/pki/tls/certs/ca-bundle.crt CApath: none
FileLoader ERROR unable to run cmd curl --retry 5 -f https://eicweb.phy.anl.gov/EIC/detectors/athena/uploads/fd64db33e1c99547ab1840d8847fd1c1/EIC_v.2.0.3_Magnetic_Field_Map_2021_09_28_radial_coords_cm_T.Bmap.txt -o fieldmaps/6567da5a4887ad9b
FileLoader ERROR check command and retry
nightly> akshaya@narval1:~/projects/def-wdconinc/akshaya/EIC$
```https://eicweb.phy.anl.gov/EIC/detectors/ip6/-/issues/55Far Backward Arb8 error2022-06-30T22:06:55ZWouter DeconinckFar Backward Arb8 errorThis seems like a degenerate arb8.
```
#25 136.2 Error in <TGeoArb8::ComputeTwist>: Shape TGeoArb8 type Arb8: Vertices must be defined clockwise in XY planes. Re-ordering...
#25 136.2 Error in <TGeoArb8::ComputeTwist>: Shape TGeoArb8 typ...This seems like a degenerate arb8.
```
#25 136.2 Error in <TGeoArb8::ComputeTwist>: Shape TGeoArb8 type Arb8: Vertices must be defined clockwise in XY planes. Re-ordering...
#25 136.2 Error in <TGeoArb8::ComputeTwist>: Shape TGeoArb8 type Arb8: Malformed polygon with crossing opposite segments
#25 144.5 *** Shape TGeoArb8: TGeoArb8 ***
#25 144.5 point #0 : x= -66.27951 y=-2399.27697 z= -20.00000
#25 144.5 point #1 : x= -9.80000 y=-1486.50000 z= -20.00000
#25 144.5 point #2 : x= -9.80000 y=-1486.50000 z= -20.00000
#25 144.5 point #3 : x= -26.28605 y=-2400.00000 z= -20.00000
#25 144.5 point #4 : x= -66.27951 y=-2399.27697 z= 20.00000
#25 144.5 point #5 : x= -9.80000 y=-1486.50000 z= 20.00000
#25 144.5 point #6 : x= -9.80000 y=-1486.50000 z= 20.00000
#25 144.5 point #7 : x= -26.28605 y=-2400.00000 z= 20.00000
#25 144.5 Bounding box:
#25 144.5 *** Shape TGeoArb8: TGeoBBox ***
#25 144.5 dX = 0.00000
#25 144.5 dY = 0.00000
#25 144.5 dZ = 0.00000
#25 144.5 origin: x= 0.00000 y= 0.00000 z= 0.00000
```Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/detectors/ecce/-/issues/12dRICH sensor material should not be AirOptical2022-07-27T04:12:53ZChristopher DilksdRICH sensor material should not be AirOpticalThis was temporary and needs to be updated
- moreover, I had to add reflectivity properties to the dRICH sensor surface from the ATHENA `master` branch, in order to see hits; not sure if this is correct (since we did not have this in the...This was temporary and needs to be updated
- moreover, I had to add reflectivity properties to the dRICH sensor surface from the ATHENA `master` branch, in order to see hits; not sure if this is correct (since we did not have this in the ATHENA proposal production)
## tasks
- [ ] check the photon count, to make sure we get the expected numbers
- [ ] check the sensor geometry and segmentation
- [ ] surface paramaterization (esp. reflectivity)
- [ ] sensor bases (how to model the back of the sensors?) - make sure there is enough space for this
## notes from ATHENA studies (from discussions with Alexander, Chandra, and Chris)
- sensor used for proposal: https://www.hamamatsu.com/us/en/product/optical-sensors/mppc/mppc_mppc-array/S13361-3050AE-08.html
- the non-negligible gas-to-resin (numerically
<5%) and resin-to-silicon (~15%) reflections must be effectively
accounted in the quoted PDEs. At all incident angles of interest
for us (up to 30 degree or so) there numbers get shared between
the polarization states, but average stays more or less the same.
- However in our case (g4dRIChOptics.hh) the sensor surface was defined as dielectric_metal with
some bogus imaginary refractive index parameters. This killed the
photons, and also I verified that the MC .root files contain only
the ones which were detected somewhere. So I temporarily changed the material to AirOptical, and added
a benign surface to the optical_metarials.xml database, and we are
back at ~10 npe. And "Chandra's number" for 350..650nm integral is
now ~72 or so.
- I think the correct way to account Cherenkov photon polarization is
to create a resin volume, and perhaps even a silicon volume inside it,
but renormalize the PDE, accounting for the normal incident losses.
- look into `G4SiPM` packageChandradoy Chatterjeechandradoy.chatterjee@ts.infn.itChandradoy Chatterjeechandradoy.chatterjee@ts.infn.ithttps://eicweb.phy.anl.gov/EIC/benchmarks/physics_benchmarks/-/issues/46Correspondence between ReconstructedParticle and MCParticles2022-07-31T00:10:47ZZhoudunming TuCorrespondence between ReconstructedParticle and MCParticlesWe would like to have a way to access a specific reconstructed particle that has matched to the MCparticle, e.g., the scattered electron, and so on. This would allow studies for smearing/resolution effect, efficiency and fake rate for tr...We would like to have a way to access a specific reconstructed particle that has matched to the MCparticle, e.g., the scattered electron, and so on. This would allow studies for smearing/resolution effect, efficiency and fake rate for tracking, user-defined quantities at reconstructed level, etc.Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/detectors/ip6/-/issues/54Job Failed #662974: increase timeout in generate_prim2022-05-26T05:06:56ZWouter DeconinckJob Failed #662974: increase timeout in generate_primJob [#662974](https://eicweb.phy.anl.gov/EIC/detectors/ip6/-/jobs/662974) failed for 8cb6388d3805dee0418c0c4f67392325ba3fe738:Job [#662974](https://eicweb.phy.anl.gov/EIC/detectors/ip6/-/jobs/662974) failed for 8cb6388d3805dee0418c0c4f67392325ba3fe738:Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/detectors/athena/-/issues/152Fix issues introduced by ip6 changes2022-05-13T19:17:51ZWouter DeconinckFix issues introduced by ip6 changesWouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/detectors/ecce/-/issues/10Re-scale the dRICH2022-05-30T17:30:37ZChristopher DilksRe-scale the dRICHRe-scale the dRICH envelope. For the first attempt, we can use the ECCE dRICH dimensions.Re-scale the dRICH envelope. For the first attempt, we can use the ECCE dRICH dimensions.Christopher DilksChristopher Dilkshttps://eicweb.phy.anl.gov/EIC/detectors/ecce/-/issues/8Use correct BaBar fieldmap2022-06-10T23:04:39ZWouter DeconinckUse correct BaBar fieldmap