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