exampleUtilCdc.cpp

Go to the documentation of this file.

/*
 * exampleCdc.cpp
 *
 *  Created on: 23 Mar 2023
 *      Author: kleinwrt
 */
 
#include "exampleUtilCdc.h"
 
using namespace Eigen;
 
namespace gbl {
 
 
GblDetectorLayer CreateWireCdc(const std::string aName, unsigned int layer,
                double xPos, double yPos, double zPos, double phi, double tanLambda,
                double stereoAngle, double uRes) {
        Vector3d aCenter(xPos, yPos, zPos);
        Vector2d aResolution(uRes, 0.);
        Vector2d aPrecision(1. / (uRes * uRes), 0.);
// track direction (N)
        const double cosLambda = 1. / sqrt(1. + tanLambda * tanLambda);
        Vector3d tDir(cos(phi) * cosLambda, sin(phi) * cosLambda,
                        tanLambda * cosLambda);
// wire direction (V)
        const double rad = sqrt(xPos * xPos + yPos * yPos);
        const double scale = sqrt(1. + stereoAngle * stereoAngle);
        Vector3d vDir(-yPos * stereoAngle / rad / scale,
                        xPos * stereoAngle / rad / scale, 1. / scale);
// measurement direction perpendicular to track and wire direction (U)
        Vector3d uDir = vDir.cross(tDir).normalized();
// normal to measurement plane
        Vector3d nDir = uDir.cross(vDir).normalized();
 
        Matrix3d measTrafo;
        measTrafo << uDir[0], uDir[1], uDir[2], vDir[0], vDir[1], vDir[2], nDir[0], nDir[1], nDir[2]; // U,V,N
        Matrix3d alignTrafo = Matrix3d::Identity();                      // I,J,K
        return GblDetectorLayer(aName, layer, 1, 0.0, aCenter, aResolution,
                        aPrecision, measTrafo, alignTrafo); // no mult. scat.
}
 
 
GblDetectorLayer CreateImpactPar(const std::string aName, unsigned int layer,
                double xPos, double yPos, double zPos, double phi, double tanLambda,
                double xRes, double yRes, double zRes) {
// track direction (N)
        const double cosLambda = 1. / sqrt(1. + tanLambda * tanLambda);
        Vector3d tDir(cos(phi) * cosLambda, sin(phi) * cosLambda,
                        tanLambda * cosLambda);
// beam direction (V) = Z
        Vector3d vDir(0., 0., 1.);
// measurement direction perpendicular to track and beam direction (U)
        Vector3d uDir = vDir.cross(tDir).normalized();
// normal to measurement plane
        Vector3d nDir = uDir.cross(vDir).normalized();
// beam variance in U direction
        const double uVar = uDir[0] * xRes * xRes * uDir[0]
                        + uDir[1] * yRes * yRes * uDir[1];
        Vector3d aCenter(xPos, yPos, zPos);
        Vector2d aResolution(sqrt(uVar), zRes);
        Vector2d aPrecision(1. / uVar, 1. / (zRes * zRes));
 
        Matrix3d measTrafo;
        measTrafo << uDir[0], uDir[1], uDir[2], vDir[0], vDir[1], vDir[2], nDir[0], nDir[1], nDir[2]; // U,V,N
        Matrix3d alignTrafo = Matrix3d::Identity();                      // I,J,K
        return GblDetectorLayer(aName, layer, 2, 0., aCenter, aResolution,
                        aPrecision, measTrafo, alignTrafo);
}
 
}