Subversion Repositories f9daq

typedef struct  {

    unsigned short offset;

    unsigned short fineadj_unipol;

    unsigned short kill;

} channel_register;

channel_register creg[144];

double gsigma=0;

double gmean=0;

double sa02calculateoffset(double x, int *offset, int * fineadj_unipol ){

    const double fscale= 2500/1023;      // 2500 mV / 1023 ch / dacch

    const double foffset = 59.0/fscale;  // 59.0 mV / dacch (coarse)

    const double ffine    = 5.4/fscale; // 5.4mV / dacch (fine)

    *offset        =  int(x/ foffset);

    x=x-(*offset)*foffset;

    *fineadj_unipol= int(x/ffine);

        x=x-(*fineadj_unipol)*ffine;

        return x;

}

int sa02fithisto(double mean0=512 , double sigma0=3){

gmean=mean0;

gsigma=sigma0;

 TH2F *h=NULL;

 if (!h) h=((TH2F * ) gDirectory->Get("h2d"));

double max= h->GetMaximum();

TF1 *f = new TF1("merfc","[2]*TMath::Erfc((x-[0])*[1])", 0 , 1000);

f->FixParameter(2,max*0.5);

 TNtuple *nt=new TNtuple("nt","nt","ch:mean:sigma");

 h->FitSlicesY(f);

 TAxis *ax= h->GetYaxis();

 double xmin=  ax->GetXmin();

 double xmax=  ax->GetXmax();

 char n0[100];

 char n1[100];

 char n2[100];

 sprintf(n0,"%s_0",h->GetName());

 sprintf(n1,"%s_1",h->GetName());

 sprintf(n2,"%s_2",h->GetName());

 TH1F *h0 = ((TH1F * ) gDirectory->Get(n0));

 TH1F *h1 = ((TH1F * ) gDirectory->Get(n1));

 TH1F *h2 = ((TH1F * ) gDirectory->Get(n2));

 //h0->Draw();

 TCanvas *c= new TCanvas("c","c",1000,1000);

 c->Divide(2,2);

 c->cd(1);

 h0->SetMaximum(xmax);

 h0->SetMinimum(xmin);

 h0->Draw();

 c->cd(3);

 h1->SetMaximum(200);

 h1->SetMinimum(0);

 sprintf(n0,"%sc",h1->GetName());

 TH1F *h1c = h1->Clone(n0);

 float fx[10];

 for (int i=1;i<h1->GetNbinsX()+1;i++){

    double first = h->GetBinContent(h->FindBin(i-1,xmin+1));

    double last  = h->GetBinContent(h->FindBin(i-1,xmax-1));

    double sigma= TMath::Abs( h1->GetBinContent(i) );

    double mean = TMath::Abs( h0->GetBinContent(i) );

    if (sigma>0) sigma=1/sigma; else sigma=0;

    fx[0]=i;

    fx[1]=mean;

    fx[2]=sigma;

        nt->Fill(fx);

        double            x = ((mean+sigma0*sigma)-mean0) ;

    int offset        = 0;

    int fineadj_unipol= 0;

        sa02calculateoffset(x, &offset, &fineadj_unipol );

        int valid = 1;

        if (last>0.5*h2d->GetMaximum())  valid = 0;

        if (first<0.5*h2d->GetMaximum()) valid = 0;

        if (valid) {

          creg[i-1].offset         = offset;

          creg[i-1].fineadj_unipol = fineadj_unipol;

          creg[i-1].kill           = 0;

        } else {

          creg[i-1].offset         = 0;

          creg[i-1].fineadj_unipol = 0;

          creg[i-1].kill           = 1;

        }

    printf("%d %f %f %f %f %d %d %f\n",i, first, last, mean,sigma,offset,fineadj_unipol ,x );

    h1c->SetBinContent(i,sigma);

    // of << ch << "\t" << (int)entries << "\t" << (((mean*2.4)-1245.2)) << "\t" << sigma*2.4 << endl;

 }

 h1c->Draw();

 c->cd(4);

 nt->Draw("sigma","sigma>0&&sigma<100");

 //h2->Draw();

 c->cd(2);

 //h->Draw("colz");

 /*

 h0->SetLineColor(kRed);

 h0->SetLineWidth(10);

 h0->Draw("same");

 */

 nt->Draw("mean","mean>0&&mean<1000");

 c->Modified();

 c->Update();

}

int sa02param( char *fname="../params/offset_orig.param", char *hname="offsettmp.param" ){

    int ndim=400;

    char line[400];

    char cmd[400];

    char sasic[400];

    char v0[400];

    char v1[400];

    int asic=0, ch=0;

    int  cval=0;

    int i;

    FILE *fout = fopen(hname,"w");

    if (!fout) {

        printf("Error! Cannot open file %s\n",hname);

        return -1;

    }

    FILE *fp = fopen(fname,"r");

    if (!fp) {

        printf("Error! Cannot open file %s\n",fname);

        return -1;

    }

        int write=1;

        fprintf(fout,"### Created from sa02fit mean=%f sigma=%f\n", gmean,gsigma);

        fprintf(fout,"### date tbi\n");

        char *null=NULL;

    while (fgets(line,ndim,fp)!=NULL) {

                write=1;

        int nb = sscanf(line,"%s%s%s%s",cmd,sasic,v0,v1);

        if (nb<1 || cmd[0]=='#') {

           fprintf(fout,line);

           continue;

        }  

        asic =   atoi (sasic);

        ch   =   atoi (v0);

        cval =   strtoul (v1,null,0);

        if (strcmp(cmd,"offset")==0)     {

                fprintf(fout,"offset %d %d %d\n", asic,ch, creg[asic*36+ch].offset);

                write=0;

        }

        if (strcmp(cmd,"fineadj_unipol")==0) {

                fprintf(fout,"fineadj_unipol %d %d %d\n", asic,ch, creg[asic*36+ch].fineadj_unipol);

                write=0;

        }

        if (strcmp(cmd,"kill")==0)       {

                fprintf(fout,"kill %d %d %d\n", asic,ch, creg[asic*36+ch].kill);

                write=0;

        }

        if (write){

           fprintf(fout,line);

        }

    }

    fclose(fp);

    fclose(fout);

return 0;

}

int sa02fit(char *hname,double mean0=400 , double sigma0=3, char *fnameparam="../../params/offset_orig.param"){

 gROOT->SetMacroPath("../macros");

 gROOT->ProcessLine(".L H2Dload.cxx");

 H2Dload(hname);

 //H2Dload(gSystem->UnixPathName(hname));

 sa02fithisto(mean0,sigma0);

 char foutparname[400];

 sprintf(foutparname,"%s.param",hname);

 sa02param(fnameparam,foutparname);

return 0;      

}