spacecharge::SpaceChargeStandard Class Reference

#include "SpaceChargeStandard.h"

Inheritance diagram for spacecharge::SpaceChargeStandard:

Public Member Functions
	SpaceChargeStandard (fhicl::ParameterSet const &pset)
	SpaceChargeStandard (SpaceChargeStandard const &)=delete
virtual	~SpaceChargeStandard ()=default
bool	Configure (fhicl::ParameterSet const &pset)
bool	Update (uint64_t ts=0)
bool	EnableSimSpatialSCE () const override
bool	EnableSimEfieldSCE () const override
bool	EnableCorrSCE () const override
	Return boolean indicating whether or not to apply SCE corrections. More...
bool	EnableCalSpatialSCE () const override
	Return boolean indicating whether or not to apply SCE corrections. More...
bool	EnableCalEfieldSCE () const override
	Return boolean indicating whether or not to apply SCE corrections. More...
geo::Vector_t	GetPosOffsets (geo::Point_t const &point) const override
geo::Vector_t	GetEfieldOffsets (geo::Point_t const &point) const override
geo::Vector_t	GetCalPosOffsets (geo::Point_t const &point, int const &TPCid) const override
geo::Vector_t	GetCalEfieldOffsets (geo::Point_t const &point, int const &TPCid) const override

Protected Member Functions
std::vector< double >	GetPosOffsetsParametric (double xVal, double yVal, double zVal) const
	Provides position offsets using a parametric representation. More...
double	GetOnePosOffsetParametric (double xVal, double yVal, double zVal, std::string axis) const
std::vector< double >	GetEfieldOffsetsParametric (double xVal, double yVal, double zVal) const
	Provides E field offsets using a parametric representation. More...
double	GetOneEfieldOffsetParametric (double xVal, double yVal, double zVal, std::string axis) const
double	TransformX (double xVal) const
	Transform X to SCE X coordinate - redefine this in experiment-specific implementation! More...
double	TransformY (double yVal) const
	Transform Y to SCE Y coordinate - redefine this in experiment-specific implementation! More...
double	TransformZ (double zVal) const
	Transform Z to SCE Z coordinate - redefine this in experiment-specific implementation! More...
bool	IsInsideBoundaries (double xVal, double yVal, double zVal) const
	Check to see if point is inside boundaries of map - redefine this in experiment-specific implementation! More...

Protected Attributes
bool	fEnableSimSpatialSCE
bool	fEnableSimEfieldSCE
bool	fEnableCalSpatialSCE
bool	fEnableCalEfieldSCE
bool	fEnableCorrSCE
std::string	fRepresentationType
std::string	fInputFilename
TGraph **	g1_x = new TGraph*[7]
TGraph **	g2_x = new TGraph*[7]
TGraph **	g3_x = new TGraph*[7]
TGraph **	g4_x = new TGraph*[7]
TGraph **	g5_x = new TGraph*[7]
TGraph **	g1_y = new TGraph*[7]
TGraph **	g2_y = new TGraph*[7]
TGraph **	g3_y = new TGraph*[7]
TGraph **	g4_y = new TGraph*[7]
TGraph **	g5_y = new TGraph*[7]
TGraph **	g6_y = new TGraph*[7]
TGraph **	g1_z = new TGraph*[7]
TGraph **	g2_z = new TGraph*[7]
TGraph **	g3_z = new TGraph*[7]
TGraph **	g4_z = new TGraph*[7]
TF1 *	f1_x = new TF1("f1_x", "pol6")
TF1 *	f2_x = new TF1("f2_x", "pol6")
TF1 *	f3_x = new TF1("f3_x", "pol6")
TF1 *	f4_x = new TF1("f4_x", "pol6")
TF1 *	f5_x = new TF1("f5_x", "pol6")
TF1 *	fFinal_x = new TF1("fFinal_x", "pol4")
TF1 *	f1_y = new TF1("f1_y", "pol5")
TF1 *	f2_y = new TF1("f2_y", "pol5")
TF1 *	f3_y = new TF1("f3_y", "pol5")
TF1 *	f4_y = new TF1("f4_y", "pol5")
TF1 *	f5_y = new TF1("f5_y", "pol5")
TF1 *	f6_y = new TF1("f6_y", "pol5")
TF1 *	fFinal_y = new TF1("fFinal_y", "pol5")
TF1 *	f1_z = new TF1("f1_z", "pol4")
TF1 *	f2_z = new TF1("f2_z", "pol4")
TF1 *	f3_z = new TF1("f3_z", "pol4")
TF1 *	f4_z = new TF1("f4_z", "pol4")
TF1 *	fFinal_z = new TF1("fFinal_z", "pol3")
TGraph **	g1_Ex = new TGraph*[7]
TGraph **	g2_Ex = new TGraph*[7]
TGraph **	g3_Ex = new TGraph*[7]
TGraph **	g4_Ex = new TGraph*[7]
TGraph **	g5_Ex = new TGraph*[7]
TGraph **	g1_Ey = new TGraph*[7]
TGraph **	g2_Ey = new TGraph*[7]
TGraph **	g3_Ey = new TGraph*[7]
TGraph **	g4_Ey = new TGraph*[7]
TGraph **	g5_Ey = new TGraph*[7]
TGraph **	g6_Ey = new TGraph*[7]
TGraph **	g1_Ez = new TGraph*[7]
TGraph **	g2_Ez = new TGraph*[7]
TGraph **	g3_Ez = new TGraph*[7]
TGraph **	g4_Ez = new TGraph*[7]
TF1 *	f1_Ex = new TF1("f1_Ex", "pol6")
TF1 *	f2_Ex = new TF1("f2_Ex", "pol6")
TF1 *	f3_Ex = new TF1("f3_Ex", "pol6")
TF1 *	f4_Ex = new TF1("f4_Ex", "pol6")
TF1 *	f5_Ex = new TF1("f5_Ex", "pol6")
TF1 *	fFinal_Ex = new TF1("fFinal_Ex", "pol4")
TF1 *	f1_Ey = new TF1("f1_Ey", "pol5")
TF1 *	f2_Ey = new TF1("f2_Ey", "pol5")
TF1 *	f3_Ey = new TF1("f3_Ey", "pol5")
TF1 *	f4_Ey = new TF1("f4_Ey", "pol5")
TF1 *	f5_Ey = new TF1("f5_Ey", "pol5")
TF1 *	f6_Ey = new TF1("f6_Ey", "pol5")
TF1 *	fFinal_Ey = new TF1("fFinal_Ey", "pol5")
TF1 *	f1_Ez = new TF1("f1_Ez", "pol4")
TF1 *	f2_Ez = new TF1("f2_Ez", "pol4")
TF1 *	f3_Ez = new TF1("f3_Ez", "pol4")
TF1 *	f4_Ez = new TF1("f4_Ez", "pol4")
TF1 *	fFinal_Ez = new TF1("fFinal_Ez", "pol3")

