EIC issueshttps://eicweb.phy.anl.gov/groups/EIC/-/issues2022-06-30T22:06:55Zhttps://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/benchmarks/reconstruction_benchmarks/-/issues/88TrackProjector benchmark2022-06-29T03:32:40ZWouter DeconinckTrackProjector benchmarkEnable the track projector benchmarkEnable the track projector benchmarkWouter 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/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/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/juggler/-/issues/90Unsigned eicd::RawCalorimeterHit amplitude causes threshold issues2022-06-21T22:43:35ZWouter DeconinckUnsigned eicd::RawCalorimeterHit amplitude causes threshold issues`(rh.getAmplitude() - m_pedMeanADC) < thresholdADC` for negative `rh.getAmplitude() - m_pedMeanADC` will results in wraparound since `unsigned int - int` is an `unsigned int`.`(rh.getAmplitude() - m_pedMeanADC) < thresholdADC` for negative `rh.getAmplitude() - m_pedMeanADC` will results in wraparound since `unsigned int - int` is an `unsigned int`.Wouter DeconinckWouter Deconinckhttps://eicweb.phy.anl.gov/EIC/benchmarks/reconstruction_benchmarks/-/issues/46Track finding efficiency2022-06-16T15:55:15ZWhitney ArmstrongTrack finding efficiencyAdd track finding efficiency analysis.
- [ ] Add basic script for CKF finding.
- [ ] Add new algorithm (Hough transform)
- [ ] Create comparison between finders.Add track finding efficiency analysis.
- [ ] Add basic script for CKF finding.
- [ ] Add new algorithm (Hough transform)
- [ ] Create comparison between finders.Complete Tracking Performance MetricsWhitney ArmstrongWhitney Armstrong2021-08-11https://eicweb.phy.anl.gov/EIC/juggler/-/issues/5Added second tracking algorithm.2022-06-16T15:26:10ZWhitney ArmstrongAdded second tracking algorithm.Whitney ArmstrongWhitney Armstronghttps://eicweb.phy.anl.gov/EIC/benchmarks/reconstruction_benchmarks/-/issues/33Clean up final artifacts/report2022-06-16T15:25:59ZWhitney ArmstrongClean up final artifacts/reportNeed some better organization for the final artifacts.
https://eicweb.phy.anl.gov/EIC/benchmarks/reconstruction_benchmarks/-/jobs/123674/artifacts/browse/results/Need some better organization for the final artifacts.
https://eicweb.phy.anl.gov/EIC/benchmarks/reconstruction_benchmarks/-/jobs/123674/artifacts/browse/results/Whitney ArmstrongWhitney Armstronghttps://eicweb.phy.anl.gov/EIC/benchmarks/reconstruction_benchmarks/-/issues/55Hough Transform Track finding2022-06-16T15:25:54ZWhitney ArmstrongHough Transform Track findingWith full search of $`\theta`$ vs $`\rho`$ 2D histogram.With full search of $`\theta`$ vs $`\rho`$ 2D histogram.Complete Tracking Performance MetricsWhitney ArmstrongWhitney Armstronghttps://eicweb.phy.anl.gov/EIC/juggler/-/issues/39Basic Vertexing algorithm2022-06-16T15:20:55ZWhitney ArmstrongBasic Vertexing algorithmAdd a basic vertexing example.Add a basic vertexing example.Whitney ArmstrongWhitney Armstronghttps://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/detectors/athena/-/issues/146DRICHHits (usually) empty after the change of `sens.setType("photoncounter")`...2022-06-12T06:15:38ZChristopher DilksDRICHHits (usually) empty after the change of `sens.setType("photoncounter")` to `tracker`After https://eicweb.phy.anl.gov/EIC/detectors/athena/-/merge_requests/330 was merged, the `DRICHHits` branch is usually empty (or there are only a few sporadic entries). If I try reverting the type back to `photoncounter`, we see the ex...After https://eicweb.phy.anl.gov/EIC/detectors/athena/-/merge_requests/330 was merged, the `DRICHHits` branch is usually empty (or there are only a few sporadic entries). If I try reverting the type back to `photoncounter`, we see the expected distribution.
Using MR https://eicweb.phy.anl.gov/EIC/detectors/athena/-/merge_requests/332 for testing: see artifact `test.png` from job `test_drich`
Steps to reproduce (executed from top-level `athena` directory):
```
npsim --runType run --compactFile compact/subsystem_views/drich_only.xml --macro test.mac --outputFile test.root --enableG4GPS
```
where `test.mac` is
```
/control/verbose 2
/run/initialize
/gps/verbose 2
/gps/particle pi+
/gps/number 1
/gps/ene/mono 20 GeV
/gps/position 0 0 0 cm
/gps/direction 0.25 0.0 1.0
/run/beamOn 50
```
Opening `test.root` and running `events->Draw("DRICHHits.position.x")` should produce a distribution that looks like the attached image, but is coming out empty (or sometimes a few sporadic hits)
It is likely the pfRICH is also affected.
![photoncounter](/uploads/b0d3bad0892f26ccb60ffe8e94b8dd27/photoncounter.png)https://eicweb.phy.anl.gov/EIC/detectors/ecce/-/issues/8Use correct BaBar fieldmap2022-06-10T23:04:39ZWouter DeconinckUse correct BaBar fieldmaphttps://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/athena/-/issues/142curl in FileLoader fails with certificate error2022-06-01T04:02:46ZWouter Deconinckcurl in FileLoader fails with certificate error```
Info in <TGeoManager::SetTopVolume>: Top volume is world_volume. Master volume is world_volume
% Total % Received % Xferd Average Speed Time Time Time Current
Dload Upload Total S...```
Info in <TGeoManager::SetTopVolume>: Top volume is world_volume. Master volume is world_volume
% Total % Received % Xferd Average Speed Time Time Time Current
Dload Upload Total Spent Left Speed
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
```Wouter 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/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 Deconinck