/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Tof/TofQElecSvc/TofQElecSvc-00-00-05/test/getConstants.cxx

Go to the documentation of this file.

void getConstants() {
  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   = 8;
  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[10]=".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( "BEastRatio.txt" );
  inRatio2.open( "BWestRatio.txt" );
  for( unsigned int i=0; i<barScin; i++ ) {
    inRatio1 >> bRatio[0][i];
    inRatio2 >> bRatio[1][i];
  }
  inRatio1.close();
  inRatio2.close();
  inRatio3.open( "EndRatio.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;
}

---