#include <SpaceChargeSBND.h>

Inheritance diagram for spacecharge::SpaceChargeSBND:

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

	SpaceChargeSBND (SpaceChargeSBND const &)=delete

virtual	~SpaceChargeSBND ()=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

bool	EnableCalEfieldSCE () const override

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=1) const override

geo::Vector_t	GetCalEfieldOffsets (geo::Point_t const &point, int const &TPCid=1) 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

double	GetOnePosOffsetParametric (double xVal, double yVal, double zVal, std::string axis) const

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

double	GetOneEfieldOffsetParametric (double xVal, double yVal, double zVal, std::string axis) const

double	TransformX (double xVal) const

double	TransformY (double yVal) const

double	TransformZ (double zVal) const

bool	IsInsideBoundaries (double xVal, double yVal, double zVal) const

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

Protected Attributes
int	initialSpatialFitPolN [3] = {3, 4, 3}

int	intermediateSpatialFitPolN [3] = {4, 4, 4}

int	initialEFieldFitPolN [3] = {3, 3, 3}

int	intermediateEFieldFitPolN [3] = {6, 4, 4}

double	DriftField = 500.0

bool	fEnableSimSpatialSCE

bool	fEnableSimEfieldSCE

bool	fEnableCalSpatialSCE

bool	fEnableCalEfieldSCE

bool	fEnableCorrSCE

std::string	fRepresentationType

std::string	fInputFilename

std::vector< TH3F * >	SCEhistograms = std::vector<TH3F*>(9)

TGraph *	gSpatialGraphX [99][99]

TF1 *	intermediateSpatialFitFunctionX [99]

TF1 *	initialSpatialFitFunctionX

TGraph *	gSpatialGraphY [99][99]

TF1 *	intermediateSpatialFitFunctionY [99]

TF1 *	initialSpatialFitFunctionY

TGraph *	gSpatialGraphZ [99][99]

TF1 *	intermediateSpatialFitFunctionZ [99]

TF1 *	initialSpatialFitFunctionZ

TGraph *	gEFieldGraphX [99][99]

TF1 *	intermediateEFieldFitFunctionX [99]

TF1 *	initialEFieldFitFunctionX

TGraph *	gEFieldGraphY [99][99]

TF1 *	intermediateEFieldFitFunctionY [99]

TF1 *	initialEFieldFitFunctionY

TGraph *	gEFieldGraphZ [99][99]

TF1 *	intermediateEFieldFitFunctionZ [99]

TF1 *	initialEFieldFitFunctionZ

Detailed Description

Definition at line 20 of file SpaceChargeSBND.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 19 of file SpaceChargeSBND.cxx.

 {
     Configure(pset);
 }

spacecharge::SpaceChargeSBND::SpaceChargeSBND ( SpaceChargeSBND const & )

delete

virtual spacecharge::SpaceChargeSBND::~SpaceChargeSBND ( )

virtualdefault

Member Function Documentation

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

