latest and greatest

Ana_tests.C

@@ -0,0 +1,213 @@
+#include "Ana_tests.h"
+#include "registers_config.h"
+#include <fstream>
+#include <cmath>
+#include <cstring>
+
+
+// CONSTRUCTOR
+Ana_tests::Ana_tests(Fpga *fpga, Salt *salt, FastComm *fastComm) {
+
+  fpga_=fpga;
+  salt_=salt;
+  fastComm_=fastComm;
+
+}
+
+void Ana_tests::Get_run(string option, int runs, bool output, string outText) {
+
+
+
+  // Define command length (will be 2 in this case)
+  uint8_t length = 90;
+  // Define command list (BXID and Sync)
+  uint8_t command[256]={0};
+  // Read out data packet
+  uint16_t length_read = 150; // number of clock cycles to read
+  //uint32_t data[5120]; // data packet
+  //uint8_t delay = 0; // clock delay
+
+  ofstream outfile;
+  outfile.open (("noise_"+option+"_"+outText+".txt").c_str());
+  float avg_noise = 0;
+  int period = length;
+  bool singleShot = true;
+  string data_string;
+  float avg_ADC[128] = {0};
+  int ADC[128] = {0};
+  int bxid = 0;
+  int parity = 0;
+  int mcm_v = 0;
+  int mcm_ch = 0;
+  int mem_space = 0;
+  int flag = 0;
+  //int runs = 1;
+  int length1= 0;
+  unsigned twelveBits;
+  //int largest=0;
+  length_read = 255;
+  float avg_chip[runs] = {0};
+  float std_dev = 0;
+
+
+  // cout << "option = " << option << endl;
+
+
+  for(int i=0; i<79; i++)
+    command[i]=0x04;
+
+  if(option == "NZS") command[79]=0x08;
+  else if(option == "Normal") command[79]=0x00;
+  else {
+    cout << "ERROR: NZS or Normal packet not properly defined." << endl;
+
+  }
+
+  double v[runs][128] = {0};
+  period = length;
+  float length_avg = 0;
+
+  for(int i = 0; i < runs; i++) {
+    fastComm_->Take_a_run(length_read, data_string, length, 0, command, period, singleShot, true );     
+    //    cout << "data string is " << data_string << endl;
+    // read 12 bits at a time
+    for(int j=0; j<data_string.length(); j+=3) {
+      twelveBits = fastComm_->read_twelveBits(data_string, j);
+      fastComm_->read_Header(twelveBits, bxid, parity, flag, length1);      
+
+
+
+      if(flag==0) {
+	// cout << "bxid = " << bxid << endl;
+	// cout << "parity = " << parity << endl;
+	// cout << "flag = " << flag << endl;
+	// cout << "length = " << dec << length1 << endl;
+	 length_avg+=length1/((float) runs);
+	if(length1==0) break;
+
+	fastComm_->read_Normal_packet(data_string, j, ADC);
+	break;
+      }
+
+
+      if(length1==6) {
+	// cout << "bxid = " << bxid << endl;
+	// cout << "parity = " << parity << endl;
+	// cout << "flag = " << flag << endl;
+	//	 cout << "length = " << dec << length1 << endl;
+	 length_avg+=length1/((float) runs);
+	fastComm_->read_NZS_packet(data_string, j, ADC, bxid, parity, mcm_v, mcm_ch, mem_space);
+	break;
+      }
+
+    }
+
+    for(int j=0; j<128; j++) {
+      v[i][j]+=ADC[j]/((float)runs);
+      avg_ADC[j]+=ADC[j]/((float)runs);
+      avg_chip[i]+=ADC[j]/(128.);
+	ADC[j]=0;
+      }     
+
+  }
+
+  cout << length_avg << endl;
+
+  outfile << "length_avg = " << dec << length_avg << endl;
+
+  for(int j = 0; j<runs; j++) avg_noise+=avg_chip[j]/((float) runs);
+
+  outfile << "runs = " << dec << runs << ", Average = " << avg_noise << ", std_dev = " << calculateSD(avg_chip)<<endl;
+
+
+
+    for(int j=0; j<128; j++) {
+      outfile << dec << avg_ADC[j] << endl;
+    }
+
+
+   outfile.close();
+
+  return;  
+}
+
+
+
+
+void Ana_tests::Check_noise() {
+
+    string option = "NZS";
+    string outText;
+    uint8_t buffer=0x00;
+    cout << "PERFORMING NOISE CHECK" << endl;
+
+    // DSP output
+    salt_->write_salt(registers::ser_source_cfg,(uint8_t) 0x20);
+    salt_->write_salt(registers::n_zs_cfg,buffer);
+
+    // read back n_zs_cfg to make proper outfile name
+    salt_->read_salt(registers::n_zs_cfg, &buffer);
+
+    if((buffer & 0b01100000) == 0) outText = "noise_masked_ADC";
+    else if((buffer & 0b01100000) == 0b0010000) outText = "noise_sync_ADC";
+    else if((buffer & 0b01100000) == 0b0100000) outText = "noise_after_ped";
+    else if((buffer & 0b01100000) == 0b0110000) outText = "noise_after_mcm";
+    else {
+
+      cout <<"ERROR::n_zs_cfg not properly defined" << endl;
+      return;
+
+    }
+
+    // Do Normal packet
+    Get_run("Normal",1000,true,outText);
+
+    // Do NZS packet
+    Get_run("NZS",1000,true,outText);
+
+
+
+}
+
+
+void Ana_tests::Check_NZS() {
+
+  cout << "DOING NZS TEST" << endl;
+
+  stringstream outText;
+
+
+  for(int i = 0; i < 32; i++) {
+
+    salt_->write_salt(registers::n_zs_cfg,(uint8_t) i);
+
+    outText << "nzs_th_" << dec << i;
+
+    cout << i << ", ";
+
+    Get_run("Normal",100,true,outText.str());
+
+    outText.str("");
+  }
+
+}
+
+float Ana_tests::calculateSD(float data[]) {
+
+  float sum = 0.0, mean, standardDeviation = 0.0;
+
+  int i;
+
+  for(i = 0; i < 10; ++i)
+    {
+      sum += data[i];
+    }
+
+  mean = sum/10;
+
+  for(i = 0; i < 10; ++i)
+    standardDeviation += pow(data[i] - mean, 2);
+
+  return sqrt(standardDeviation / 10);
+}
+

Ana_tests.h

@@ -0,0 +1,34 @@
+#include <stddef.h>
+#include <stdio.h>
+#include <cstring>
+#include <cstdlib>
+#include <ctime>
+#include <sstream>
+#include "Salt.h"
+#include "Fpga.h"
+#include "fastComm.h"
+#include <fstream>
+
+using namespace std;
+
+class Ana_tests {
+
+ public:
+
+  Ana_tests(){};
+  Ana_tests(Fpga*, Salt*, FastComm*);
+  ~Ana_tests(){};
+
+  void Get_run(string option, int runs, bool output, string outText);
+  void Trim_dac_scan();
+  void Check_noise();
+  void Check_NZS();
+  float calculateSD(float data[]);
+
+
+  private:
+  Fpga* fpga_;
+  Salt* salt_;
+  FastComm* fastComm_;
+
+};