Detailed Description

Definition at line 32 of file SpaceChargeStandard.h.

Constructor & Destructor Documentation

spacecharge::SpaceChargeStandard::SpaceChargeStandard ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 26 of file SpaceChargeStandard.cxx.

References Configure().

{
  Configure(pset);
}

spacecharge::SpaceChargeStandard::SpaceChargeStandard ( SpaceChargeStandard const &  )

delete

virtual spacecharge::SpaceChargeStandard::~SpaceChargeStandard ( )

virtualdefault

Member Function Documentation

bool spacecharge::SpaceChargeStandard::Configure

(

fhicl::ParameterSet const &

pset

)

Definition at line 32 of file SpaceChargeStandard.cxx.

References art::errors::Configuration, fEnableCalEfieldSCE, fEnableCalSpatialSCE, fEnableCorrSCE, fEnableSimEfieldSCE, fEnableSimSpatialSCE, fInputFilename, fRepresentationType, g1_Ex, g1_Ey, g1_Ez, g1_x, g1_y, g1_z, g2_Ex, g2_Ey, g2_Ez, g2_x, g2_y, g2_z, g3_Ex, g3_Ey, g3_Ez, g3_x, g3_y, g3_z, g4_Ex, g4_Ey, g4_Ez, g4_x, g4_y, g4_z, g5_Ex, g5_Ey, g5_x, g5_y, g6_Ey, g6_y, fhicl::ParameterSet::get(), and fhicl::ParameterSet::has_key().

Referenced by spacecharge::SpaceChargeServiceStandard::reconfigure(), and SpaceChargeStandard().

