EXP1803

#include "TLine.h"

#include "TPad.h"

#include "TCutG.h"

void showT1Sim()

{

        TFile *fc = new TFile("cut.root", "READ");

        TCutG *cHe = (TCutG*)fc->Get("cut3He");

        TFile *fcPos = new TFile("cutsPos.root", "READ");

        TCutG *c3He = (TCutG*)fcPos->Get("c3HePos");

        TCutG *c6He = (TCutG*)fcPos->Get("c6HePos");

        TFile *fr = new TFile("sim_digi.root", "READ");

        TTree *tr = (TTree*)fr->Get("er");

        TCanvas *c3 = new TCanvas("c3", "Telescope info", 1200, 800);

        c3->Divide(3,2);

        c3->cd(1);

//        tr->Draw("(MCTrack.fEnergy - MCTrack.fMass)*1000", "MCTrack.fPdgCode == 1000020060");

        tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000",

                        "ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000<6 && ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000<70","");

        cHe->Draw("same");

        c3->cd(2);

//        tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000",

//                        "ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000<6 && ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000<70","");

//        cHe->Draw("same");

        tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "cut3He", "col");

        c3->cd(3);

        tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "", "col");

        c3He->SetLineWidth(3);

        c3He->SetLineColor(kRed);

        c3He->Draw("same");

        c6He->Draw("same");

        c6He->SetLineWidth(3);

        c6He->SetLineColor(kGreen);

        c3->cd(4);

        tr->Draw("(MCTrack.fEnergy - MCTrack.fMass)*1000", "MCTrack.fPdgCode == 1000010010");

        tr->Draw("(MCTrack.fEnergy - MCTrack.fMass)*1000", "MCTrack.fPdgCode == 2212");

        tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000",

                        "ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000<6 && ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000<70","");

        cHe->Draw("same");

        tr->SetMarkerColor(kRed);

        tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000",

                        "ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000<6 "

                        "&& ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000<70 "

                        "&& cut3He","same");

        c3->cd(5);

        tr->Draw("ERQTelescopeSiDigi_DoubleSi_SD1_XY_0_X.fEdep*1000:ERQTelescopeSiDigi_DoubleSi_SD2_XY_1_X.fEdep*1000","","col");

        c3->cd(6);

        tr->Draw("ERQTelescopeSiPoint_DoubleSi_SD2_XY_1_X.fY:ERQTelescopeSiPoint_DoubleSi_SD2_XY_1_X.fX", "", "col");

        return;

        TCanvas *c1 = new TCanvas("c1", "Telescope info", 800, 800);

//        c1->Divide(3,2);

        std::cout << tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "", "col") << std::endl;

        c3He->Draw("same");

//        std::cout << tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "c3HePos", "col") << std::endl;

//        std::cout << tr->Draw("ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeY:ERQTelescopeTrack_DoubleSi_SD2_XY_1.fTelescopeX", "c6HePos", "col") << std::endl;

}