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

Go to the documentation of this file.

#include <TROOT.h>
#include <TRint.h>
#include <TCanvas.h>
#include <TTree.h>
#include <TFile.h>
#include <TSystem.h>
#include <TGeoManager.h>
#include <TGeoMatrix.h>
#include <iostream>
#include <TreeClass.hh>
#include <sstream>
#include "Log.hh"

#include "Run.hh"
#include "Chamber.hh"
#include "Event.hh"
#include "Detector.hh"
#include "Centaure.hh"
#include "DiracLabview.hh"
#include "DiracChamber1.hh"
#include "GassiplexChamber1.hh"
#include "GassiplexChamber4.hh"
#include "HardRockChamber6.hh"
#include "XMLTool.hh"
#include "MicroException.hh"


#include <TApplication.h>
#include <TGClient.h>
//#include <TGButton.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 != 4  ) {
   FILE_LOG(logERROR)  << "usage: EventDisplay streeFile rootFile eventNumber" << endl;
   exit(1);
  }



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

  string rootName;
  rootName.assign(argv[2]);

  int event = atoi(argv[3]);

  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;

    Run run(detector);
    if  ( !steerDesc.setRun(run))
    {
        exit(-1);
    }
    fclose(file);

    FILE_LOG(logINFO) << run << endl;



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



  gSystem->Load("libMicro.so");
  TFile *f = new TFile(rootName.c_str());
  if ( f->IsZombie()) {
    FILE_LOG(logERROR) << "Root file ["<< rootName.c_str() << "] does not exist" << endl;
    exit(-1);
  }
  TTree *tree = (TTree *) gROOT->FindObject("t1");


gSystem->Load("libGeom");
TGeoManager *frame =new TGeoManager("Frame","mM event displayer");

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

//TGeoVolume *world = frame->MakeBox("World", med, xLength, yLength, zLength); // x, y , z size unit cm

// -----  create World volume who contains all other.
TGeoVolume *world = frame->MakeBox("World", med, 1, 1, 1); // x, y , z size unit cm

FILE_LOG(logINFO) << "Number of reoncstruct events: " << tree->GetEntries() << endl;
MTEvent *evt = new MTEvent();
TBranch *branch= tree->GetBranch("MTEvent");
branch->SetAddress(&evt);

tree->GetEntry(event);

Int_t NC=evt->GetNchannel();
frame->SetMaxVisNodes(NC);

  float x,y,z;
  for(int i=0;i<NC;i++)
  {
      int nodes = world->CountNodes();

      if ( world->CountNodes() < 1000 )
      {
         //the only known way to get the information from the channels number i
         MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(i);

         float xLength = .098 * 4; //cm
         float yLength = .098 * 4; // cm
         float zLength = .12; // cm

         // ChannelPad name
         string padID="";
         stringstream currentpad;// (stringstream::in | stringstream::out);
         currentpad<<i;
         padID="pad_"+currentpad.str();

         TGeoVolume *pad = frame->MakeBox(padID.c_str(),med,xLength,yLength,zLength );
         z=channel->GetZ();
         x=channel->GetX();//+xOffset;
         y=channel->GetY();//+yOffset;
         //set now the color of the pad
         int color;
         int sMax = 30;//set to limit of saturation
         stringstream ch;
         ch<<i;
         //string command="fChannels.fValue["+ch.str()+"]";

         color = channel->GetValue();
         color=60+40.*min(1.,(1.*color-20)/sMax);
        //world->AddNode(pad,1,new TGeoTranslation(z,x,y));
        world->AddNode(pad,1,new TGeoTranslation(x,y,z));
        pad->SetLineColor(color);
      }
      else
      {
          FILE_LOG(logERROR) << "Too Mnay Nodes !" << endl;

      }

  }
  frame->SetTopVolume(world);
  frame->CloseGeometry();

  world->Draw();

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

   return 0;
}

---