histCorr.C

void histCorr()
{
	gSystem->Load("../libData.so");
	TFile *f = new TFile("../data/dataDSR4/analysis_07_8.root");
	TTree *t = (TTree*)f->Get("atree");
	TCanvas *c1 = new TCanvas("canvas", "Correlation pictures");
	c1->Divide(2,2);

	TH2F *h1 = new TH2F("h1", "Time difference between two channels vs an integral of the first channel",500,10,1e-10,500,-300,300);
	TH2F *h2 = new TH2F("h2", "Time difference between two channels vs an integral of the second channel",100,10,1e-10,100,-300,300);
	TH2F *h3 = new TH2F("h3", "Integral of the first channel vs integral of the second channel",100,10,1e-10,100,-300,300);
	TH2F *h4 = new TH2F("h4", "Time difference between two channels vs amplitude of the first channel",100,10,1e-10,100,-300,300);

	c1->cd(1);
	t->Draw("( Ach0.fTimeCFD - Ach2.fTimeCFD ) : Ach0.fChargeCFD >> h1","","col");

	c1->cd(2);
	t->Draw("( Ach0.fTimeCFD - Ach2.fTimeCFD ) : Ach2.fChargeCFD >> h2","","col");

	c1->cd(3);
	t->Draw("Ach0.fChargeCFD:Ach2.fChargeCFD >> h3","","col");

	c1->cd(4);

	//TCanvas *c2 = new TCanvas("Correlations 2");
	//c2->Divide(1,2);

}
/*Time difference between two channels vs an integral of the first channel,
2. Time difference between two channels vs an integral of the second channel,
3. Integral of the first channel vs integral of the second channel,
4. Time difference between two channels vs amplitude of the first channel,
5. Time difference between two channels vs amplitude of the second channel,
6. Amplitude of the first channel vs amplitude of the second channel
*/