//Created by Laszo Varga & Jan Glorius
//last update: 20.11.2020
//run in root:
//root -l
//.L script.c+
//run()
////////////////////////
// INPUT FILES
#define INPUT1 "/data.local3/test_data_2020/Si2_run002.lmd.root"
//OUTPUT FILE
#define OUTPUT "Si2_run002.ana.root"
//Energy thresholds for the silicon:
#define Si_Thrs_LOW 100
#define Si_Thrs_HIGH 8000
////////////////////////
#include <cmath>
#include <string>
#include <cstdlib>
#include <cstdio>
#include <ctime>
#include <fstream>
#include <iostream>
#include <stdint.h>
#include "TROOT.h"
#include "TAttText.h"
#include "TAxis.h"
#include "TCanvas.h"
#include "TChain.h"
#include "TCut.h"
#include "TF1.h"
#include "TFile.h"
#include "TGraph.h"
#include "TGraphAsymmErrors.h"
#include "TGraphErrors.h"
#include "TH1.h"
#include "TH2.h"
#include "THistPainter.h"
#include "TKey.h"
#include "TLatex.h"
#include "TLegend.h"
#include "TMath.h"
#include "TMatrixD.h"
#include "TMinuit.h"
#include "TMultiGraph.h"
#include "TNtuple.h"
#include "TPave.h"
#include "TPaveText.h"
#include "TPoint.h"
#include "TRandom.h"
#include "TRint.h"
#include "TStyle.h"
#include "TString.h"
#include "TTree.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TSystem.h"
#include "TProfile.h"
#include "TVirtualFitter.h"
#include "TCanvas.h"
#include "TLegend.h"
#include "TColor.h"
using namespace std;
void loop(TChain *fChain){
#include "/data.local3/test_data_2020/e127.h"
//creating histos
//counting of variables starts from 0 --> eg.: first Si strip is 0.
TH2D *h_pos_Si_xy_raw=new TH2D("h_pos_si_xy_raw", "h_pos_si_xy_raw",16,-0.5,15.5,16,-0.5,15.5);
TH2D *h_pos_Si_xy_noDC=new TH2D("h_pos_si_xy_noDC", "h_pos_si_xy_noDC",16,-0.5,15.5,16,-0.5,15.5);
TH2D *h_E_Si_all=new TH2D("h_E_Si_all", "E vs. strip Nr.; strip nr.; energy ",32,-0.5,31.5,8192,0.5,8191.5);
TH1D *h_E_Si_x[16];
TH1D *h_E_Si_y[16];
TH2D *h_E_ev_x[16];
TH2D *h_E_ev_y[16];
for(int a=0;a<16;a++){
h_E_Si_x[a]=new TH1D(Form("h_E_Si_x%d",a),Form("h_E_Si_x%d",a),8192,0.5,8191.5);
h_E_Si_y[a]=new TH1D(Form("h_E_Si_y%d",a),Form("h_E_Si_y%d",a),8192,0.5,8191.5);
h_E_ev_x[a]=new TH2D(Form("h_E_ev_x%d",a),Form("h_E_ev_x%d;ev;E",a),1024,0,55e6,4096,0.5,8191.5);
h_E_ev_y[a]=new TH2D(Form("h_E_ev_y%d",a),Form("h_E_ev_y%d;ev;E",a),1024,0,55e6,4096,0.5,8191.5);
}
////////////////////////////////////////////////////////////////////////
Long64_t nentries = fChain->GetEntries();
Long64_t nbytes = 0;
int dc_y_max=-999;
int dc_x_max=-999;
int dc_Ey_max=0;
int dc_Ex_max=0;
int dc_ty_max=0;
int dc_tx_max=0;
for (Long64_t i=0; i<nentries;i++){
//for (Long64_t i=0; i<100000;i++){
nbytes += fChain->GetEntry(i);
// event countdown
if ((float(i)/100000.)==int(i/100000)){cout << "event: " << i << " \tof " << nentries << endl;}
if(trigger==1){
for (int i_x=0;i_x<16;i_x++){
for (int i_y=0;i_y<16;i_y++){
if(int(E_Si_x[i_x])> 2000 && int(E_Si_y[i_y])> 2000)
{
h_pos_Si_xy_raw->Fill(i_x,i_y);
}
if(int(E_Si_x[i_x])>dc_Ex_max && int(E_Si_y[i_y])>dc_Ey_max){
dc_y_max = i_y;
dc_x_max = i_x;
dc_Ey_max = E_Si_y[i_y];
dc_Ex_max = E_Si_x[i_x];
dc_ty_max = t_Si_y[i_y];
dc_tx_max = t_Si_x[i_x];
}
}
//x-loop is missused for y also
if(t_Si_x[i_x] > 0){
h_E_Si_all->Fill(i_x,E_Si_x[i_x]);
h_E_Si_x[i_x]->Fill(E_Si_x[i_x]);
}
if(t_Si_y[i_x] > 0){
h_E_Si_all->Fill(i_x + 16,E_Si_y[i_x]);
h_E_Si_y[i_x]->Fill(E_Si_y[i_x]);
}
}
if(dc_x_max!=-999 && dc_y_max!=-999){
h_pos_Si_xy_noDC->Fill(dc_x_max,dc_y_max);
h_E_ev_x[dc_x_max]->Fill(i,dc_Ex_max);
h_E_ev_y[dc_y_max]->Fill(i,dc_Ey_max);
}
dc_y_max=-999;
dc_x_max=-999;
dc_Ex_max=0;
dc_Ey_max=0;
dc_tx_max=0;
dc_ty_max=0;
}
}//event loop
TFile *graphfile = TFile::Open(OUTPUT, "RECREATE");
//graphfile->mkdir("Si_map");
//graphfile->cd("Si_map");
h_pos_Si_xy_raw->Write();
h_pos_Si_xy_noDC->Write();
h_E_Si_all->Write();
//graphfile->mkdir("Si_E_histos");
//graphfile->cd("Si_E_histos");
for(int a=0;a<16;a++){
h_E_Si_x[a]->Write();
h_E_Si_y[a]->Write();
h_E_ev_x[a]->Write();
h_E_ev_y[a]->Write();
}
graphfile->Close();
cout << "\033[0;32m" << OUTPUT << " is created!\033[0m" << endl;
} // loop
void run(){
const char *command = new char[1000];
char filename[100];
TChain *fChain = new TChain("h101");
sprintf(filename,INPUT1);
cout<<"\033[0;37m//loading run: "<<filename << "\033[0m" <<endl;
fChain->Add(filename);
#ifdef INPUT2 // optional second input file
sprintf(filename,INPUT2);
cout<<"\033[0;37m//loading run: "<<filename << "\033[0m" <<endl;
fChain->Add(filename);
#endif
//sprintf(filename,OUTPUT);
// TFile *outfile = new TFile(filename,"RECREATE");
loop(fChain);
// outfile->Write();
command = "rm *.so";
gSystem->Exec(command);
command = "rm *.d";
gSystem->Exec(command);
command = "rm *.pcm";
gSystem->Exec(command);
return;
}
int main(){
run();
return(0);
}
|
//Declaration of leaves types
UInt_t tpat;
UInt_t trigger;
UInt_t Si_voltage;
UInt_t Si_current;
UInt_t E_Si_y[16];
UInt_t E_Si_x[16];
UInt_t E_Xray[3];
UInt_t E_BaF[6];
UInt_t t_Si_y[16];
UInt_t t_Si_x[16];
UInt_t t_Xray[3];
UInt_t t_BaF[6];
UInt_t v_mega_clock;
UInt_t v_Si_X;
UInt_t v_Si_Y;
UInt_t v_Xray[3];
UInt_t v_BaF[6];
UInt_t v_TRAFO;
UInt_t v_I_cool;
UInt_t v_U_cool;
UInt_t v_jet_S[1];
UInt_t bdt_si_ON;
UInt_t bdt_si_OFF;
UInt_t bdt_xray_ON;
UInt_t bdt_xray_OFF;
UInt_t bdt_baf_ON;
UInt_t bdt_baf_OFF;
UInt_t adt_si_ON;
UInt_t adt_si_OFF;
UInt_t adt_xray_ON;
UInt_t adt_xray_OFF;
UInt_t adt_baf_ON;
UInt_t adt_baf_OFF;
UInt_t red_si_ON;
UInt_t red_si_OFF;
UInt_t red_xray_ON;
UInt_t red_xray_OFF;
UInt_t red_baf_ON;
UInt_t red_baf_OFF;
// Set branch addresses.
fChain->SetBranchAddress("tpat",&tpat);
fChain->SetBranchAddress("TRIGGER",&trigger);
fChain->SetBranchAddress("Si_voltage",&Si_voltage);
fChain->SetBranchAddress("Si_current",&Si_current);
fChain->SetBranchAddress("E_Si_y1",&E_Si_y[0]);
fChain->SetBranchAddress("E_Si_y2",&E_Si_y[1]);
fChain->SetBranchAddress("E_Si_y3",&E_Si_y[2]);
fChain->SetBranchAddress("E_Si_y4",&E_Si_y[3]);
fChain->SetBranchAddress("E_Si_y5",&E_Si_y[4]);
fChain->SetBranchAddress("E_Si_y6",&E_Si_y[5]);
fChain->SetBranchAddress("E_Si_y7",&E_Si_y[6]);
fChain->SetBranchAddress("E_Si_y8",&E_Si_y[7]);
fChain->SetBranchAddress("E_Si_y9",&E_Si_y[8]);
fChain->SetBranchAddress("E_Si_y10",&E_Si_y[9]);
fChain->SetBranchAddress("E_Si_y11",&E_Si_y[10]);
fChain->SetBranchAddress("E_Si_y12",&E_Si_y[11]);
fChain->SetBranchAddress("E_Si_y13",&E_Si_y[12]);
fChain->SetBranchAddress("E_Si_y14",&E_Si_y[13]);
fChain->SetBranchAddress("E_Si_y15",&E_Si_y[14]);
fChain->SetBranchAddress("E_Si_y16",&E_Si_y[15]);
fChain->SetBranchAddress("E_Si_x1",&E_Si_x[0]);
fChain->SetBranchAddress("E_Si_x2",&E_Si_x[1]);
fChain->SetBranchAddress("E_Si_x3",&E_Si_x[2]);
fChain->SetBranchAddress("E_Si_x4",&E_Si_x[3]);
fChain->SetBranchAddress("E_Si_x5",&E_Si_x[4]);
fChain->SetBranchAddress("E_Si_x6",&E_Si_x[5]);
fChain->SetBranchAddress("E_Si_x7",&E_Si_x[6]);
fChain->SetBranchAddress("E_Si_x8",&E_Si_x[7]);
fChain->SetBranchAddress("E_Si_x9",&E_Si_x[8]);
fChain->SetBranchAddress("E_Si_x10",&E_Si_x[9]);
fChain->SetBranchAddress("E_Si_x11",&E_Si_x[10]);
fChain->SetBranchAddress("E_Si_x12",&E_Si_x[11]);
fChain->SetBranchAddress("E_Si_x13",&E_Si_x[12]);
fChain->SetBranchAddress("E_Si_x14",&E_Si_x[13]);
fChain->SetBranchAddress("E_Si_x15",&E_Si_x[14]);
fChain->SetBranchAddress("E_Si_x16",&E_Si_x[15]);
fChain->SetBranchAddress("E_Xray1",&E_Xray[0]);
fChain->SetBranchAddress("E_Xray2",&E_Xray[1]);
fChain->SetBranchAddress("E_Xray3",&E_Xray[2]);
fChain->SetBranchAddress("E_BaF1",&E_BaF[0]);
fChain->SetBranchAddress("E_BaF2",&E_BaF[1]);
fChain->SetBranchAddress("E_BaF3",&E_BaF[2]);
fChain->SetBranchAddress("E_BaF4",&E_BaF[3]);
fChain->SetBranchAddress("E_BaF5",&E_BaF[4]);
fChain->SetBranchAddress("E_BaF6",&E_BaF[5]);
fChain->SetBranchAddress("t_Si_y1",&t_Si_y[0]);
fChain->SetBranchAddress("t_Si_y2",&t_Si_y[1]);
fChain->SetBranchAddress("t_Si_y3",&t_Si_y[2]);
fChain->SetBranchAddress("t_Si_y4",&t_Si_y[3]);
fChain->SetBranchAddress("t_Si_y5",&t_Si_y[4]);
fChain->SetBranchAddress("t_Si_y6",&t_Si_y[5]);
fChain->SetBranchAddress("t_Si_y7",&t_Si_y[6]);
fChain->SetBranchAddress("t_Si_y8",&t_Si_y[7]);
fChain->SetBranchAddress("t_Si_y9",&t_Si_y[8]);
fChain->SetBranchAddress("t_Si_y10",&t_Si_y[9]);
fChain->SetBranchAddress("t_Si_y11",&t_Si_y[10]);
fChain->SetBranchAddress("t_Si_y12",&t_Si_y[11]);
fChain->SetBranchAddress("t_Si_y13",&t_Si_y[12]);
fChain->SetBranchAddress("t_Si_y14",&t_Si_y[13]);
fChain->SetBranchAddress("t_Si_y15",&t_Si_y[14]);
fChain->SetBranchAddress("t_Si_y16",&t_Si_y[15]);
fChain->SetBranchAddress("t_Si_x1",&t_Si_x[0]);
fChain->SetBranchAddress("t_Si_x2",&t_Si_x[1]);
fChain->SetBranchAddress("t_Si_x3",&t_Si_x[2]);
fChain->SetBranchAddress("t_Si_x4",&t_Si_x[3]);
fChain->SetBranchAddress("t_Si_x5",&t_Si_x[4]);
fChain->SetBranchAddress("t_Si_x6",&t_Si_x[5]);
fChain->SetBranchAddress("t_Si_x7",&t_Si_x[6]);
fChain->SetBranchAddress("t_Si_x8",&t_Si_x[7]);
fChain->SetBranchAddress("t_Si_x9",&t_Si_x[8]);
fChain->SetBranchAddress("t_Si_x10",&t_Si_x[9]);
fChain->SetBranchAddress("t_Si_x11",&t_Si_x[10]);
fChain->SetBranchAddress("t_Si_x12",&t_Si_x[11]);
fChain->SetBranchAddress("t_Si_x13",&t_Si_x[12]);
fChain->SetBranchAddress("t_Si_x14",&t_Si_x[13]);
fChain->SetBranchAddress("t_Si_x15",&t_Si_x[14]);
fChain->SetBranchAddress("t_Si_x16",&t_Si_x[15]);
fChain->SetBranchAddress("t_Xray1",&t_Xray[0]);
fChain->SetBranchAddress("t_Xray2",&t_Xray[1]);
fChain->SetBranchAddress("t_Xray3",&t_Xray[2]);
fChain->SetBranchAddress("t_BaF1",&t_BaF[0]);
fChain->SetBranchAddress("t_BaF2",&t_BaF[1]);
fChain->SetBranchAddress("t_BaF3",&t_BaF[2]);
fChain->SetBranchAddress("t_BaF4",&t_BaF[3]);
fChain->SetBranchAddress("t_BaF5",&t_BaF[4]);
fChain->SetBranchAddress("t_BaF6",&t_BaF[5]);
fChain->SetBranchAddress("v_mega_clock",&v_mega_clock);
fChain->SetBranchAddress("v_Si_X",&v_Si_X);
fChain->SetBranchAddress("v_Si_Y",&v_Si_Y);
fChain->SetBranchAddress("v_Xray1",&v_Xray[0]);
fChain->SetBranchAddress("v_Xray2",&v_Xray[1]);
fChain->SetBranchAddress("v_Xray3",&v_Xray[2]);
fChain->SetBranchAddress("v_BaF1",&v_BaF[0]);
fChain->SetBranchAddress("v_BaF2",&v_BaF[1]);
fChain->SetBranchAddress("v_BaF3",&v_BaF[2]);
fChain->SetBranchAddress("v_BaF4",&v_BaF[3]);
fChain->SetBranchAddress("v_BaF5",&v_BaF[4]);
fChain->SetBranchAddress("v_BaF6",&v_BaF[5]);
fChain->SetBranchAddress("v_TRAFO",&v_TRAFO);
fChain->SetBranchAddress("v_I_cool",&v_I_cool);
fChain->SetBranchAddress("v_U_cool",&v_U_cool);
fChain->SetBranchAddress("v_jet_S1",&v_jet_S[0]);
fChain->SetBranchAddress("bdt_si_ON",&bdt_si_ON);
fChain->SetBranchAddress("bdt_si_OFF",&bdt_si_OFF);
fChain->SetBranchAddress("bdt_xray_ON",&bdt_xray_ON);
fChain->SetBranchAddress("bdt_xray_OFF",&bdt_xray_OFF);
fChain->SetBranchAddress("bdt_baf_ON",&bdt_baf_ON);
fChain->SetBranchAddress("bdt_baf_OFF",&bdt_baf_OFF);
fChain->SetBranchAddress("adt_si_ON",&adt_si_ON);
fChain->SetBranchAddress("adt_si_OFF",&adt_si_OFF);
fChain->SetBranchAddress("adt_xray_ON",&adt_xray_ON);
fChain->SetBranchAddress("adt_xray_OFF",&adt_xray_OFF);
fChain->SetBranchAddress("adt_baf_ON",&adt_baf_ON);
fChain->SetBranchAddress("adt_baf_OFF",&adt_baf_OFF);
fChain->SetBranchAddress("red_si_ON",&red_si_ON);
fChain->SetBranchAddress("red_si_OFF",&red_si_OFF);
fChain->SetBranchAddress("red_xray_ON",&red_xray_ON);
fChain->SetBranchAddress("red_xray_OFF",&red_xray_OFF);
fChain->SetBranchAddress("red_baf_ON",&red_baf_ON);
fChain->SetBranchAddress("red_baf_OFF",&red_baf_OFF);
/*
fChain->SetBranchAddress("tpat",&tpat);
fChain->SetBranchAddress("trigger",&trigger);
fChain->SetBranchAddress("timestamp_lsb",×tamp_lsb);
fChain->SetBranchAddress("timestamp_msb",×tamp_msb);
fChain->SetBranchAddress("Tsc_before_lmu1",&Tsc_before_lmu[1]);
fChain->SetBranchAddress("Tsc_before_lmu1",&Tsc_before_lmu1);
fChain->SetBranchAddress("Tsc_before_lmu2",&Tsc_before_lmu2);
fChain->SetBranchAddress("Tsc_before_lmu3",&Tsc_before_lmu3);
fChain->SetBranchAddress("Tsc_before_lmu4",&Tsc_before_lmu4);
fChain->SetBranchAddress("Tsc_before_lmu5",&Tsc_before_lmu5);
fChain->SetBranchAddress("Tsc_before_lmu6",&Tsc_before_lmu6);
fChain->SetBranchAddress("Tsc_before_lmu7",&Tsc_before_lmu7);
fChain->SetBranchAddress("Tsc_before_lmu8",&Tsc_before_lmu8);
fChain->SetBranchAddress("Tsc_before_lmu9",&Tsc_before_lmu9);
fChain->SetBranchAddress("Tsc_before_lmu10",&Tsc_before_lmu10);
fChain->SetBranchAddress("Tsc_before_lmu11",&Tsc_before_lmu11);
fChain->SetBranchAddress("Tsc_before_lmu12",&Tsc_before_lmu12);
fChain->SetBranchAddress("Tsc_before_lmu13",&Tsc_before_lmu13);
fChain->SetBranchAddress("Tsc_before_lmu14",&Tsc_before_lmu14);
fChain->SetBranchAddress("Tsc_before_lmu15",&Tsc_before_lmu15);
fChain->SetBranchAddress("Tsc_before_lmu16",&Tsc_before_lmu16);
fChain->SetBranchAddress("Tsc_before_dt1",&Tsc_before_dt1);
fChain->SetBranchAddress("Tsc_before_dt2",&Tsc_before_dt2);
fChain->SetBranchAddress("Tsc_before_dt3",&Tsc_before_dt3);
fChain->SetBranchAddress("Tsc_before_dt4",&Tsc_before_dt4);
fChain->SetBranchAddress("Tsc_before_dt5",&Tsc_before_dt5);
fChain->SetBranchAddress("Tsc_before_dt6",&Tsc_before_dt6);
fChain->SetBranchAddress("Tsc_before_dt7",&Tsc_before_dt7);
fChain->SetBranchAddress("Tsc_before_dt8",&Tsc_before_dt8);
fChain->SetBranchAddress("Tsc_before_dt9",&Tsc_before_dt9);
fChain->SetBranchAddress("Tsc_before_dt10",&Tsc_before_dt10);
fChain->SetBranchAddress("Tsc_before_dt11",&Tsc_before_dt11);
fChain->SetBranchAddress("Tsc_before_dt12",&Tsc_before_dt12);
fChain->SetBranchAddress("Tsc_before_dt13",&Tsc_before_dt13);
fChain->SetBranchAddress("Tsc_before_dt14",&Tsc_before_dt14);
fChain->SetBranchAddress("Tsc_before_dt15",&Tsc_before_dt15);
fChain->SetBranchAddress("Tsc_before_dt16",&Tsc_before_dt16);
fChain->SetBranchAddress("Tsc_after_dt1",&Tsc_after_dt1);
fChain->SetBranchAddress("Tsc_after_dt2",&Tsc_after_dt2);
fChain->SetBranchAddress("Tsc_after_dt3",&Tsc_after_dt3);
fChain->SetBranchAddress("Tsc_after_dt4",&Tsc_after_dt4);
fChain->SetBranchAddress("Tsc_after_dt5",&Tsc_after_dt5);
fChain->SetBranchAddress("Tsc_after_dt6",&Tsc_after_dt6);
fChain->SetBranchAddress("Tsc_after_dt7",&Tsc_after_dt7);
fChain->SetBranchAddress("Tsc_after_dt8",&Tsc_after_dt8);
fChain->SetBranchAddress("Tsc_after_dt9",&Tsc_after_dt9);
fChain->SetBranchAddress("Tsc_after_dt10",&Tsc_after_dt10);
fChain->SetBranchAddress("Tsc_after_dt11",&Tsc_after_dt11);
fChain->SetBranchAddress("Tsc_after_dt12",&Tsc_after_dt12);
fChain->SetBranchAddress("Tsc_after_dt13",&Tsc_after_dt13);
fChain->SetBranchAddress("Tsc_after_dt14",&Tsc_after_dt14);
fChain->SetBranchAddress("Tsc_after_dt15",&Tsc_after_dt15);
fChain->SetBranchAddress("Tsc_after_dt16",&Tsc_after_dt16);
fChain->SetBranchAddress("E_Xray1",&E_Xray1);
fChain->SetBranchAddress("E_Xray2",&E_Xray2);
fChain->SetBranchAddress("E_Xray3",&E_Xray3);
fChain->SetBranchAddress("E_Xray4",&E_Xray4);
fChain->SetBranchAddress("E_Xray5",&E_Xray5);
fChain->SetBranchAddress("E_si_y1",&E_si_y[1]);
fChain->SetBranchAddress("E_si_y2",&E_si_y[2]);
fChain->SetBranchAddress("E_si_y3",&E_si_y[3]);
fChain->SetBranchAddress("E_si_y4",&E_si_y[4]);
fChain->SetBranchAddress("E_si_y5",&E_si_y[5]);
fChain->SetBranchAddress("E_si_y6",&E_si_y[6]);
fChain->SetBranchAddress("E_si_y7",&E_si_y[7]);
fChain->SetBranchAddress("E_si_y8",&E_si_y[8]);
fChain->SetBranchAddress("E_si_y9",&E_si_y[9]);
fChain->SetBranchAddress("E_si_y10",&E_si_y[10]);
fChain->SetBranchAddress("E_si_y11",&E_si_y[11]);
fChain->SetBranchAddress("E_si_y12",&E_si_y[12]);
fChain->SetBranchAddress("E_si_y13",&E_si_y[13]);
fChain->SetBranchAddress("E_si_y14",&E_si_y[14]);
fChain->SetBranchAddress("E_si_y15",&E_si_y[15]);
fChain->SetBranchAddress("E_si_y16",&E_si_y[16]);
fChain->SetBranchAddress("E_si_x1",&E_si_x[1]);
fChain->SetBranchAddress("E_si_x2",&E_si_x[2]);
fChain->SetBranchAddress("E_si_x3",&E_si_x[3]);
fChain->SetBranchAddress("E_si_x4",&E_si_x[4]);
fChain->SetBranchAddress("E_si_x5",&E_si_x[5]);
fChain->SetBranchAddress("E_si_x6",&E_si_x[6]);
fChain->SetBranchAddress("E_si_x7",&E_si_x[7]);
fChain->SetBranchAddress("E_si_x8",&E_si_x[8]);
fChain->SetBranchAddress("E_si_x9",&E_si_x[9]);
fChain->SetBranchAddress("E_si_x10",&E_si_x[10]);
fChain->SetBranchAddress("E_si_x11",&E_si_x[11]);
fChain->SetBranchAddress("E_si_x12",&E_si_x[12]);
fChain->SetBranchAddress("E_si_x13",&E_si_x[13]);
fChain->SetBranchAddress("E_si_x14",&E_si_x[14]);
fChain->SetBranchAddress("E_si_x15",&E_si_x[15]);
fChain->SetBranchAddress("E_si_x16",&E_si_x[16]);
fChain->SetBranchAddress("t_Xray1",&t_Xray1);
fChain->SetBranchAddress("t_Xray2",&t_Xray2);
fChain->SetBranchAddress("t_Xray3",&t_Xray3);
fChain->SetBranchAddress("t_Xray4",&t_Xray4);
fChain->SetBranchAddress("t_Xray5",&t_Xray5);
fChain->SetBranchAddress("t_MWPC_anode",&t_MWPC_anode);
fChain->SetBranchAddress("t_MWPC_x1",&t_MWPC_x1);
fChain->SetBranchAddress("t_MWPC_x2",&t_MWPC_x2);
fChain->SetBranchAddress("t_MWPC_y1",&t_MWPC_y1);
fChain->SetBranchAddress("t_MWPC_y2",&t_MWPC_y2);
fChain->SetBranchAddress("t_si_y1",&t_si_y[1]);
fChain->SetBranchAddress("t_si_y2",&t_si_y[2]);
fChain->SetBranchAddress("t_si_y3",&t_si_y[3]);
fChain->SetBranchAddress("t_si_y4",&t_si_y[4]);
fChain->SetBranchAddress("t_si_y5",&t_si_y[5]);
fChain->SetBranchAddress("t_si_y6",&t_si_y[6]);
fChain->SetBranchAddress("t_si_y7",&t_si_y[7]);
fChain->SetBranchAddress("t_si_y8",&t_si_y[8]);
fChain->SetBranchAddress("t_si_y9",&t_si_y[9]);
fChain->SetBranchAddress("t_si_y10",&t_si_y[10]);
fChain->SetBranchAddress("t_si_y11",&t_si_y[11]);
fChain->SetBranchAddress("t_si_y12",&t_si_y[12]);
fChain->SetBranchAddress("t_si_y13",&t_si_y[13]);
fChain->SetBranchAddress("t_si_y14",&t_si_y[14]);
fChain->SetBranchAddress("t_si_y15",&t_si_y[15]);
fChain->SetBranchAddress("t_si_y16",&t_si_y[16]);
fChain->SetBranchAddress("t_si_x1",&t_si_x[1]);
fChain->SetBranchAddress("t_si_x2",&t_si_x[2]);
fChain->SetBranchAddress("t_si_x3",&t_si_x[3]);
fChain->SetBranchAddress("t_si_x4",&t_si_x[4]);
fChain->SetBranchAddress("t_si_x5",&t_si_x[5]);
fChain->SetBranchAddress("t_si_x6",&t_si_x[6]);
fChain->SetBranchAddress("t_si_x7",&t_si_x[7]);
fChain->SetBranchAddress("t_si_x8",&t_si_x[8]);
fChain->SetBranchAddress("t_si_x9",&t_si_x[9]);
fChain->SetBranchAddress("t_si_x10",&t_si_x[10]);
fChain->SetBranchAddress("t_si_x11",&t_si_x[11]);
fChain->SetBranchAddress("t_si_x12",&t_si_x[12]);
fChain->SetBranchAddress("t_si_x13",&t_si_x[13]);
fChain->SetBranchAddress("t_si_x14",&t_si_x[14]);
fChain->SetBranchAddress("t_si_x15",&t_si_x[15]);
fChain->SetBranchAddress("t_si_x16",&t_si_x[16]);
fChain->SetBranchAddress("sc_TRIDI_clk",&sc[1]);
fChain->SetBranchAddress("sc_ext_clk",&sc[2]);
fChain->SetBranchAddress("SC3",&sc[3]);
fChain->SetBranchAddress("SC4",&sc[4]);
fChain->SetBranchAddress("SC5",&sc[5]);
fChain->SetBranchAddress("SC6",&sc[6]);
fChain->SetBranchAddress("SC7",&sc[7]);
fChain->SetBranchAddress("SC8",&sc[8]);
fChain->SetBranchAddress("sc_esr_trafo",&sc[9]);
fChain->SetBranchAddress("sc_I_cooler",&sc[10]);
fChain->SetBranchAddress("sc_U_cooler",&sc[11]);
fChain->SetBranchAddress("sc_jet_E1",&sc[12]);
fChain->SetBranchAddress("sc_jet_S1",&sc[13]);
fChain->SetBranchAddress("sc_jet_S2",&sc[14]);
fChain->SetBranchAddress("sc_jet_S3",&sc[15]);
fChain->SetBranchAddress("sc_jet_S4",&sc[16]);
fChain->SetBranchAddress("sc_si_y",&sc[17]);
fChain->SetBranchAddress("sc_si_x",&sc[18]);
fChain->SetBranchAddress("sc_xray1",&sc[19]);
fChain->SetBranchAddress("sc_xray2",&sc[20]);
fChain->SetBranchAddress("sc_xray3",&sc[21]);
fChain->SetBranchAddress("sc_xray4",&sc[22]);
fChain->SetBranchAddress("sc_mwpc",&sc[23]);
fChain->SetBranchAddress("sc_xray5",&sc[24]);
fChain->SetBranchAddress("sc_mwpcCFD",&sc[25]);
fChain->SetBranchAddress("SC26",&sc[26]);
fChain->SetBranchAddress("SC27",&sc[27]);
fChain->SetBranchAddress("SC28",&sc[28]);
fChain->SetBranchAddress("SC29",&sc[29]);
fChain->SetBranchAddress("SC30",&sc[30]);
fChain->SetBranchAddress("SC31",&sc[31]);
fChain->SetBranchAddress("SC32",&sc[32]);
*/
|