FMD/scripts/VA1Train.C

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