Plotter.cxx

Go to the documentation of this file.

// Plotter.cxx

#include <iostream>
#include "PulserCalibration/Plotter.h"
#include "TROOT.h"
#include "TPad.h"
#include "TCanvas.h"
#include "TGraphErrors.h"
#include "TGraph.h"
#include "TAxis.h"
#include "TF1.h"

using std::cout;
using std::endl;

//ClassImp(Plotter)
  //...

Plotter::Plotter(VldContext vld) : fPinRP(vld), fNearRP(vld), fFarRP(vld), fPlex(vld)
{
  fVld = vld;
  fRes = 0;
  cout<<"Got Tables for "<<vld<<endl;
  cout<<"Pin: "<<fPinRP.GetNumRows()<<" rows"<<endl;
  cout<<"Near: "<<fNearRP.GetNumRows()<<" rows"<<endl;
  cout<<"Far: "<<fFarRP.GetNumRows()<<" rows"<<endl;
}

void Plotter::BasicPlot(int plane, int strip)
{
  // Make Ids...
  PlexStripEndId seid(Detector::kCalDet, plane, strip, StripEnd::kEast);
  PlexStripEndId other(Detector::kCalDet, plane, strip, StripEnd::kWest);

  PlexLedId led = myGetLedId(fPlex,seid);
  PlexLedId otherled = myGetLedId(fPlex,other);

  std::pair<PlexPinDiodeId,PlexPinDiodeId> pins = myGetPinDiodeIds(fPlex,led);
  std::pair<PlexPinDiodeId,PlexPinDiodeId> otherpins = myGetPinDiodeIds(fPlex,otherled);

  const PulserRawGainPin* nph = fPinRP.GetRowByIndex(pins.first.GetEncoded());
  const PulserRawGainPin* npl = fPinRP.GetRowByIndex(pins.second.GetEncoded());
  const PulserRawGainPin* fph = fPinRP.GetRowByIndex(otherpins.first.GetEncoded());
  const PulserRawGainPin* fpl = fPinRP.GetRowByIndex(otherpins.second.GetEncoded());
  cout<<"PIN Ids: "<<pins.first.GetEncoded()<<" "<<pins.second.GetEncoded()
      <<" "<<otherpins.first.GetEncoded()<<" "<<otherpins.second.GetEncoded()
      <<endl;
  TVirtualPad *parent = gPad;
  if (!parent) {
    //TCanvas *canvas = 
    new TCanvas("plotter","Plotter",0,0,700,700);
    parent = gPad;
  }
  parent->Clear();
  parent->Divide(2,2,0,0,0);
  

  const PulserRawGain* near = fNearRP.GetRowByIndex(seid.BuildPlnStripEndKey());
  const PulserRawGainF* far = fFarRP.GetRowByIndex(other.BuildPlnStripEndKey());
  if (near && far) {
    parent->cd(1);
    TGraphErrors *g1h = 
      new TGraphErrors(near->GetNumPoints(),
                       nph->GetMean(),near->GetMean(),
                       nph->GetError(),near->GetError());
    g1h->SetMarkerStyle(2);
    g1h->SetMarkerColor(2);
    g1h->Draw("AP");
    TGraphErrors *g1l = 
      new TGraphErrors(near->GetNumPoints(),
                       npl->GetMean(),near->GetMean(),
                       npl->GetError(),near->GetError());
    g1l->SetMarkerStyle(2);
    g1l->SetMarkerColor(4);
    g1l->Draw("P");

    parent->cd(2);
    TGraphErrors *g2h = 
      new TGraphErrors(far->GetNumPoints(),
                       nph->GetMean(),far->GetMean(),
                       nph->GetError(),far->GetError());
    g2h->SetMarkerStyle(2);
    g2h->SetMarkerColor(2);
    g2h->Draw("AP");
    TGraphErrors *g2l = 
      new TGraphErrors(far->GetNumPoints(),
                       npl->GetMean(),far->GetMean(),
                       npl->GetError(),far->GetError());
    g2l->SetMarkerStyle(2);
    g2l->SetMarkerColor(4);
    g2l->Draw("P");
  }    
  near = fNearRP.GetRowByIndex(other.BuildPlnStripEndKey());
  far = fFarRP.GetRowByIndex(seid.BuildPlnStripEndKey());
  if (near && far) {
    parent->cd(4);
    TGraphErrors *g4h = 
      new TGraphErrors(near->GetNumPoints(),
                       fph->GetMean(),near->GetMean(),
                       fph->GetError(),near->GetError());
    g4h->SetMarkerStyle(2);
    g4h->SetMarkerColor(2);
    g4h->Draw("AP");
    TGraphErrors *g4l = 
      new TGraphErrors(near->GetNumPoints(),
                       fpl->GetMean(),near->GetMean(),
                       fpl->GetError(),near->GetError());
    g4l->SetMarkerStyle(2);
    g4l->SetMarkerColor(4);
    g4l->Draw("P");

    parent->cd(3);
    TGraphErrors *g3h = 
      new TGraphErrors(far->GetNumPoints(),
                       fph->GetMean(),far->GetMean(),
                       fph->GetError(),far->GetError());
    g3h->SetMarkerStyle(2);
    g3h->SetMarkerColor(2);
    g3h->Draw("AP");
    TGraphErrors *g3l = 
      new TGraphErrors(far->GetNumPoints(),
                       fpl->GetMean(),far->GetMean(),
                       fpl->GetError(),far->GetError());
    g3l->SetMarkerStyle(2);
    g3l->SetMarkerColor(4);
    g3l->Draw("P");
  }    
                       
  parent->cd(0);
  

}



