|
|
|
#include <LASBarrelAlgorithm.h>
Public Member Functions | |
| LASBarrelAlignmentParameterSet | CalculateParameters (LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &) |
| void | Dump (void) |
| LASBarrelAlgorithm () | |
Private Member Functions | |
| void | ReadMisalignmentFromFile (const char *, LASGlobalData< LASCoordinateSet > &, LASGlobalData< LASCoordinateSet > &) |
| void | ReadStartParametersFromFile (const char *, float[52]) |
Private Attributes | |
| TMinuit * | minuit |
implementation of the alignment tube algorithm
Definition at line 23 of file LASBarrelAlgorithm.h.
| LASBarrelAlgorithm::LASBarrelAlgorithm | ( | ) |
Definition at line 18 of file LASBarrelAlgorithm.cc.
References minuit.
{
minuit = new TMinuit( 52 );
}
| LASBarrelAlignmentParameterSet LASBarrelAlgorithm::CalculateParameters | ( | LASGlobalData< LASCoordinateSet > & | measuredCoordinates, |
| LASGlobalData< LASCoordinateSet > & | nominalCoordinates | ||
| ) |
The minimization of the equation system for the barrel. For documentation, please refer to the TkLasATModel TWiki page: TWiki > CMS Web > CMSTrackerLaserAlignmenSystem > TkLasBarrelAlgorithm > TkLasATModel
Definition at line 33 of file LASBarrelAlgorithm.cc.
References fitWZ::arglist, gather_cfg::cout, fcn(), LASBarrelAlignmentParameterSet::GetParameter(), minuit, and Gflash::par.
Referenced by LaserAlignment::endRun().
{
std::cout << " [LASBarrelAlgorithm::CalculateParameters] -- Starting." << std::endl;
// for testing..
//ReadMisalignmentFromFile( "misalign-var.txt", measuredCoordinates, nominalCoordinates );
// statics for minuit
aMeasuredCoordinates = &measuredCoordinates;
aNominalCoordinates = &nominalCoordinates;
// minimizer and variables for it
minuit->SetFCN( fcn );
double arglist[10];
int _ierflg = 0;
// toggle minuit blabla
arglist[0] = -1;
minuit->mnexcm( "SET PRI", arglist, 1, _ierflg );
// set par errors
arglist[0] = 1;
minuit->mnexcm( "SET ERR", arglist ,1, _ierflg );
//
// define 52 parameters
//
// start values: to be evacuated to cfg
static float _vstart[52] = {
0.00, 0.00, 0.0, 0.0, 0.0, 0.0, // subdet for TIB+
0.00, 0.00, 0.0, 0.0, 0.0, 0.0, // subdet for TIB-
0.00, 0.00, 0.0, 0.0, 0.0, 0.0, // subdet for TOB+
0.00, 0.00, 0.0, 0.0, 0.0, 0.0, // subdet for TOB-
0.00, 0.00, 0.0, 0.0, 0.0, 0.0, // subdet for TEC+
0.00, 0.00, 0.0, 0.0, 0.0, 0.0, // subdet for TEC-
0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, // beams 0-3
0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 // beams 4-7
};
// ReadStartParametersFromFile( "startParameters.txt", _vstart ); // debug
// step sizes: to be tuned, to be evacuated to cfg
static float _vstep[52] = {
0.001, 0.001, 0.1, 0.1, 0.1, 0.1, // subdet for TIB+
0.001, 0.001, 0.1, 0.1, 0.1, 0.1, // subdet for TIB-
0.001, 0.001, 0.1, 0.1, 0.1, 0.1, // subdet for TOB+
0.001, 0.001, 0.1, 0.1, 0.1, 0.1, // subdet for TOB-
0.001, 0.001, 0.1, 0.1, 0.1, 0.1, // subdet for TEC+
0.001, 0.001, 0.1, 0.1, 0.1, 0.1, // subdet for TEC-
0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, // beams 0-3
0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001 // beams 4-7
};
// subdetector parameters for TIB+:
// rotation around z of first end face
minuit->mnparm( 0, "subRot1TIB+", _vstart[0], _vstep[0], 0, 0, _ierflg );
// rotation around z of second end face
minuit->mnparm( 1, "subRot2TIB+", _vstart[1], _vstep[1], 0, 0, _ierflg );
// translation along x of first end face
minuit->mnparm( 2, "subTransX1TIB+", _vstart[2], _vstep[2], 0, 0, _ierflg );
// translation along x of second end face
minuit->mnparm( 3, "subTransX2TIB+", _vstart[3], _vstep[3], 0, 0, _ierflg );
// translation along y of first end face
minuit->mnparm( 4, "subTransY1TIB+", _vstart[4], _vstep[4], 0, 0, _ierflg );
// translation along y of second end face
minuit->mnparm( 5, "subTransY2TIB+", _vstart[5], _vstep[5], 0, 0, _ierflg );
// subdetector parameters for TIB-:
// rotation around z of first end face
minuit->mnparm( 6, "subRot1TIB-", _vstart[6], _vstep[6], 0, 0, _ierflg );
// rotation around z of second end face
minuit->mnparm( 7, "subRot2TIB-", _vstart[7], _vstep[7], 0, 0, _ierflg );
// translation along x of first end face
minuit->mnparm( 8, "subTransX1TIB-", _vstart[8], _vstep[8], 0, 0, _ierflg );
// translation along x of second end face
minuit->mnparm( 9, "subTransX2TIB-", _vstart[9], _vstep[9], 0, 0, _ierflg );
// translation along y of first end face
minuit->mnparm( 10, "subTransY1TIB-", _vstart[10], _vstep[10], 0, 0, _ierflg );
// translation along y of second end face
minuit->mnparm( 11, "subTransY2TIB-", _vstart[11], _vstep[11], 0, 0, _ierflg );
// subdetector parameters for TOB+:
// rotation around z of first end face
minuit->mnparm( 12, "subRot1TOB+", _vstart[12], _vstep[12], 0, 0, _ierflg );
// rotation around z of second end face
minuit->mnparm( 13, "subRot2TOB+", _vstart[13], _vstep[13], 0, 0, _ierflg );
// translation along x of first end face
minuit->mnparm( 14, "subTransX1TOB+", _vstart[14], _vstep[14], 0, 0, _ierflg );
// translation along x of second end face
minuit->mnparm( 15, "subTransX2TOB+", _vstart[15], _vstep[15], 0, 0, _ierflg );
// translation along y of first end face
minuit->mnparm( 16, "subTransY1TOB+", _vstart[16], _vstep[16], 0, 0, _ierflg );
// translation along y of second end face
minuit->mnparm( 17, "subTransY2TOB+", _vstart[17], _vstep[17], 0, 0, _ierflg );
// subdetector parameters for TOB-:
// rotation around z of first end face
minuit->mnparm( 18, "subRot1TOB-", _vstart[18], _vstep[18], 0, 0, _ierflg );
// rotation around z of second end face
minuit->mnparm( 19, "subRot2TOB-", _vstart[19], _vstep[19], 0, 0, _ierflg );
// translation along x of first end face
minuit->mnparm( 20, "subTransX1TOB-", _vstart[20], _vstep[20], 0, 0, _ierflg );
// translation along x of second end face
minuit->mnparm( 21, "subTransX2TOB-", _vstart[21], _vstep[21], 0, 0, _ierflg );
// translation along y of first end face
minuit->mnparm( 22, "subTransY1TOB-", _vstart[22], _vstep[22], 0, 0, _ierflg );
// translation along y of second end face
minuit->mnparm( 23, "subTransY2TOB-", _vstart[23], _vstep[23], 0, 0, _ierflg );
// subdetector parameters for TEC+:
// rotation around z of first end face
minuit->mnparm( 24, "subRot1TEC+", _vstart[24], _vstep[24], 0, 0, _ierflg );
// rotation around z of second end face
minuit->mnparm( 25, "subRot2TEC+", _vstart[25], _vstep[25], 0, 0, _ierflg );
// translation along x of first end face
minuit->mnparm( 26, "subTransX1TEC+", _vstart[26], _vstep[26], 0, 0, _ierflg );
// translation along x of second end face
minuit->mnparm( 27, "subTransX2TEC+", _vstart[27], _vstep[27], 0, 0, _ierflg );
// translation along y of first end face
minuit->mnparm( 28, "subTransY1TEC+", _vstart[28], _vstep[28], 0, 0, _ierflg );
// translation along y of second end face
minuit->mnparm( 29, "subTransY2TEC+", _vstart[29], _vstep[29], 0, 0, _ierflg );
// subdetector parameters for TEC-:
// rotation around z of first end face
minuit->mnparm( 30, "subRot1TEC-", _vstart[30], _vstep[30], 0, 0, _ierflg );
// rotation around z of second end face
minuit->mnparm( 31, "subRot2TEC-", _vstart[31], _vstep[31], 0, 0, _ierflg );
// translation along x of first end face
minuit->mnparm( 32, "subTransX1TEC-", _vstart[32], _vstep[32], 0, 0, _ierflg );
// translation along x of second end face
minuit->mnparm( 33, "subTransX2TEC-", _vstart[33], _vstep[33], 0, 0, _ierflg );
// translation along y of first end face
minuit->mnparm( 34, "subTransY1TEC-", _vstart[34], _vstep[34], 0, 0, _ierflg );
// translation along y of second end face
minuit->mnparm( 35, "subTransY2TEC-", _vstart[35], _vstep[35], 0, 0, _ierflg );
// beam parameters (+-z pairs, duplicated for beams 0-7):
// rotation around z at zt1
minuit->mnparm( 36, "beamRot1Beam0", _vstart[36], _vstep[36], 0, 0, _ierflg );
// rotation around z at zt2
minuit->mnparm( 37, "beamRot2Beam0", _vstart[37], _vstep[37], 0, 0, _ierflg );
// rotation around z at zt1
minuit->mnparm( 38, "beamRot1Beam1", _vstart[38], _vstep[38], 0, 0, _ierflg );
// rotation around z at zt2
minuit->mnparm( 39, "beamRot2Beam1", _vstart[39], _vstep[39], 0, 0, _ierflg );
// rotation around z at zt1
minuit->mnparm( 40, "beamRot1Beam2", _vstart[40], _vstep[40], 0, 0, _ierflg );
// rotation around z at zt2
minuit->mnparm( 41, "beamRot2Beam2", _vstart[41], _vstep[41], 0, 0, _ierflg );
// rotation around z at zt1
minuit->mnparm( 42, "beamRot1Beam3", _vstart[42], _vstep[42], 0, 0, _ierflg );
// rotation around z at zt2
minuit->mnparm( 43, "beamRot2Beam3", _vstart[43], _vstep[43], 0, 0, _ierflg );
// rotation around z at zt1
minuit->mnparm( 44, "beamRot1Beam4", _vstart[44], _vstep[44], 0, 0, _ierflg );
// rotation around z at zt2
minuit->mnparm( 45, "beamRot2Beam4", _vstart[45], _vstep[45], 0, 0, _ierflg );
// rotation around z at zt1
minuit->mnparm( 46, "beamRot1Beam5", _vstart[46], _vstep[46], 0, 0, _ierflg );
// rotation around z at zt2
minuit->mnparm( 47, "beamRot2Beam5", _vstart[47], _vstep[47], 0, 0, _ierflg );
// rotation around z at zt1
minuit->mnparm( 48, "beamRot1Beam6", _vstart[48], _vstep[48], 0, 0, _ierflg );
// rotation around z at zt2
minuit->mnparm( 49, "beamRot2Beam6", _vstart[49], _vstep[49], 0, 0, _ierflg );
// rotation around z at zt1
minuit->mnparm( 50, "beamRot1Beam7", _vstart[50], _vstep[50], 0, 0, _ierflg );
// rotation around z at zt2
minuit->mnparm( 51, "beamRot2Beam7", _vstart[51], _vstep[51], 0, 0, _ierflg );
// we fix the respective outer disks 9 of each endcap
// as a reference system (pars 25,27,29 & 30,32,34)
// note: minuit numbering is fortran style...
arglist[0] = 26; arglist[1] = 28; arglist[2] = 30;
// minuit->mnexcm( "FIX", arglist ,3, _ierflg ); // TEC+
arglist[0] = 31; arglist[1] = 33; arglist[2] = 35;
// minuit->mnexcm( "FIX", arglist ,3, _ierflg ); // TEC-
// DEBUG: FIX BEAM PARAMETERS /////////////////////////////////////////////////////////////////////
double parlist[16];
for( int par = 37; par <= 52; ++par ) parlist[par-37] = par;
minuit->mnexcm( "FIX", parlist ,16, _ierflg );
// DEBUG: FIX ALGN PARAMETERS /////////////////////////////////////////////////////////////////////
// double parlist[36];
// for( int par = 1; par <= 36; ++par ) parlist[par-1] = par;
// minuit->mnexcm( "FIX", parlist ,36, _ierflg );
// now ready for minimization step
arglist[0] = 10000;
arglist[1] = 0.1;
minuit->mnexcm( "MIGRAD", arglist , 2, _ierflg ); // minimizer
// minuit->mnexcm( "MINOS", arglist , 1, _ierflg ); // error recalculation
// now fill the result vector.
// turned out that the parameter numbering is stupid, change this later..
LASBarrelAlignmentParameterSet theResult;
double par = 0., parError = 0.;
// TEC+ rot
minuit->GetParameter( 24, par, parError ); theResult.GetParameter( 0, 0, 0 ).first = par; theResult.GetParameter( 0, 0, 0 ).second = parError;
minuit->GetParameter( 25, par, parError ); theResult.GetParameter( 0, 1, 0 ).first = par; theResult.GetParameter( 0, 1, 0 ).second = parError;
// TEC+ x
minuit->GetParameter( 26, par, parError ); theResult.GetParameter( 0, 0, 1 ).first = par; theResult.GetParameter( 0, 0, 1 ).second = parError;
minuit->GetParameter( 27, par, parError ); theResult.GetParameter( 0, 1, 1 ).first = par; theResult.GetParameter( 0, 1, 1 ).second = parError;
// TEC+ x
minuit->GetParameter( 28, par, parError ); theResult.GetParameter( 0, 0, 2 ).first = par; theResult.GetParameter( 0, 0, 2 ).second = parError;
minuit->GetParameter( 29, par, parError ); theResult.GetParameter( 0, 1, 2 ).first = par; theResult.GetParameter( 0, 1, 2 ).second = parError;
// TEC- rot
minuit->GetParameter( 30, par, parError ); theResult.GetParameter( 1, 0, 0 ).first = par; theResult.GetParameter( 1, 0, 0 ).second = parError;
minuit->GetParameter( 31, par, parError ); theResult.GetParameter( 1, 1, 0 ).first = par; theResult.GetParameter( 1, 1, 0 ).second = parError;
// TEC- x
minuit->GetParameter( 32, par, parError ); theResult.GetParameter( 1, 0, 1 ).first = par; theResult.GetParameter( 1, 0, 1 ).second = parError;
minuit->GetParameter( 33, par, parError ); theResult.GetParameter( 1, 1, 1 ).first = par; theResult.GetParameter( 1, 1, 1 ).second = parError;
// TEC- x
minuit->GetParameter( 34, par, parError ); theResult.GetParameter( 1, 0, 2 ).first = par; theResult.GetParameter( 1, 0, 2 ).second = parError;
minuit->GetParameter( 35, par, parError ); theResult.GetParameter( 1, 1, 2 ).first = par; theResult.GetParameter( 1, 1, 2 ).second = parError;
// TIB+ rot
minuit->GetParameter( 0, par, parError ); theResult.GetParameter( 2, 0, 0 ).first = par; theResult.GetParameter( 2, 0, 0 ).second = parError;
minuit->GetParameter( 1, par, parError ); theResult.GetParameter( 2, 1, 0 ).first = par; theResult.GetParameter( 2, 1, 0 ).second = parError;
// TIB+ x
minuit->GetParameter( 2, par, parError ); theResult.GetParameter( 2, 0, 1 ).first = par; theResult.GetParameter( 2, 0, 1 ).second = parError;
minuit->GetParameter( 3, par, parError ); theResult.GetParameter( 2, 1, 1 ).first = par; theResult.GetParameter( 2, 1, 1 ).second = parError;
// TIB+ x
minuit->GetParameter( 4, par, parError ); theResult.GetParameter( 2, 0, 2 ).first = par; theResult.GetParameter( 2, 0, 2 ).second = parError;
minuit->GetParameter( 5, par, parError ); theResult.GetParameter( 2, 1, 2 ).first = par; theResult.GetParameter( 2, 1, 2 ).second = parError;
// TIB- rot
minuit->GetParameter( 6, par, parError ); theResult.GetParameter( 3, 0, 0 ).first = par; theResult.GetParameter( 3, 0, 0 ).second = parError;
minuit->GetParameter( 7, par, parError ); theResult.GetParameter( 3, 1, 0 ).first = par; theResult.GetParameter( 3, 1, 0 ).second = parError;
// TIB- x
minuit->GetParameter( 8, par, parError ); theResult.GetParameter( 3, 0, 1 ).first = par; theResult.GetParameter( 3, 0, 1 ).second = parError;
minuit->GetParameter( 9, par, parError ); theResult.GetParameter( 3, 1, 1 ).first = par; theResult.GetParameter( 3, 1, 1 ).second = parError;
// TIB- x
minuit->GetParameter( 10, par, parError ); theResult.GetParameter( 3, 0, 2 ).first = par; theResult.GetParameter( 3, 0, 2 ).second = parError;
minuit->GetParameter( 11, par, parError ); theResult.GetParameter( 3, 1, 2 ).first = par; theResult.GetParameter( 3, 1, 2 ).second = parError;
// TOB+ rot
minuit->GetParameter( 12, par, parError ); theResult.GetParameter( 4, 0, 0 ).first = par; theResult.GetParameter( 4, 0, 0 ).second = parError;
minuit->GetParameter( 13, par, parError ); theResult.GetParameter( 4, 1, 0 ).first = par; theResult.GetParameter( 4, 1, 0 ).second = parError;
// TOB+ x
minuit->GetParameter( 14, par, parError ); theResult.GetParameter( 4, 0, 1 ).first = par; theResult.GetParameter( 4, 0, 1 ).second = parError;
minuit->GetParameter( 15, par, parError ); theResult.GetParameter( 4, 1, 1 ).first = par; theResult.GetParameter( 4, 1, 1 ).second = parError;
// TOB+ x
minuit->GetParameter( 16, par, parError ); theResult.GetParameter( 4, 0, 2 ).first = par; theResult.GetParameter( 4, 0, 2 ).second = parError;
minuit->GetParameter( 17, par, parError ); theResult.GetParameter( 4, 1, 2 ).first = par; theResult.GetParameter( 4, 1, 2 ).second = parError;
// TOB- rot
minuit->GetParameter( 18, par, parError ); theResult.GetParameter( 5, 0, 0 ).first = par; theResult.GetParameter( 5, 0, 0 ).second = parError;
minuit->GetParameter( 19, par, parError ); theResult.GetParameter( 5, 1, 0 ).first = par; theResult.GetParameter( 5, 1, 0 ).second = parError;
// TOB- x
minuit->GetParameter( 20, par, parError ); theResult.GetParameter( 5, 0, 1 ).first = par; theResult.GetParameter( 5, 0, 1 ).second = parError;
minuit->GetParameter( 21, par, parError ); theResult.GetParameter( 5, 1, 1 ).first = par; theResult.GetParameter( 5, 1, 1 ).second = parError;
// TOB- x
minuit->GetParameter( 22, par, parError ); theResult.GetParameter( 5, 0, 2 ).first = par; theResult.GetParameter( 5, 0, 2 ).second = parError;
minuit->GetParameter( 23, par, parError ); theResult.GetParameter( 5, 1, 2 ).first = par; theResult.GetParameter( 5, 1, 2 ).second = parError;
std::cout << " [LASBarrelAlgorithm::CalculateParameters] -- Done." << std::endl;
return theResult;
}
| void LASBarrelAlgorithm::Dump | ( | void | ) |
Print resulting parameters to stdout and to a file - for debugging only
Definition at line 558 of file LASBarrelAlgorithm.cc.
References ExpressReco_HICollisions_FallBack::cerr, gather_cfg::cout, error, dbtoweb::file, minuit, Gflash::par, relativeConstraints::value, and z.
{
if( !minuit ) {
std::cerr << " [LASBarrelAlgorithm::Dump] ** WARNING: minimizer object uninitialized." << std::endl;
return;
}
std::cout << std::endl << " [LASBarrelAlgorithm::Dump] -- Parameter dump: " << std::endl;
const int subdetParMap[6] = { 24, 30, 0, 6, 12, 18 }; // map to one-dim array
const std::string subdetNames[6] = { " TEC+ ", " TEC- ", " TIB+ ", " TIB- ", " TOB+ ", " TOB- " };
double value, error;
std::cout << " Detector parameters: " << std::endl;
std::cout << " -------------------" << std::endl;
std::cout << " Values: PHI1 X1 Y1 PHI2 X2 Y2 " << std::endl;
for( int subdet = 0; subdet < 6; ++subdet ) {
std::cout <<subdetNames[subdet];
for( int par = subdetParMap[subdet]; par <= subdetParMap[subdet] + 4; par += 2 ) {
minuit->GetParameter( par, value, error );
std::cout << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << value;
}
for( int par = subdetParMap[subdet] + 1; par <= subdetParMap[subdet] + 5; par += 2 ) {
minuit->GetParameter( par, value, error );
std::cout << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << value;
}
std::cout << std::endl;
}
std::cout << " Errors: PHI1 X1 Y1 PHI2 X2 Y2 " << std::endl;
for( int subdet = 0; subdet < 6; ++subdet ) {
std::cout <<subdetNames[subdet];
for( int par = subdetParMap[subdet]; par <= subdetParMap[subdet] + 4; par += 2 ) {
minuit->GetParameter( par, value, error );
std::cout << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << error;
}
for( int par = subdetParMap[subdet] + 1; par <= subdetParMap[subdet] + 5; par += 2 ) {
minuit->GetParameter( par, value, error );
std::cout << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << error;
}
std::cout << std::endl;
}
std::cout << std::endl;
std::cout << " Beam parameters: " << std::endl;
std::cout << " ---------------" << std::endl;
std::cout << " Values: PHI1 PHI2" << std::endl;
for( int beam = 0; beam < 8; ++beam ) {
std::cout << " " << beam << " ";
for( int z = 0; z < 2; ++z ) {
minuit->GetParameter( 36 + 2 * beam + z, value, error );
std::cout << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << value;
}
std::cout << std::endl;
}
std::cout << " Errors: PHI1 PHI2" << std::endl;
for( int beam = 0; beam < 8; ++beam ) {
std::cout << " " << beam << " ";
for( int z = 0; z < 2; ++z ) {
minuit->GetParameter( 36 + 2 * beam + z, value, error );
std::cout << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << error;
}
std::cout << std::endl;
}
// det parameters once again without leading column (for easy read-in), into a file
ofstream file( "/afs/cern.ch/user/o/olzem/public/parameters_det.txt" );
for( int subdet = 0; subdet < 6; ++subdet ) {
for( int par = subdetParMap[subdet]; par <= subdetParMap[subdet] + 4; par += 2 ) {
minuit->GetParameter( par, value, error );
file << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << value;
}
for( int par = subdetParMap[subdet] + 1; par <= subdetParMap[subdet] + 5; par += 2 ) {
minuit->GetParameter( par, value, error );
file << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << value;
}
file << std::endl;
}
file.close();
// same for beam parameters
file.open( "/afs/cern.ch/user/o/olzem/public/parameters_beam.txt" );
for( int beam = 0; beam < 8; ++beam ) {
for( int z = 0; z < 2; ++z ) {
minuit->GetParameter( 36 + 2 * beam + z, value, error );
file << std::setw( 12 ) << std::setprecision( 6 ) << std::fixed << value;
}
file << std::endl;
}
file.close();
std:: cout << " [LASBarrelAlgorithm::Dump] -- End parameter dump." << std::endl;
std::cout << std::endl;
}
| void LASBarrelAlgorithm::ReadMisalignmentFromFile | ( | const char * | filename, |
| LASGlobalData< LASCoordinateSet > & | measuredCoordinates, | ||
| LASGlobalData< LASCoordinateSet > & | nominalCoordinates | ||
| ) | [private] |
allows to push in a simple simulated misalignment for quick internal testing purposes; overwrites LASGlobalData<LASCoordinateSet>& measuredCoordinates; call at beginning of LASBarrelAlgorithm::CalculateParameters method
one line per module, format for TEC: det ring beam disk phi phiErr format for TEC(at) & TIBTOB: det beam z "-1" phi phiErr
Definition at line 671 of file LASBarrelAlgorithm.cc.
References ExpressReco_HICollisions_FallBack::cerr, dbtoweb::file, LASGlobalData< T >::GetTEC2TECEntry(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), phi, pos, relativeConstraints::ring, LASCoordinateSet::SetPhi(), LASCoordinateSet::SetPhiError(), LASGlobalLoop::TEC2TECLoop(), LASGlobalLoop::TECLoop(), LASGlobalLoop::TIBTOBLoop(), and z.
{
ifstream file( filename );
if( file.bad() ) {
std::cerr << " [LASBarrelAlgorithm::ReadMisalignmentFromFile] ** ERROR: cannot open file \"" << filename << "\"." << std::endl;
return;
}
std::cerr << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << std::endl;
std::cerr << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << std::endl;
std::cerr << " [LASBarrelAlgorithm::ReadMisalignmentFromFile] ** WARNING: you are reading a fake measurement from a file!" << std::endl;
std::cerr << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << std::endl;
std::cerr << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << std::endl;
// the measured coordinates will finally be overwritten;
// first, set them to the nominal values
measuredCoordinates = nominalCoordinates;
// and put large errors on all values;
{
LASGlobalLoop moduleLoop;
int det, ring, beam, disk, pos;
det = 0; ring = 0; beam = 0; disk = 0;
do {
measuredCoordinates.GetTECEntry( det, ring, beam, disk ).SetPhiError( 1000. );
} while( moduleLoop.TECLoop( det, ring, beam, disk ) );
det = 2; beam = 0; pos = 0;
do {
measuredCoordinates.GetTIBTOBEntry( det, beam, pos ).SetPhiError( 1000. );
} while( moduleLoop.TIBTOBLoop( det, beam, pos ) );
det = 0; beam = 0; disk = 0;
do {
measuredCoordinates.GetTEC2TECEntry( det, beam, disk ).SetPhiError( 1000. );
} while( moduleLoop.TEC2TECLoop( det, beam, disk ) );
}
// buffers for read-in
int det, beam, z, ring;
double phi, phiError;
while( !file.eof() ) {
file >> det;
if( file.eof() ) break; // do not read the last line twice, do not fill trash if file empty
file >> beam; file >> z; file >> ring;
file >> phi; file >> phiError;
if( det > 1 ) { // TIB/TOB
measuredCoordinates.GetTIBTOBEntry( det, beam, z ).SetPhi( phi );
measuredCoordinates.GetTIBTOBEntry( det, beam, z ).SetPhiError( phiError );
} else { // TEC or TEC(at)
if( ring > -1 ) { // TEC
measuredCoordinates.GetTECEntry( det, ring, beam, z ).SetPhi( phi );
measuredCoordinates.GetTECEntry( det, ring, beam, z ).SetPhiError( phiError );
}
else { // TEC(at)
measuredCoordinates.GetTEC2TECEntry( det, beam, z ).SetPhi( phi );
measuredCoordinates.GetTEC2TECEntry( det, beam, z ).SetPhiError( phiError );
}
}
}
file.close();
}
| void LASBarrelAlgorithm::ReadStartParametersFromFile | ( | const char * | filename, |
| float | values[52] | ||
| ) | [private] |
this function is here only for debugging, don't use. file format: <phi1> <x1> <y1> <phi2> <x2> <y2> // for TEC* " " " " " " // TEC- .. then for TIB+, TIB-, TOB+, TOB- index 1 if for lower z, 2 for higher z
Definition at line 758 of file LASBarrelAlgorithm.cc.
References ExpressReco_HICollisions_FallBack::cerr, and dbtoweb::file.
{
ifstream file( filename );
if( file.bad() ) {
std::cerr << " [LASBarrelAlgorithm::ReadStartParametersFromFile] ** ERROR: cannot open file \"" << filename << "\"." << std::endl;
return;
}
std::cerr << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << std::endl;
std::cerr << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << std::endl;
std::cerr << " [LASBarrelAlgorithm::ReadStartParametersFrom File] ** WARNING: you are reading parameter start values from a file!" << std::endl;
std::cerr << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << std::endl;
std::cerr << " @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" << std::endl;
// map to the minuit par array
const int subdetParMap[6] = { 24, 30, 0, 6, 12, 18 };
for( int det = 0; det < 6; ++det ) {
file >> values[subdetParMap[det]]; // phi1
file >> values[subdetParMap[det]+2]; // x1
file >> values[subdetParMap[det]+4]; // y1
file >> values[subdetParMap[det]+1]; // phi2
file >> values[subdetParMap[det]+3]; // x2
file >> values[subdetParMap[det]+5]; // y2
}
}
TMinuit* LASBarrelAlgorithm::minuit [private] |
Definition at line 33 of file LASBarrelAlgorithm.h.
Referenced by CalculateParameters(), Dump(), and LASBarrelAlgorithm().
1.7.3