{TFile f("pythia_053.root");
TH1F *h10 = new TH1F("hpz",   "pz", 50, -10., 10.);
TH1F *h11 = new TH1F("hpzabs","abs(pz)", 50, 0., 10.);
TH1F *h12 = new TH1F("hx1",   "x1", 50, 0., 0.5);
TH2F *h13 = new TH2F("hx1pz",  "x1 vs pz",            50, 0., 10., 50, 0., 0.5);
TH2F *h14 = new TH2F("hx1pzprof","x1 vs pz, profile", 50, 0., 10., 50, 0., 0.5);

h1000->Draw("Pz1>>hpz","","goff");
h1000->Draw("abs(Pz1)>>hpzabs","Pz1>1.5","goff");
h1000->Draw("x1>>hx1","","goff");
h1000->Draw("x1:abs(Pz1)>>hx1pz","Pz1>1.5","goff");

h1000.Draw("x1:abs(Pz1)>>hx1pzprof","Pz1>1.5&&Pz1<10.0","profsgoff");
hx1pzprof.SetMarkerStyle(21);
hx1pzprof.SetMarkerSize(0.5);
hx1pzprof.Fit("pol1","","samegoff");

TCanvas *c1 = new TCanvas("c1","run 53, charm->single mu",480,600); 
c1->Divide(2,3);
c1_1->cd();
hpz->Draw();


c1_2->cd();
hpzabs->Draw();

c1_3->cd();
c1_3->SetLogy();
hx1->Draw();

c1_4->cd();
hx1pz->Draw();


c1_5->cd();
hx1pzprof->Draw();



}