#include "TAxis.h"
#include "TH2.h"
#include "TStyle.h"
#include "TCanvas.h"
#include "TNtupleD.h"
#include <stdio.h>
#include <stdlib.h>
#include <vector>

void script(){



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


   
  TNtuple* ntuple = (TNtuple*)*file->Get("ntuple_muon");
   
 
  
  TH2F* hist = new TH2F("hist","veritices distribution",40,-31.,31.,40,-11.,11.);
  hist->GetXaxis()->SetTitle("z (mm)");
  hist->GetYaxis()->SetTitle("x (mm)");
  ntuple->Draw("x0:z0>>hist","1.>0.","CONTZ");

  // CHECK HIST HAS ENTRIES
  cout << hist->Integral() << endl;

  c=new TCanvas();
  hist->Draw("CONTZ");

  c=new TCanvas();
  hist->Draw();

  c=new TCanvas();
  hist->Draw("SCAT");
  
  


  file->Close(); 
 

}