#include <SpaceCharge3x1x1dphase.h>

Inheritance diagram for spacecharge::SpaceCharge3x1x1dphase:

Public Member Functions
	SpaceCharge3x1x1dphase (fhicl::ParameterSet const &pset)

	SpaceCharge3x1x1dphase (SpaceCharge3x1x1dphase const &)=delete

virtual	~SpaceCharge3x1x1dphase ()=default

bool	Configure (fhicl::ParameterSet const &pset)

bool	Update (uint64_t ts=0)

bool	EnableSimSpatialSCE () const override

bool	EnableSimEfieldSCE () const override

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...

bool	EnableCorrSCE () const override

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

Public Member Functions inherited from spacecharge::SpaceCharge
	SpaceCharge (const SpaceCharge &)=delete

	SpaceCharge (SpaceCharge &&)=delete

SpaceCharge &	operator= (const SpaceCharge &)=delete

SpaceCharge &	operator= (SpaceCharge &&)=delete

virtual	~SpaceCharge ()=default

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: [0.0,3.6] –> [0.0,3.6]. More...

double	TransformY (double yVal) const
	Transform Y to SCE Y coordinate: [0.0,6.08] –> [0.0,6.0]. More...

double	TransformZ (double zVal) const
	Transform Z to SCE Z coordinate: [0.0,6.97] –> [0.0,7.2]. More...

bool	IsInsideBoundaries (double xVal, double yVal, double zVal) const
	Check to see if point is inside boundaries of map (allow to go slightly out of range) More...

Protected Member Functions inherited from spacecharge::SpaceCharge
	SpaceCharge ()=default

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 31 of file SpaceCharge3x1x1dphase.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 26 of file SpaceCharge3x1x1dphase.cxx.

 {
   Configure(pset);
 }

spacecharge::SpaceCharge3x1x1dphase::SpaceCharge3x1x1dphase ( SpaceCharge3x1x1dphase const & )

delete

virtual spacecharge::SpaceCharge3x1x1dphase::~SpaceCharge3x1x1dphase ( )

virtualdefault

Member Function Documentation

bool spacecharge::SpaceCharge3x1x1dphase::Configure ( fhicl::ParameterSet const & pset )

