#include "TProfile.h"

Float_t x[10000], y[10000], ex[10000], ey[10000];
Int_t   nentries;
TGraphErrors* plot;

void plot_zee_mass(const char* filename) {

   c1 = new TCanvas("zee_mass","Zee mass vs run number",200,10,700,500);

   c1->GetFrame()->SetBorderSize(12);

  TFile* f = new TFile(filename);

  TTree* tree = (TTree*) f->Get("physmon");

  TPhysMon*  pm = new TPhysMon();

  tree->SetBranchAddress("PhysMon",&pm);

  nentries = tree->GetEntries();

  int nb;

  int npt = 0;
  for (Int_t jentry=0; jentry<nentries; jentry++) {
    nb = tree->GetEntry(jentry);
    //    printf(" cal.time, njobs: %i  %i \n",tq->fCalendarTime,tq->fNJobs);

    if (pm->fMZeeFitErr < 100) {
      x [npt] = pm->fRunNumber;
      y [npt] = pm->fMZee;
      ey[npt] = pm->fMZeeFitErr;
      ex[npt] = 0;
      npt++;
    }
  }

  plot = new TGraphErrors(npt,x,y,ex,ey);
  
  plot->SetTitle("Zee mass vs run number");
  plot->SetMarkerStyle(20);
  plot->SetMarkerSize(1.);
  plot->Draw("AP");
  plot->GetHistogram()->GetXaxis()->SetTitle("run number");
  plot->Draw("AP");
  c1->Update();
}