structure including all variables of interest with uncertainties from boot-strap More...

Collaboration diagram for UnknownPars:

[legend]

Public Member Functions
void	add (SpotParam sPar, double SizeX, double SizeY, double SizeXY, double SizeZ)
	add next boot-strap replica of the BS parameters

void	printStat ()
	Print interesting statistics from boot-strap.

void	getOutput (std::vector< VectorXd > &vtxPos, std::vector< MatrixXd > &vtxErr, MatrixXd &sizeMat)
	get output in Belle2-like format

void	setBranchVal (TTree T, std::vector< double > vec, TString n)
	save bootstrap variable to TTree

void	setBranchSpline (TTree T, Spline spl, TString n)
	save bootstrap spline to TTree

void	save2tree (TString fName)
	save everything to TTree

Public Attributes
UnknowSpline	x
	BS position (x coordinate)

UnknowSpline	y
	BS position (y coordinate)

UnknowSpline	kX
	BS angle in xz plane.

UnknowSpline	kY
	BS angle in yz plane.

UnknowSpline	z
	BS position (z coordinate)

UnknowVar	sizeX
	BS size in x direction.

UnknowVar	sizeY
	BS size in y direction.

UnknowVar	sizeXY
	off-diagonal component of BS size cov matrix in frame where z' is aligned with z

UnknowVar	sizeMin
	smallest eigenvalue of the BS size cov matrix (similar to sizeY)

UnknowVar	sizeMax
	middle eigenvalue of the BS size cov matrix (similar to sizeX)

UnknowVar	xyAngle
	angle of the BS in xy plane when z' is aligned with z

UnknowVar	sizeZ
	BS size in z direction.

UnknowVar	crossAngle
	derived value of the crossing angle of the HER & LER beams

UnknowVar	matXX
	XX element of BS size cov matrix.

UnknowVar	matYY
	YY element of BS size cov matrix.

UnknowVar	matZZ
	ZZ element of BS size cov matrix.

UnknowVar	matXY
	XY element of BS size cov matrix.

UnknowVar	matXZ
	XZ element of BS size cov matrix.

UnknowVar	matYZ
	YZ element of BS size cov matrix.

Detailed Description

structure including all variables of interest with uncertainties from boot-strap

Definition at line 333 of file BeamSpotStandAlone.cc.

Member Function Documentation

◆ add()

void add	(	SpotParam	sPar,
		double	SizeX,
		double	SizeY,
		double	SizeXY,
		double	SizeZ
	)

inline

add next boot-strap replica of the BS parameters

Definition at line 359 of file BeamSpotStandAlone.cc.

    {
      x.add(sPar.x);
      y.add(sPar.y);
      kX.add(sPar.kX);
      kY.add(sPar.kY);
      z.add(sPar.z);
 
 
      sizeX.add(SizeX);
      sizeY.add(SizeY);
      sizeXY.add(SizeXY);
      sizeZ.add(SizeZ);
 
 
      // Calculate the eigen-values
      double SizeMin, SizeMax;
      std::tie(SizeMin, SizeMax) = getSizeMinMax(SizeX, SizeY, SizeXY);
 
      sizeMin.add(SizeMin);
      sizeMax.add(SizeMax);
 
      // and angle in mrad
      double angle = 1e3 * getAngle(SizeX, SizeY, SizeXY);
 
 
      xyAngle.add(angle);
 
      //Get whole cov matrix
      MatrixXd matSize = getRotatedSizeMatrix({sqrS(SizeX), sqrS(SizeY), sqrS(SizeXY)}, sqrS(SizeZ), sPar.kX.val(sPar.kX.center()),
                                              sPar.kY.val(sPar.kY.center()));
 
      // Store elements in [um]
      matXX.add(sqrtS(matSize(0, 0)));
      matYY.add(sqrtS(matSize(1, 1)));
      matZZ.add(sqrtS(matSize(2, 2)));
      matXY.add(sqrtS(matSize(0, 1)));
 
      matXZ.add(sqrtS(matSize(0, 2)));
      matYZ.add(sqrtS(matSize(1, 2)));
 
 
      // crossing-angle in mrad
      double crossAngleVal = 1e3 * 2 * sqrtS(matSize(0, 0)) / sqrtS(matSize(2, 2));
      crossAngle.add(crossAngleVal);
    }

◆ getOutput()

void getOutput	(	std::vector< VectorXd > &	vtxPos,
		std::vector< MatrixXd > &	vtxErr,
		MatrixXd &	sizeMat
	)

inline

get output in Belle2-like format

Definition at line 436 of file BeamSpotStandAlone.cc.

    {
      //Store the vertex position
      int nVals = x.spls[0].vals.size();
 
      vtxPos.clear();
      vtxErr.clear();
 
      const double toCm = 1e-4;
 
      for (int i = 0; i < nVals; ++i) {
        //vertex position
        Vector3d vtx(x.spls[0].vals[i]*toCm, y.spls[0].vals[i]*toCm, z.spls[0].vals[i]*toCm);
 
        //vertex error matrix (symmetric)
        Matrix3d mS = Matrix3d::Zero();
        mS(0, 0) = sqrS(x.spls[0].errs[i] * toCm);
        mS(1, 1) = sqrS(y.spls[0].errs[i] * toCm);
        mS(2, 2) = sqrS(z.spls[0].errs[i] * toCm);
 
        vtxPos.push_back(vtx);
        vtxErr.push_back(mS);
      }
 
      //BeamSpot size matrix (from boot-strap iteration 0)
 
      sizeMat.resize(3, 3);
      sizeMat(0, 0) = sqrS(matXX.vars[0] * toCm);
      sizeMat(1, 1) = sqrS(matYY.vars[0] * toCm);
      sizeMat(2, 2) = sqrS(matZZ.vars[0] * toCm);
 
      sizeMat(0, 1) = sqrS(matXY.vars[0] * toCm);
      sizeMat(0, 2) = sqrS(matXZ.vars[0] * toCm);
      sizeMat(1, 2) = sqrS(matYZ.vars[0] * toCm);
 
      sizeMat(1, 0) = sizeMat(0, 1);
      sizeMat(2, 0) = sizeMat(0, 2);
      sizeMat(2, 1) = sizeMat(1, 2);
    }

