#include <iostream>
using namespace std;

void correlationplot()
{
	TStopwatch t;
	TFile *file = new TFile ("cal0061.root");
	TTree *tree = (TTree*)file->Get("cal");
	TChain *ch1 = new TChain("cal");
	ch1 -> Add ("cal0061.root");
	
//	cout<< ch1->GetEntries()<<endl;
	ch1 -> Add ("cal0062.root");
	ch1 -> Add ("cal0063.root");
	ch1 -> Add ("cal0064.root");
	ch1 -> Add ("cal0065.root");
	ch1 -> Add ("cal0066.root");
	ch1 -> Add ("cal0067.root");
	ch1 -> Add ("cal0068.root");
	ch1 -> Add ("cal0069.root");
	ch1 -> Add ("cal0070.root");
	ch1 -> Add ("cal0071.root");
	ch1 -> Add ("cal0072.root");
	ch1 -> Add ("cal0073.root");
	ch1 -> Add ("cal0074.root");

	cout<< ch1->GetEntries()<<endl;

    TCanvas *c1 = new TCanvas("c1", "Tx2 vs Tx1 correlation plot"); 
    
    for(int ippac=0; ippac<7; ippac++){
        if (ippac == 5) {
            continue;
        } //to exclude the default detector ippac == 5; 
        for (int j=0; j<4; j++){
           
	        c1->Divide(ippac,j); 
	        c1->cd(ippac+j*6+1);
	        gPad-> SetLogz();
	        gPad->SetGrid(1,1);
            gStyle->SetPalette(kRainBow);

            ch1->Draw(TString::Format("PPAC_Tx2[%d][%d] : PPAC_Tx1[%d][%d]", ippac, j, ippac, j), TString::Format("PPAC_Tx2[%d][%d]>-9999 && PPAC_Tx1[%d][%d]>-9999", ippac, j , ippac, j),"colz");
            



        }
    }
	t.Print();

}