Subversion Repositories f9daq

Compare Revisions

Ignore whitespace Rev 316 → Rev 315

/drs/tree.C
File deleted
/drs/process_results.cxx
File deleted
\ No newline at end of file
/drs/runroot.bat
File deleted
/drs/tree.h
File deleted
/drs/AnalyseWaveForms.cxx
File deleted
\ No newline at end of file
/drs/DrsChannel.h
8,7 → 8,6
float vcut[2];
int edge;
float offset;
float ymin,ymax;
} DrsChannel;
 
#endif
/drs/drs.C
3,7 → 3,6
#include "TH2.h"
#include "TFile.h"
#include "TGraph.h"
#include "TMultiGraph.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TLegend.h"
42,7 → 41,6
TH1F *h1_qsumy;
TGraph *gr_wfm[NCH];
TLegend *leg;
TMultiGraph *wfm;
 
int drs::histoinit( int nch, DrsChannel* d, int trgch, double ymin, double ymax){
char histname[128];
50,16 → 48,15
c_wfm->Divide(1,1);
c_wfm->cd(1);
const double maxt=200;
const double maxt=1800;
leg = new TLegend(0.8, 0.8, 1, 1);
wfm = new TMultiGraph();
for (int i=0; i<nch; i++) {
sprintf(histname,"chmax_%02d",i);
h1_chmax[i]=new TH1F(histname,Form("%d"),1500,d[i].ymin,d[i].ymax);
h1_chmax[i]=new TH1F(histname,histname,1500,ymin,ymax);
sprintf(histname,"chmaxt_%02d",i);
h1_chmaxt[i]=new TH1F(histname,histname,DWIDTH,-0.1,maxt);
sprintf(histname,"chmin_%02d",i);
h1_chmin[i]=new TH1F(histname,histname,1500,d[i].ymin,d[i].ymax);
h1_chmin[i]=new TH1F(histname,histname,1500,ymin,ymax);
sprintf(histname,"chmint_%02d",i);
h1_chmint[i]=new TH1F(histname,histname,DWIDTH,-0.1,maxt);
sprintf(histname,"letime_%02d",i);
69,45 → 66,25
sprintf(histname,"cftdif_%02d",i);
h1_cftdif[i]=new TH1F(histname,histname,2000,d[i].twin[0]-d[trgch].twin[0]-40,d[i].twin[1]-d[trgch].twin[0]+20);
sprintf(histname,"qdc_%02d",i);
float dx = 7*(d[i].ymax - d[i].ymin);
double ymin, ymax;
if(d[i].edge){
ymin=d[i].ymin*8;
ymax=d[i].ymax*8;
}
else{
ymin=-d[i].ymax*8;
ymax=-d[i].ymin*8;
}
h1_qdc[i]=new TH1F(histname,histname,800, ymin , ymax);
h1_qdc[i]=new TH1F(histname,histname,500,-4.5,4.5);
sprintf(histname,"dpo_%02d",i);
h2_dpo[i]=new TH2F(histname,histname,DWIDTH,-0.1,maxt,1500,d[i].ymin,d[i].ymax);
h2_dpo[i]=new TH2F(histname,histname,DWIDTH,-0.1,maxt,1500,ymin,ymax);
sprintf(histname,"tdcadc_%02d",i);
h2_tdcadc[i]=new TH2F(histname,histname,250,d[i].ymin,d[i].ymax,400,d[i].twin[0]-d[trgch].twin[0],d[i].twin[1]-d[trgch].twin[0]);
h2_tdcadc[i]=new TH2F(histname,histname,250,ymin,ymax,400,d[i].twin[0]-d[trgch].twin[0],d[i].twin[1]-d[trgch].twin[0]);
sprintf(histname,"ctdcadc_%02d",i);
h2_ctdcadc[i]=new TH2F(histname,histname,250,d[i].ymin,d[i].ymax,400,d[i].twin[0]-d[trgch].twin[0],d[i].twin[1]-d[trgch].twin[0]);
h2_ctdcadc[i]=new TH2F(histname,histname,250,ymin,ymax,400,d[i].twin[0]-d[trgch].twin[0],d[i].twin[1]-d[trgch].twin[0]);
gr_wfm[i]=new TGraph(DWIDTH);
//gr_wfm[i]->SetMinimum(d[i].ymin);
//gr_wfm[i]->SetMaximum(d[i].ymax);
gr_wfm[i]->SetMinimum(ymin);
gr_wfm[i]->SetMaximum(ymax);
gr_wfm[i]->SetLineColor(i+1);
gr_wfm[i]->SetMarkerColor(i+1);
sprintf(histname,"ch %d",i);
leg->AddEntry(gr_wfm[i], histname, "l");
wfm->Add(gr_wfm[i]);
//if (i) gr_wfm[i]->Draw("LP"); else gr_wfm[i]->Draw("ALP");
//gr_wfm[i]->GetHistogram()->SetTitle("WaveForms");
//gr_wfm[i]->GetHistogram()->SetDirectory(0);
if (i) gr_wfm[i]->Draw("LP"); else gr_wfm[i]->Draw("ALP");
gr_wfm[i]->GetHistogram()->SetTitle("WaveForms");
gr_wfm[i]->GetHistogram()->SetDirectory(0);
}
wfm->SetMinimum(ymin);
wfm->SetMaximum(ymax);
wfm->Draw("ALP");
//wfm->SetTitle("WaveForms");
//wfm->SetDirectory(0);
 
leg->Draw();
h2_ctalk=new TH2F("ctalk","ctalk",250,0.,0.5,250,0.,0.5);
h2_pos=new TH2F("cpos","cpos",40,-1,1,40,-1,1);
164,9 → 141,8
histoinit(nch, d, trgch, ymin,ymax);
char rootname[0xff];
sprintf(rootname, "%s_x%d_y%d.root", FileName, 0, 0);
froot = new TFile(rootname,"RECREATE");
 
fin=fopen(FileName,"rb");
if (fin==NULL) {
260,6 → 236,7
break;
}
if (iev<first) continue;
 
if (recid[0] == 'E') {
for ( i=0; i<nch; i++) {
if (!lwfm[i]) continue;
287,7 → 264,6
double ph1 = (d[1].edge)? chmax[1][0]: -chmin[1][0];
h2_ctalk->Fill(ph0,ph1);
}
if (event<1) continue;
if (recid[0]!='C') continue;
if (ich>=nch) continue;
388,7 → 364,7
h1_letime[ich]->Fill(letime[ich]);
h1_cftime[ich]->Fill(cftime[ich]);
time(&t0);
 
if ((t0!=told) || ((updfrq>0)&&(!(iev%updfrq) ) ) ) {
c_wfm->Update();
for (i=0; i<DWIDTH; i++) {
403,27 → 379,25
printf("number of events: %d\n",iev);
fclose(fin);
TString pdfname = Form("%s.pdf",FileName);
gStyle->SetOptStat(0);
gStyle->SetOptStat(1111111);
TCanvas* c_dpo= new TCanvas("c_dpo","DPO",750,50,700,700);
c_dpo->Divide(2,2);
 
for (i=0; i<nch; i++) {c_dpo->cd(i+1)->SetLogz(1); h2_dpo[i]->DrawCopy("colz");}
 
gStyle->SetOptStat(1111111);
TCanvas* c_cftdif= new TCanvas("c_cftdif","Constant fraction time difference",750,50,700,700);
 
TCanvas* c_cftdif= new TCanvas("c_cftdif","DPO",750,50,700,700);
c_cftdif->Divide(2,2);
for (i=0; i<nch; i++) {c_cftdif->cd(i+1); h1_cftdif[i]->DrawCopy();}
/*
TCanvas* c_cftcorr= new TCanvas("c_cftcorr","pulseheight vs timing",750,50,700,700);
c_cftcorr->Divide(2,2);
for (i=0; i<nch; i++) {c_cftcorr->cd(i+1); h2_tdcadc[i]->DrawCopy("colz");}
*/
TCanvas* c_adc= new TCanvas("c_adc","adc",750,50,700,700);
c_adc->Divide(2,2);
for (i=0; i<nch; i++) {c_adc->cd(i+1); h1_qdc[i]->DrawCopy();}
/*
TCanvas* c_qsum= new TCanvas("c_qsum","qsum",750,50,700,700);
c_qsum->Divide(1,2);
c_qsum->cd(1); h1_qsumx->DrawCopy();
434,18 → 408,15
c_pos->cd();
h2_pos->DrawCopy("colz");
*/
//c_pos->Print(pdfname+"(","pdf");
//c_wfm->Print(pdfname,"pdf");
c_adc->Print(pdfname+"(","pdf");
//c_qsum->Print(pdfname,"pdf");
c_pos->Print(pdfname+"(","pdf");
c_wfm->Print(pdfname,"pdf");
c_adc->Print(pdfname,"pdf");
c_qsum->Print(pdfname,"pdf");
c_dpo->Print(pdfname,"pdf");
c_cftdif->Print(pdfname + ")","pdf");
c_dpo->Write();
c_adc->Write();
c_cftdif->Write();
froot->Write();
froot->Close();
 
501,21 → 472,20
 
c[1].edge = edge[1];
c[1].offset = offset[1];
drs * d = new drs(FileName, 2, c, trgch, 0, nev, updfrq, -0.55,0.55 );
drs * d = new drs(FileName, 2, c, trgch, nev, updfrq, -0.55,0.55 );
}
 
 
 
void sensl(int nev=1000, int updfrq=20, char *FileName="./data/x1.dat", int trgch=3){
void sensl(int nev=1000, int updfrq=20, char *FileName="./data/x1.dat", int trgch=0){
// nev = -1 all events
double threshold[4]={-0.005,-0.05,-0.005,-0.05}; // threshold for determination of timing
float twin[4][2]={{20,40.},{980.,1150.},{100.,160.},{100.,160.}}; // window for calculation of the timing
float adcgate[4][2]={{20.,30.}, {980.,1150.},{112.,123.},{150.,163.}}; //integration window
double threshold[4]={-0.005,-0.005,-0.005,-0.005}; // threshold for determination of timing
float twin[4][2]={{980.,1150.},{980.,1150.},{980.,1150.},{980,1150.}}; // window for calculation of the timing
float adcgate[4][2]={{980.,1150.},{980.,1150.},{980.,1150.},{980.,1150.}}; //integration window
double vcut[4][2]={{-1,1},{-1,1},{-1,1},{-1,1}}; // cut on the pulseheight for filling the timing histograms
int edge[4]={1,0,1,0};// 1 -- rising edge, 0 -- falling edge
double offset[4]={0.25,0.25, 0.025,0.25};// qdc offset
double ymin[4]={-0.1,-0.5,-0.01,-0.6};// ymin offset
double ymax[4]={0.5,0.5,0.07,0.1};// ymin offset
int edge[4]={0,0,0,0};// 1 -- rising edge, 0 -- falling edge
float offset[4]={0.25,0.25,0.25,0.25};// qdc offset
DrsChannel c[NCH];
for (int k=0;k<NCH;k++){
c[k].cfrac = CFRAC;
526,42 → 496,7
c[k].vcut[j] = vcut[k][j];
c[k].edge = edge[k];
c[k].offset = offset[k];
c[k].ymin = ymin[k];
c[k].ymax = ymax[k];
}
}
drs * d = new drs(FileName, NCH, c, trgch, 0, nev, updfrq , -1,1);
drs * d = new drs(FileName, NCH, c, trgch, nev, updfrq , -0.1,0.12);
}
 
 
 
void ana(TString FileName="./data/x1.dat"){
int nev=-1;
int updfrq=1000;
int trgch=0;
// nev = -1 all events
double threshold[4]={-0.005,-0.05,-0.005,-0.05}; // threshold for determination of timing
float twin[4][2]={{20,40.},{980.,1150.},{100.,160.},{100.,160.}}; // window for calculation of the timing
float adcgate[4][2]={{20.,30.}, {980.,1150.},{112.,123.},{150.,163.}}; //integration window
double vcut[4][2]={{-1,1},{-1,1},{-1,1},{-1,1}}; // cut on the pulseheight for filling the timing histograms
int edge[4]={1,0,1,0};// 1 -- rising edge, 0 -- falling edge
double offset[4]={0.25,0.25, 0.025,0.25};// qdc offset
double ymin[4]={-0.1,-0.5,-0.01,-0.6};// ymin offset
double ymax[4]={0.5,0.5,0.07,0.1};// ymin offset
DrsChannel c[NCH];
for (int k=0;k<NCH;k++){
c[k].cfrac = CFRAC;
c[k].threshold = threshold[k];
for (int j=0;j<2;j++){
c[k].twin[j] = twin[k][j];
c[k].adcgate[j] = adcgate[k][j];
c[k].vcut[j] = vcut[k][j];
c[k].edge = edge[k];
c[k].offset = offset[k];
c[k].ymin = ymin[k];
c[k].ymax = ymax[k];
}
}
cout << "DRS4 " << FileName << endl;
drs * d = new drs(FileName.Data(), NCH, c, trgch, 0, nev, updfrq , -0.5,0.5);
}
/drs/README.txt
1,30 → 1,3
Analiza podatkov SiPM, zajetih z DRS4 samplerjem:
root
.L drs.C+
.L process_results.cxx+
// sprocesiraj vse .dat file v direktoriju datasipm (in poddirektorijih), ki se nimajo ob sebi root fileov
// pozene drs analizo s prednastavljenimi parametri na vsakem fileu in sporducira .root in .pdf
process_results("datasipm",".dat", 0);
 
// Zberi vse histograme v skupnem root fileu in jim dodaj threshold scan grafe sipm.root:
process_results("datasipm",".dat", 1);
 
 
// analiza iz root tree v sipm.root
// output v results (extension je drug parameter ),
// ena stran za vsak sipm type in sample z matriko po temperaturah in napetostih
.L tree.C+
tree t;
tree.Loop(0,"pdf"); //qdc2
tree.Loop(1,"png"); //qdc3
tree.Loop(2,"png"); //dpo2
tree.Loop(3,"png"); //dpo3
tree.Loop(4,"png"); //cftdif2
tree.Loop(5,"png"); //cftdif3
tree.Loop(6,"png"); //thscan
 
__________________________________________________________________________
 
Analysis of the data acquired with the readout chip
 
https://www.psi.ch/drs/drs-chip