EXP1803

{

        TFile* f = new TFile("sim_digi.root");

           TTree* t = (TTree*)f->Get("er");

           t->AddFriend("er","reco.root");

        TCanvas *c1 = new TCanvas("c1", "Telescope Coordinate reconstruction", 1400, 1000);

        c1->Divide(3,2);

        c1->cd(1);

        gPad->SetLogy();

        TH1F *T1X = new TH1F("T1X","T1 Xreco-Xsim",100,-1,1);

        T1X->GetXaxis()->SetTitle("Xreco-Xsim[cm]");

        //T1deE->GetYaxis()->SetTitle("de");

        t->Draw("ERQTelescopeTrack_T1_DoubleSi_SD2_XY_0.fTelescopeGlobalX-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fX>>T1X","","gof");

        T1X->Draw();

        c1->cd(2);

        gPad->SetLogy();

        TH1F *T1Y = new TH1F("T1Y","T1 Yreco-Ysim",100,-1,1);

        T1Y->GetXaxis()->SetTitle("Yreco-Ysim[cm]");

        //T1deE->GetYaxis()->SetTitle("de");

        t->Draw("ERQTelescopeTrack_T1_DoubleSi_SD2_XY_0.fTelescopeGlobalY-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fY>>T1Y","","gof");

        T1Y->Draw();

        c1->cd(3);

        gPad->SetLogy();

        TH1F *T1Z = new TH1F("T1Z","T1 Zreco-Zsim",100,-1,1);

        T1Z->GetXaxis()->SetTitle("Zreco-Zsim[cm]");

        //T1deE->GetYaxis()->SetTitle("de");

        t->Draw("ERQTelescopeTrack_T1_DoubleSi_SD2_XY_0.fTelescopeGlobalZ-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fZ>>T1Z","","gof");

        T1Z->Draw();

        c1->cd(4);

        gPad->SetLogy();

        TH1F *T2X = new TH1F("T2X","T2 Xreco-Xsim",100,-1,1);

        T2X->GetXaxis()->SetTitle("Xreco-Xsim[cm]");

        //T1deE->GetYaxis()->SetTitle("de");

        t->Draw("ERQTelescopeTrack_T2_DoubleSi_SD2_XY_1.fTelescopeGlobalX-ERQTelescopeSiPoint_T2_DoubleSi_SD2_XY_1_YX.fX>>T2X","","gof");

        T2X->Draw();

        c1->cd(5);

        gPad->SetLogy();

        TH1F *T2Y = new TH1F("T2Y","T2 Yreco-Ysim",100,-1,1);

        T2Y->GetXaxis()->SetTitle("Yreco-Ysim[cm]");

        //T1deE->GetYaxis()->SetTitle("de");

        t->Draw("ERQTelescopeTrack_T2_DoubleSi_SD2_XY_1.fTelescopeGlobalY-ERQTelescopeSiPoint_T2_DoubleSi_SD2_XY_1_Y.fY>>T2Y","","gof");

        T2Y->Draw();

        c1->cd(6);

        gPad->SetLogy();

        TH1F *T2Z = new TH1F("T2Z","T2 Zreco-Zsim",100,-1,1);

        T2Z->GetXaxis()->SetTitle("Zreco-Zsim[cm]");

        //T1deE->GetYaxis()->SetTitle("de");

        t->Draw("ERQTelescopeTrack_T2_DoubleSi_SD2_XY_1.fTelescopeGlobalZ-ERQTelescopeSiPoint_T2_DoubleSi_SD2_XY_1_Y.fZ>>T2Z","","gof");

        T2Z->Draw();

}