{
  fEnableSimSpatialSCE = pset.get<bool>("EnableSimSpatialSCE");
  fEnableSimEfieldSCE = pset.get<bool>("EnableSimEfieldSCE");
  fEnableCalSpatialSCE = pset.get<bool>("EnableCalSpatialSCE");
  fEnableCalEfieldSCE = pset.get<bool>("EnableCalEfieldSCE");
  fEnableCorrSCE = pset.get<bool>("EnableCorrSCE");

  // check that the old obsoleted parameter is not in configuration:
  if (pset.has_key("EnableSimulationSCE")) {
    throw art::Exception(art::errors::Configuration)
      << "Configuration parameter 'EnableSimulationSCE' has been replaced by "
         "'EnableSimSpatialSCE'.\n";
  }

  if ((fEnableSimSpatialSCE == true) | (fEnableSimEfieldSCE == true)) {
    fRepresentationType = pset.get<std::string>("RepresentationType");
    fInputFilename = pset.get<std::string>("InputFilename");

    std::string fname;
    cet::search_path sp("FW_SEARCH_PATH");
    sp.find_file(fInputFilename, fname);

    auto infile = std::make_unique<TFile>(fname.c_str(), "READ");
    if (!infile->IsOpen())
      throw art::Exception(art::errors::Configuration)
        << "Could not find the space charge effect file '" << fname << "'!\n";

    if (fRepresentationType == "Parametric") {
      for (int i = 0; i < 5; i++) {
        g1_x[i] = (TGraph*)infile->Get(Form("deltaX/g1_%d", i));
        g2_x[i] = (TGraph*)infile->Get(Form("deltaX/g2_%d", i));
        g3_x[i] = (TGraph*)infile->Get(Form("deltaX/g3_%d", i));
        g4_x[i] = (TGraph*)infile->Get(Form("deltaX/g4_%d", i));
        g5_x[i] = (TGraph*)infile->Get(Form("deltaX/g5_%d", i));

        g1_y[i] = (TGraph*)infile->Get(Form("deltaY/g1_%d", i));
        g2_y[i] = (TGraph*)infile->Get(Form("deltaY/g2_%d", i));
        g3_y[i] = (TGraph*)infile->Get(Form("deltaY/g3_%d", i));
        g4_y[i] = (TGraph*)infile->Get(Form("deltaY/g4_%d", i));
        g5_y[i] = (TGraph*)infile->Get(Form("deltaY/g5_%d", i));
        g6_y[i] = (TGraph*)infile->Get(Form("deltaY/g6_%d", i));

        g1_z[i] = (TGraph*)infile->Get(Form("deltaZ/g1_%d", i));
        g2_z[i] = (TGraph*)infile->Get(Form("deltaZ/g2_%d", i));
        g3_z[i] = (TGraph*)infile->Get(Form("deltaZ/g3_%d", i));
        g4_z[i] = (TGraph*)infile->Get(Form("deltaZ/g4_%d", i));

        g1_Ex[i] = (TGraph*)infile->Get(Form("deltaExOverE/g1_%d", i));
        g2_Ex[i] = (TGraph*)infile->Get(Form("deltaExOverE/g2_%d", i));
        g3_Ex[i] = (TGraph*)infile->Get(Form("deltaExOverE/g3_%d", i));
        g4_Ex[i] = (TGraph*)infile->Get(Form("deltaExOverE/g4_%d", i));
        g5_Ex[i] = (TGraph*)infile->Get(Form("deltaExOverE/g5_%d", i));

        g1_Ey[i] = (TGraph*)infile->Get(Form("deltaEyOverE/g1_%d", i));
        g2_Ey[i] = (TGraph*)infile->Get(Form("deltaEyOverE/g2_%d", i));
        g3_Ey[i] = (TGraph*)infile->Get(Form("deltaEyOverE/g3_%d", i));
        g4_Ey[i] = (TGraph*)infile->Get(Form("deltaEyOverE/g4_%d", i));
        g5_Ey[i] = (TGraph*)infile->Get(Form("deltaEyOverE/g5_%d", i));
        g6_Ey[i] = (TGraph*)infile->Get(Form("deltaEyOverE/g6_%d", i));

        g1_Ez[i] = (TGraph*)infile->Get(Form("deltaEzOverE/g1_%d", i));
        g2_Ez[i] = (TGraph*)infile->Get(Form("deltaEzOverE/g2_%d", i));
        g3_Ez[i] = (TGraph*)infile->Get(Form("deltaEzOverE/g3_%d", i));
        g4_Ez[i] = (TGraph*)infile->Get(Form("deltaEzOverE/g4_%d", i));
      }

      g1_x[5] = (TGraph*)infile->Get("deltaX/g1_5");
      g2_x[5] = (TGraph*)infile->Get("deltaX/g2_5");
      g3_x[5] = (TGraph*)infile->Get("deltaX/g3_5");
      g4_x[5] = (TGraph*)infile->Get("deltaX/g4_5");
      g5_x[5] = (TGraph*)infile->Get("deltaX/g5_5");

      g1_y[5] = (TGraph*)infile->Get("deltaY/g1_5");
      g2_y[5] = (TGraph*)infile->Get("deltaY/g2_5");
      g3_y[5] = (TGraph*)infile->Get("deltaY/g3_5");
      g4_y[5] = (TGraph*)infile->Get("deltaY/g4_5");
      g5_y[5] = (TGraph*)infile->Get("deltaY/g5_5");
      g6_y[5] = (TGraph*)infile->Get("deltaY/g6_5");

      g1_x[6] = (TGraph*)infile->Get("deltaX/g1_6");
      g2_x[6] = (TGraph*)infile->Get("deltaX/g2_6");
      g3_x[6] = (TGraph*)infile->Get("deltaX/g3_6");
      g4_x[6] = (TGraph*)infile->Get("deltaX/g4_6");
      g5_x[6] = (TGraph*)infile->Get("deltaX/g5_6");

      g1_Ex[5] = (TGraph*)infile->Get("deltaExOverE/g1_5");
      g2_Ex[5] = (TGraph*)infile->Get("deltaExOverE/g2_5");
      g3_Ex[5] = (TGraph*)infile->Get("deltaExOverE/g3_5");
      g4_Ex[5] = (TGraph*)infile->Get("deltaExOverE/g4_5");
      g5_Ex[5] = (TGraph*)infile->Get("deltaExOverE/g5_5");

      g1_Ey[5] = (TGraph*)infile->Get("deltaEyOverE/g1_5");
      g2_Ey[5] = (TGraph*)infile->Get("deltaEyOverE/g2_5");
      g3_Ey[5] = (TGraph*)infile->Get("deltaEyOverE/g3_5");
      g4_Ey[5] = (TGraph*)infile->Get("deltaEyOverE/g4_5");
      g5_Ey[5] = (TGraph*)infile->Get("deltaEyOverE/g5_5");
      g6_Ey[5] = (TGraph*)infile->Get("deltaEyOverE/g6_5");

      g1_Ex[6] = (TGraph*)infile->Get("deltaExOverE/g1_6");
      g2_Ex[6] = (TGraph*)infile->Get("deltaExOverE/g2_6");
      g3_Ex[6] = (TGraph*)infile->Get("deltaExOverE/g3_6");
      g4_Ex[6] = (TGraph*)infile->Get("deltaExOverE/g4_6");
      g5_Ex[6] = (TGraph*)infile->Get("deltaExOverE/g5_6");
    }

    infile->Close();
  }

  if (fEnableCorrSCE == true) {
    // Grab other parameters from pset
  }

  return true;
}

