EXPERT ROOT

void show_sim()

{

        TFile *fr = new TFile("sim.root");

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

        TCanvas *c1 = new TCanvas();

        tr->Draw("RTelescope1Si1Point.fELoss+RTelescope1Si2Point.fELoss", "RTelescope1Si1Point.fPID<3000 && RTelescope1Si2Point.fPID<3000");

        c1->Divide(3,2);

        c1->cd(1);

        tr->Draw("RTelescope1Si1Point.fY:RTelescope1Si1Point.fX", "", "box");

//        tr->Draw("RTelescope1Si1Point.fX", "", "box");

        c1->cd(3);

        tr->Draw("RTelescope1Si1Point.fELoss+RTelescope1Si2Point.fELoss+RTelescope1CsIPoint.fELoss",

                        "RTelescope1Si1Point.fPID<3000 && RTelescope1Si2Point.fPID<3000 && RTelescope1CsIPoint.fPID<3000");

        c1->cd(2);

        tr->Draw("RTelescope1Si1Point.fPID", "RTelescope1Si1Point.fPID>800000000");

        c1->cd(3);

        tr->Draw("RTelescope1Si1Point.fPID", "RTelescope1Si1Point.fPID<800000000");

        c1->cd(4);

        tr->Draw("Sum$(RTelescope1Si1Point.fELoss)+Sum$(RTelescope1Si2Point.fELoss)+Sum$(RTelescope1CsIPoint.fELoss)");

        c1->cd(5);

                tr->Draw("RTelescope1Si1Point.fELoss+RTelescope1Si2Point.fELoss",

                                "RTelescope1Si1Point.fPID<3000 && RTelescope1Si2Point.fPID<3000");

        c1->cd(6);

        tr->Draw("RTelescope1Si1Point.fELoss+RTelescope1Si2Point.fELoss+RTelescope1CsIPoint.fELoss",

                        "RTelescope1Si1Point.fPID>3000 && RTelescope1Si2Point.fPID>3000 && RTelescope1CsIPoint.fPID>3000");

}