Definition at line 24 of file SpaceChargeSBND.cxx.

 {
     fEnableSimSpatialSCE = pset.get<bool>("EnableSimSpatialSCE");
     fEnableSimEfieldSCE = pset.get<bool>("EnableSimEfieldSCE");
     //fEnableCorrSCE = pset.get<bool>("EnableCorrSCE");
     fEnableCalSpatialSCE = pset.get<bool>("EnableCalSpatialSCE");
     fEnableCalEfieldSCE = pset.get<bool>("EnableCalEfieldSCE");
 
     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("SpaceChargeSBND") << "Could not find the space charge effect file '" << fname << "'!\n";
                 }
 
             if(fRepresentationType == "Voxelized_TH3"){
               std::cout << "begin loading voxelized TH3s..." << std::endl;
 
               //Load in histograms
               TH3F* hTrueFwdX = (TH3F*) infile->Get("TrueFwd_Displacement_X");
               TH3F* hTrueFwdY = (TH3F*) infile->Get("TrueFwd_Displacement_Y");
               TH3F* hTrueFwdZ = (TH3F*) infile->Get("TrueFwd_Displacement_Z");
               TH3F* hTrueBkwdX = (TH3F*) infile->Get("TrueBkwd_Displacement_X");
               TH3F* hTrueBkwdY = (TH3F*) infile->Get("TrueBkwd_Displacement_Y");
               TH3F* hTrueBkwdZ = (TH3F*) infile->Get("TrueBkwd_Displacement_Z");
               TH3F* hTrueEFieldX = (TH3F*) infile->Get("True_ElecField_X");
               TH3F* hTrueEFieldY = (TH3F*) infile->Get("True_ElecField_Y");
               TH3F* hTrueEFieldZ = (TH3F*) infile->Get("True_ElecField_Z");
 
               //https://root.cern.ch/doc/master/classTH1.html#a0367fe04ae8709fd4b82795d0a5462c3
               //Set hist directories so they can be referenced elsewhere
               //This needs to be done because they were read in from ext file
               //Note this is not a property of the TH3F, so does't survive copying
               hTrueFwdX->SetDirectory(0);
               hTrueFwdY->SetDirectory(0);
               hTrueFwdZ->SetDirectory(0);
               hTrueBkwdX->SetDirectory(0);
               hTrueBkwdY->SetDirectory(0);
               hTrueBkwdZ->SetDirectory(0);
               hTrueEFieldX->SetDirectory(0);
               hTrueEFieldY->SetDirectory(0);
               hTrueEFieldZ->SetDirectory(0);
 
               //SCEhistograms can be accessed globally in this script
               SCEhistograms = {hTrueFwdX, hTrueFwdY, hTrueFwdZ,
                                hTrueBkwdX, hTrueBkwdY, hTrueBkwdZ,
                                hTrueEFieldX, hTrueEFieldY, hTrueEFieldZ};
 
 
               std::cout << "...finished loading TH3s" << std::endl;
             }else if(fRepresentationType == "Parametric")
                 {
                     for(int i = 0; i < initialSpatialFitPolN[0] + 1; i++)
                         {
                             for(int j = 0; j < intermediateSpatialFitPolN[0] + 1; j++)
                                 {
                                     gSpatialGraphX[i][j] = (TGraph*)infile->Get(Form("deltaX/g%i_%i", i, j));
 
                                 }
                             intermediateSpatialFitFunctionX[i] = new TF1(Form("intermediateSpatialFitFunctionX_%i", i), Form("pol%i", intermediateSpatialFitPolN[0]));
                         }
                     for(int i = 0; i < initialSpatialFitPolN[1] + 1; i++)
                         {
                             for(int j = 0; j < intermediateSpatialFitPolN[1] + 1; j++)
                                 {
                                     gSpatialGraphY[i][j] = (TGraph*)infile->Get(Form("deltaY/g%i_%i", i, j));
 
                                 }
                             intermediateSpatialFitFunctionY[i] = new TF1(Form("intermediateSpatialFitFunctionY_%i", i), Form("pol%i", intermediateSpatialFitPolN[1]));
                         }
                     for(int i = 0; i < initialSpatialFitPolN[2] + 1; i++)
                         {
                             for(int j = 0; j < intermediateSpatialFitPolN[2] + 1; j++)
                                 {
                                     gSpatialGraphZ[i][j] = (TGraph*)infile->Get(Form("deltaZ/g%i_%i", i, j));
 
                                 }
                             intermediateSpatialFitFunctionZ[i] = new TF1(Form("intermediateSpatialFitFunctionZ_%i", i), Form("pol%i", intermediateSpatialFitPolN[2]));
                         }
 
                     initialSpatialFitFunctionX =  new TF1("initialSpatialFitFunctionX", Form("pol%i", initialSpatialFitPolN[0]));
                     initialSpatialFitFunctionY =  new TF1("initialSpatialFitFunctionY", Form("pol%i", initialSpatialFitPolN[1]));
                     initialSpatialFitFunctionZ =  new TF1("initialSpatialFitFunctionZ", Form("pol%i", initialSpatialFitPolN[2]));
 
                     for(int i = 0; i < initialEFieldFitPolN[0] + 1; i++)
                         {
                             for(int j = 0; j < intermediateEFieldFitPolN[0] + 1; j++)
                                 {
                                     gEFieldGraphX[i][j] = (TGraph*)infile->Get(Form("deltaEx/g%i_%i", i, j));
 
                                 }
                             intermediateEFieldFitFunctionX[i] = new TF1(Form("intermediateEFieldFitFunctionX_%i", i), Form("pol%i", intermediateEFieldFitPolN[0]));
                         }
                     for(int i = 0; i < initialEFieldFitPolN[1] + 1; i++)
                         {
                             for(int j = 0; j < intermediateEFieldFitPolN[1] + 1; j++)
                                 {
                                     gEFieldGraphY[i][j] = (TGraph*)infile->Get(Form("deltaEy/g%i_%i", i, j));
 
                                 }
                             intermediateEFieldFitFunctionY[i] = new TF1(Form("intermediateEFieldFitFunctionY_%i", i), Form("pol%i", intermediateEFieldFitPolN[1]));
                         }
                     for(int i = 0; i < initialEFieldFitPolN[2] + 1; i++)
                         {
                             for(int j = 0; j < intermediateEFieldFitPolN[2] + 1; j++)
                                 {
                                     gEFieldGraphZ[i][j] = (TGraph*)infile->Get(Form("deltaEz/g%i_%i", i, j));
 
                                 }
                             intermediateEFieldFitFunctionZ[i] = new TF1(Form("intermediateEFieldFitFunctionZ_%i", i), Form("pol%i", intermediateEFieldFitPolN[2]));
                         }
 
                     initialEFieldFitFunctionX =  new TF1("initialEFieldFitFunctionX", Form("pol%i", initialEFieldFitPolN[0]));
                     initialEFieldFitFunctionY =  new TF1("initialEFieldFitFunctionY", Form("pol%i", initialEFieldFitPolN[1]));
                     initialEFieldFitFunctionZ =  new TF1("initialEFieldFitFunctionZ", Form("pol%i", initialEFieldFitPolN[2]));
                 }else{
                   std::cout << "fRepresentationType not known!!!" << std::endl;
                 }
             infile->Close();
         }
 
     if(fEnableCorrSCE == true)
         {
             // Grab other parameters from pset
         }
     return true;
 }