TGraphErrors * Plotter::Plot(int yplane, int ystrip, StripEnd::StripEnd_t yend, 
                             int xplane, int xstrip, StripEnd::StripEnd_t xend,
                             int led, char *draw_opt)
{
 TPad *parent = dynamic_cast<TPad *>(gPad);
  if (!parent) {
    //TCanvas *canvas = 
    new TCanvas("plotter","Plotter",0,0,700,700);
    parent = dynamic_cast<TPad *>(gPad);
  }


  const float *x=0;
  const float *ex = 0;
  const float *y = 0;
  const float *ey = 0;
  int num;

  // Get X
  PlexStripEndId xid(Detector::kCalDet, xplane, xstrip, xend);
 
  PlexLedId tryled = myGetLedId(fPlex, xid);
  if (tryled.GetLedInBox() == led) {
    // Near end
    const PulserRawGain *xn = fNearRP.GetRowByIndex(xid.BuildPlnStripEndKey());
    if (!xn) cout<<" Oh Bugger!"<<endl;
    x = xn->GetMean();
    ex = xn->GetError();
    num = xn->GetNumPoints();
  }
  else if (myGetLedId(fPlex,xid.GetOppositeSEId()).GetLedInBox() == led) {
    // Far end
    const PulserRawGainF *xn = fFarRP.GetRowByIndex(xid.BuildPlnStripEndKey());
    if (!xn) cout<<" Oh Bugger!"<<endl;
    x = xn->GetMean();
    ex = xn->GetError();
    num = xn->GetNumPoints();
  }
  else {
    cout<<" No X data for this LED"<<endl;
    return 0;
  }

  // Get Y
  PlexStripEndId yid(Detector::kCalDet, yplane, ystrip, yend);
 
  
  if (myGetLedId(fPlex,yid).GetLedInBox() == led) {
    // Near end
    const PulserRawGain *yn = fNearRP.GetRowByIndex(yid.BuildPlnStripEndKey());
    if (!yn) cout<<" Oh Bugger!"<<endl;
    y = yn->GetMean();
    ey = yn->GetError();
  }
  else if (myGetLedId(fPlex,yid.GetOppositeSEId()).GetLedInBox() == led) {
    // Far end
    const PulserRawGainF *yn = fFarRP.GetRowByIndex(yid.BuildPlnStripEndKey());
    if (!yn) cout<<" Oh Bugger!"<<endl;
    y = yn->GetMean();
    ey = yn->GetError();
  }
  else {
    cout<<" No Y data for this LED"<<endl;
    return 0;
  }

  TGraphErrors *g = new TGraphErrors(num,x,y,ex,ey);
  g->SetMarkerStyle(2);
  g->SetMarkerColor(2);
  g->SetTitle("");
  g->Draw(draw_opt);
  char name[100];
  sprintf(name,"Plane %d Strip %d %s",xplane,xstrip,StripEnd::AsString(xend));
  if (g->GetXaxis()) g->GetXaxis()->SetTitle(name);
  sprintf(name,"Plane %d Strip %d %s",yplane,ystrip,StripEnd::AsString(yend));
  if (g->GetXaxis()) g->GetYaxis()->SetTitle(name);
  g->Draw(draw_opt);
  return g;
}

