exampleSit.cpp File Reference

Example silicon tracker application. More...

#include <time.h>
#include "exampleSit.h"
#include "GblTrajectory.h"

Go to the source code of this file.

Namespaces
namespace	gbl
	Namespace for the general broken lines package.

Functions
void	exampleSit ()
	Silicon tracker example. More...
GblDetectorLayer	gbl::CreateLayerSit (const std::string aName, unsigned int layer, double xPos, double yPos, double zPos, double thickness, double uAngle, double uRes)
	Create a silicon layer with 1D measurement. More...
GblDetectorLayer	gbl::CreateLayerSit (const std::string aName, unsigned int layer, double xPos, double yPos, double zPos, double thickness, double uAngle, double uRes, double vAngle, double vRes)
	Create a silicon layer with 2D measurement. More...

Detailed Description

Example silicon tracker application.

Author

Claus Kleinwort, DESY, 2018 (Claus.nosp@m..Kle.nosp@m.inwor.nosp@m.t@de.nosp@m.sy.de)

Copyright

Copyright (c) 2018-2021 Deutsches Elektronen-Synchroton, Member of the Helmholtz Association, (DESY), HAMBURG, GERMANY

This library is free software; you can redistribute it and/or modify it under the terms of the GNU Library General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.

This library is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public License for more details.

You should have received a copy of the GNU Library General Public License along with this program (see the file COPYING.LIB for more details); if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.

Definition in file exampleSit.cpp.

Function Documentation

exampleSit()

void exampleSit

(

)

Silicon tracker example.

Simulate and reconstruct helical tracks in silicon pixel and (1D or 2D) strip detectors.

Create points on initial trajectory, create trajectory from points, fit and write trajectory to MP-II binary file (for rigid body alignment).

Setup:

• Beam (mainly) in X direction
• Constant magnetic field in Z direction
• Silicon sensors measuring in YZ plane, orthogonal (pixel) or non-orthogonal (stereo strips) measurement systems
• Multiple scattering in sensors (air in between ignored)
• Curvilinear system (T,U,V) as local coordinate system and (q/p, slopes, offsets) as local track parameters

Remarks

To exercise (mis)alignment different sets of layers (with different geometry) for simulation and reconstruction can be used.

Example steering file for Millepede-II (B=0):

Cfiles
milleBinaryISN.dat
 
method inversion 3 0.1
chiscut 30. 6.
printcounts
! fix first pixel and last stereo layer as reference
parameter
  1  0.  -1.
  2  0.  -1.
  3  0.  -1.
  4  0.  -1.
  5  0.  -1.
  6  0.  -1.
 61  0.  -1.
 62  0.  -1.
 63  0.  -1.
 64  0.  -1.
 65  0.  -1.
 66  0.  -1.
end

Definition at line 79 of file exampleSit.cpp.

References gbl::GblPoint::addGlobals(), gbl::GblPoint::addMeasurement(), gbl::GblPoint::addScatterer(), gbl::CreateLayerSit(), gbl::GblTrajectory::fit(), gbl::gblMultipleScatteringError(), gbl::gblSimpleJacobian(), gbl::GblHelixPrediction::getArcLength(), gbl::GblHelixPrediction::getCosIncidence(), gbl::GblHelixPrediction::getCurvilinearDirs(), gbl::GblHelixPrediction::getDirection(), gbl::GblDetectorLayer::getLayerID(), gbl::GblHelixPrediction::getMeasPred(), gbl::GblDetectorLayer::getMeasSystemDirs(), gbl::GblHelixPrediction::getPosition(), gbl::GblDetectorLayer::getPrecision(), gbl::GblDetectorLayer::getRadiationLength(), gbl::GblDetectorLayer::getRigidBodyDerLocal(), gbl::GblDetectorLayer::intersectWithHelix(), gbl::GblTrajectory::milleOut(), gbl::GblSimpleHelix::moveToXY(), gbl::unif(), and gbl::unrm().