LASBarrelAlgorithm Class Reference

#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

Detailed Description

implementation of the alignment tube algorithm

Definition at line 22 of file LASBarrelAlgorithm.h.

Constructor & Destructor Documentation

LASBarrelAlgorithm()

LASBarrelAlgorithm::LASBarrelAlgorithm

(

)

Definition at line 13 of file LASBarrelAlgorithm.cc.

{ minuit = new TMinuit(52); }

References minuit.

Member Function Documentation

CalculateParameters()

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 20 of file LASBarrelAlgorithm.cc.

                                                                                                               {
  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;
}

References aMeasuredCoordinates, aNominalCoordinates, fitWZ::arglist, gather_cfg::cout, fcn(), LASBarrelAlignmentParameterSet::GetParameter(), and minuit.

Referenced by LaserAlignment::endRunProduce().

Dump()

void LASBarrelAlgorithm::Dump

(

void

)

Print resulting parameters to stdout and to a file - for debugging only

Definition at line 603 of file LASBarrelAlgorithm.cc.

                                  {
  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
  std::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;
}

References EcalCondDBWriter_cfi::beam, beam_dqm_sourceclient-live_cfg::cerr, gather_cfg::cout, relativeConstraints::error, FrontierConditions_GlobalTag_cff::file, alignBH_cfg::fixed, minuit, AlCaHLTBitMon_QueryRunRegistry::string, subdetNames, relativeConstraints::value, and z.

ReadMisalignmentFromFile()

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 707 of file LASBarrelAlgorithm.cc.

                                                                                                       {
  std::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();
}

References EcalCondDBWriter_cfi::beam, beam_dqm_sourceclient-live_cfg::cerr, FrontierConditions_GlobalTag_cff::file, corrVsCorr::filename, LASGlobalData< T >::GetTEC2TECEntry(), LASGlobalData< T >::GetTECEntry(), LASGlobalData< T >::GetTIBTOBEntry(), phi, muonME0PseudoDigis_cfi::phiError, relativeConstraints::ring, LASCoordinateSet::SetPhi(), LASCoordinateSet::SetPhiError(), LASGlobalLoop::TEC2TECLoop(), LASGlobalLoop::TECLoop(), LASGlobalLoop::TIBTOBLoop(), and z.

ReadStartParametersFromFile()

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 805 of file LASBarrelAlgorithm.cc.

                                                                                           {
  std::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
  }
}

References beam_dqm_sourceclient-live_cfg::cerr, FrontierConditions_GlobalTag_cff::file, corrVsCorr::filename, and contentValuesCheck::values.

Member Data Documentation

minuit

TMinuit* LASBarrelAlgorithm::minuit

private

Definition at line 32 of file LASBarrelAlgorithm.h.

Referenced by CalculateParameters(), Dump(), and LASBarrelAlgorithm().