Subversion Repositories f9daq

#include "TROOT.h"

#include "TFile.h"

#include "TBenchmark.h"

#include "TH1F.h"

#include "TH2F.h"

#include "TCanvas.h"

#include "TStyle.h"

#include "TPad.h"

#include "TPaveText.h"

//#include "TLabel.h"

#include "TF1.h"

#include "TGraph.h"

#include "TSpectrum.h"

#include "stdio.h"

#include "include/RTUtil.h"

int focus(char filename[256] = "test", int chX=3, int chY=2, double parameter1=16.90, double parameter2=17.20, bool debug = false)

{

  const int c_nChannels = 64;

  int map[8][8]={{32,34,53,55,40,42,61,63},

                 {48,50,37,39,56,58,45,47},

                 {33,35,52,54,41,43,60,62},

                 {49,51,36,38,57,59,44,46},

                 {17,19,4,6,25,27,12,14},

                 {1,3,20,22,9,11,28,30},

                 {16,18,5,7,24,26,13,15},

                 {0,2,21,23,8,10,29,31}

                };

  char fnameroot[256];

  TFile* rootfile;

	sprintf(fnameroot, "root/%s.root", filename);

	rootfile = (TFile *) gROOT->FindObject(filename);

	if(rootfile==NULL) rootfile = new TFile(fnameroot);

	if(rootfile==NULL) {

	  printf("Cannot open root file %s!!!\n",fnameroot);

	  return(0);

	}

	// set draw style

	RTSetStyle(gStyle);

  /** Draws the beam profile and fits it with error function

    * on some background function

    */

  TCanvas* canvas9 = new TCanvas("canvas9","canvas9", 500,500);

  canvas9->cd();

  char hname[128];

  int channel = map[chX][chY];

  sprintf(hname, "hnhitsy%d", channel);

  TH1F* h_laser = (TH1F*)rootfile->Get(hname);

  h_laser->Draw();

  h_laser->SetStats(1);

  //TF1* err = new TF1("err","[0]+[1]*TMath::Erf((x-[2])/[3])",parameter1,parameter2);

  TF1* err = new TF1("err","[0]+[1]*TMath::Erf((x-[2])/[3])",h_laser->GetBinCenter(3), h_laser->GetBinCenter(h_laser->GetMaximumBin()));

  err->SetParameter(0, h_laser->GetMinimum());

  err->SetParameter(1, h_laser->GetMaximum() - h_laser->GetMinimum());

  err->SetParameter(2, h_laser->GetBinCenter(h_laser->GetMaximumBin()));

  err->SetParameter(3, 0.025);

  h_laser->Fit(err,"qr");

  h_laser->Fit(err,"lr");

  double sigma = err->GetParameter(3);

  printf("sigma = %2.0f um,     FWHM = %2.0f um\n", sigma*1000, 2.35*sigma*1000);

  return 0;

  }