◆ printStat()

void printStat ( )

inline

Print interesting statistics from boot-strap.

Definition at line 407 of file BeamSpotStandAlone.cc.

    {
      x.getMeanSigma().print("x");
      y.getMeanSigma().print("y");
      kX.getMeanSigma().print("kX");
      kY.getMeanSigma().print("kY");
      z.getMeanSigma().print("z");
 
      sizeX.printStat("sizeX");
      sizeY.printStat("sizeY");
      sizeXY.printStat("sizeXY");
 
 
      sizeMin.printStat("sizeMin");
      sizeMax.printStat("sizeMax");
      xyAngle.printStat("xyAngle");
      sizeZ.printStat("sizeZ");
      crossAngle.printStat("crossAngle");
 
 
      matXX.printStat("matXX");
      matYY.printStat("matYY");
      matZZ.printStat("matZZ");
      matXY.printStat("matXY");
      matXZ.printStat("matXZ");
      matYZ.printStat("matYZ");
    }

◆ save2tree()

void save2tree ( TString fName )

inline

save everything to TTree

Definition at line 496 of file BeamSpotStandAlone.cc.

    {
 
      TTree* T = new TTree("runs", "beam conditions of runs");
 
      int run = -99; //currently dummy
      T->Branch("run", &run, "run/I");
 
      Spline xAvg = x.getMeanSigma();
      setBranchSpline(T, &xAvg, "x");
      Spline yAvg = y.getMeanSigma();
      setBranchSpline(T, &yAvg, "y");
      Spline zAvg = z.getMeanSigma();
      setBranchSpline(T, &zAvg, "z");
 
      Spline kxAvg = kX.getMeanSigma();
      setBranchSpline(T, &kxAvg, "kX");
      Spline kyAvg = kY.getMeanSigma();
      setBranchSpline(T, &kyAvg, "kY");
 
      std::vector<double> sizeXVar = sizeX.getStats();
      setBranchVal(T, &sizeXVar, "sizeX");
      std::vector<double> sizeYVar = sizeY.getStats();
      setBranchVal(T, &sizeYVar, "sizeY");
      std::vector<double> sizeXYVar = sizeXY.getStats();
      setBranchVal(T, &sizeXYVar, "sizeXY");
      std::vector<double> sizeZVar = sizeZ.getStats();
      setBranchVal(T, &sizeZVar, "sizeZ");
 
      std::vector<double> sizeMinVar = sizeMin.getStats();
      setBranchVal(T, &sizeMinVar, "sizeMin");
      std::vector<double> sizeMaxVar = sizeMax.getStats();
      setBranchVal(T, &sizeMaxVar, "sizeMax");
      std::vector<double> xyAngleVar = xyAngle.getStats();
      setBranchVal(T, &xyAngleVar, "xyAngle");
 
      std::vector<double> crossAngleVar = crossAngle.getStats();
      setBranchVal(T, &crossAngleVar, "crossAngle");
 
 
      std::vector<double> matXXVar = matXX.getStats();
      std::vector<double> matYYVar = matYY.getStats();
      std::vector<double> matZZVar = matZZ.getStats();
      std::vector<double> matXYVar = matXY.getStats();
      std::vector<double> matXZVar = matXZ.getStats();
      std::vector<double> matYZVar = matYZ.getStats();
 
      setBranchVal(T, &matXXVar, "matXX");
      setBranchVal(T, &matYYVar, "matYY");
      setBranchVal(T, &matZZVar, "matZZ");
 
      setBranchVal(T, &matXYVar, "matXY");
      setBranchVal(T, &matXZVar, "matXZ");
      setBranchVal(T, &matYZVar, "matYZ");
 
 
      T->Fill();
      T->SaveAs(fName);
    }

◆ setBranchSpline()

void setBranchSpline	(	TTree *	T,
		Spline *	spl,
		TString	n
	)

inline

save bootstrap spline to TTree

Definition at line 486 of file BeamSpotStandAlone.cc.

    {
      T->Branch(n + "_nodes", &spl->nodes);
      T->Branch(n + "_vals",  &spl->vals);
      T->Branch(n + "_errs",  &spl->errs);
    }

◆ setBranchVal()

void setBranchVal	(	TTree *	T,
		std::vector< double > *	vec,
		TString	n
	)

inline

save bootstrap variable to TTree

Definition at line 478 of file BeamSpotStandAlone.cc.

    {
      T->Branch(n, &vec->at(0), n + "/D");
      T->Branch(n + "_low", &vec->at(1), n + "_low/D");
      T->Branch(n + "_high", &vec->at(2), n + "_high/D");
    }