Commit 580c6137 authored by Muzalevsky I.A's avatar Muzalevsky I.A

script for painting summ of amplitudes of every channel added

parent d5758504
void analyseRawData()
{
gSystem->Load("../libData.so");
TFile *f = new TFile("../data/rawDataTektronix/exp2.root");
TTree *tr = (TTree*)f->Get("rtree");
const Int_t noBranches = 4;
TString bName;
RawEvent *revent[noBranches]; // pointer to the array (of RawEvent class) in which raw data for each channel will be put
for (Int_t j = 0; j<noBranches; j++) {
revent[j] = new RawEvent(); //each raw event element is of class RawEvent()
bName.Form("ch%d.", j);
tr->SetBranchAddress(bName.Data(), &revent[j]); //read the tree tr with raw data and fill array revent with raw data
}
TH1F *hist1 = new TH1F("hist1", "h1 title", 1000, 0, 100);
TH1F *hist2 = new TH1F("hist2", "h2 title", 1000, 0, 100);
TH1F *hist3 = new TH1F("hist3", "h3 title", 1000, 0, 100);
TH1F *hist4 = new TH1F("hist4", "h4 title", 1000, 0, 100);
TF1 *fit1 = new TF1("fit1","-[0]*exp(-x*[1])");
fit1->SetRange(34,50);
fit1->SetParName(1,"tD");
Long64_t nEntries = tr->GetEntries();
//loop over events
for (j = 0; j < nEntries; j++) {
tr->GetEntry(j);
for(Int_t i =0;i<1000;i++){
//hist->SetBinContent(i,revent[0]->GetAmp(i));
hist1->AddBinContent(i,revent[0]->GetAmp(i));
hist2->AddBinContent(i,revent[1]->GetAmp(i));
hist3->AddBinContent(i,revent[2]->GetAmp(i));
hist4->AddBinContent(i,revent[3]->GetAmp(i));
}
}
TCanvas *c1 = new TCanvas("c1","test",10,10,1000,600);
c1->Divide(2,2);
c1->cd(1);
hist1->Draw();
c1->cd(2);
hist2->Draw();
hist2->Fit("fit1","R","");
c1->cd(3);
hist3->Draw();
c1->cd(4);
hist4->Draw();
}
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