bool spacecharge::SpaceChargeSBND::EnableCalEfieldSCE ( ) const

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 195 of file SpaceChargeSBND.cxx.

 {
   return fEnableCalEfieldSCE;
 }

bool spacecharge::SpaceChargeSBND::EnableCalSpatialSCE ( ) const

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 189 of file SpaceChargeSBND.cxx.

 {
   return fEnableCalSpatialSCE;
 }

bool spacecharge::SpaceChargeSBND::EnableCorrSCE ( ) const

inlineoverridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 36 of file SpaceChargeSBND.h.

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

spacecharge::SpaceChargeSBND::EnableCalEfieldSCE

bool EnableCalEfieldSCE() const override

Definition: SpaceChargeSBND.cxx:195

spacecharge::SpaceChargeSBND::EnableCalSpatialSCE

bool EnableCalSpatialSCE() const override

Definition: SpaceChargeSBND.cxx:189

bool spacecharge::SpaceChargeSBND::EnableSimEfieldSCE ( ) const

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 177 of file SpaceChargeSBND.cxx.

 {
     return fEnableSimEfieldSCE;
 }

bool spacecharge::SpaceChargeSBND::EnableSimSpatialSCE ( ) const

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 171 of file SpaceChargeSBND.cxx.

 {
     return fEnableSimSpatialSCE;
 }

geo::Vector_t spacecharge::SpaceChargeSBND::GetCalEfieldOffsets	(	geo::Point_t const &	point,
		int const &	TPCid = `1`
	)		const

inlineoverridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 41 of file SpaceChargeSBND.h.

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

