/Users/jacquem/NetBeansProjects/MicromegasFramework/trunk/src/analyse/geometryDisplay.C

Go to the documentation of this file.

#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 <TreeClass.hh>
#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 "HardRockChamber6.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);
      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="";
             stringstream currentpad;// (stringstream::in | stringstream::out);
             currentpad<<channel.getSoftId();
             padID="pad_"+currentpad.str();

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


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

             color = 2;

             if ( channel.getHardId() % 96 == 0  )
             { color = 4; }

             color =  channel.getHardId() % 96;
             cout << color << endl;


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

  world->Draw();


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

   return 0;
}

---