#define invm_cxx
#include "invm.h"
#include "TH2.h"
#include "TStyle.h"
#include "TCanvas.h"
#include <iostream>
#include "TCanvas.h"   //<===========================
using namespace std;
#include "TMatrixD.h"

void invm::Loop()
{// invm loop starts 

  Dca_phi_1 = new TProfile("Dca_phi_1","dca vs phi(gtr401):pt>1,-40<z<-30",70,-0.1,7.0,-0.3,0.3);
  Dca_phi_2 = new TProfile("Dca_phi_2","dca vs phi(gtr401):pt>1,-30<z<-20",70,-0.1,7.0,-0.3,0.3);
  Dca_phi_3 = new TProfile("Dca_phi_3","dca vs phi(gtr401):pt>1,-20<z<-10",70,-0.1,7.0,-0.3,0.3);
  
  Double_t PI = 3.1428571428571 ; cout << "  PI  " << PI << endl;

  if (fChain == 0) return;
  Int_t nentries = Int_t(fChain->GetEntries()); Int_t nbytes = 0; Int_t nb = 0;

  for(Int_t ev=0 ; ev < nentries ; ev++ )
    {//first  event loop start with   ev = 0  
      if(LoadTree(ev) < 0) break;  nb = fChain->GetEntry(ev); nbytes += nb;

      for(Int_t l=0; l<Gtr401atDCA_ntrk401DCA; l++)
        {//gtr401 event loop 
	  
	  Double_t dca = Gtr401atDCA_par0[l];
	  Double_t phi = Gtr401atDCA_par2[l];
	  if(Gtr401atDCA_par2[l]<0.0)phi = 2*PI + Gtr401atDCA_par2[l];
	  
          if(Gtr401atDCA_nhit[l]>10 && fabs(1/Gtr401atDCA_par4[l])>1.0)
	    {//nhits > 10 && pt>1

	      if(Gtr401atDCA_par1[l]>-40.0 && Gtr401atDCA_par1[l]<-30.0){Dca_phi_1->Fill(phi,dca);}
	      if(Gtr401atDCA_par1[l]>-30.0 && Gtr401atDCA_par1[l]<-20.0){Dca_phi_2->Fill(phi,dca);}
	      if(Gtr401atDCA_par1[l]>-20.0 && Gtr401atDCA_par1[l]<-10.0){Dca_phi_3->Fill(phi,dca);}
	    }
     	}
    }
  
  myfit(); // do for Dca_phi_1 ??  && calculate  parA and par B
  myfit(); // do for Dca_phi_2 ??  && calculate  parA and par B
  myfit(); // do for Dca_phi_3 ??  && calculate  parA and par B

}