void p2n_makeplots() {

TCut ptot_gt_25 = "ptot>2.5";
TCut zcut_gt_1.0 = "(z12*sqrt(x9^2+y9^2)-z9*sqrt(x12^2+y12^2))/(sqrt(x9^2+y9^2)-sqrt(x12^2+y12^2)) > 1.0";

TCanvas *c1 = new TCanvas("c1","c->D->mu",800,800); 
c1->Divide(2,2); 
 
  c1_1->cd();
  TH1F *a1 = new TH1F("a1","z vertex for p>2.5", 50, -2.5, 2.5);
  p2n->Draw("(z12*sqrt(x9^2+y9^2)-z9*sqrt(x12^2+y12^2))/(sqrt(x9^2+y9^2)-sqrt(x12^2+y12^2)) >> a1", ptot_gt_25);
  c1_1->SetLogy();
  a1->Draw();

  c1_2->cd();
  TH1F *a2 = new TH1F("a2"," ptot", 50, 0.0, 10.0);
  p2n->Draw("ptot>>a2");
  c1_2->SetLogy();
  a2->Draw();

    c1_3->cd();
  TH1F *a3 = new TH1F("a3"," ptot, z>1cm", 50, 0.0, 10.0);
  p2n->Draw("ptot>>a3","ptot>2.5 && (z12*sqrt(x9^2+y9^2)-z9*sqrt(x12^2+y12^2))/(sqrt(x9^2+y9^2)-sqrt(x12^2+y12^2)) >1.0");
  c1_3->SetLogy();
  a3->Draw();

  }