JInterDomCal_IO.hh File Reference

#include <iostream>
#include <tuple>
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JModuleLocation.hh"
#include "JDetector/JPMTRouter.hh"
#include "JLang/JObjectID.hh"
#include "Detector.hh"
#include "Control_utils.hh"
#include "JDetector/JModuleRouter.hh"
#include "JDetector/JDetectorToolkit.hh"

Go to the source code of this file.

Classes
struct	IO
	Structure to store the different command line arguments for JInterDomCal. More...

Functions
int	read_options (IO &options, int argc, char **argv)
	Parses the command line options and fills an IO structure with them. More...
void	write_output_det (string filename, const JDetector &detector)
	Writes a .detx file with a JDetector. More...
void	write_output_compare (string filename, const JDetector &old_detector, const JDetector &new_detector, int string_number, int ref_pmt, int tgt_pmt)
	Writes a .txt with the comparison of the t0s in a given DU for two JDetectors. More...
void	write_output_checks (string filename, NBRun &run)
	Writes a .root file with some objects that can be used to check the calibration. More...

Function Documentation

read_options()

int read_options ( IO &  options, int  argc, char **  argv  )

inline

Parses the command line options and fills an IO structure with them.

Parameters

options	an option structure
argc	the number of command line arguments
argv	the command line arguments

Returns

1 if works 2 if it doesn't work

Definition at line 82 of file JInterDomCal_IO.hh.

                                                           {
 
  int a = 0 ;
 
  try {
 
    JParser<string> zap;
 
    zap["c"] = make_field(options.ofname_checks) = "" ;
 
    zap["x"] = make_field(options.ofname_detx) = "out.detx" ;
 
    zap["t"] = make_field(options.ofname_txt) = "" ;
 
    zap["f"] = make_field(options.ifname) ; 
 
    zap["a"] = make_field(options.detector_file) ;
 
    zap["s"] = make_field(options.string_number) ;
 
    zap["u"] = make_field(options.up_pmts) = 1 ;
 
    zap["d"] = make_field(options.down_pmts) = 4 ;
 
    zap["n"] = make_field(options.number_neighbors) = 2 ;
 
    zap["l"] = make_field(options.analysis_level) = 0 ;
 
    zap["rpm"] = make_field(options.ref_pmt) = 0 ;
 
    zap["tpm"] = make_field(options.tgt_pmt) = 22 ;
 
    if (zap.read(argc, argv) != 0)
 
      a = 1;
 
  }
 
  catch(const exception &error) {
 
    ERROR(error.what() << endl);
 
    a = 2;
 
  }
 
  return a ;
 
}

write_output_det()

void write_output_det ( string  filename, const JDetector &  detector  )

inline

Writes a .detx file with a JDetector.

Parameters

filename	path to the name of the .detx
detector	the detector to be stored

Definition at line 139 of file JInterDomCal_IO.hh.

                                                                         {
 
  store(filename , detector) ;
 
  cout << endl << "New detector stored in file: " << endl ;
 
  cout << filename << endl ;
 
  cout << "------------------------------------" << endl << endl ;
 
}

write_output_compare()

void write_output_compare ( string  filename, const JDetector &  old_detector, const JDetector &  new_detector, int  string_number, int  ref_pmt, int  tgt_pmt  )

inline

Writes a .txt with the comparison of the t0s in a given DU for two JDetectors.

Parameters

filename	path to the name of the .txt
old_detector	one JDetector
new_detector	the other JDetector
string_number	the number of the DU to compare
ref_pmt	pmt in upper hemisphere used for the comparison
tgt_pmt	pmt in lower hemisphere used for the comparsion