bool spacecharge::SpaceChargeStandard::EnableCalEfieldSCE ( ) const

overridevirtual

Return boolean indicating whether or not to apply SCE corrections.

Implements spacecharge::SpaceCharge.

Definition at line 186 of file SpaceChargeStandard.cxx.

References fEnableCalEfieldSCE.

{
  return fEnableCalEfieldSCE;
}

bool spacecharge::SpaceChargeStandard::EnableCalSpatialSCE ( ) const

overridevirtual

Return boolean indicating whether or not to apply SCE corrections.

Implements spacecharge::SpaceCharge.

Definition at line 180 of file SpaceChargeStandard.cxx.

References fEnableCalSpatialSCE.

{
  return fEnableCalSpatialSCE;
}

bool spacecharge::SpaceChargeStandard::EnableCorrSCE ( ) const

overridevirtual

Return boolean indicating whether or not to apply SCE corrections.

Implements spacecharge::SpaceCharge.

Definition at line 174 of file SpaceChargeStandard.cxx.

References fEnableCorrSCE.

{
  return fEnableCorrSCE;
}

bool spacecharge::SpaceChargeStandard::EnableSimEfieldSCE ( ) const

overridevirtual

Return boolean indicating whether or not to turn simulation of SCE on for E field distortions

Implements spacecharge::SpaceCharge.

Definition at line 167 of file SpaceChargeStandard.cxx.

References fEnableSimEfieldSCE.

{
  return fEnableSimEfieldSCE;
}

bool spacecharge::SpaceChargeStandard::EnableSimSpatialSCE ( ) const

overridevirtual

Return boolean indicating whether or not to turn simulation of SCE on for spatial distortions

Implements spacecharge::SpaceCharge.

Definition at line 159 of file SpaceChargeStandard.cxx.

References fEnableSimSpatialSCE.

{
  return fEnableSimSpatialSCE;
}

geo::Vector_t spacecharge::SpaceChargeStandard::GetCalEfieldOffsets ( geo::Point_t const &  point, int const &  TPCid  ) const

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 363 of file SpaceChargeStandard.cxx.

{

  return {0., 0., 0.};
}

geo::Vector_t spacecharge::SpaceChargeStandard::GetCalPosOffsets ( geo::Point_t const &  point, int const &  TPCid  ) const

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 211 of file SpaceChargeStandard.cxx.

{

  return {0., 0., 0.};
}

geo::Vector_t spacecharge::SpaceChargeStandard::GetEfieldOffsets ( geo::Point_t const &  point ) const

overridevirtual

Primary working method of service that provides E field offsets to be used in charge/light yield calculation (e.g.)

Implements spacecharge::SpaceCharge.

Definition at line 351 of file SpaceChargeStandard.cxx.

References fRepresentationType, and GetEfieldOffsetsParametric().

{
  std::vector<double> theEfieldOffsets;

  if (fRepresentationType == "Parametric")
    theEfieldOffsets = GetEfieldOffsetsParametric(point.X(), point.Y(), point.Z());
  else
    theEfieldOffsets.resize(3, 0.0);

  return {-theEfieldOffsets[0], -theEfieldOffsets[1], -theEfieldOffsets[2]};
}

std::vector< double > spacecharge::SpaceChargeStandard::GetEfieldOffsetsParametric ( double  xVal, double  yVal, double  zVal  ) const

protected

Provides E field offsets using a parametric representation.

Definition at line 372 of file SpaceChargeStandard.cxx.

References GetOneEfieldOffsetParametric(), TransformX(), TransformY(), and TransformZ().

Referenced by GetEfieldOffsets().

{
  std::vector<double> theEfieldOffsetsParametric;

  double xValNew = TransformX(xVal);
  double yValNew = TransformY(yVal);
  double zValNew = TransformZ(zVal);

  theEfieldOffsetsParametric.push_back(
    GetOneEfieldOffsetParametric(xValNew, yValNew, zValNew, "X"));
  theEfieldOffsetsParametric.push_back(
    GetOneEfieldOffsetParametric(xValNew, yValNew, zValNew, "Y"));
  theEfieldOffsetsParametric.push_back(
    GetOneEfieldOffsetParametric(xValNew, yValNew, zValNew, "Z"));

  return theEfieldOffsetsParametric;
}

double spacecharge::SpaceChargeStandard::GetOneEfieldOffsetParametric ( double  xValNew, double  yValNew, double  zValNew, std::string  axis  ) const

