KVTGIDGrid.cpp

/*
$Id: KVTGIDGrid.cpp,v 1.2 2009/03/03 14:27:15 franklan Exp $
$Revision: 1.2 $
$Date: 2009/03/03 14:27:15 $
*/
 
//Created by KVClassFactory on Fri Apr 25 15:55:08 2008
//Author: franklan
 
#include "KVTGIDGrid.h"
#include "KVTGID.h"
#include "KVIDGridManager.h"
 
using namespace std;
 
ClassImp(KVTGIDGrid)
 
 
 
 
KVTGIDGrid::KVTGIDGrid()
{
   // Default constructor
   fTGID = 0;
}
 
 
 
 
KVTGIDGrid::KVTGIDGrid(KVTGID* tgid, KVIDZAGrid* original)
   : KVIDZAGrid()
{
   // Create a grid in order to visualize the results of a fit of an identification grid.
   // We keep any "cuts" from the original grid and add them to the fitted grid
 
   fTGID = tgid;
   SetVarX(fTGID->GetVarX());
   SetVarY(fTGID->GetVarY());
   SetRuns(fTGID->GetValidRuns());
 
//    TList   *idtels = (TList *)fTGID->GetIDTelescopes();
//    TObject *idt    = NULL;
//     TIter next_idt( idtels );
//    while( (idt = next_idt()) ) fTelescopes.Add( idt );
//     delete idtels;
 
   TString name(fTGID->GetName());
   Int_t ver = 1;
   if (gIDGridManager->GetGrid(name.Data())) {
      do {
         name.Form("%s_%d", fTGID->GetName(), ver);
         ver++;
      }
      while (gIDGridManager->GetGrid(name.Data()));
   }
   SetName(name.Data());
   if (original) original->GetCuts()->Copy(fCuts);
 
   GetParameters()->SetValue("HasTGID", 1);
   FindAxisLimits();
}
 
 
 
 
KVTGIDGrid::~KVTGIDGrid()
{
   // Destructor
}
 
 
 
 
 
void KVTGIDGrid::WriteToAsciiFile(ofstream& gridfile)
{
   // Write parameters of LTG fit used to generate grid in gridfile.
   KVIDGraph::WriteToAsciiFile(gridfile);
   if (fTGID) fTGID->WriteToAsciiFile(gridfile);
}
 
 
 
 
 
void KVTGIDGrid::ReadFromAsciiFile(ifstream& gridfile)
{
   // Read grid and parameters of LTG fit used to generate grid in gridfile.
   KVIDGraph::ReadFromAsciiFile(gridfile);
   KVString line;
   if (fPar.GetIntValue("HasTGID")) {
      line.ReadLine(gridfile);
      if (line.BeginsWith("++KVTGID"))  fTGID = KVTGID::ReadFromAsciiFile(GetName(), gridfile);
   }
   FindAxisLimits();
}
 
 
 
 
 
void KVTGIDGrid::Generate(Double_t xmax, Double_t xmin, Int_t ID_min, Int_t ID_max, Int_t npoints, Bool_t logscale)
{
   //(Re)generate ID grid from the functional using current values of parameters.
   //
   //      xmin, xmax - min and max values of 'x' coordinate used in corresponding
   //                                      identification map
   //      ID_min, ID_max - min and max 'ID' of lines in grid
   //      npoints    - number of points in each line
   //
   //DEFAULT ARGUMENT VALUES:
   //=========================
   //   xmin = 0.
   // If the grid gets "ugly" for small x, you might try setting xmin = fitted X-pedestal
   // This is the value of parameter 7 in most cases....
   //
   //  ID_min, ID_max = min and max ID of lines on which functional was fitted
   //  npoints = 100
   //
   // if logscale=kTRUE (default is kFALSE) lines are generated with more points at
   // the beginning of the lines than at the end.
 
   if (!GetTGID()) {
      Error("Generate", "No parameter stored to initialize a functional !");
      return;
   }
 
   ID_min = (ID_min ? ID_min : (Int_t) fTGID->GetIDmin());
   ID_max = (ID_max ? ID_max : (Int_t) fTGID->GetIDmax());
 
   // clear existing lines
   Bool_t was_drawn = kFALSE;
   TVirtualPad* last_pad = fPad;
   if (fPad) {
      UnDraw();
      was_drawn = kTRUE;
   }
 
   // delete all previous identification lines
   fIdentifiers.Delete();
   fXmin = fYmin = fXmax = fYmax = 0;
 
   for (Int_t ID = ID_min; ID <= ID_max; ID++) {
      fTGID->AddLineToGrid(this, ID, npoints, xmin, xmax, logscale);
   }
   if (was_drawn) {
      last_pad->cd();
      Draw();
   }
   FindAxisLimits();
}
 
 
 
 
const KVTGID* KVTGIDGrid::GetTGID() const
{
   // a completer mais j'ai la fleme...
   return fTGID;
}
 
 
 
 
void KVTGIDGrid::AddIdentifier(KVIDentifier* id)
{
   id->SetLineColor(kGray + 2);
   id->SetLineStyle(2);
   KVIDGraph::AddIdentifier(id);
}