geo::Vector_t spacecharge::SpaceChargeSBND::GetCalPosOffsets	(	geo::Point_t const &	point,
		int const &	TPCid = `1`
	)		const

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 240 of file SpaceChargeSBND.cxx.

 {
   std::vector<double> theCalPosOffsets;
   double xx=point.X(), yy=point.Y(), zz=point.Z();
 
   if(fRepresentationType == "Voxelized_TH3"){
     //handle OOAV by projecting edge cases
     if(xx<-199.999){xx=-199.999;}
     else if(xx>199.999){xx=199.999;}
     if(yy<-199.999){yy=-199.999;}
     else if(yy>199.999){yy=199.999;}
     if(zz<0.001){zz=0.001;}
     else if(zz>499.999){zz=499.999;}
     //correct for charge drifted across cathode
     if ((TPCid == 0) and (xx > -2.5)) { xx = -2.5; }
     if ((TPCid == 1) and (xx < 2.5)) { xx = 2.5; }
     double offset_x=0., offset_y=0., offset_z=0.;
     offset_x = SCEhistograms.at(3)->Interpolate(xx,yy,zz);
     offset_y = SCEhistograms.at(4)->Interpolate(xx,yy,zz);
     offset_z = SCEhistograms.at(5)->Interpolate(xx,yy,zz);
     theCalPosOffsets = {offset_x, offset_y, offset_z};
     
   }else if(fRepresentationType == "Parametric"){     
     //this is not supported for parametric
     std::cout << "Change Representation Type to Voxelized TH3 if you want to use the backward offset function" << std::endl;
     theCalPosOffsets.resize(3, 0.0);
   }
   
   return { theCalPosOffsets[0], theCalPosOffsets[1], theCalPosOffsets[2] };
 }

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

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 383 of file SpaceChargeSBND.cxx.

 {
     std::vector<double> theEfieldOffsets;
     double xx=point.X(), yy=point.Y(), zz=point.Z();
     double offset_x=0., offset_y=0., offset_z=0.;
 
     if(fRepresentationType == "Voxelized_TH3"){
       //handle OOAV by projecting edge cases
       if(xx<-199.999){xx=-199.999;}
       else if(xx>199.999){xx=199.999;}
       if(yy<-199.999){yy=-199.999;}
       else if(yy>199.999){yy=199.999;}
       if(zz<0.001){zz=0.001;}
       else if(zz>499.999){zz=499.999;}
       offset_x = SCEhistograms.at(6)->Interpolate(xx, yy, zz);
       offset_y = SCEhistograms.at(7)->Interpolate(xx, yy, zz);
       offset_z = SCEhistograms.at(8)->Interpolate(xx, yy, zz);
 
       theEfieldOffsets = {offset_x, offset_y, offset_z};
       
     }else if(fRepresentationType == "Parametric"){
 
       if(IsInsideBoundaries(point.X(), point.Y(), point.Z()) == false){
         theEfieldOffsets.resize(3, 0.0);
       }
       else
         {
           theEfieldOffsets = GetEfieldOffsetsParametric(point.X(), point.Y(), point.Z());
 
           // GetOneEfieldOffsetParametric returns V/m
           // The E-field offsets are returned as -dEx/|E_nominal|, -dEy/|E_nominal|, and -dEz/|E_nominal| where |E_nominal| is DriftField
           theEfieldOffsets[0] = -1.0 * theEfieldOffsets[0] / (100.0 * DriftField);
           theEfieldOffsets[1] = -1.0 * theEfieldOffsets[1] / (100.0 * DriftField);
           theEfieldOffsets[2] = -1.0 * theEfieldOffsets[2] / (100.0 * DriftField);
         }
     }
 
     return { theEfieldOffsets[0], theEfieldOffsets[1], theEfieldOffsets[2] };
 }

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

protected

Definition at line 424 of file SpaceChargeSBND.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::SpaceChargeSBND::GetOneEfieldOffsetParametric	(	double	xVal,
		double	yVal,
		double	zVal,
		std::string	axis
	)		const

protected

