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...
geo::Vector_t	GetPosOffsets (geo::Point_t const &point) const override
geo::Vector_t	GetEfieldOffsets (geo::Point_t const &point) 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	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 31 of file SpaceChargeStandard.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 25 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 33 of file SpaceChargeStandard.cxx.

References art::errors::Configuration, 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 SpaceChargeStandard().

{  
  fEnableSimSpatialSCE = pset.get<bool>("EnableSimSpatialSCE");
  fEnableSimEfieldSCE = pset.get<bool>("EnableSimEfieldSCE");
  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);

    std::unique_ptr<TFile> infile(new 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::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::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 341 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 355 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 373 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 221 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 182 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 203 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 522 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 492 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 502 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 512 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.

{
  if (ts == 0) return false;

  return true;
}

Member Data Documentation

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

protected

Definition at line 124 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 131 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 139 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 85 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 92 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 100 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 125 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 132 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 140 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 86 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 93 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 101 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 126 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 133 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 141 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 87 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 94 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 102 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 127 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 134 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 142 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 88 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 95 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 103 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 128 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 135 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 89 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 96 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 136 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 97 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

bool spacecharge::SpaceChargeStandard::fEnableCorrSCE

protected

Definition at line 62 of file SpaceChargeStandard.h.

Referenced by Configure(), and EnableCorrSCE().

bool spacecharge::SpaceChargeStandard::fEnableSimEfieldSCE

protected

Definition at line 61 of file SpaceChargeStandard.h.

Referenced by Configure(), and EnableSimEfieldSCE().

bool spacecharge::SpaceChargeStandard::fEnableSimSpatialSCE

protected

Definition at line 60 of file SpaceChargeStandard.h.

Referenced by Configure(), and EnableSimSpatialSCE().

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

protected

Definition at line 129 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 137 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 143 of file SpaceChargeStandard.h.

Referenced by GetOneEfieldOffsetParametric().

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

protected

Definition at line 90 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 98 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

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

protected

Definition at line 104 of file SpaceChargeStandard.h.

Referenced by GetOnePosOffsetParametric().

std::string spacecharge::SpaceChargeStandard::fInputFilename

protected

Definition at line 65 of file SpaceChargeStandard.h.

Referenced by Configure().

std::string spacecharge::SpaceChargeStandard::fRepresentationType

protected

Definition at line 64 of file SpaceChargeStandard.h.

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

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

protected

Definition at line 106 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 112 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 119 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 67 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 73 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 80 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 107 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 113 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 120 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 68 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 74 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 81 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 108 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 114 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 121 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 69 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 75 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 82 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 109 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 115 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 122 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 70 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 76 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 83 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 110 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 116 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 71 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 77 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

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

protected

Definition at line 117 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOneEfieldOffsetParametric().

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

protected

Definition at line 78 of file SpaceChargeStandard.h.

Referenced by Configure(), and GetOnePosOffsetParametric().

---

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

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