FMD/scripts/VA1Train.C

   1 //
   2 // $Id$
   3 //
   4 // Small script that shows a signal train from a VA1 pre-amp.
   5 //
   6 void
   7 VA1Train()
   8 {
   9   gStyle->SetOptTitle(0);
  10   gStyle->SetOptFit(0);
  11   gStyle->SetOptStat(0);
  12
  13   TCanvas* c = new TCanvas("c", "C", 800, 400);
  14   c->SetFillColor(0);
  15   c->SetBorderMode(0);
  16   c->SetBorderSize(0);
  17
  18
  19   TArrayF measurements(6);
  20   std::cout << "Measurements are: " << std::flush;
  21   for (Int_t i = 0; i < measurements.fN; i++) {
  22     measurements[i] = gRandom->Uniform(0,1);
  23     std::cout << measurements[i] << " " << std::flush;
  24   }
  25   std::cout << std::endl;
  26
  27   Float_t last = 0;
  28   Float_t B    = 5;
  29   TH2* frame = new TH2F("frame", "Frame", measurements.fN, 0,
  30                         measurements.fN, 10, 0, 1.1);
  31   frame->Draw();
  32   for (Int_t i = 0; i < measurements.fN; i++) {
  33     TF1* f = new TF1("f", "[2] + exp(-[1] * (x - [3])) * ([0] - [2])",
  34                      i, i + 1);
  35     f->SetParameter(3, i);
  36     f->SetParameter(1, B);
  37     f->SetParameter(2, measurements[i]);
  38     f->SetParameter(0, last);
  39
  40     if (measurements[i] > last) {
  41       // f = new TF1("f", "[0] * (1 - exp(-[1] * (x - [2]))) + [3]", i, i+1);
  42       // f->SetParameters(measurements[i] - last, B, i, last);
  43       // f->SetParameter(0, measurements[i]);
  44       // f->SetParameter(2, last);
  45     }
  46     else {
  47       // f->SetParameter(2, measurements[i]);
  48       // f->SetParameter(0, last);
  49       // f = new TF1("f", "[0] * (exp(-[1] * (x - [2]))) + [3]", i, i+1);
  50       // f->SetParameters(last - measurements[i], B, i, measurements[i]);
  51     }
  52     f->Draw("same");
  53     last = measurements[i];
  54   }
  55 }
  56 //____________________________________________________________________
  57 //
  58 // EOF
  59 //