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

Go to the documentation of this file.

/* @version $Revision: 1566 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-03-08 14:00:10 +0100 (Thu, 08 Mar 2012) $ */
#include <TROOT.h>
#include <TPaveLabel.h>
#include <TRint.h>
#include <TCanvas.h>
#include <TTree.h>
#include <TFile.h>
#include <TSystem.h>
#include <TGeoManager.h>
#include <TGeoMatrix.h>
#include <TCanvas.h>
#include <iostream>
#include <sstream>
#include <Toolbox.hh>

#include "event/Run.hh"
#include "event/Event.hh"

#include "geometry/Board.hh"
#include "geometry/Chip.hh"
#include "geometry/Dif.hh"
#include "geometry/Chamber.hh"
#include "geometry/Detector.hh"
#include "geometry/GassiplexChamber1.hh"
#include "geometry/GassiplexChamber4.hh"
#include "geometry/Hardroc1Chamber.hh"

#include "xml/XMLTool.hh"

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


#include <TApplication.h>
#include <TGClient.h>
#include <TGListBox.h>
#include <TList.h>



using namespace std;


//______________________________________________________________________________
int main(int argc, char **argv)

{  //un-named script to be run in ROOT framework

  FILELog::ReportingLevel() = FILELog::FromString(INFO);

  //----- Control usage and option
  if ( argc != 3 ) {
   FILE_LOG(logERROR)  << "usage: channelDisplay streeFile chamberId " << endl;
   exit(1);
  }

  string steerName;
  steerName.assign(argv[1]);

  string chamberId;
  chamberId.assign(argv[2]);
  ui32 chamId = atoi(chamberId.c_str());

  TApplication *theApp;

  argc = 1;
  theApp = new TRint("App", &argc, argv);


  /******************************************/
  /******************************************/
  /******************************************/

// test if sterr 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);
  }
// create Run, set detector and steerDesc from steer XML file
    SteerDesc steerDesc;
    XMLTool xml;
    xml.parse(steerDesc, steerName);

    Detector detector;
    detector.build(steerDesc);

    Run run(detector);

    steerDesc.setRun(run);

    fclose(file);


  /******************************************/
  /******************************************/
  /******************************************/
  /******************************************/
  /******************************************/
  /******************************************/


   string txt = "";

   int xLength, yLength = 0;

   try {
    const Chamber &chamber = detector.getChamberById(chamId);
    xLength = int(chamber.getXLength());// + chamber.getX();
    yLength = int(chamber.getYLength());// + chamber.getY();
    if ( chamber.getZrotation() == 90 )
     {
         int foo = xLength;
         xLength = yLength;
         yLength = foo;
     }



//   TCanvas *c1 = new TCanvas("c1","Analyze.mac",yLength*20,xLength*20);
     int hMax = 1000;
     int vMax = 1000;
     int hMin = 100;
     int vMin = 100;
     int hCanvas,vCanvas;

                 if ( yLength > xLength)
     {
       vCanvas = vMax;
       float coef = (float) yLength / (float)xLength;
                         cout << "--yLength["<< yLength << "]" << endl;
                         cout << "--xLength["<< xLength << "]" << endl;
                         cout << "--coef["<< coef << "]" << endl;
       hCanvas = int(hMax / coef);
       if  ( hCanvas < hMin ) { hCanvas = hMin; }

     }
     else
     {
       hCanvas = hMax;
       float coef = (float)xLength / (float)yLength;
                         cout << "--yLength["<< yLength << "]" << endl;
                         cout << "--xLength["<< xLength << "]" << endl;
                         cout << "--coef["<< coef << "]" << endl;
       vCanvas = int(vMax / coef);
       if  ( vCanvas < vMin ) { vCanvas = vMin; }
     }
     //hCanvas = 1600;
     //vCanvas = 1600;


     FILE_LOG(logINFO) << "Canvas size hCanvas[" << hCanvas << "], vCanvas[" << vCanvas << "]" << endl;
     TCanvas *c1 = new TCanvas("c1","Analyze.mac",vCanvas,hCanvas);


   FILE_LOG(logINFO) << "Chamber size xLength[" << xLength << "], yLenght[" << yLength << "]" << endl;

    c1->Range(0,0,yLength,xLength);// (xmin,ymin,xmax,ymax
    int nbChannel = 0;
    int color = 0;
    const BoardMap_t &boards = chamber.getBoards();
    for (BoardMap_t::const_iterator boardIt = boards.begin(); boardIt != boards.end(); ++boardIt)
    {
      Board& board = *(boardIt->second);
      const ChipMap_t &chips = board.getChips();
      for (ChipMap_t::const_iterator chipIt = chips.begin(); chipIt != chips.end(); ++chipIt)
      {
        Chip& chip = *(chipIt->second);
     
                                        ////cout << " chipId[" << chip.getId() << "] getX()[" << chip.getX() << "]" << endl;
                                        //cout << " chipId[" << chip.getId() << "] getY()[" << chip.getY() << "]" << endl;
           float x=chip.getLeftBottomX();
           float y=chip.getLeftBottomY();
           float z=chamber.getZ() ;

           stringstream currentpad;
           currentpad<<chip.getId();

           float x2=chip.getRightUpperX();
           float y2=chip.getRightUpperY();


           FILE_LOG(logDEBUG1) << "-----chip id[" << chip.getId()<< "]" << endl;
           FILE_LOG(logDEBUG1) << "x[" << x << "], x2[" << x2 << "]" << endl;
           FILE_LOG(logDEBUG1) << "y[" << y << "], y2[" << y2 << "]" << endl;

           // On inverse x et y pour repere micromegas
           TPaveLabel *pl1= new TPaveLabel(y+yLength/2, x+xLength/2,y2+yLength/2,  x2+xLength/2,currentpad.str().c_str()); // x1 x2 y1 y2
           if ( chamber.getType().find("STRIP") != string::npos)
           { 
            pl1->SetTextSize(1);
//            pl1->SetTextAngle(90);
           }
           //TPaveLabel *pl1= new TPaveLabel(y, x,y2,  x2,""); // x1 x2 y1 y2
           if ( chamber.getType().find("GASSIPLEX") != string::npos)
           {
               // GASSIPLEX
                 pl1->SetFillColor(board.getId() + 3);
           }
           else {
               //HARDROC
               pl1->SetFillColor(board.getId() + chip.getId() + 2);
           }
           //pl1->SetFillColor(3);
           pl1->Draw();
         }
       }

    Toolbox::printRepere(chamber.getXrotation(),chamber.getYrotation(),chamber.getZrotation());


     theApp->Run();
     delete theApp;
     // run ROOT application
   } //end try
   catch ( MicroException e )
   {
   FILE_LOG(logERROR) << e.getMessage() << endl;
   }

   return 0;
}

---