EXP1803

void showBeamDetElossSim()

{

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

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

        TFile *freco = new TFile("reco.root", "READ");

        TTree *treco = (TTree*)freco->Get("er");

        TCanvas *c1 = new TCanvas("c1", "MWPC profiles", 1200, 800);

        c1->Divide(3,2);

        c1->cd(1);

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

        treco->Draw("BeamDetParticle.fTofState.fE*1000 - 5606.556", "");

        c1->cd(2);

        treco->Draw("BeamDetParticle.fTargetState.fE*1000 - 5606.556", "");

        c1->cd(3);

        treco->Draw("BeamDetParticle.fTargetState.fE*1000 - 5606.556", "BeamDetTrack.fXt && BeamDetTrack.fYt");

        c1->cd(4);

        treco->Draw("BeamDetParticle.fTofState.fE*1000 - 5606.556", "BeamDetParticle.fPID == 1000020060");

        c1->cd(5);

        treco->Draw("BeamDetParticle.fTargetState.fE*1000 - 5606.556", "TMath::Abs(BeamDetTrack.fXt)<0.1 && TMath::Abs(BeamDetTrack.fYt)<0.1");

        c1->cd(6);

}