/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/TB2010/Testbeam_june_10/Complete_detector/Analysis_tb_m2.C

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <map>
#include "TKey.h"
#include <TROOT.h>
#include <TRint.h>
#include <TKey.h>
#include <TTree.h>
#include <TFile.h>
#include <TSystem.h>
#include <TreeClass.hh>
#include <TF1.h>
#include <TH1I.h>
#include <TH2I.h>
#include <TH3F.h>
#include <TH2F.h>
#include <TGraph.h>
#include <sstream>
#include <vector>
#include "Log.hh"

#include "MicroException.hh"


using namespace std;


int main(int argc, char**argv){


  bool verbose_steerfile = 0;
  bool verbose_runloop = 0;
  bool verbose_ped_file = 0;
  bool verbose_offset = 1;

  int nsigma_adc_cut = 2;



  if (argc!=2){

    cout<<endl;
    cout<<"Number of argument invalid"<<endl;
    cout<<" 1 -> steerfile"<<endl;
    cout<<endl;}

  else{



    /***************Read steerfile***************************************/

    TString ped_gas_string;
    TString ped_m2_string;
    TString data_string;

    int i = 0;
    int j = 0;
    int run = 0;
    int nrun = 0;

    vector <TString> ped_gas_file_vector;
    vector <TString> ped_m2_file_vector;
    vector <TString> data_file_vector;

    ifstream steerfile(argv[1]);

    while(steerfile>>ped_gas_string>>ped_m2_string>>data_string){

      if (ped_gas_string=="end_fo_file"){break;}

      else{

        nrun++;
        ped_gas_file_vector.push_back(ped_gas_string);
        ped_m2_file_vector.push_back(ped_m2_string);
        data_file_vector.push_back(data_string);}}

    if (verbose_steerfile){

      cout<<endl;
      cout<<"Number of runs = "<<nrun<<endl;
      cout<<endl;
      for (i=0;i<nrun;i++){
        cout<<"Pedestal gassiplex file    = "<<ped_gas_file_vector[i]<<endl;
        cout<<"Pedestal m2 prototype file = "<<ped_m2_file_vector[i]<<endl;
        cout<<"ROOT data file             = "<<data_file_vector[i]<<endl;
        cout<<endl;}
      cout<<endl;}




    /*Pedestal file variables*/

    int chn,chb,xpos,ypos;
    float ped,ped_err,width,width_err,gain;

    /*M2 chip settings*/

    int nchn_m2 = 9216;
    float ped_m2_tab[nchn_m2];
    float ped_err_m2_tab[nchn_m2];
    float width_m2_tab[nchn_m2];
    float width_err_m2_tab[nchn_m2];
    float gain_m2_tab[nchn_m2];

    /*Gassiplex chip settings*/

    int nchb_gas = 4;
    int nchn_gas = 384;
    float ped_gas_tab[nchb_gas][nchn_gas];
    float ped_err_gas_tab[nchb_gas][nchn_gas];
    float width_gas_tab[nchb_gas][nchn_gas];
    float width_err_gas_tab[nchb_gas][nchn_gas];
    float thr_gas_tab[nchb_gas][nchn_gas];
    float gain_gas_tab[nchb_gas][nchn_gas];





    /*******************Loop over runs******************************************/

    for (run=0;run<1;run++){

      if (verbose_runloop){

        cout<<endl;
        cout<<"Analysis of run "<<run+1<<endl;
        cout<<endl;
        cout<<"Pedestal gassiplex file    = "<<ped_gas_file_vector[run]<<endl;
        cout<<"Pedestal m2 prototype file = "<<ped_m2_file_vector[run]<<endl;
        cout<<"ROOT data file             = "<<data_file_vector[run]<<endl;
        cout<<endl;}




      /*************************Read pedestal files***************************/

      /*Gassiplex chip settings*/

      for (i=0;i<nchb_gas;i++){
        for (j=0;j<nchn_gas;j++){

          ped_gas_tab[i][j] = 0;
          ped_err_gas_tab[i][j] = 0;
          width_gas_tab[i][j] = 0;
          width_err_gas_tab[i][j] = 0;
          gain_gas_tab[i][j] = 0;
          thr_gas_tab[i][j] = 0;}}

      /*M2 chip settings*/

      for (i=0;i<nchn_m2;i++){

        ped_m2_tab[i] = 0;
        ped_err_m2_tab[i] = 0;
        width_m2_tab[i] = 0;
        width_err_m2_tab[i] = 0;
        gain_m2_tab[i] = 0;}

      ifstream ped_gas_file(ped_gas_file_vector[run].Data());
      ifstream ped_m2_file(ped_m2_file_vector[run].Data());

      if(verbose_ped_file){cout<<"Read gassiplex pedestal file "<<ped_gas_file_vector[run]<<endl;}
      while(ped_gas_file>>chn>>xpos>>ypos>>chb>>ped>>ped_err>>width>>width_err>>gain){

        if (gain!=0){

          if (verbose_ped_file){cout<<chn<<"\t"<<xpos<<"\t"<<ypos<<"\t"<<chb<<"\t"<<ped<<"\t"<<ped_err<<"\t"<<width<<"\t"<<width_err<<"\t"<<gain<<"\r"<<flush;}

          ped_gas_tab[chb-1][chn] = ped;
          ped_err_gas_tab[chb-1][chn] = ped_err;
          width_gas_tab[chb-1][chn] = width;
          width_err_gas_tab[chb-1][chn] = width_err;
          gain_gas_tab[chb-1][chn] = gain;
          thr_gas_tab[chb-1][chn] = floor(ped + nsigma_adc_cut * width + 1);}}

      if(verbose_ped_file){cout<<endl;cout<<"Read m2 prototype pedestal file "<<ped_m2_file_vector[run]<<endl;}
      while(ped_m2_file>>chn>>xpos>>ypos>>chb>>ped>>ped_err>>width>>width_err>>gain){

        if (gain!=0){

          if (verbose_ped_file){cout<<chn<<"\t"<<xpos<<"\t"<<ypos<<"\t"<<chb<<"\t"<<ped<<"\t"<<ped_err<<"\t"<<width<<"\t"<<width_err<<"\t"<<gain<<"\r"<<flush;}

          ped_m2_tab[chn] = ped;
          ped_err_m2_tab[chn] = ped_err;
          width_m2_tab[chn] = width;
          width_err_m2_tab[chn] = width_err;
          gain_m2_tab[chn] = gain;}}

      if(verbose_ped_file){cout<<endl;}






      /*******************Measure offset between telescope and m2 prototype**********************/

      int ntrigger = 0;

      TFile * rootfile = new TFile(data_file_vector[run]);

      if (verbose_offset){cout<<endl;cout<<"Reading ROOT file "<<data_file_vector[run]<<endl;}

      TIter nextkey(rootfile->GetListOfKeys());
      TKey *key;

      while (key = (TKey*)nextkey()) {
    
        TTree *tree = (TTree*)key->ReadObj();                
        MTRun* run = (MTRun*)tree->GetUserInfo()->FindObject("MTRun");
        
        ntrigger = tree->GetEntries();

        MTEvent *evt = new MTEvent();
        TBranch *branch= tree->GetBranch("MTEvent");
        branch->SetAddress(&evt);

        if (verbose_offset){cout<<"Number of recorded triggers = "<<ntrigger<<endl;cout<<endl;}
      }








    }












  }


return 0;}

---