Definition at line 34 of file SpaceCharge3x1x1dphase.cxx.

 {  
   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");
 
   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 cet::exception("SpaceCharge3x1x1dphase") << "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::SpaceCharge3x1x1dphase::EnableCalEfieldSCE ( ) const

overridevirtual

Return boolean indicating whether or not to apply SCE corrections.

Implements spacecharge::SpaceCharge.

Definition at line 183 of file SpaceCharge3x1x1dphase.cxx.

 {
   return fEnableCalEfieldSCE;
 }

bool spacecharge::SpaceCharge3x1x1dphase::EnableCalSpatialSCE ( ) const

overridevirtual

Return boolean indicating whether or not to apply SCE corrections.

Return boolean indicating whether or not to apply SCE corrections

Implements spacecharge::SpaceCharge.

Definition at line 177 of file SpaceCharge3x1x1dphase.cxx.

 {
   return fEnableCalSpatialSCE;
 }

bool spacecharge::SpaceCharge3x1x1dphase::EnableCorrSCE ( ) const

inlineoverridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 47 of file SpaceCharge3x1x1dphase.h.

47 {return (EnableCalSpatialSCE()||EnableCalEfieldSCE()) ;}

spacecharge::SpaceCharge3x1x1dphase::EnableCalSpatialSCE

bool EnableCalSpatialSCE() const override

Return boolean indicating whether or not to apply SCE corrections.

Definition: SpaceCharge3x1x1dphase.cxx:177

spacecharge::SpaceCharge3x1x1dphase::EnableCalEfieldSCE

bool EnableCalEfieldSCE() const override

Return boolean indicating whether or not to apply SCE corrections.

Definition: SpaceCharge3x1x1dphase.cxx:183

bool spacecharge::SpaceCharge3x1x1dphase::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 164 of file SpaceCharge3x1x1dphase.cxx.

 {
   return fEnableSimEfieldSCE;
 }

bool spacecharge::SpaceCharge3x1x1dphase::EnableSimSpatialSCE ( ) const

overridevirtual

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

Implements spacecharge::SpaceCharge.

Definition at line 156 of file SpaceCharge3x1x1dphase.cxx.

 {
   return fEnableSimSpatialSCE;
 }

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

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 374 of file SpaceCharge3x1x1dphase.cxx.

 {
   return { 0.0, 0.0, 0.0 };
 }

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

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 210 of file SpaceCharge3x1x1dphase.cxx.

 {
   return { 0.0, 0.0, 0.0 };
 }

geo::Vector_t spacecharge::SpaceCharge3x1x1dphase::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 355 of file SpaceCharge3x1x1dphase.cxx.

 {
   std::vector<double> theEfieldOffsets;
 
   if(IsInsideBoundaries(point.X(), point.Y(), point.Z()) == false)
   {
     theEfieldOffsets.resize(3,0.0);
   }
   else
   {
     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::SpaceCharge3x1x1dphase::GetEfieldOffsetsParametric	(	double	xVal,
		double	yVal,
		double	zVal
	)		const

protected

Provides E field offsets using a parametric representation.

Definition at line 381 of file SpaceCharge3x1x1dphase.cxx.

 {
   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::SpaceCharge3x1x1dphase::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

Definition at line 399 of file SpaceCharge3x1x1dphase.cxx.

 {      
   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::SpaceCharge3x1x1dphase::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 235 of file SpaceCharge3x1x1dphase.cxx.

 {      
   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 = 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 = 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 = fFinal_z->Eval(bValNew);
   }
    
   return offsetValNew;
 }

geo::Vector_t spacecharge::SpaceCharge3x1x1dphase::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 191 of file SpaceCharge3x1x1dphase.cxx.

 {
   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::SpaceCharge3x1x1dphase::GetPosOffsetsParametric	(	double	xVal,
		double	yVal,
		double	zVal
	)		const

protected

Provides position offsets using a parametric representation.

Definition at line 217 of file SpaceCharge3x1x1dphase.cxx.

 {
   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::SpaceCharge3x1x1dphase::IsInsideBoundaries	(	double	xVal,
		double	yVal,
		double	zVal
	)		const

protected

Check to see if point is inside boundaries of map (allow to go slightly out of range)

Definition at line 545 of file SpaceCharge3x1x1dphase.cxx.

 {
   bool isInside = true;
 
   if((xVal < -50.0) || (xVal > 50.0) || (yVal < -50.0) || (yVal > 50.0) || (zVal < 0.0) || (zVal > 300.0))
   {
     isInside = false;
   }
 
   return isInside;
 }

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

protected

Transform X to SCE X coordinate: [0.0,3.6] –> [0.0,3.6].

Definition at line 518 of file SpaceCharge3x1x1dphase.cxx.

 {
   double xValNew = xVal;
 
   return xValNew;
 }

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

protected

Transform Y to SCE Y coordinate: [0.0,6.08] –> [0.0,6.0].

Definition at line 527 of file SpaceCharge3x1x1dphase.cxx.

 {
   double yValNew =  yVal;
 
   return yValNew;
 }

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

protected

Transform Z to SCE Z coordinate: [0.0,6.97] –> [0.0,7.2].

Definition at line 536 of file SpaceCharge3x1x1dphase.cxx.

 {
   double zValNew = zVal;
 
   return zValNew;
 }

bool spacecharge::SpaceCharge3x1x1dphase::Update ( uint64_t ts = 0 )

Definition at line 146 of file SpaceCharge3x1x1dphase.cxx.

 {
   if (ts == 0) return false;
 
   return true;
 }

Member Data Documentation

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

protected

Definition at line 132 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 139 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 147 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 93 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 100 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 108 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 133 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 140 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 148 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 94 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 101 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 109 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 134 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 141 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 149 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 95 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 102 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 110 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 135 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 142 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 150 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 96 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 103 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 111 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 136 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 143 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 97 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 104 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 144 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 105 of file SpaceCharge3x1x1dphase.h.

bool spacecharge::SpaceCharge3x1x1dphase::fEnableCalEfieldSCE

protected

Definition at line 69 of file SpaceCharge3x1x1dphase.h.

bool spacecharge::SpaceCharge3x1x1dphase::fEnableCalSpatialSCE

protected

Definition at line 68 of file SpaceCharge3x1x1dphase.h.

bool spacecharge::SpaceCharge3x1x1dphase::fEnableCorrSCE

protected

Definition at line 70 of file SpaceCharge3x1x1dphase.h.

bool spacecharge::SpaceCharge3x1x1dphase::fEnableSimEfieldSCE

protected

Definition at line 67 of file SpaceCharge3x1x1dphase.h.

bool spacecharge::SpaceCharge3x1x1dphase::fEnableSimSpatialSCE

protected

Definition at line 66 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 137 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 145 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 151 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 98 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 106 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 112 of file SpaceCharge3x1x1dphase.h.

std::string spacecharge::SpaceCharge3x1x1dphase::fInputFilename

protected

Definition at line 73 of file SpaceCharge3x1x1dphase.h.

std::string spacecharge::SpaceCharge3x1x1dphase::fRepresentationType

protected

Definition at line 72 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 114 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 120 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 127 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 75 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 81 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 88 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 115 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 121 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 128 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 76 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 82 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 89 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 116 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 122 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 129 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 77 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 83 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 90 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 117 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 123 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 130 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 78 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 84 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 91 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 118 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 124 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 79 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 85 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 125 of file SpaceCharge3x1x1dphase.h.

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

protected

Definition at line 86 of file SpaceCharge3x1x1dphase.h.

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

dunesim/dunesim/SpaceCharge/SpaceCharge3x1x1dphase.h
dunesim/dunesim/SpaceCharge/SpaceCharge3x1x1dphase.cxx

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation