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

Go to the documentation of this file.

/* @version $Revision: 1678 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-04-26 10:59:36 +0200 (Thu, 26 Apr 2012) $ */
#include <TROOT.h>
#include <TRint.h>
#include <TKey.h>
#include <TCanvas.h>
#include <TTree.h>
#include <TFile.h>
#include <TSystem.h>
#include <TGeoManager.h>
#include <TGeoMatrix.h>
#include <iostream>
#include <sstream>
#include "Log.hh"

//#include "Run.hh"
#include "Chamber.hh"
//#include "Event.hh"
#include "Detector.hh"
#include "DiracChamber1.hh"
#include "GassiplexChamber1.hh"
#include "GassiplexChamber4.hh"
#include "Hardroc1Chamber.hh"
//#include "XMLTool.hh"
#include "MicroException.hh"
#include <iostream>


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

#include "root/MTRun.hh"
#include "root/MTEvent.hh"
#include "root/MTChannel.hh"


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: EventDisplay  rootFile eventNum" << endl;
   exit(1);
  }



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

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

  int eventNum = atoi(argv[2]);

  TApplication *theApp;

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


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


  TFile *f = new TFile(rootName.c_str());
  if ( f->IsZombie()) {
    FILE_LOG(logERROR) << "Root file ["<< rootName.c_str() << "] does not exist" << endl;
    exit(-1);
  }


gSystem->Load("libGeom");
TGeoManager *frame =new TGeoManager("Frame","mM event displayer");
frame->SetMaxVisNodes(40000);
//cout << "--------max vis node["  <<frame->GetMaxVisNodes() << "]" << endl;

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

map<string,bool> node;

MTEvent *evt = new MTEvent();
TIter nextkey(f->GetListOfKeys());
TKey *key;
int padNumber = 0;
while ((key = (TKey*)nextkey()))
    {
    TTree *tree = (TTree*)key->ReadObj();
    MTRun* run = (MTRun*)tree->GetUserInfo()->FindObject("MTRun");

//    FILE_LOG(logINFO) << "Number of reoncstruct events: " << tree->GetEntries() << endl;

    TBranch *branch= tree->GetBranch("MTEvent");
    branch->SetAddress(&evt);

                Int_t globalTriggerCounter = -1;
                Int_t nbEvent = tree->GetEntries();
                bool find = false;
    tree->GetEntry(eventNum);
    std::vector<UInt_t>& dtValid = evt->GetDtCandidate();
      
    UInt_t dt = 0;
    if (dtValid.size() >0) 
    {
      dt = dtValid[0];
      cout << "valid dt [" << dt << "]" << endl;
    }
    
    Int_t NC=evt->GetNchannel();
    frame->SetMaxVisNodes(NC);

    float x,y,z;//                   694 node MaxMax
    for(int i=0;i<NC  ;i++)
    //for(int i=0;i<10  ;i++)
    {
      //if ( padNumber < 693)
      //cout << "pad number]" << padNumber << "]" << endl;
      if ( padNumber < 169300)
      {
        int nodes = world->CountNodes();
        //cout << "nodes]" << nodes << "]" << endl;

        //the only known way to get the information from the channels number i
        MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(i);
        UInt_t t2 = channel->GetBcId_Dif();
        UInt_t t3 = channel->GetBcId_Hit();
        UInt_t dtcurr = t2 - t3;

        if ( dtcurr == dt )
        {

          float xLength = .098 * 4; //cm
          float yLength = .098 * 4; // cm
          float zLength = .12; // cm
          // ChannelPad name
          string padID="";
          stringstream currentpad;// (stringstream::in | stringstream::out);
          padNumber++;
          currentpad<<padNumber;
          padID="pad_"+currentpad.str();

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

          color = channel->GetAnalogValue()+channel->GetDigitalValue()*10;
          ostringstream ox,oy,oz;
          ox << x; oy << y; oz <<z ;
          string key = "x[" + ox.str() + "]y[" + oy.str() + "]z[" + oz.str() +"]" ;
          //cout << "key[" << key <<"]" <<endl;
          if ( node[key]==false)
          {
            world->AddNode(pad,1,new TGeoTranslation(x,y,z));
            node[key] = true;
            color = color * 2;
          }
          pad->SetLineColor(color);

        }
      }
        else
        {
          FILE_LOG(logERROR)   << "------------------Too Many Nodes !" << endl<< endl << endl ;
          break;
        }
      }
    }
    frame->SetTopVolume(world);
    frame->CloseGeometry();

    world->Draw();

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

    return 0;
    }

---