TGraphErrors * Plotter::Plot(int plane, int strip, StripEnd::StripEnd_t end,
                             float *x, float *ex, char *xtitle,
                             int led, char *draw_opt)
{
 TPad *parent = dynamic_cast<TPad *>(gPad);
  if (!parent) {
    //TCanvas *canvas = 
    new TCanvas("plotter","Plotter",0,0,700,700);
    parent = dynamic_cast<TPad *>(gPad);
  }
  const float *y = 0;
  const float *ey = 0;
  int num;

  // Get Y
  PlexStripEndId yid(Detector::kCalDet, plane, strip, end);
 
  
  if (myGetLedId(fPlex,yid).GetLedInBox() == led) {
    // Near end
    const PulserRawGain *yn = fNearRP.GetRowByIndex(yid.BuildPlnStripEndKey());
    if (!yn) cout<<" Oh Bugger!"<<endl;
    y = yn->GetMean();
    ey = yn->GetError();
    num = yn->GetNumPoints();
  }
  else if (myGetLedId(fPlex,yid.GetOppositeSEId()).GetLedInBox() == led) {
    // Far end
    const PulserRawGainF *yn = fFarRP.GetRowByIndex(yid.BuildPlnStripEndKey());
    if (!yn) cout<<" Oh Bugger!"<<endl;
    y = yn->GetMean();
    ey = yn->GetError();
    num = yn->GetNumPoints();
  }
  else {
    cout<<" No Y data for this LED"<<endl;
    return 0;
  }
 TGraphErrors *g = new TGraphErrors(num,x,y,ex,ey);
  g->SetMarkerStyle(2);
  g->SetMarkerColor(2);
  g->SetTitle("");
  g->Draw(draw_opt);
  char name[100];
  if (g->GetXaxis()) g->GetXaxis()->SetTitle(xtitle);
  sprintf(name,"Plane %d Strip %d %s",plane,strip,StripEnd::AsString(end));
  if (g->GetXaxis()) g->GetYaxis()->SetTitle(name);
  g->Draw(draw_opt);
  return g;

}

TGraphErrors * Plotter::Plot(int plane, int strip, StripEnd::StripEnd_t end,
                             int gain,
                             int led, char *draw_opt)
{
  const float *x = 0;
  const float *ex = 0;
  const float *y = 0;
  const float *ey = 0;
  int num;
  // Get Y
  PlexStripEndId yid(Detector::kCalDet, plane, strip, end);
 
  
  if (myGetLedId(fPlex,yid).GetLedInBox() == led) {
    // Near end
    const PulserRawGain *yn = fNearRP.GetRowByIndex(yid.BuildPlnStripEndKey());
    if (!yn) cout<<" Oh Bugger!"<<endl;
    y = yn->GetMean();
    ey = yn->GetError();
    num = yn->GetNumPoints();
  }
  else if (myGetLedId(fPlex,yid.GetOppositeSEId()).GetLedInBox() == led) {
    // Far end
    const PulserRawGainF *yn = fFarRP.GetRowByIndex(yid.BuildPlnStripEndKey());
    if (!yn) cout<<" Oh Bugger!"<<endl;
    y = yn->GetMean();
    ey = yn->GetError();
    num = yn->GetNumPoints();
  }
  else {
    cout<<" No Y data for this LED"<<endl;
    return 0;
  }

  // Get Pin (X)
  PlexLedId ledid(Detector::kCalDet, 0, led);
  cout<<ledid<<endl;
  std::pair<PlexPinDiodeId,PlexPinDiodeId> pinids = myGetPinDiodeIds(fPlex,ledid);
  cout<<pinids.first<<" "<<pinids.second<<endl;
  PlexPinDiodeId pinid;
  if (gain==1) {
    pinid = pinids.first; // high
  } else {
    pinid = pinids.second; // low
  }

  const PulserRawGainPin *thepin = fPinRP.GetRowByIndex(pinid.GetEncoded());
  cout<<"PIN Id: "<<pinid.GetEncoded()<<endl;
  if (!thepin) cout<<"Oh Bugger"<<endl;
  x = thepin->GetMean();
  ex = thepin->GetError();

  TGraphErrors *g = new TGraphErrors(num,x,y,ex,ey);
  g->SetMarkerStyle(2);
  g->SetMarkerColor(2);
  g->SetTitle("");
  g->Draw(draw_opt);
  char name[100];
  sprintf(name,"PIN Diode, LED %d gain %d",led,gain);
  if (g->GetXaxis()) g->GetXaxis()->SetTitle(name);
  sprintf(name,"Plane %d Strip %d %s",plane,strip,StripEnd::AsString(end));
  if (g->GetXaxis()) g->GetYaxis()->SetTitle(name);
  g->Draw(draw_opt);
  return g;

}

