/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/root/slowcontrol.cpp

Go to the documentation of this file.

/* @version $Revision: 1337 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2011-10-06 10:48:29 +0200 (Thu, 06 Oct 2011) $ */
#include "event/Run.hh"

#include "xml/XMLTool.hh"

#include "tools/Log.hh"
#include "tools/MicroException.hh"

#include "slowControl/SlowControlManager.hh"
#include "slowControl/SlowControlEntry.hh"
#include "slowControl/BoardSlowControl.hh"
#include "slowControl/PMSlowControl.hh"

#include <iostream>
 #include <time.h>



#include "TROOT.h"
#include "TGraph.h"
#include "TCanvas.h"
#include "TAxis.h"
#include "TTimeStamp.h"
#include "TFrame.h"
#include "TApplication.h"
#include <TRint.h>




using namespace std;

int main(int argc, char **argv)
{
string steerName;


  TApplication *theApp  = NULL;
  // set verbose level
  FILELog::ReportingLevel() = FILELog::FromString(INFO);

  ui32 slowcontrolNotFound = 0;

//----- Control usage and option
  if ( argc != 2  ) {
  }

  // Get command parameter : here steering file
   steerName.assign(argv[1]);
    
  // test if steer file exist
  FILE *file = fopen(steerName.c_str(), "r");
  if(file == NULL){
    FILE_LOG(logERROR) << "Steer file ["<< steerName.c_str() << "] does not exist" << endl;
    exit(0);
    fclose(file);
  }
  else
  {


    //theApp  = new TRint("App", &argc, argv);
    theApp  = new TRint("App",NULL,NULL );
    SteerDesc steerDesc;
    XMLTool xml;
    xml.parse(steerDesc, steerName); // parse steer file and fill steerDesc
    fclose(file);

    string slowControlPath = steerDesc.getSlowControlPath();
    SlowControlManager sc;
    sc.build(steerDesc);
    sc.readDB(slowControlPath);

    FILE_LOG(logINFO) <<  "slowControl entries:[" << sc.getNbEntries() << "]" << endl;

    const Int_t n = sc.getNbEntries();
    Double_t timestamp[n] ;
    Double_t temperature[n] ;
    Double_t pressure[n] ;
    Double_t rapport[n] ;
    Double_t voltage[6][n] ;
    TTimeStamp timeformat[n];

//    Double_t x[n] , y[n];


   UInt_t index = 0;
    for (  std::set<SlowControlEntry>::const_iterator iter = sc.begin(); iter!=sc.end(); iter++)
    {
      const SlowControlEntry& entry = *iter;
      timestamp[index] = entry.getTimeStamp();
      temperature[index] = entry.getTemperature();
      pressure[index] = entry.getPressure();
      time_t val = time_t(( entry.getTimeStamp() / 1000));
      timeformat[index] = TTimeStamp( val,0);
      rapport[index] = entry.getPressure()/entry.getTemperature();

      try 
      {
        voltage[0][index] = entry.getBoard(1,84,13).getDrift().voltage;
        voltage[1][index] = entry.getBoard(1,84,13).getMesh().voltage;
        voltage[2][index] = entry.getPM(1).getParams().v1;
        voltage[3][index] = entry.getPM(1).getParams().v2;
        voltage[4][index] = entry.getBoard(22,0,1).getDrift().voltage;
        voltage[5][index] = entry.getBoard(22,0,1).getMesh().voltage;

//cout << "timestamp     :"<< timestamp[index] << endl;
//cout << "temperature   :"<< temperature[index] << endl;
//cout << "Drift ASU 13  :"<< voltage[0][index] << endl;
//cout << "Mesh  ASU 13  :"<< voltage[1][index] << endl;
 const PMSlowControl &pm = entry.getPM(1);
 const PMParams_t &param = pm.getParams();
//cout << "PM1           :"<< param.v1 << endl << endl;
      }
      catch ( MicroException &e)
      {
        cout << e.getMessage() << endl;
        exit(-1);
      }
      index++;
    }
      


  // Fill Root Graph and display it
  TCanvas *c1 = new TCanvas("c1","A Simple Graph Example",200,10,700,500);

   c1->SetFillColor(42);
   c1->SetGrid();
   c1->Divide(2,3);


 TGraph *gr_temp = new TGraph(n,timestamp,temperature);
 TGraph *gr_pressure_temp = new TGraph(n,timestamp,rapport);
 TGraph *gr_vol_1_84_13 = new TGraph(n,timestamp,voltage[0]);
 TGraph *gr_vol_1_84_14 = new TGraph(n,timestamp,voltage[1]);
 TGraph *gr_vol_1_85_15 = new TGraph(n,timestamp,voltage[2]);
 TGraph *gr_vol_1_85_16 = new TGraph(n,timestamp,voltage[3]);
 TGraph *gr_vol_1_88_17 = new TGraph(n,timestamp,voltage[4]);
 TGraph *gr_vol_1_88_18 = new TGraph(n,timestamp,voltage[5]);

  for (UInt_t i=1;i<=index;i=i+1000) {
    const char * foo =  timeformat[i].AsString("s");
  }

/*

   TGraph *gr_pressure = new TGraph(n,timestamp,pressure);
   gr_temp->SetTitle("SPS H4 August temperature form slow control");
   gr_temp->GetXaxis()->SetTitle("timestamp");
   gr_temp->GetYaxis()->SetTitle("temperature");
   gr_temp->Draw("ACP");

   gr_pressure->SetTitle("SPS H4 August pressure form slow control");
   gr_pressure->GetXaxis()->SetTitle("timestamp");
   gr_pressure->GetYaxis()->SetTitle("pressure");
   gr_pressure->Draw("ACP");
*/


   c1->cd(1);
   gr_vol_1_84_13->SetTitle("Drift voltage ASU 13");
   gr_vol_1_84_13->GetXaxis()->SetTitle("timestamp");
   gr_vol_1_84_13->GetYaxis()->SetTitle("voltage");
   gr_vol_1_84_13->Draw("AL");

   c1->cd(2);
   gr_vol_1_84_14->SetTitle("Mesh voltage ASU 13");
   gr_vol_1_84_14->GetXaxis()->SetTitle("timestamp");
   gr_vol_1_84_14->GetYaxis()->SetTitle("voltage");
   gr_vol_1_84_14->Draw("AL");

   c1->cd(3);
   gr_vol_1_85_15->SetTitle("PM1 v1");
   gr_vol_1_85_15->GetXaxis()->SetTitle("timestamp");
   gr_vol_1_85_15->GetYaxis()->SetTitle("voltage");
   gr_vol_1_85_15->Draw("AL");

   c1->cd(4);
   gr_vol_1_85_16->SetTitle("PM1 v2");
   gr_vol_1_85_16->GetXaxis()->SetTitle("timestamp");
   gr_vol_1_85_16->GetYaxis()->SetTitle("voltage");
   gr_vol_1_85_16->Draw("AL");

   c1->cd(5);
   gr_vol_1_88_17->SetTitle("Strip 1 Drift");
   gr_vol_1_88_17->GetXaxis()->SetTitle("timestamp");
   gr_vol_1_88_17->GetYaxis()->SetTitle("voltage");
   gr_vol_1_88_17->Draw("AL");

   c1->cd(6);
   gr_vol_1_88_18->SetTitle("Strip 1 Mesh");
   gr_vol_1_88_18->GetXaxis()->SetTitle("timestamp");
   gr_vol_1_88_18->GetYaxis()->SetTitle("voltage");
   gr_vol_1_88_18->Draw("AL");

 c1->Update();


  }   

   theApp->Run();
   delete theApp;
  return 0;
}

---