Merge branch 'master' of http://er.jinr.ru:90/vratislav.chudoba/NeuRad_tests

macros/integralFormSignal.cpp

+void integralFormSignal()
+{
+	gSystem->Load("../libData.so");
+
+	TFile *f = new TFile("../data/dataTektronix/exp2.root");
+	TTree *tr = (TTree*)f->Get("atree");
+
+	const Int_t noBranches = 4;
+	Double_t ZeroTime[4];
+	Int_t iZero[4],PosZero[4],deltaT[4]; 
+	
+	TString bName;
+	AEvent *aevent[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++) {
+		aevent[j] = new AEvent();	//each raw event element is of class RawEvent()
+		bName.Form("Ach%d.", j);
+		tr->SetBranchAddress(bName.Data(), &aevent[j]);	//read the tree tr with raw data and fill array revent with raw data
+	}
+
+	TH1F *hist0 = new TH1F("hist1", "h1 title", 2000, -100, 100);
+	TH1F *hist1 = new TH1F("hist2", "h2 title", 2000, -100, 100);
+	TH1F *hist2 = new TH1F("hist3", "h3 title", 2000, -100, 100);
+	TH1F *hist3 = new TH1F("hist4", "h4 title", 2000, -100, 100);
+
+	TF1 *fit1 = new TF1("fit1","-[0]*exp(-x*[1])");
+	fit1->SetRange(-15,-5);
+	fit1->SetParName(1,"tD");
+
+	Long64_t nEntries = tr->GetEntries();
+
+	//loop over events
+	for (j = 0; j < nEntries; j++) {
+		tr->GetEntry(j);
+		for(Int_t k=0;k<4;k++) { // loop for channels
+			ZeroTime[k] = aevent[k]->GetfCFD(); // get fCFD for channel №k
+			iZero[k] = ZeroTime[k]*10; 
+			if(j==0) { PosZero[k] = iZero[k]; }
+			deltaT[k] = iZero[k] - PosZero[k];
+		}
+		
+
+		//for(k=0;k<4;k++) {
+			for(Int_t i = 0; i<1000; i++){
+				if( ((i+deltaT[0])>-1) && ((i+deltaT[0])<1000))	hist0->AddBinContent(i+500,aevent[0]->GetOnefAmpPos(i+deltaT[0]));	
+				if( ((i+deltaT[1])>-1) && ((i+deltaT[1])<1000))	hist1->AddBinContent(i+500,aevent[1]->GetOnefAmpPos(i+deltaT[1]));
+				hist2->AddBinContent(i+500,aevent[0]->GetOnefAmpPos(i));
+				hist3->AddBinContent(i+500,aevent[1]->GetOnefAmpPos(i));
+			}
+		//}
+	}
+
+	TCanvas *c1 = new TCanvas("c1","test",10,10,1000,600);
+	c1->Divide(2,2);
+
+	c1->cd(1);
+	hist0->Draw();
+	//hist1->Fit("fit1","R","");
+
+	c1->cd(2);
+	//hist2->Draw();
+	//hist2->Fit("fit1","R","");
+	//tr->Draw("Ach1.fCFD");
+	hist1->Draw();
+
+	c1->cd(3);
+	hist2->Draw();
+
+	c1->cd(4);
+	hist3->Draw();
+
+}