#include <iostream>
#include <iomanip>
#include <vector>
#include <fstream>

using namespace std;

int Resolution_formula()
{
TFile *CRC = new TFile("CRC-Run065_Icomp100_Th15_MSPL1_lat32-32_I781uA_I771uA_asynchrous-121114_075331-0_ZeroOffSet.root","READ");
TFile *RECO = new TFile("RECO-CRC-Run065_Icomp100_Th15_MSPL1_lat32-32_I781uA_I771uA_asynchrous-121114_075331-0_ZeroOffSet.root","READ");

TTree *rd51tbgeo = (TTree*)CRC->Get("rd51tbgeo");
Float_t zpos_cms, xpos_cms,g1xcl_geoposX;

// populate the defined variables
rd51tbgeo->SetBranchAddress("g1xcl.geoposX",g1xcl_geoposX);
Long64_t nevent=rd51tbgeo->GetEntries();

cout<<"####### n = "<<nevent<<endl;

for (Int_t i = 0; i<nevent; i++){
	rd51tbgeo->GetEntry(i);
	cout<<"xpos_cms = "<<g1xcl_geoposX[i]<<endl;
}

cout<<"####### n = "<<nevent<<endl;

return 0;
}