protected

Provides one E field offset using a parametric representation, for a given axis, with returned E field offsets normalized to nominal drift E field

Definition at line 395 of file SpaceChargeStandard.cxx.

References f1_Ex, f1_Ey, f1_Ez, f2_Ex, f2_Ey, f2_Ez, f3_Ex, f3_Ey, f3_Ez, f4_Ex, f4_Ey, f4_Ez, f5_Ex, f5_Ey, f6_Ey, fFinal_Ex, fFinal_Ey, fFinal_Ez, g1_Ex, g1_Ey, g1_Ez, g2_Ex, g2_Ey, g2_Ez, g3_Ex, g3_Ey, g3_Ez, g4_Ex, g4_Ey, g4_Ez, g5_Ex, g5_Ey, and g6_Ey.

Referenced by GetEfieldOffsetsParametric().

{
  double parA[6][7];
  double parB[6];

  for (int j = 0; j < 6; j++) {
    for (int i = 0; i < 7; i++)
      parA[j][i] = 0.0;

    parB[j] = 0.0;
  }

  if (axis == "X") {
    for (int j = 0; j < 7; j++) {
      parA[0][j] = g1_Ex[j]->Eval(zValNew);
      parA[1][j] = g2_Ex[j]->Eval(zValNew);
      parA[2][j] = g3_Ex[j]->Eval(zValNew);
      parA[3][j] = g4_Ex[j]->Eval(zValNew);
      parA[4][j] = g5_Ex[j]->Eval(zValNew);
    }

    f1_Ex->SetParameters(parA[0]);
    f2_Ex->SetParameters(parA[1]);
    f3_Ex->SetParameters(parA[2]);
    f4_Ex->SetParameters(parA[3]);
    f5_Ex->SetParameters(parA[4]);
  }
  else if (axis == "Y") {
    for (int j = 0; j < 6; j++) {
      parA[0][j] = g1_Ey[j]->Eval(zValNew);
      parA[1][j] = g2_Ey[j]->Eval(zValNew);
      parA[2][j] = g3_Ey[j]->Eval(zValNew);
      parA[3][j] = g4_Ey[j]->Eval(zValNew);
      parA[4][j] = g5_Ey[j]->Eval(zValNew);
      parA[5][j] = g6_Ey[j]->Eval(zValNew);
    }

    f1_Ey->SetParameters(parA[0]);
    f2_Ey->SetParameters(parA[1]);
    f3_Ey->SetParameters(parA[2]);
    f4_Ey->SetParameters(parA[3]);
    f5_Ey->SetParameters(parA[4]);
    f6_Ey->SetParameters(parA[5]);
  }
  else if (axis == "Z") {
    for (int j = 0; j < 5; j++) {
      parA[0][j] = g1_Ez[j]->Eval(zValNew);
      parA[1][j] = g2_Ez[j]->Eval(zValNew);
      parA[2][j] = g3_Ez[j]->Eval(zValNew);
      parA[3][j] = g4_Ez[j]->Eval(zValNew);
    }

    f1_Ez->SetParameters(parA[0]);
    f2_Ez->SetParameters(parA[1]);
    f3_Ez->SetParameters(parA[2]);
    f4_Ez->SetParameters(parA[3]);
  }

  double aValNew;
  double bValNew;

  if (axis == "Y") {
    aValNew = xValNew;
    bValNew = yValNew;
  }
  else {
    aValNew = yValNew;
    bValNew = xValNew;
  }

  double offsetValNew = 0.0;
  if (axis == "X") {
    parB[0] = f1_Ex->Eval(aValNew);
    parB[1] = f2_Ex->Eval(aValNew);
    parB[2] = f3_Ex->Eval(aValNew);
    parB[3] = f4_Ex->Eval(aValNew);
    parB[4] = f5_Ex->Eval(aValNew);

    fFinal_Ex->SetParameters(parB);
    offsetValNew = fFinal_Ex->Eval(bValNew);
  }
  else if (axis == "Y") {
    parB[0] = f1_Ey->Eval(aValNew);
    parB[1] = f2_Ey->Eval(aValNew);
    parB[2] = f3_Ey->Eval(aValNew);
    parB[3] = f4_Ey->Eval(aValNew);
    parB[4] = f5_Ey->Eval(aValNew);
    parB[5] = f6_Ey->Eval(aValNew);

    fFinal_Ey->SetParameters(parB);
    offsetValNew = fFinal_Ey->Eval(bValNew);
  }
  else if (axis == "Z") {
    parB[0] = f1_Ez->Eval(aValNew);
    parB[1] = f2_Ez->Eval(aValNew);
    parB[2] = f3_Ez->Eval(aValNew);
    parB[3] = f4_Ez->Eval(aValNew);

    fFinal_Ez->SetParameters(parB);
    offsetValNew = fFinal_Ez->Eval(bValNew);
  }

  return offsetValNew;
}

double spacecharge::SpaceChargeStandard::GetOnePosOffsetParametric ( double  xValNew, double  yValNew, double  zValNew, std::string  axis  ) const

protected

Provides one position offset using a parametric representation, for a given axis

Definition at line 240 of file SpaceChargeStandard.cxx.

References f1_x, f1_y, f1_z, f2_x, f2_y, f2_z, f3_x, f3_y, f3_z, f4_x, f4_y, f4_z, f5_x, f5_y, f6_y, fFinal_x, fFinal_y, fFinal_z, g1_x, g1_y, g1_z, g2_x, g2_y, g2_z, g3_x, g3_y, g3_z, g4_x, g4_y, g4_z, g5_x, g5_y, and g6_y.

Referenced by GetPosOffsetsParametric().

{
  double parA[6][7];
  double parB[6];

  for (int j = 0; j < 6; j++) {
    for (int i = 0; i < 7; i++)
      parA[j][i] = 0.0;

    parB[j] = 0.0;
  }

  if (axis == "X") {
    for (int j = 0; j < 7; j++) {
      parA[0][j] = g1_x[j]->Eval(zValNew);
      parA[1][j] = g2_x[j]->Eval(zValNew);
      parA[2][j] = g3_x[j]->Eval(zValNew);
      parA[3][j] = g4_x[j]->Eval(zValNew);
      parA[4][j] = g5_x[j]->Eval(zValNew);
    }

    f1_x->SetParameters(parA[0]);
    f2_x->SetParameters(parA[1]);
    f3_x->SetParameters(parA[2]);
    f4_x->SetParameters(parA[3]);
    f5_x->SetParameters(parA[4]);
  }
  else if (axis == "Y") {
    for (int j = 0; j < 6; j++) {
      parA[0][j] = g1_y[j]->Eval(zValNew);
      parA[1][j] = g2_y[j]->Eval(zValNew);
      parA[2][j] = g3_y[j]->Eval(zValNew);
      parA[3][j] = g4_y[j]->Eval(zValNew);
      parA[4][j] = g5_y[j]->Eval(zValNew);
      parA[5][j] = g6_y[j]->Eval(zValNew);
    }

    f1_y->SetParameters(parA[0]);
    f2_y->SetParameters(parA[1]);
    f3_y->SetParameters(parA[2]);
    f4_y->SetParameters(parA[3]);
    f5_y->SetParameters(parA[4]);
    f6_y->SetParameters(parA[5]);
  }
  else if (axis == "Z") {
    for (int j = 0; j < 5; j++) {
      parA[0][j] = g1_z[j]->Eval(zValNew);
      parA[1][j] = g2_z[j]->Eval(zValNew);
      parA[2][j] = g3_z[j]->Eval(zValNew);
      parA[3][j] = g4_z[j]->Eval(zValNew);
    }

    f1_z->SetParameters(parA[0]);
    f2_z->SetParameters(parA[1]);
    f3_z->SetParameters(parA[2]);
    f4_z->SetParameters(parA[3]);
  }

  double aValNew;
  double bValNew;

  if (axis == "Y") {
    aValNew = xValNew;
    bValNew = yValNew;
  }
  else {
    aValNew = yValNew;
    bValNew = xValNew;
  }

  double offsetValNew = 0.0;
  if (axis == "X") {
    parB[0] = f1_x->Eval(aValNew);
    parB[1] = f2_x->Eval(aValNew);
    parB[2] = f3_x->Eval(aValNew);
    parB[3] = f4_x->Eval(aValNew);
    parB[4] = f5_x->Eval(aValNew);

    fFinal_x->SetParameters(parB);
    offsetValNew = 100.0 * fFinal_x->Eval(bValNew);
  }
  else if (axis == "Y") {
    parB[0] = f1_y->Eval(aValNew);
    parB[1] = f2_y->Eval(aValNew);
    parB[2] = f3_y->Eval(aValNew);
    parB[3] = f4_y->Eval(aValNew);
    parB[4] = f5_y->Eval(aValNew);
    parB[5] = f6_y->Eval(aValNew);

    fFinal_y->SetParameters(parB);
    offsetValNew = 100.0 * fFinal_y->Eval(bValNew);
  }
  else if (axis == "Z") {
    parB[0] = f1_z->Eval(aValNew);
    parB[1] = f2_z->Eval(aValNew);
    parB[2] = f3_z->Eval(aValNew);
    parB[3] = f4_z->Eval(aValNew);

    fFinal_z->SetParameters(parB);
    offsetValNew = 100.0 * fFinal_z->Eval(bValNew);
  }

  return offsetValNew;
}

geo::Vector_t spacecharge::SpaceChargeStandard::GetPosOffsets ( geo::Point_t const &  point ) const

overridevirtual

Primary working method of service that provides position offsets to be used in ionization electron drift

Implements spacecharge::SpaceCharge.

Definition at line 194 of file SpaceChargeStandard.cxx.

References fRepresentationType, GetPosOffsetsParametric(), and IsInsideBoundaries().

{
  std::vector<double> thePosOffsets;

  if (IsInsideBoundaries(point.X(), point.Y(), point.Z()) == false) {
    thePosOffsets.resize(3, 0.0);
  }
  else {
    if (fRepresentationType == "Parametric")
      thePosOffsets = GetPosOffsetsParametric(point.X(), point.Y(), point.Z());
    else
      thePosOffsets.resize(3, 0.0);
  }

  return {thePosOffsets[0], thePosOffsets[1], thePosOffsets[2]};
}

