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