Definition at line 440 of file SpaceChargeSBND.cxx.

 {
     double parA[99][99];
     double parB[99];
 
     for(int i = 0; i < 99; i++)
         {
             for(int j = 0; j < 99; j++)
                 {
                     parA[i][j] = 0.0;
                 }
             parB[i] = 0.0;
         }
     if(axis == "X")
         {
             for(int i = 0; i < initialEFieldFitPolN[0] + 1; i++)
                 {
                     for(int j = 0; j < intermediateEFieldFitPolN[0] + 1; j++)
                         {
                             parA[i][j] = gEFieldGraphX[i][j]->Eval(zValNew);
                         }
                     intermediateEFieldFitFunctionX[i]->SetParameters(parA[i]);
                 }
         }
     else if(axis == "Y")
         {
             for(int i = 0; i < initialEFieldFitPolN[1] + 1; i++)
                 {
                     for(int j = 0; j < intermediateEFieldFitPolN[1] + 1; j++)
                         {
                             parA[i][j] = gEFieldGraphY[i][j]->Eval(zValNew);
                         }
                     intermediateEFieldFitFunctionY[i]->SetParameters(parA[i]);
                 }
         }
     else if(axis == "Z")
         {
             for(int i = 0; i < initialEFieldFitPolN[2] + 1; i++)
                 {
                     for(int j = 0; j < intermediateEFieldFitPolN[2] + 1; j++)
                         {
                             parA[i][j] = gEFieldGraphZ[i][j]->Eval(zValNew);
                         }
                     intermediateEFieldFitFunctionZ[i]->SetParameters(parA[i]);
                 }
         }
 
     double aValNew;
     double bValNew;
 
     if(axis == "Y")
         {
             aValNew = xValNew;
             bValNew = yValNew;
         }
     else
         {
             aValNew = yValNew;
             bValNew = xValNew;
         }
 
     double offsetValNew = 0.0;
     if(axis == "X")
         {
             for(int i = 0; i < initialEFieldFitPolN[0] + 1; i++)
                 {
                     parB[i] = intermediateEFieldFitFunctionX[i]->Eval(aValNew);
                 }
             initialEFieldFitFunctionX->SetParameters(parB);
             offsetValNew = initialEFieldFitFunctionX->Eval(bValNew);
         }
     else if(axis == "Y")
         {
             for(int i = 0; i < initialEFieldFitPolN[1] + 1; i++)
                 {
                     parB[i] = intermediateEFieldFitFunctionY[i]->Eval(aValNew);
                 }
             initialEFieldFitFunctionY->SetParameters(parB);
             offsetValNew = initialEFieldFitFunctionY->Eval(bValNew);
         }
     else if(axis == "Z")
         {
             for(int i = 0; i < initialEFieldFitPolN[2] + 1; i++)
                 {
                     parB[i] = intermediateEFieldFitFunctionZ[i]->Eval(aValNew);
                 }
             initialEFieldFitFunctionZ->SetParameters(parB);
             offsetValNew = initialEFieldFitFunctionZ->Eval(bValNew);
         }
 
     return offsetValNew;
 }

double spacecharge::SpaceChargeSBND::GetOnePosOffsetParametric	(	double	xVal,
		double	yVal,
		double	zVal,
		std::string	axis
	)		const

protected

Definition at line 290 of file SpaceChargeSBND.cxx.

 {
     double parA[99][99];
     double parB[99];
 
     for(int i = 0; i < 99; i++)
         {
             for(int j = 0; j < 99; j++)
                 {
                     parA[i][j] = 0.0;
                 }
             parB[i] = 0.0;
         }
     if(axis == "X")
         {
             for(int i = 0; i < initialSpatialFitPolN[0] + 1; i++)
                 {
                     for(int j = 0; j < intermediateSpatialFitPolN[0] + 1; j++)
                         {
                             parA[i][j] = gSpatialGraphX[i][j]->Eval(zValNew);
                         }
                     intermediateSpatialFitFunctionX[i]->SetParameters(parA[i]);
                 }
         }
     else if(axis == "Y")
         {
             for(int i = 0; i < initialSpatialFitPolN[1] + 1; i++)
                 {
                     for(int j = 0; j < intermediateSpatialFitPolN[1] + 1; j++)
                         {
                             parA[i][j] = gSpatialGraphY[i][j]->Eval(zValNew);
                         }
                     intermediateSpatialFitFunctionY[i]->SetParameters(parA[i]);
                 }
         }
     else if(axis == "Z")
         {
             for(int i = 0; i < initialSpatialFitPolN[2] + 1; i++)
                 {
                     for(int j = 0; j < intermediateSpatialFitPolN[2] + 1; j++)
                         {
                             parA[i][j] = gSpatialGraphZ[i][j]->Eval(zValNew);
                         }
                     intermediateSpatialFitFunctionZ[i]->SetParameters(parA[i]);
                 }
         }
 
     double aValNew;
     double bValNew;
 
     if(axis == "Y")
         {
             aValNew = xValNew;
             bValNew = yValNew;
         }
     else
         {
             aValNew = yValNew;
             bValNew = xValNew;
         }
     double offsetValNew = 0.0;
     if(axis == "X")
         {
             for(int i = 0; i < initialSpatialFitPolN[0] + 1; i++)
                 {
                     parB[i] = intermediateSpatialFitFunctionX[i]->Eval(aValNew);
                 }
             initialSpatialFitFunctionX->SetParameters(parB);
             offsetValNew = initialSpatialFitFunctionX->Eval(bValNew);
         }
     else if(axis == "Y")
         {
             for(int i = 0; i < initialSpatialFitPolN[1] + 1; i++)
                 {
                     parB[i] = intermediateSpatialFitFunctionY[i]->Eval(aValNew);
                 }
             initialSpatialFitFunctionY->SetParameters(parB);
             offsetValNew = initialSpatialFitFunctionY->Eval(bValNew);
         }
     else if(axis == "Z")
         {
             for(int i = 0; i < initialSpatialFitPolN[2] + 1; i++)
                 {
                     parB[i] = intermediateSpatialFitFunctionZ[i]->Eval(aValNew);
                 }
             initialSpatialFitFunctionZ->SetParameters(parB);
             offsetValNew = initialSpatialFitFunctionZ->Eval(bValNew);
         }
 
     return offsetValNew;
 }

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

