JVoltageOptimizer_utils.hh

Go to the documentation of this file.

#ifndef __OPT_UTILS__
#define __OPT_UTILS__
 
//c++ standard library
#include <iostream>
 
//root
#include "TH2D.h"
 
//Nanobeacons
#include "NBRun.hh"
 
// namespaces
using namespace std;
 
 
/**
 * \author rgruiz
 */
 
 
 
/**
 * Structure to store the ToT mean and standard deviation of the hits produced by a nanobeacon in a source and target modules. 
 */
struct event{
 
  double V = 999;
 
  double meanToT_ref = 999;
 
  double meanToT_tgt = 999 ;
 
  double sigmaToT_ref = 999 ;
 
  double sigmaToT_tgt = 999 ;
 
};
 
 
 
/**
 * Optimizes the voltage of each nanobeacon 
 *
 * \param max_tot maximum allowed ToT value
 * \param evt_table table of event structures for each nanobeacon and run
 * \return vector of pairs with optimal voltage and corresponding run number
 */
inline vector < pair < double , int > > optimize_voltage (vector < vector < event > > evt_table , double max_tot){
 
  int Nbeacons = evt_table.size() ;
 
  int Nruns = evt_table[0].size() ;
 
  vector < pair < double , int > > optimal_voltages (Nbeacons) ;
 
  for (int i=0 ; i<Nbeacons ; i++){
    
    for(int j=Nruns-1 ; j>=0 ; j--){
 
      if(evt_table[i][j].meanToT_ref < max_tot){
        
        cout << evt_table[i][j].meanToT_ref << "\t " << evt_table[i][j].V << endl ;
        
        optimal_voltages[i] = make_pair(evt_table[i][j].V , j);
 
        break ;
 
      }
      
    }
    
  }
 
  return optimal_voltages ;
  
}
 
 
#endif