TGraphErrors * Plotter::Plot(int gain,
                             float *x, float *ex, char *xtitle,
                             int led, char *draw_opt)
{

  const float *y = 0;
  const float *ey = 0;
  int num;
  // Get Pin (Y)
  PlexLedId ledid(Detector::kCalDet, 0, led);
  cout<<ledid<<endl;
  std::pair<PlexPinDiodeId,PlexPinDiodeId> pinids = myGetPinDiodeIds(fPlex,ledid);
  cout<<pinids.first<<" "<<pinids.second<<endl;
  PlexPinDiodeId pinid;
  if (gain==1) {
    pinid = pinids.first; // high
  } else {
    pinid = pinids.second; // low
  }

  const PulserRawGainPin *thepin = fPinRP.GetRowByIndex(pinid.GetEncoded());
  cout<<"PIN Id: "<<pinid.GetEncoded()<<endl;
  if (!thepin) cout<<"Oh Bugger"<<endl;
  y = thepin->GetMean();
  ey = thepin->GetError();
  num = thepin->GetNumPoints();

  TGraphErrors *g = new TGraphErrors(num,x,y,ex,ey);
  g->SetMarkerStyle(2);
  g->SetMarkerColor(2);
  g->SetTitle("");
  g->Draw(draw_opt);
  char name[100];
  sprintf(name,"PIN Diode, LED %d gain %d",led,gain);
  if (g->GetYaxis()) g->GetYaxis()->SetTitle(name);
 
  if (g->GetXaxis()) g->GetXaxis()->SetTitle(xtitle);
  g->Draw(draw_opt);
  return g;


}

TGraphErrors * Plotter::Plot(int gain,
                             int gain2,
                             int led, char *draw_opt)
{
 const float *x = 0;
  const float *ex = 0;
  const float *y = 0;
  const float *ey = 0;
  int num;

  // Get Pin (Y)
  PlexLedId ledid(Detector::kCalDet, 0, led);
  cout<<ledid<<endl;
  std::pair<PlexPinDiodeId,PlexPinDiodeId> pinids = myGetPinDiodeIds(fPlex,ledid);
  cout<<pinids.first<<" "<<pinids.second<<endl;
  PlexPinDiodeId pinid;
  if (gain==1) {
    pinid = pinids.first; // high
  } else {
    pinid = pinids.second; // low
  }

  const PulserRawGainPin *thepin = fPinRP.GetRowByIndex(pinid.GetEncoded());
  cout<<"PIN Id: "<<pinid.GetEncoded()<<endl;
  if (!thepin) cout<<"Oh Bugger"<<endl;
  y = thepin->GetMean();
  ey = thepin->GetError();
  num = thepin->GetNumPoints();

 // Get Pin (X)
  {
    PlexLedId ledid(Detector::kCalDet, 0, led);
    cout<<ledid<<endl;
    std::pair<PlexPinDiodeId,PlexPinDiodeId> pinids = myGetPinDiodeIds(fPlex,ledid);
    cout<<pinids.first<<" "<<pinids.second<<endl;
    PlexPinDiodeId pinid;
    if (gain2==1) {
      pinid = pinids.first; // high
    } else {
      pinid = pinids.second; // low
    }
    
    const PulserRawGainPin *thepin = fPinRP.GetRowByIndex(pinid.GetEncoded());
    cout<<"PIN Id: "<<pinid.GetEncoded()<<endl;
    if (!thepin) cout<<"Oh Bugger"<<endl;
    x = thepin->GetMean();
    ex = thepin->GetError();
  }
  TGraphErrors *g = new TGraphErrors(num,x,y,ex,ey);
  g->SetMarkerStyle(2);
  g->SetMarkerColor(2);
  g->SetTitle("");
  g->Draw(draw_opt);
  char name[100];
  sprintf(name,"PIN Diode, LED %d gain %d",led,gain);
  if (g->GetYaxis()) g->GetYaxis()->SetTitle(name);
  sprintf(name,"PIN Diode, LED %d gain %d",led,gain2);
  if (g->GetXaxis()) g->GetXaxis()->SetTitle(name);
  g->Draw(draw_opt);
  return g;



}