overridevirtual

Implements spacecharge::SpaceCharge.

Definition at line 201 of file SpaceChargeSBND.cxx.

 {
     std::vector<double> thePosOffsets;
     double xx=point.X(), yy=point.Y(), zz=point.Z();
 
     if(fRepresentationType == "Voxelized_TH3"){
       //handle OOAV by projecting edge cases
       if(xx<-199.999){xx=-199.999;}
       else if(xx>199.999){xx=199.999;}
       if(yy<-199.999){yy=-199.999;}
       else if(yy>199.999){yy=199.999;}
       if(zz<0.001){zz=0.001;}
       else if(zz>499.999){zz=499.999;}
       //larsim requires negative sign in TPC 0
       int corr = 1;
       if (xx < 0) { corr = -1; }
       double offset_x=0., offset_y=0., offset_z=0.;
       offset_x = corr*SCEhistograms.at(0)->Interpolate(xx,yy,zz);
       offset_y = SCEhistograms.at(1)->Interpolate(xx,yy,zz);
       offset_z = SCEhistograms.at(2)->Interpolate(xx,yy,zz);
       thePosOffsets = {offset_x, offset_y, offset_z};
 
     }else if(fRepresentationType == "Parametric"){
       if(IsInsideBoundaries(point.X(), point.Y(), point.Z()) == false){
         thePosOffsets.resize(3, 0.0);
       }else{
         // GetPosOffsetsParametric returns m; the PosOffsets should be in cm
         thePosOffsets = GetPosOffsetsParametric(xx, yy, zz);
         for(int i=0; i<3; i++){thePosOffsets[i]=100.*thePosOffsets[i];}
       }
     }
     else{
       thePosOffsets.resize(3, 0.0);
     }
 
     return { thePosOffsets[0], thePosOffsets[1], thePosOffsets[2] };
 }

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

protected

Definition at line 274 of file SpaceChargeSBND.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::SpaceChargeSBND::IsInsideBoundaries	(	double	xVal,
		double	yVal,
		double	zVal
	)		const

protected

Definition at line 565 of file SpaceChargeSBND.cxx.

 {
     bool isInside = true;
 
     if((xVal < -196.5) || (xVal > 196.5) || (yVal < -200.0) || (yVal > 200.0) || (zVal < 0) || (zVal > 500.0))
         {
             isInside = false;
         }
 
     return isInside;
 }

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

protected

Definition at line 535 of file SpaceChargeSBND.cxx.

 {
     xVal = xVal / 100.0;
 
     // We use the same map twice; 0 is cathod
     // Map [0, 196.5] to [0, 2.0]
     double xValNew = (2.0 / 1.965) * fabs(xVal);
 
     return xValNew;
 }

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

protected

Definition at line 548 of file SpaceChargeSBND.cxx.

 {
     yVal = yVal / 100.0;
 
     return (yVal + 2.0);
 }

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

protected