std::vector< double > spacecharge::SpaceChargeStandard::GetPosOffsetsParametric ( double  xVal, double  yVal, double  zVal  ) const

protected

Provides position offsets using a parametric representation.

Definition at line 220 of file SpaceChargeStandard.cxx.

References GetOnePosOffsetParametric(), TransformX(), TransformY(), and TransformZ().

Referenced by GetPosOffsets().

{
  std::vector<double> thePosOffsetsParametric;

  double xValNew = TransformX(xVal);
  double yValNew = TransformY(yVal);
  double zValNew = TransformZ(zVal);

  thePosOffsetsParametric.push_back(GetOnePosOffsetParametric(xValNew, yValNew, zValNew, "X"));
  thePosOffsetsParametric.push_back(GetOnePosOffsetParametric(xValNew, yValNew, zValNew, "Y"));
  thePosOffsetsParametric.push_back(GetOnePosOffsetParametric(xValNew, yValNew, zValNew, "Z"));

  return thePosOffsetsParametric;
}

bool spacecharge::SpaceChargeStandard::IsInsideBoundaries ( double  xVal, double  yVal, double  zVal  ) const

protected

Check to see if point is inside boundaries of map - redefine this in experiment-specific implementation!

Definition at line 535 of file SpaceChargeStandard.cxx.

Referenced by GetPosOffsets().

{
  bool isInside = false;

  return isInside;
}

double spacecharge::SpaceChargeStandard::TransformX ( double  xVal ) const

protected

Transform X to SCE X coordinate - redefine this in experiment-specific implementation!

Definition at line 505 of file SpaceChargeStandard.cxx.

Referenced by GetEfieldOffsetsParametric(), and GetPosOffsetsParametric().

{
  double xValNew;
  xValNew = xVal / 100.0;

  return xValNew;
}

double spacecharge::SpaceChargeStandard::TransformY ( double  yVal ) const

protected

Transform Y to SCE Y coordinate - redefine this in experiment-specific implementation!

Definition at line 515 of file SpaceChargeStandard.cxx.

Referenced by GetEfieldOffsetsParametric(), and GetPosOffsetsParametric().

{
  double yValNew;
  yValNew = yVal / 100.0;

  return yValNew;
}

double spacecharge::SpaceChargeStandard::TransformZ ( double  zVal ) const

protected

Transform Z to SCE Z coordinate - redefine this in experiment-specific implementation!

Definition at line 525 of file SpaceChargeStandard.cxx.

Referenced by GetEfieldOffsetsParametric(), and GetPosOffsetsParametric().

{
  double zValNew;
  zValNew = zVal / 100.0;

  return zValNew;
}

bool spacecharge::SpaceChargeStandard::Update

(

uint64_t

ts = 0

)

Definition at line 149 of file SpaceChargeStandard.cxx.

Referenced by spacecharge::SpaceChargeServiceStandard::preBeginRun().

{
  if (ts == 0) return false;

  return true;
}

Member Data Documentation

TF1* spacecharge::SpaceChargeStandard::f1_Ex = new TF1("f1_Ex", "pol6")

protected

Definition at line 133 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f1_Ey = new TF1("f1_Ey", "pol5")

protected

Definition at line 140 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f1_Ez = new TF1("f1_Ez", "pol4")

protected

Definition at line 148 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f1_x = new TF1("f1_x", "pol6")

protected

Definition at line 94 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f1_y = new TF1("f1_y", "pol5")

protected

Definition at line 101 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f1_z = new TF1("f1_z", "pol4")

protected

Definition at line 109 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f2_Ex = new TF1("f2_Ex", "pol6")

protected

Definition at line 134 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f2_Ey = new TF1("f2_Ey", "pol5")

protected

Definition at line 141 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f2_Ez = new TF1("f2_Ez", "pol4")

protected

Definition at line 149 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f2_x = new TF1("f2_x", "pol6")

protected

Definition at line 95 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f2_y = new TF1("f2_y", "pol5")

protected

Definition at line 102 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f2_z = new TF1("f2_z", "pol4")

protected

Definition at line 110 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f3_Ex = new TF1("f3_Ex", "pol6")

protected

Definition at line 135 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f3_Ey = new TF1("f3_Ey", "pol5")

protected

Definition at line 142 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f3_Ez = new TF1("f3_Ez", "pol4")

protected

Definition at line 150 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f3_x = new TF1("f3_x", "pol6")

protected

Definition at line 96 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f3_y = new TF1("f3_y", "pol5")

protected

Definition at line 103 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f3_z = new TF1("f3_z", "pol4")

protected

Definition at line 111 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f4_Ex = new TF1("f4_Ex", "pol6")

protected

Definition at line 136 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f4_Ey = new TF1("f4_Ey", "pol5")

protected

Definition at line 143 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f4_Ez = new TF1("f4_Ez", "pol4")

protected

Definition at line 151 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f4_x = new TF1("f4_x", "pol6")

protected

Definition at line 97 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f4_y = new TF1("f4_y", "pol5")

protected

Definition at line 104 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f4_z = new TF1("f4_z", "pol4")