std::pair<PlexPinDiodeId,PlexPinDiodeId> Plotter::myGetPinDiodeIds(PlexHandle plex, PlexLedId led)
{
   Int_t highpinchads[6]={4274,4146,4402,4786,4914,5042};
   Int_t lowpinchads[6]={4242,4114,4370,4626,4882,5010};
   PlexPinDiodeId high = 
     plex.GetPinDiodeId(RawChannelId(Detector::kCalDet,ElecType::kVA,
                                     0,highpinchads[led.GetLedInBox()-1]));

   PlexPinDiodeId low = 
     plex.GetPinDiodeId(RawChannelId(Detector::kCalDet,ElecType::kVA,
                                     0,lowpinchads[led.GetLedInBox()-1]));

    return std::make_pair(high,low);
}

PlexLedId Plotter::myGetLedId(PlexHandle plex, PlexStripEndId seid)
{
  int led = getLed(plex.GetRawChannelId(seid).GetCrate(),seid.GetStrip(),seid.GetPlane());
  return PlexLedId(Detector::kCalDet, 0 , led);

}

int Plotter::getLed(int crate,int strip, int plane)
{
  switch (crate)
    {
    case 0:
      if (strip<5) {
        if (plane==40) return 2;
        if ((plane%4)==0||(plane%4)==1) return 1;
        return 2;
      } else if (strip<15) {
        if (plane==50) return 3;
        if (plane==40) return 1;
        if ((plane%4)==0||(plane%4)==1) return 2;
        return 1;
      } else {
        if (plane==50) return 1;
        return 3;
      }
    case 1:
      if (strip<11) {
        if (plane==56||plane==58) return 4;
        return 6;
      } else if (strip<21) {
        if (plane==58) return 6;
        if ((plane%4)==0||(plane%4)==1) return 5;
        return 4;
      } else {
        if ((plane%4)==0||(plane%4)==1) return 4;
        return 5;
      }
    default:
      return 0;
    }

}

void Plotter::Residuals(TGraphErrors *g)
{
  // Canvas...
  //  if (!fRes) 
fRes = new TCanvas("residuals","Residuals",600,0,600,600);
//fRes->cd();

  double *x = g->GetX();
  double *y = g->GetY();
  int num = g->GetN();

  // Don't fit above 8000 ADCs
  int pos; 
  for (pos = num-1;;pos--) {
    if (y[pos]<8000) break;
  }
  for (;;pos--) {
    if (x[pos]<8000) break;
  }
  double xmax = x[pos];
  TF1 *func = new TF1("func","[0] + [1]*x",0,xmax);
  g->Fit("func","R");

  // residuals.
  double *resx = new double[num];
  double *resy = new double[num];

  for (int i=0;i<num;i++) {
    resx[i] = x[i];
    resy[i] = y[i] - func->GetParameter(0) - x[i]*func->GetParameter(1);
    if (y[i]>0.1) {
      resy[i]/=y[i];
    } else {
      resy[i] = 0;
    }
  }

  TGraph *gres = new TGraph(num,resx,resy);
  gres->SetMarkerStyle(2);
  gres->SetMarkerColor(4);
  gres->SetMaximum(0.1);
  gres->SetMinimum(-0.1);
  gres->SetTitle("Residuals");
  gres->Draw("AP");
  if (g->GetXaxis()) {
    gres->GetXaxis()->SetTitle(g->GetXaxis()->GetTitle());
    gres->GetYaxis()->SetTitle(g->GetYaxis()->GetTitle());
    gres->Draw("AP");
  }

  delete [] resx;
  delete [] resy;
}

---