diff --git a/convertTektronix/IntegralCFD.c b/convertTektronix/IntegralCFD.c index b7e60da71bc8b5252efc0416ca7122c5060c77cb..fa85964cdc080b0ec7464f74a50f7562ff4ab78c 100644 --- a/convertTektronix/IntegralCFD.c +++ b/convertTektronix/IntegralCFD.c @@ -6,7 +6,7 @@ void IntegralCFD(){ Int_t j,i,nentry,NumM,imax,imin; Double_t A1[1000],A2[1000],A3[1000],A4[1000]; /// амплитуды с сырых данных (отрицательные) Double_t T[1000]; - Double_t Am[1000]; // амплитуды с перевёрнутых пиков + Double_t Am[1000],Ahist[1000]; // амплитуды с перевёрнутых пиков,.... амплитуда в гистограмме Double_t Tim[1000]; Double_t Acfd[1000]; Double_t sum,zeroLEVEL1,zeroLEVEL2,tm,t1_15,m,m1_15,t2_15,m2_15,tL_15,mL_15,Integrall,LinPar,ChiSquare,t1_90,t2_90,minCFD,maxCFD,zCFD,Tcfd,deltaZero; @@ -66,6 +66,7 @@ void IntegralCFD(){ ///////// заполнение дерева TCanvas *c1 = new TCanvas("c","zero line"); j=0; + TH1F *hA = new TH1F("hA","full amlitude",1000,5e-04,100+5e-04); // гистограмма для заполнения данными из триггерного канала for(j=0;jGetEntries();j++) // j - номер события , i - внутри! { tree1->GetEntry(j); @@ -89,21 +90,22 @@ void IntegralCFD(){ if(Acfd[i]Draw("Al*"); - c1->Update(); - cout<(tm-3)) && (Tim[i]<(tm+17)) ) {Integrall = Integrall + Am[i];} } - //cout<Draw("Al*"); + + //hA = hA + (gr1->GetHistogram()); + //hA->Add(gr1->GetHistogram()); + //hA = gr1->GetHistogram(); + //hA->Draw(); + //c1->Update(); + //gr1->Draw("Al*"); gr1->Fit("fit1","R","goff"); zeroLEVEL1 = fit1->GetParameter(0); c1->Update(); @@ -143,7 +145,10 @@ void IntegralCFD(){ LinPar = fit3->GetParameter(0); theTree->Fill(); - } + } + c1->cd(); + hA->Draw(); + c1->Update(); ////////// theTree->Write(); f1->Close(); diff --git a/convertTektronix/exp2.root b/convertTektronix/exp2.root index 862ca8338b506ba016b1f4d7bba94a7b118cdf54..3940361a479d9b54b32a1cbad59c8175ae77e561 100644 Binary files a/convertTektronix/exp2.root and b/convertTektronix/exp2.root differ diff --git a/convertTektronix/read1.cpp b/convertTektronix/read1.cpp index c56e33be550c7d8edffd3d66227163686b853f7d..e70ab099374ff0a31adb7efdd0d86a9c2af72cd5 100644 --- a/convertTektronix/read1.cpp +++ b/convertTektronix/read1.cpp @@ -1,6 +1,8 @@ -void read1() { +void read1(Int_t input=0) { +cout<Load("../libData.so"); - + Double_t A1[1000],A2[1000],A3[1000],A4[1000]; Double_t T[1000]; Int_t i,j,k,n,p,Nchannel,Nlines; @@ -55,7 +57,7 @@ void read1() { for(i=0;;i++){ fileName.ReadLine(myfile); NameLines[i] = fileName; - if ( fileName.IsNull() ) break; + if ( fileName.IsNull() ) {break;} } myfile.close(); diff --git a/dataClasses/AEvent.cpp b/dataClasses/AEvent.cpp index 362afefd7a34801d8b4e8a4b6b06a4372a05a4bf..ca56e016fcb09de1c868b6b0a8f2ade98814d7b8 100644 --- a/dataClasses/AEvent.cpp +++ b/dataClasses/AEvent.cpp @@ -168,7 +168,16 @@ void AEvent::SetCFD() { } } - return; +} +Double_t AEvent::GetfCFD() { + return fCFD; +} + +Double_t AEvent::GetOnefTime(Int_t i) { + return fTime[i]; +} +Double_t AEvent::GetOnefAmpPos(Int_t i) { + return fAmpPos[i]; } diff --git a/dataClasses/AEvent.h b/dataClasses/AEvent.h index 145155835643c8896a402e8a98cdb1d35658e317..5898cec50f929a4772af85323069ce4b79959b4a 100644 --- a/dataClasses/AEvent.h +++ b/dataClasses/AEvent.h @@ -54,7 +54,9 @@ public: AEvent(const Int_t npoints); virtual ~AEvent(); ClassDef(AEvent,1); - + Double_t GetfCFD(); + Double_t GetOnefTime(Int_t i); + Double_t GetOnefAmpPos(Int_t i); void SetRawDataFile(const char* inprawfile, const char* treename); void ProcessEvent(); // void Integral()