FitNdNonlinQuad.h

Go to the documentation of this file.

//  FitNdNonlinQuad.h
//
//  Nathaniel Tagg, Feb 2006
//  tagg@minos.phy.tufts.edu
//
//  Code to read the PulserGain and PulserGainPin tables,
//  perform a quadratic fit to the function, then convert it to
//  an ADC-centric model, create an ntuple for analysis, make
//  plots of each channel for validation, and store the 
//  output in CalPulserFits tables.
//
//  Creates:
//   A root file called ndlin_fits.root, with histograms and trees
//   Fills the ndfits/ directory with plots (if turned on)
//   Fills the CalPulserFits tables with task 44.
//
//  Code is crude and badly documented, but it works.

#include "Validity/VldContext.h"
#include "Plex/PlexStripEndId.h"
#include "TClonesArray.h"

// The (prototype) of the call to actually do stuff.
// To generate DB, "writeToDb" must be true.
// Be default, low gain PINs are used
// Drawing channel plots takes a lot of time
void FitNdNonlinQuad(
                     VldContext inputContext,
                     const char* outputDirectory = "ndfits/",
                     bool writeToDb = false,
                     bool useHighGainPins = false,
                     bool drawChannelPlots = false
                     );


class NDQStripPulsePoint : public TObject
{
 public:
  NDQStripPulsePoint() {};
  
  Int_t   i;
  Float_t x;  // pin
  Float_t dx; // pin err
  Float_t z;  // pmt
  Float_t dz; // pmterrr
  Float_t zerocorr; // zero correction factor on the pmt
  Float_t f;
  Float_t y;
  Float_t dy;
  Float_t caly;
  Float_t dcaly;
  Float_t res;
  Float_t fres;
  Float_t xres;
  Float_t fxres;
  Float_t gain;                                 
  Float_t pe;

  Float_t ph; // pulse height setting, if available.
  Float_t philph; // Phil's measurement, if available.
  ClassDef(NDQStripPulsePoint,1);
};

class NDQStripfit : public TObject
{
 public:
 NDQStripfit(int N=40) : data(0) {
    Reset(N);
  };
  
  void Reset(int N=40) {
    if(!data) data = new TClonesArray("NDQStripPulsePoint",N);
    data->Clear();
    data->ExpandCreate(N);
    fndata = N;
  }
  
  ~NDQStripfit() { delete data; };
  
  // Id.
  PlexStripEndId seid;
  Int_t          plane;
  Int_t          strip;
  Int_t          end;
  Int_t          pb;
  Int_t          led;
  Int_t          pin;
  Int_t          pmt;
  
  // Has data.
  Int_t          ok;

  // Fit.
  Float_t  a;
  Float_t  da;
  Float_t  b;
  Float_t  db;
  Float_t  c;
  Float_t  dc;
  Float_t  d;
  Float_t  dd;
  Float_t  e;
  Float_t  de;
  Float_t  beta;
  Float_t  dbeta;
  
  // Data.
  Int_t n; // true number of points
  Int_t dof; // number used in fit.
  Int_t fndata; // size of data array (may be bigger)
  TClonesArray* data;
  
  // Evaluation.
  Float_t  chisq;

  Float_t  maxres;
  Float_t  maxresf;
  Float_t  maxres2500;
  Float_t  maxresf2500;
  Float_t  meanres2500;
  Float_t  adcmax;

  Int_t    bad;
  ClassDef(NDQStripfit,1);
};


                     

---