protected

Definition at line 112 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f5_Ex = new TF1("f5_Ex", "pol6")

protected

Definition at line 137 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f5_Ey = new TF1("f5_Ey", "pol5")

protected

Definition at line 144 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f5_x = new TF1("f5_x", "pol6")

protected

Definition at line 98 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f5_y = new TF1("f5_y", "pol5")

protected

Definition at line 105 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f6_Ey = new TF1("f6_Ey", "pol5")

protected

Definition at line 145 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::f6_y = new TF1("f6_y", "pol5")

protected

Definition at line 106 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

bool spacecharge::SpaceChargeStandard::fEnableCalEfieldSCE

protected

Definition at line 70 of file SpaceChargeStandard.h.

Referenced by Configure(), and EnableCalEfieldSCE().

bool spacecharge::SpaceChargeStandard::fEnableCalSpatialSCE

protected

Definition at line 69 of file SpaceChargeStandard.h.

Referenced by Configure(), and EnableCalSpatialSCE().

bool spacecharge::SpaceChargeStandard::fEnableCorrSCE

protected

Definition at line 71 of file SpaceChargeStandard.h.

Referenced by Configure(), and EnableCorrSCE().

bool spacecharge::SpaceChargeStandard::fEnableSimEfieldSCE

protected

Definition at line 68 of file SpaceChargeStandard.h.

Referenced by Configure(), and EnableSimEfieldSCE().

bool spacecharge::SpaceChargeStandard::fEnableSimSpatialSCE

protected

Definition at line 67 of file SpaceChargeStandard.h.

Referenced by Configure(), and EnableSimSpatialSCE().

TF1* spacecharge::SpaceChargeStandard::fFinal_Ex = new TF1("fFinal_Ex", "pol4")

protected

Definition at line 138 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::fFinal_Ey = new TF1("fFinal_Ey", "pol5")

protected

Definition at line 146 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::fFinal_Ez = new TF1("fFinal_Ez", "pol3")

protected

Definition at line 152 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::fFinal_x = new TF1("fFinal_x", "pol4")

protected

Definition at line 99 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::fFinal_y = new TF1("fFinal_y", "pol5")

protected

Definition at line 107 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

TF1* spacecharge::SpaceChargeStandard::fFinal_z = new TF1("fFinal_z", "pol3")

protected

Definition at line 113 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

std::string spacecharge::SpaceChargeStandard::fInputFilename

protected

Definition at line 74 of file SpaceChargeStandard.h.

Referenced by Configure().

std::string spacecharge::SpaceChargeStandard::fRepresentationType

protected

Definition at line 73 of file SpaceChargeStandard.h.

Referenced by Configure(), GetEfieldOffsets(), and GetPosOffsets().

TGraph** spacecharge::SpaceChargeStandard::g1_Ex = new TGraph*[7]

protected

Definition at line 115 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g1_Ey = new TGraph*[7]

protected

Definition at line 121 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g1_Ez = new TGraph*[7]

protected

Definition at line 128 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g1_x = new TGraph*[7]

protected

Definition at line 76 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g1_y = new TGraph*[7]

protected

Definition at line 82 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g1_z = new TGraph*[7]

protected

Definition at line 89 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g2_Ex = new TGraph*[7]

protected

Definition at line 116 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g2_Ey = new TGraph*[7]

protected

Definition at line 122 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g2_Ez = new TGraph*[7]

protected

Definition at line 129 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g2_x = new TGraph*[7]

protected

Definition at line 77 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g2_y = new TGraph*[7]

protected

Definition at line 83 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g2_z = new TGraph*[7]

protected

Definition at line 90 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g3_Ex = new TGraph*[7]

protected

Definition at line 117 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g3_Ey = new TGraph*[7]

protected

Definition at line 123 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g3_Ez = new TGraph*[7]

protected

Definition at line 130 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g3_x = new TGraph*[7]

protected

Definition at line 78 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g3_y = new TGraph*[7]

protected

Definition at line 84 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g3_z = new TGraph*[7]

protected

Definition at line 91 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g4_Ex = new TGraph*[7]

protected

Definition at line 118 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g4_Ey = new TGraph*[7]

protected

Definition at line 124 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g4_Ez = new TGraph*[7]

protected

Definition at line 131 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g4_x = new TGraph*[7]

protected

Definition at line 79 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g4_y = new TGraph*[7]

protected

Definition at line 85 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g4_z = new TGraph*[7]

protected

Definition at line 92 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g5_Ex = new TGraph*[7]

protected

Definition at line 119 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g5_Ey = new TGraph*[7]

protected

Definition at line 125 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g5_x = new TGraph*[7]

protected

Definition at line 80 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g5_y = new TGraph*[7]

protected

Definition at line 86 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g6_Ey = new TGraph*[7]

protected

Definition at line 126 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

TGraph** spacecharge::SpaceChargeStandard::g6_y = new TGraph*[7]

protected

Definition at line 87 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

---

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

• larevt/v09_10_03/source/larevt/SpaceCharge/SpaceChargeStandard.h
• larevt/v09_10_03/source/larevt/SpaceCharge/SpaceChargeStandard.cxx