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

Go to the documentation of this file.

/* @version $Revision: 1602 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-03-26 17:21:50 +0200 (Mon, 26 Mar 2012) $ */

#include <iostream>
#include <sstream>
#include "Log.hh"

#include "Run.hh"
#include "Board.hh"
#include "Chip.hh"
#include "Dif.hh"
#include "Chamber.hh"
#include "Event.hh"
#include "Detector.hh"
#include "Centaure.hh"
#include "GassiplexChamber1.hh"
#include "GassiplexChamber4.hh"
#include "Hardroc1Chamber.hh"
#include "XMLTool.hh"
#include "Toolbox.hh"
#include "MicroException.hh"


#include <TApplication.h>
#include <TGClient.h>
#include <TGListBox.h>
#include <TList.h>
#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>


/*
 Le fait d'échanger les X et les Y par rapport a ce root attends, fabrique la translation desiree pour micromegas
 *soit Vue de PAD en direct et axe X horizontal
 **/

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 != 4 ) {
   FILE_LOG(logERROR)  << "usage: channelDisplay streeFile chamberId boardId" << endl;
   exit(1);
  }

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

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


  string sboardId;
  sboardId.assign(argv[3]);
  ui32 boardId = atoi(sboardId.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;
    if ( detector.build(steerDesc) == steerDesc.chambers.size() )
    {


         Run run(detector);

         steerDesc.setRun(run);

         fclose(file);


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


       string txt = "";

       int xLength, yLength = 0;

       try {
         const Chamber &chamber = detector.getChamberById(chamId);
         const BoardMap_t &boards = chamber.getBoards();


         xLength = int((*(boards.find(boardId)->second)).getXLength()) ;//+ chamber.getX();
         yLength = int((*(boards.find(boardId)->second)).getYLength()) ;//+ chamber.getX();


         if ( chamber.getZrotation() == 90 ) 
         {
             int foo = xLength;
             xLength = yLength;
             yLength = foo; 
         }

        
         //xLength = int(chamber.getXLength()) ;//+ chamber.getX();
         //yLength = int(chamber.getYLength()) ;//+ chamber.getY();

         int hMax = 1600;
         int vMax = 1600;
         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);
         }
         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);
         }

         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;
         for (BoardMap_t::const_iterator boardIt = boards.begin(); boardIt != boards.end(); ++boardIt) 
         {
           Board& board = *(boardIt->second);
         if ( board.getId()==boardId)
           {
                 const ChipMap_t &chips = board.getChips();
                 for (ChipMap_t::const_iterator chipIt = chips.begin(); chipIt != chips.end(); ++chipIt)
                 {
                   Chip& chip = *(chipIt->second);

                   cout << "board [" << board.getId() << "] chip[" << chip.getId() << "]" <<  endl;
                   const ChannelMap_t &channels = chip.getChannels();
                   for (ChannelMap_t::const_iterator channelIt = channels.begin(); channelIt != channels.end(); ++channelIt)
                   {
                     Channel& channel = *(channelIt->second);
{

                      
                   cout << "chip[" << chip.getId() << "] channel [" << channel.getHardId() << "] Xpos[" << channel.getX() << "] yPos["  << channel.getY() << "]" <<   endl;
  
                     float x=channel.getLeftBottomX();
                     float y=channel.getLeftBottomY();
                     float z=channel.getZ() ;


                     stringstream currentpad;
                     currentpad<<channel.getHardId();

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

                     TPaveLabel *pl1= new TPaveLabel(y+yLength/2-board.getY(), x+xLength/2-board.getX(),y2+yLength/2-board.getY(),  x2+xLength/2-board.getX() ,currentpad.str().c_str()); // x1 x2 y1 y2
                     if ( chamber.getType().find("GASSIPLEX") != string::npos)
                     {
                         // GASSIPLEX
                        //   pl1->SetFillColor(board.getId() + 3);
      cout << "------DEBUG GASSIPLERX" << endl;
                     }
                     else {
                         //HARDROC Or  MICROROC
                    //      pl1->SetFillColor(board.getId() - 11);
                         //pl1->SetFillColor( 3);
                     }
                     //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;
}

---