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);


}