getSimConstants.cxx

Go to the documentation of this file.

void getSimConstants() {
  gROOT->Reset();
  gROOT->SetStyle( "Plain" );
  gStyle->SetOptStat( 00000 );
 
  const unsigned int barScin = 176;
  const unsigned int endScin = 96;
  const unsigned int barEnd  = 2;
 
  const unsigned int crateNum        = 2;
  const unsigned int crate[crateNum] = { 64, 65 };
  const unsigned int boardNum        = 14;
  unsigned int       board[boardNum];
  for ( unsigned int i = 0; i < boardNum; i++ ) { board[i] = i + 1; }
  const unsigned int channelNum = 16;
  const unsigned int paramNum   = 10;
  double             param[crateNum][boardNum][channelNum][paramNum];
 
  unsigned int ifee[crateNum][boardNum][channelNum];
  unsigned int icrate[crateNum][boardNum][channelNum];
  unsigned int iboard[crateNum][boardNum][channelNum];
  unsigned int ichannel[crateNum][boardNum][channelNum];
  ifstream     inData[crateNum][boardNum];
  char         inFile[crateNum][boardNum][1000];
 
  //  char kind[100]=".txt";
  char kind[100] = "-Curve.txt";
 
  ofstream outData1, outData2, outData3;
  char     outFile1[800], outFile2[800], outFile3[800];
  sprintf( outFile1, "ElecBarEast%s", kind );
  sprintf( outFile2, "ElecBarWest%s", kind );
  sprintf( outFile3, "ElecEndcap%s", kind );
 
  for ( unsigned int i = 0; i < crateNum; i++ )
  {
    for ( unsigned int j = 0; j < boardNum; j++ )
    {
 
      sprintf( inFile[i][j],
               "/ihepbatch/besd09/sunss/Electronics/09summer/Data/Crate%i/Board%i/"
               "Crate%i-Board%i%s",
               crate[i], board[j], crate[i], board[j], kind );
      cout << inFile[i][j] << endl;
      inData[i][j].open( inFile[i][j] );
      for ( unsigned int k = 0; k < channelNum; k++ )
      {
        inData[i][j] >> ifee[i][j][k] >> icrate[i][j][k] >> iboard[i][j][k] >>
            ichannel[i][j][k];
        for ( unsigned int l = 0; l < paramNum; l++ ) { inData[i][j] >> param[i][j][k][l]; }
      }
      inData[i][j].close();
    }
  }
 
  int be[crateNum][boardNum][channelNum];
  int counter[crateNum][boardNum][channelNum];
  int end[crateNum][boardNum][channelNum];
  for ( unsigned int i = 0; i < crateNum; i++ )
  {
    for ( unsigned int j = 0; j < boardNum; j++ )
    {
      for ( unsigned int k = 0; k < channelNum; k++ )
      {
        if ( j < 3 )
        {
          be[i][j][k]  = 0;
          end[i][j][k] = 0;
          if ( i == 0 )
          {
            if ( j == 0 ) { counter[i][j][k] = k; }
            else if ( j == 1 )
            {
              if ( k < channelNum / 2 )
              {                                    // 8
                counter[i][j][k] = k + channelNum; // 16
              }
              else
              {
                counter[i][j][k] = k - channelNum / 2 + 48; // 8
              }
            }
            else if ( j == 2 )
            {
              counter[i][j][k] = k + channelNum / 2 + 48; // 8
            }
          }
          else if ( i == 1 )
          {
            if ( j == 0 )
            {
              counter[i][j][k] = k + channelNum + channelNum / 2; // 24
            }
            else if ( j == 1 )
            {
              if ( k < channelNum / 2 )
              {                                                         // 8
                counter[i][j][k] = k + channelNum * 2 + channelNum / 2; // 40
              }
              else
              {
                counter[i][j][k] = k + channelNum + 48; // 16
              }
            }
            else if ( j == 2 )
            {
              counter[i][j][k] = k + channelNum * 2 + 48; // 32
            }
          }
        }
        else
        {
          be[i][j][k] = 1;
          if ( ( k % 2 ) == 0 )
          {
            end[i][j][k] = 0; // east
          }
          else
          {
            end[i][j][k] = 1; // west
          }
          if ( i == 0 )
          {
            if ( ( k % 4 ) < 2 )
            {
              counter[i][j][k] = ( j - 3 ) * channelNum / 4 + k / 4; // (j-3)*4+1/4
            }
            else
            {
              counter[i][j][k] = ( j - 3 ) * channelNum / 4 + k / 4 + 88; // (j-3)*4+1/4
            }
          }
          else if ( i == 1 )
          {
            if ( ( k % 4 ) < 2 )
            { counter[i][j][k] = ( boardNum - 3 + j - 3 ) * channelNum / 4 + k / 4; }
            else { counter[i][j][k] = ( boardNum - 3 + j - 3 ) * channelNum / 4 + k / 4 + 88; }
          }
        }
        //  cout << i << " " << j << " " << k << " " << be[i][j][k] << " " <<
        // counter[i][j][k] << " " << end[i][j][k] << endl;
      }
    }
  }
 
  for ( unsigned int i = 0; i < crateNum; i++ )
  {
    for ( unsigned int j = 0; j < boardNum; j++ )
    {
      for ( unsigned int k = 0; k < channelNum; k++ )
      {
        if ( be[i][j][k] == 0 )
        {
          counter[i][j][k] = counter[i][j][k] - 12;
          if ( counter[i][j][k] < 0 || ( counter[i][j][k] >= 36 && counter[i][j][k] < 48 ) )
          { counter[i][j][k] = counter[i][j][k] + 48; }
        }
        else
        {
          counter[i][j][k] = counter[i][j][k] - 11;
          if ( counter[i][j][k] < 0 || ( counter[i][j][k] >= 77 && counter[i][j][k] < 88 ) )
          { counter[i][j][k] = counter[i][j][k] + 88; }
        }
        //  cout << i << " " << j << " " << k << " " << be[i][j][k] << " " <<
        // counter[i][j][k] << " " << end[i][j][k] << endl;
      }
    }
  }
 
  for ( unsigned int i = 0; i < crateNum; i++ )
  {
    for ( unsigned int j = 0; j < boardNum; j++ )
    {
      for ( unsigned int k = 0; k < channelNum; k++ )
      {
        if ( i == 0 ) { icrate[i][j][k] = 64; }
        else { icrate[i][j][k] = 65; }
        iboard[i][j][k]   = j + 1;
        ichannel[i][j][k] = k + 1;
      }
    }
  }
 
  unsigned int Fee[barScin][barEnd], ecFee[endScin];
  unsigned int Crate[barScin][barEnd], ecCrate[endScin];
  unsigned int Board[barScin][barEnd], ecBoard[endScin];
  unsigned int Channel[barScin][barEnd], ecChannel[endScin];
  double       pbeast[barScin][paramNum], pbwest[barScin][paramNum];
  double       pendcap[endScin][paramNum], ;
  for ( unsigned int i = 0; i < crateNum; i++ )
  {
    for ( unsigned int j = 0; j < boardNum; j++ )
    {
      for ( unsigned int k = 0; k < channelNum; k++ )
      {
        if ( be[i][j][k] == 1 )
        {
          if ( end[i][j][k] == 0 )
          {
            Fee[counter[i][j][k]][0]     = ifee[i][j][k];
            Crate[counter[i][j][k]][0]   = icrate[i][j][k];
            Board[counter[i][j][k]][0]   = iboard[i][j][k];
            Channel[counter[i][j][k]][0] = ichannel[i][j][k];
            for ( unsigned int l = 0; l < paramNum; l++ )
            { pbeast[counter[i][j][k]][l] = param[i][j][k][l]; }
          }
          else
          {
            Fee[counter[i][j][k]][1]     = ifee[i][j][k];
            Crate[counter[i][j][k]][1]   = icrate[i][j][k];
            Board[counter[i][j][k]][1]   = iboard[i][j][k];
            Channel[counter[i][j][k]][1] = ichannel[i][j][k];
            for ( unsigned int l = 0; l < paramNum; l++ )
            { pbwest[counter[i][j][k]][l] = param[i][j][k][l]; }
          }
        }
        else
        {
          ecFee[counter[i][j][k]]     = ifee[i][j][k];
          ecCrate[counter[i][j][k]]   = icrate[i][j][k];
          ecBoard[counter[i][j][k]]   = iboard[i][j][k];
          ecChannel[counter[i][j][k]] = ichannel[i][j][k];
          for ( unsigned int l = 0; l < paramNum; l++ )
          { pendcap[counter[i][j][k]][l] = param[i][j][k][l]; }
        }
      }
    }
  }
 
  double   bRatio[2][barScin], eRatio[endScin];
  ifstream inRatio1, inRatio2, inRatio3;
  inRatio1.open( "BEastSimRatio.txt" );
  inRatio2.open( "BWestSimRatio.txt" );
  for ( unsigned int i = 0; i < barScin; i++ )
  {
    inRatio1 >> bRatio[0][i];
    inRatio2 >> bRatio[1][i];
  }
  inRatio1.close();
  inRatio2.close();
  inRatio3.open( "EndSimRatio.txt" );
  for ( unsigned int i = 0; i < endScin; i++ ) { inRatio3 >> eRatio[i]; }
  inRatio3.close();
 
  outData1.open( outFile1 );
  outData2.open( outFile2 );
  for ( unsigned int i = 0; i < barScin; i++ )
  {
    if ( paramNum == 8 )
    {
      outData1 << Fee[i][0] << " " << Crate[i][0] << " " << Board[i][0] << " " << Channel[i][0]
               << " ";
      outData2 << Fee[i][1] << " " << Crate[i][1] << " " << Board[i][1] << " " << Channel[i][1]
               << " ";
    }
    outData1 << bRatio[0][i];
    outData2 << bRatio[1][i];
    for ( unsigned int l = 0; l < paramNum; l++ )
    {
      outData1 << " " << pbeast[i][l];
      outData2 << " " << pbwest[i][l];
    }
    outData1 << endl;
    outData2 << endl;
  }
  outData1.close();
  outData2.close();
 
  outData3.open( outFile3 );
  for ( unsigned int i = 0; i < endScin; i++ )
  {
    if ( paramNum == 8 )
    {
      outData3 << ecFee[i] << " " << ecCrate[i] << " " << ecBoard[i] << " " << ecChannel[i]
               << " ";
    }
    outData3 << eRatio[i];
    for ( unsigned int l = 0; l < paramNum; l++ ) { outData3 << " " << pendcap[i][l]; }
    outData3 << endl;
  }
  outData3.close();
 
  return;
}