Definition at line 163 of file JInterDomCal_IO.hh.

                                                                                                                                                                {
 
  if (filename.length()==0){
 
    cout << "Omitting comparison with old detector" << endl ;
 
    return ;
 
  }else{
 
    cout << "Writing comparison between new and old detx in file: " << endl ;
 
    cout << filename << endl ;
 
    cout << "------------------------------------" << endl << endl ;
 
    vector<double> t0_diff_new = getDetFile_t0_differences(string_number , new_detector , ref_pmt , tgt_pmt) ;
    
    vector<double> t0_diff_old = getDetFile_t0_differences(string_number , old_detector , ref_pmt , tgt_pmt) ;
 
    vector < tuple <double , double , double> > t0_vs_depth_new = get_t0_offsets_vs_depth(string_number , new_detector) ;
 
    vector < tuple <double , double , double> > t0_vs_depth_old = get_t0_offsets_vs_depth(string_number , old_detector) ;
    
    ofstream txtfile;
 
    txtfile.open (filename);
 
    txtfile << "# Comparison between detector files: column 1 = Delta_t0 new detx , column 2 = Delta_t0 old detx , column 3 = column1 - column 2" << endl;
  
    for(int i=0 ; i<(int)t0_diff_new.size() ; i++){
      
      txtfile << t0_diff_new[i] << "\t" << t0_diff_old[i] << "\t" << t0_diff_new[i] - t0_diff_old[i] << endl;
      
    }
  
    txtfile << "\n# t0 vs depth: column 1 = floor , columns 2,3,4 = depth_new, mean_centered_t0_new , error_new , columns 5,6,7 = depth_old, mean_centered_t0_old , error_old , column 8 = t0_old - t0_new " << endl ;
 
    for (int i=0 ; i < (int)t0_vs_depth_new.size() ; i++){
 
      txtfile << i+1 << "\t" 
              << get<0>(t0_vs_depth_new[i]) << "\t" 
              << get<1>(t0_vs_depth_new[i]) << "\t" 
              << get<2>(t0_vs_depth_new[i]) << "\t" 
              << get<0>(t0_vs_depth_old[i]) << "\t" 
              << get<1>(t0_vs_depth_old[i]) << "\t" 
              << get<2>(t0_vs_depth_old[i]) << "\t"
              << get<1>(t0_vs_depth_new[i]) - get<1>(t0_vs_depth_old[i]) << endl ; 
    }   
 
  }
 
}

write_output_checks()

void write_output_checks ( string  filename, NBRun &  run  )

inline

Writes a .root file with some objects that can be used to check the calibration.

Parameters

filename	path to the name of the .root
run	The nanobeacon calibration run

Definition at line 225 of file JInterDomCal_IO.hh.

                                                             {
 
  if (filename.length()==0){
 
    cout << "Omitting calibration checks..." << endl ;
 
    return ;
 
  }else{
 
    cout << "Writing calibration checks in file: " << endl ;
 
    cout << filename << endl ;
 
    cout << "------------------------------------" << endl << endl ;
 
    TFile outfile(filename.c_str() , "recreate") ;
    
    outfile.cd() ;
    
    TH2D* srcs = srcs_th2(&run) ;
    TH2D* tgts = tgts_th2(&run) ;
    TH2D* refs = good_refs_th2(&run) ;
    TH2D* tgt_pmts = good_tgt_pmts_th2(&run) ;
 
    srcs->Write() ;
    tgts->Write() ;
    refs->Write() ;
    tgt_pmts->Write();
 
    outfile.mkdir("REF/Good") ;
    outfile.mkdir("REF/Weak") ;
    outfile.mkdir("REF/Saturated") ;
    
    outfile.mkdir("TGT/Good") ;
    outfile.mkdir("TGT/Weak") ;
    outfile.mkdir("TGT/Saturated") ;
 
    vector <SuperModule*> SuperMods = run.getSuperModules() ;
 
      
    for (auto & sm : SuperMods){
 
      for (auto & spm : sm->get_ref_pmts()){
   
        if (spm->getNBPulse()->IsGood()==true) outfile.cd("REF/Good") ;
      
        if (spm->getNBPulse()->IsSaturatedHit()==true) outfile.cd("REF/Saturated") ;
      
        if (spm->getNBPulse()->IsWeak()==true) outfile.cd("REF/Weak") ;
      
        spm->getNBPulse()->getHtime_full()->Write() ;
      
        if (spm->getNBPulse()->IsFitted()==true){
        
          RooWorkspace w = spm->getNBPulse()->getWorkspace() ;
        
          TH1D* h = spm->getNBPulse()->getHtime_full() ;
        
          TCanvas* c1 = RooCanvas(w , h) ;
        
          c1->Write() ;
        
          w.Write() ;
        
        }
      
        outfile.cd() ;
      
      }
             
      for (auto & sms : sm->get_sources()){
 
        for (auto & spm : sms.second){
   
          if (spm->getNBPulse()->IsGood()==true)  outfile.cd("TGT/Good") ;
            
          if (spm->getNBPulse()->IsSaturatedHit()==true) outfile.cd("TGT/Saturated") ;
            
          if (spm->getNBPulse()->IsWeak()==true) outfile.cd("TGT/Weak") ;
        
          spm->getNBPulse()->getHtime_full()->Write() ;
        
          if (spm->getNBPulse()->IsFitted()==true){
          
            RooWorkspace w = spm->getNBPulse()->getWorkspace() ;
 
            TH1D* h = spm->getNBPulse()->getHtime_full() ;
 
            TCanvas* c1 = RooCanvas(w , h) ;
 
            c1->Write() ;
 
            w.Write() ;
 
          }
        
          outfile.cd() ;
 
        }
 
      }
 
    }
 
    outfile.Close() ;
 
  }
 
}