/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/root/geometryDisplay.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 "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 "MicroException.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 != 2  ) {
   FILE_LOG(logERROR)  << "usage: geometryDispay streeFile" << endl;
   exit(1);
  }



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

  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);

    FILE_LOG(logINFO) << run << endl;



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


   TCanvas *c1 = new TCanvas("c1","Analyze.mac",620,790);
   c1->Range(-1,0,19,30);
   string txt = "";


TGeoManager *frame =new TGeoManager("Frame","mM event displayer");

TGeoMaterial *mat = new TGeoMaterial("SomeMaterial",1,1,1);
TGeoMedium *med = new TGeoMedium("SomeMedium",1,mat);

// -----  create World volume who contains all other.
TGeoVolume *world = frame->MakeBox("World", med, 1, 1, 1); // x, y , z size unit cm
frame->SetMaxVisNodes(9216 * detector.getNumberOfChambers());  // 10 chambers de 9216 channel

int nbChannel = 0;

    try { 
    const ChamberMap_t &chambers = detector.getChambers();
    int color = 0;
    for (ChamberMap_t::const_iterator chamberIt = chambers.begin(); chamberIt != chambers.end(); ++chamberIt) {
      Chamber& chamber = *(chamberIt->second);

      cout << "========= Cahmber [" << chamber.getId() << "]" <<endl;
//cout << " getXLength[ "<<chamber.getXLength() <<"]" << endl;
//cout << " getYLength[ "<<chamber.getYLength() <<"]" << endl;
//cout << " getMaxX[ "<<chamber.getMaxX() <<"]" << endl;
//cout << " getMaxY[ "<<chamber.getMaxY() <<"]" << endl;
                FILE_LOG(logINFO) <<    endl << "Chamber ID[" << chamber.getId() << "]" << endl;
                string result = (chamber.getXrotation() == false) ? "No rotation over X" : "There is a rotation over X";
    FILE_LOG(logINFO) << result    << endl;
    result = (chamber.getYrotation() == false) ? "No rotation over Y" : "There is a rotation over Y";
    FILE_LOG(logINFO) << result    << endl;


      const BoardMap_t &boards = chamber.getBoards();

        int color = 3;
      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);
          const ChannelMap_t &channels = chip.getChannels();

          for (ChannelMap_t::const_iterator channelIt = channels.begin(); channelIt != channels.end(); ++channelIt) 
                {

             Channel& channel = *(channelIt->second);
             //the only known way to get the information from the channels number i
             //MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(i);
     
             float xHalfLength = .098; //cm
             float yHalfLength = .098; // cm
             float zHalfLength = .12; // cm

             xHalfLength = (float) channel.getXLength() / 2 - .1; //cm
             yHalfLength = (float) channel.getYLength() / 2 - .1 ; // cm
             zHalfLength = .1; // cm


             float x=channel.getX() ;//+xOffset;
             float y=channel.getY() ;//+yOffset;
             float z=channel.getZ() ;


      
             // ChannelPad name
                                                 string padID;
                                                 //padID = channel.getStringSoftId();
                                                 padID = channel.getSoftId().toString();
    
             //TGeoVolume *pad = frame->MakeBox(padID.c_str(),med,zHalfLength,xHalfLength,yHalfLength );
             TGeoVolume *pad = frame->MakeBox(padID.c_str(),med,zHalfLength,yHalfLength,xHalfLength );
             //set now the color of the pad
             stringstream ch;
             //ch<<channel.getSoftId();


             world->AddNode(pad,1,new TGeoTranslation(z,-y,x));
             nbChannel++;

             //color = ( channel.getChip().getBoard().getId()  + 2 + channel.getChip().getId())  ;
             color = ( channel.getChip().getBoard().getId()) + 1  ;




             pad->SetLineColor( color);
             pad->SetFillStyle(1001);
          }
        }
      }
    }
    }
    catch ( MicroException e )
    {
     cout << e.getMessage();
    }
  frame->SetTopVolume(world);
  frame->CloseGeometry();

  world->Draw();


          FILE_LOG(logINFO) << "" << endl;
    FILE_LOG(logINFO) << "                                            X" << endl;
    FILE_LOG(logINFO) << "                                             " << endl;
    FILE_LOG(logINFO) << "                                            /\\" << endl;
    FILE_LOG(logINFO) << "                                             |" << endl;
    FILE_LOG(logINFO) << "                                             |" << endl;
    FILE_LOG(logINFO) << "                                             |" << endl;
    FILE_LOG(logINFO) << "SYSTEM de COORDONNEES: (0,0) en bas a gauche X--->  Y    " << endl;
    FILE_LOG(logINFO) << "  Z sortant." << endl;


 theApp->Run();
 delete theApp;
   // run ROOT application

   return 0;
}

---