#define Correlator_cxx
#include "Correlator.h"
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>

void Correlator::Loop()
{

 
 // Main method.

 // a) Measure total execution time; 
 // b) Book all objects;
 // c) Calculate Q-vectors for available angles and weights;
 // d) Calculate n-particle correlations from Q-vectors (using standalone formulas);
 // e) Calculate n-particle correlations from Q-vectors (using recursion);
 // f) Calculate n-particle correlations with nested loops (brute force);
 // g) Final printout with comparison of results of d), e) and f). 

 // a) Measure total execution time: 
 TStopwatch timerStandalone; 
 
 // b) Book all objects:
 Cosmetics();

 // c) Calculate Q-vectors for available angles and weights;
 CalculateQvectors();

 // d) Calculate n-particle correlations from Q-vectors (using standalone formulas):
 //  2-p correlations:
 cout<<endl;
 timerStandalone.Start();
 cout<<" => Calculating 2-p correlations (using standalone formulas)...\r"<<flush;
 TComplex two = Two(h1,h2)/Two(0,0).Re();
 Double_t wTwo = Two(0,0).Re(); // Weight is 'number of combinations' by default
 
 
 correlations[0][0]->Fill(0.5,two.Re(),wTwo); // <<cos(h1*phi1+h2*phi2)>>
 correlations[1][0]->Fill(0.5,two.Im(),wTwo); // <<sin(h1*phi1+h2*phi2)>>
 
  printf("\n\n"); timerStandalone.Print(); printf("\n");

  
 
// e) Final printout: 
 printf("\n cosine:");
 printf("\n  2-p => Standalone:   %.12f",correlations[0][0]->GetBinContent(1));
 printf("\n\n sinus:");
 printf("\n  2-p => Standalone:   %.12f",correlations[1][0]->GetBinContent(1));
 printf("\n\n"); 
      
}