#include <SpaceChargeDUNE35t.h>

Inheritance diagram for spacecharge::SpaceChargeDUNE35t:

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

	SpaceChargeDUNE35t (SpaceChargeDUNE35t const &)=delete

virtual	~SpaceChargeDUNE35t ()=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: [2.275,-0.3] –> [0.0,2.25]. More...

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

double	TransformZ (double zVal) const
	Transform Z to SCE Z coordinate: [0.0,1.6] –> [0.0,1.6]. 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 30 of file SpaceChargeDUNE35t.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 25 of file SpaceChargeDUNE35t.cxx.

 {
   Configure(pset);
 }

spacecharge::SpaceChargeDUNE35t::SpaceChargeDUNE35t ( SpaceChargeDUNE35t const & )

delete

virtual spacecharge::SpaceChargeDUNE35t::~SpaceChargeDUNE35t ( )

virtualdefault

Member Function Documentation

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

Definition at line 33 of file SpaceChargeDUNE35t.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("SpaceChargeDUNE35t") << "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::SpaceChargeDUNE35t::EnableCalEfieldSCE ( ) const

overridevirtual

Return boolean indicating whether or not to apply SCE corrections.

Implements spacecharge::SpaceCharge.

Definition at line 182 of file SpaceChargeDUNE35t.cxx.

 {
   return fEnableCalEfieldSCE;
 }

bool spacecharge::SpaceChargeDUNE35t::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 176 of file SpaceChargeDUNE35t.cxx.

 {
   return fEnableCalSpatialSCE;
 }

bool spacecharge::SpaceChargeDUNE35t::EnableCorrSCE ( ) const

inlineoverridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 46 of file SpaceChargeDUNE35t.h.

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

spacecharge::SpaceChargeDUNE35t::EnableCalSpatialSCE

bool EnableCalSpatialSCE() const override

Return boolean indicating whether or not to apply SCE corrections.

Definition: SpaceChargeDUNE35t.cxx:176

spacecharge::SpaceChargeDUNE35t::EnableCalEfieldSCE

bool EnableCalEfieldSCE() const override

Return boolean indicating whether or not to apply SCE corrections.

Definition: SpaceChargeDUNE35t.cxx:182

bool spacecharge::SpaceChargeDUNE35t::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 163 of file SpaceChargeDUNE35t.cxx.

 {
   return fEnableSimEfieldSCE;
 }

bool spacecharge::SpaceChargeDUNE35t::EnableSimSpatialSCE ( ) const

overridevirtual

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

Implements spacecharge::SpaceCharge.

Definition at line 155 of file SpaceChargeDUNE35t.cxx.

 {
   return fEnableSimSpatialSCE;
 }

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

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 374 of file SpaceChargeDUNE35t.cxx.

 {
   return { 0.0, 0.0, 0.0 };
 }

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

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 209 of file SpaceChargeDUNE35t.cxx.

 {
 
   return { 0.0, 0.0, 0.0 };
 }

geo::Vector_t spacecharge::SpaceChargeDUNE35t::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 SpaceChargeDUNE35t.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::SpaceChargeDUNE35t::GetEfieldOffsetsParametric	(	double	xVal,
		double	yVal,
		double	zVal
	)		const

protected

Provides E field offsets using a parametric representation.

Definition at line 381 of file SpaceChargeDUNE35t.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::SpaceChargeDUNE35t::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 SpaceChargeDUNE35t.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::SpaceChargeDUNE35t::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 SpaceChargeDUNE35t.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 = 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::SpaceChargeDUNE35t::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 190 of file SpaceChargeDUNE35t.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::SpaceChargeDUNE35t::GetPosOffsetsParametric	(	double	xVal,
		double	yVal,
		double	zVal
	)		const

protected

Provides position offsets using a parametric representation.

Definition at line 217 of file SpaceChargeDUNE35t.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::SpaceChargeDUNE35t::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 550 of file SpaceChargeDUNE35t.cxx.

 {
   bool isInside = true;
 
   if((xVal < 2.5) || (xVal > 230.0) || (yVal < -5.0) || (yVal > 205.0) || (zVal < -5.0) || (zVal > 165.0))
   {
     isInside = false;
   }
 
   return isInside;
 }

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

protected

Transform X to SCE X coordinate: [2.275,-0.3] –> [0.0,2.25].

Definition at line 518 of file SpaceChargeDUNE35t.cxx.

 {
   double xValNew;
   xValNew = 2.25 - (2.25/2.25)*((xVal-2.5)/100.0);
   xValNew -= 1.125;
 
   return xValNew;
 }

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

protected

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

Definition at line 529 of file SpaceChargeDUNE35t.cxx.

 {
   double yValNew;
   yValNew = (2.0/2.0)*(yVal/100.0);
   yValNew -= 1.0;
 
   return yValNew;
 }

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

protected

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

Definition at line 540 of file SpaceChargeDUNE35t.cxx.

 {
   double zValNew;
   zValNew = (1.6/1.6)*(zVal/100.0);
 
   return zValNew;
 }

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

Definition at line 145 of file SpaceChargeDUNE35t.cxx.

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

Member Data Documentation

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

protected

Definition at line 131 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 138 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 146 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 92 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 99 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 107 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 132 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 139 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 147 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 93 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 100 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 108 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 133 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 140 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 148 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 94 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 101 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 109 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 134 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 141 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 149 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 95 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 102 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 110 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 135 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 142 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 96 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 103 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 143 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 104 of file SpaceChargeDUNE35t.h.

bool spacecharge::SpaceChargeDUNE35t::fEnableCalEfieldSCE

protected

Definition at line 68 of file SpaceChargeDUNE35t.h.

bool spacecharge::SpaceChargeDUNE35t::fEnableCalSpatialSCE

protected

Definition at line 67 of file SpaceChargeDUNE35t.h.

bool spacecharge::SpaceChargeDUNE35t::fEnableCorrSCE

protected

Definition at line 69 of file SpaceChargeDUNE35t.h.

bool spacecharge::SpaceChargeDUNE35t::fEnableSimEfieldSCE

protected

Definition at line 66 of file SpaceChargeDUNE35t.h.

bool spacecharge::SpaceChargeDUNE35t::fEnableSimSpatialSCE

protected

Definition at line 65 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 136 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 144 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 150 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 97 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 105 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 111 of file SpaceChargeDUNE35t.h.

std::string spacecharge::SpaceChargeDUNE35t::fInputFilename

protected

Definition at line 72 of file SpaceChargeDUNE35t.h.

std::string spacecharge::SpaceChargeDUNE35t::fRepresentationType

protected

Definition at line 71 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 113 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 119 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 126 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 74 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 80 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 87 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 114 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 120 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 127 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 75 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 81 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 88 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 115 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 121 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 128 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 76 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 82 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 89 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 116 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 122 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 129 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 77 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 83 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 90 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 117 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 123 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 78 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 84 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 124 of file SpaceChargeDUNE35t.h.

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

protected

Definition at line 85 of file SpaceChargeDUNE35t.h.

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

dunesim/dunesim/SpaceCharge/SpaceChargeDUNE35t.h
dunesim/dunesim/SpaceCharge/SpaceChargeDUNE35t.cxx

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation