82 auto const* SCE = lar::providerFrom<spacecharge::SpaceChargeService>();
83 cout <<
"Is Spatial SCE enabled? " <<
bool(SCE->EnableSimSpatialSCE()) << endl;
84 cout <<
"Is E-field SCE enabled? " <<
bool(SCE->EnableSimEfieldSCE()) << endl;
87 for(
int iX = xMin; iX <= xMax - nSkip; iX++)
90 for(
int iY = yMin; iY <= yMax - nSkip; iY++)
93 for(
int iZ = zMin; iZ <= zMax - nSkip; iZ++)
96 cout << iX <<
", " << iY <<
", " << iZ << endl;
98 geo::Point_t point = {double(iX), double(iY), double(iZ)};
101 if(!((spatialOffsets.X() == spatialOffsets.X()) &&
102 (spatialOffsets.Y() == spatialOffsets.Y()) &&
103 (spatialOffsets.Z() == 0.0)))
105 hDx->Fill(spatialOffsets.X());
106 hDy->Fill(spatialOffsets.Y());
107 hDz->Fill(spatialOffsets.Z());
111 if(!((efieldOffsets.X() == efieldOffsets.X()) &&
112 (efieldOffsets.Y() == efieldOffsets.Y()) &&
113 (efieldOffsets.Z() == 0.0)))
115 hEx->Fill(efieldOffsets.X());
116 hEy->Fill(efieldOffsets.Y());
117 hEz->Fill(efieldOffsets.Z());
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Vector_t
Type for representation of momenta in 3D space.
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Point_t
Type for representation of position in physical 3D space.
BEGIN_PROLOG could also be cout