#define normist_cxx
#include "normist.h"
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <string>
#include <math.h>
#include <iostream>


void normist::Loop()
{
  if (fChain == 0) return;
  TFile *outfile= new TFile("Dphi.root", "RECREATE");
  TH1F * hdphi1  = new TH1F("hdphi1","delta phi 1",100,0,6.28);
  float deltaphi1;
  float deltaphi2;
  Long64_t nentries = fChain->GetEntriesFast();
  Long64_t nbytes = 0, nb = 0;
  
  for (Long64_t jentry=0; jentry<nentries;jentry++) 
    {
      Long64_t ientry = LoadTree(jentry);
      if (ientry < 0) break;
      nb = fChain->GetEntry(jentry);   nbytes += nb;
      
      deltaphi1 = fabs(mu1phi-mu2phi);
      if(deltaphi1>3.14159264) deltaphi1=2*3.14159264-deltaphi1;
      hdphi1->Fill(deltaphi1);
    }
  
  hdphi1->Write("dphi"); // write the histo in the output file
  outfile->Write();
  outfile->Close();  
  /*
  TCanvas * c1 = new TCanvas("c1", "phi1");
  TCanvas * c2 = new TCanvas("c2", "phi2");
  Double_t scale=1/hdphi1->Integral();
  hdphi1->Scale(scale);
  c1->cd();
  hdphi1->Draw();
  hdphi1->SetLineWidth(2);
  Double_t scale=1/hdphi2->Integral();
  hdphi2->Scale(scale);
  c2->cd();
  hdphi2->Draw();
  hdphi2->SetLineColor(2);
  hdphi2->SetLineWidth(2);
  */
   
}