Definition at line 557 of file SpaceChargeSBND.cxx.

 {
     return (zVal / 100.0);
 }

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

Definition at line 160 of file SpaceChargeSBND.cxx.

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

Member Data Documentation

double spacecharge::SpaceChargeSBND::DriftField = 500.0

protected

Definition at line 50 of file SpaceChargeSBND.h.

bool spacecharge::SpaceChargeSBND::fEnableCalEfieldSCE

protected

Definition at line 55 of file SpaceChargeSBND.h.

bool spacecharge::SpaceChargeSBND::fEnableCalSpatialSCE

protected

Definition at line 54 of file SpaceChargeSBND.h.

bool spacecharge::SpaceChargeSBND::fEnableCorrSCE

protected

Definition at line 56 of file SpaceChargeSBND.h.

bool spacecharge::SpaceChargeSBND::fEnableSimEfieldSCE

protected

Definition at line 53 of file SpaceChargeSBND.h.

bool spacecharge::SpaceChargeSBND::fEnableSimSpatialSCE

protected

Definition at line 52 of file SpaceChargeSBND.h.

std::string spacecharge::SpaceChargeSBND::fInputFilename

protected

Definition at line 59 of file SpaceChargeSBND.h.

std::string spacecharge::SpaceChargeSBND::fRepresentationType

protected

Definition at line 58 of file SpaceChargeSBND.h.

TGraph* spacecharge::SpaceChargeSBND::gEFieldGraphX[99][99]

protected

Definition at line 83 of file SpaceChargeSBND.h.

TGraph* spacecharge::SpaceChargeSBND::gEFieldGraphY[99][99]

protected

Definition at line 86 of file SpaceChargeSBND.h.

TGraph* spacecharge::SpaceChargeSBND::gEFieldGraphZ[99][99]

protected

Definition at line 89 of file SpaceChargeSBND.h.

TGraph* spacecharge::SpaceChargeSBND::gSpatialGraphX[99][99]

protected

Definition at line 73 of file SpaceChargeSBND.h.

TGraph* spacecharge::SpaceChargeSBND::gSpatialGraphY[99][99]

protected

Definition at line 76 of file SpaceChargeSBND.h.

TGraph* spacecharge::SpaceChargeSBND::gSpatialGraphZ[99][99]

protected

Definition at line 79 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::initialEFieldFitFunctionX

protected

Definition at line 85 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::initialEFieldFitFunctionY

protected

Definition at line 88 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::initialEFieldFitFunctionZ

protected

Definition at line 91 of file SpaceChargeSBND.h.

int spacecharge::SpaceChargeSBND::initialEFieldFitPolN[3] = {3, 3, 3}

protected

Definition at line 48 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::initialSpatialFitFunctionX

protected

Definition at line 75 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::initialSpatialFitFunctionY

protected

Definition at line 78 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::initialSpatialFitFunctionZ

protected

Definition at line 81 of file SpaceChargeSBND.h.

int spacecharge::SpaceChargeSBND::initialSpatialFitPolN[3] = {3, 4, 3}

protected

Definition at line 46 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::intermediateEFieldFitFunctionX[99]

protected

Definition at line 84 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::intermediateEFieldFitFunctionY[99]

protected

Definition at line 87 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::intermediateEFieldFitFunctionZ[99]

protected

Definition at line 90 of file SpaceChargeSBND.h.

int spacecharge::SpaceChargeSBND::intermediateEFieldFitPolN[3] = {6, 4, 4}

protected

Definition at line 49 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::intermediateSpatialFitFunctionX[99]

protected

Definition at line 74 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::intermediateSpatialFitFunctionY[99]

protected

Definition at line 77 of file SpaceChargeSBND.h.

TF1* spacecharge::SpaceChargeSBND::intermediateSpatialFitFunctionZ[99]

protected

Definition at line 80 of file SpaceChargeSBND.h.

int spacecharge::SpaceChargeSBND::intermediateSpatialFitPolN[3] = {4, 4, 4}

protected

Definition at line 47 of file SpaceChargeSBND.h.

std::vector<TH3F*> spacecharge::SpaceChargeSBND::SCEhistograms = std::vector<TH3F*>(9)

protected

Definition at line 71 of file SpaceChargeSBND.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation