Commit b6b64e7c authored by Vratislav Chudoba's avatar Vratislav Chudoba

AEvent sligthly modified. Script "histCorr.C" modified.

parent 20aa21b3
......@@ -60,18 +60,7 @@ void AEvent::ProcessEvent() {
fAmpPos[j] = fAmpPos[j] - fZeroLevel;
}
Double_t maxAmp = 0.;
Double_t maxAmpT = 0.;
maxAmp = fAmpPos[0];
for(Int_t j=0; j < fNPoints; j++) {
if(fAmpPos[j] > maxAmp) {
maxAmp = fAmpPos[j];
maxAmpT = fTime[j];
}
}
fAmpMax = maxAmp;
fTimeAmpMax = maxAmpT;
SetMaxAmplitudes();
SetGraphs();
......@@ -227,3 +216,18 @@ Double_t AEvent::GetOnefTime(Int_t i) {
Double_t AEvent::GetOnefAmpPos(Int_t i) {
return fAmpPos[i];
}
void AEvent::SetMaxAmplitudes() {
Double_t maxAmp = 0.;
Double_t maxAmpT = 0.;
maxAmp = fAmpPos[0];
for(Int_t j=0; j < fNPoints; j++) {
if(fAmpPos[j] > maxAmp) {
maxAmp = fAmpPos[j];
maxAmpT = fTime[j];
}
}
fAmpMax = maxAmp;
fTimeAmpMax = maxAmpT;
}
......@@ -102,6 +102,7 @@ public:
private:
void Init();
void SetMaxAmplitudes();
void SetGraphs();
void SetCFD();
};
......
//#include "TCanvas.h"
void histCorr()
{
gSystem->Load("../libData.so");
......@@ -6,19 +8,38 @@ void histCorr()
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 *h1 = new TH2F("h1", "Time difference between two channels vs an integral of the first channel",500,-150.,150.,500,0.,5e-9);
TH2F *h1 = new TH2F("h1", "Time difference between two channels vs an integral of the first channel",125,-50.,50.,125,0.,5e-9);
TH2F *h2 = new TH2F("h2", "Time difference between two channels vs an integral of the second channel",125,-50.,50.,125,0e-9,5e-9);
// TH2F *h3 = new TH2F("h3", "Integral of the first channel vs integral of the second channel",500,0.,5e-9,500,0.,5e-9);
TH2F *h3 = new TH2F("h3", "Integral of the first channel vs integral of the second channel",500,0.,2e-9,500,0.,2e-9);
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");
t->Draw("Ach0.fChargeCFD : ( Ach0.fTimeCFD - Ach2.fTimeCFD ) >> h1","","col");
//x
h1->GetXaxis()->SetTitle("\\Delta \\tau_{0-2} [ns]");
h1->GetXaxis()->CenterTitle();
//y
h1->GetYaxis()->SetTitle("Q_{0} [C]");
h1->GetYaxis()->CenterTitle();
c1->cd(2);
t->Draw("( Ach0.fTimeCFD - Ach2.fTimeCFD ) : Ach2.fChargeCFD >> h2","","col");
t->Draw("Ach2.fChargeCFD : ( Ach0.fTimeCFD - Ach2.fTimeCFD ) >> h2","","col");
h2->GetXaxis()->SetTitle("\\Delta \\tau_{0-2} [ns]");
h2->GetXaxis()->CenterTitle();
//y
h2->GetYaxis()->SetTitle("Q_{2} [C]");
h2->GetYaxis()->CenterTitle();
c1->cd(3);
t->Draw("Ach0.fChargeCFD:Ach2.fChargeCFD >> h3","","col");
//x
h3->GetXaxis()->SetTitle("Q_{2} [C]");
h3->GetXaxis()->CenterTitle();
//y
h3->GetYaxis()->SetTitle("Q_{0} [C]");
h3->GetYaxis()->CenterTitle();
c1->cd(4);
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment