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

Go to the documentation of this file.

/* @version $Revision: 1779 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-06-15 10:48:21 +0200 (Fri, 15 Jun 2012) $ */
#include <map>
#include "TKey.h"
#include "TList.h"
#include <TROOT.h>
#include <TRint.h>
#include <TKey.h>
#include <TList.h>
#include <TTree.h>
#include <TFile.h>
#include <stdlib.h>
#include <TSystem.h>
#include <iostream>
#include <sstream>
#include "Log.hh"

#include "MicroException.hh"
#include <string>

#include "mysql/Mysql.hh"


#include "root/MTRun.hh"
#include "root/MTTrack.hh"
#include "root/MTChamber.hh"
#include "root/MTDetector.hh"
#include "root/MTBoard.hh"
#include "root/MTChip.hh"
#include "root/MTChannel.hh"
#include "root/MTEvent.hh"
#include "root/MTDif.hh"
#include "root/MTHardroc1Chip.hh"
#include "root/MTHardroc2Chip.hh"
#include "root/MTDiracChip.hh"
#include "root/MTMicrorocChip.hh"


using namespace std;


int main(int argc, char**argv){

   if ( argc !=2  ) {
   FILE_LOG(logERROR)  << "usage: example rootFile " << endl;
   exit(1);
  }

  string rootName;
  rootName.assign(argv[1]);
  int nbHit = 0;
  int nbEvt = 0;

  TFile f(rootName.c_str());
  f.ls();



  TIter nextkey(f.GetListOfKeys());
  TKey *key;
  while (key = (TKey*)nextkey()) 
  {
    TTree *tree = (TTree*)key->ReadObj();                


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

  MTChannel* channel =NULL;
  int nbEntries = tree->GetEntries();

  for ( int evtNum = 0; evtNum < 2 ; evtNum++)
  {
    tree->GetEntry(evtNum);
    nbEvt++;
     bool show = false;




      UInt_t evtMin = 0xffffffff;
      UInt_t evtMax = 0x0;

      int nbChannel = evt->GetNchannel();
      //evt->GetChannels()->Sort();
      for(int i=0;i<nbChannel  ;i++)
      {
        nbHit++;
        channel = (MTChannel*)evt->GetChannels()->UncheckedAt(i);


        i32 deltaT = channel->GetBcId_Dif() - channel->GetBcId_Hit();
        ui32 absTime = channel->GetAbsTime();
        if (deltaT < 0  )
        {
          cout << "deltaT < 0"<<endl;
        }
        show = true;
        cout << "**** Event Information "<< dec << evtNum+1 << "/" << nbEntries<< " ****" <<  endl;
          cout << "      channel num :           " << i <<endl;
          cout << "      deltaT :           " << deltaT <<endl;
          cout << "      bcId dif:           " << channel->GetBcId_Dif() <<endl;
          cout << "      bcId hit:           " << channel->GetBcId_Hit() <<endl;
          cout << "    event id:" << evt->GetEventId() <<endl;

          cout << "      chamber id:         " << channel->GetChamberId() <<endl;
          cout << "      dif id:             " << channel->GetDifId() <<endl;
          cout << "      boardId id:         " << channel->GetBoardId() <<endl;
          cout << "      chip id(format dec):" << dec << channel->GetChipId() <<endl;
          cout << "      hard id    :        " << dec << channel->GetHardId() <<endl;
          cout << "      bcId abs:           " << channel->GetBcId_Abs() <<endl;
          cout << "      absolute time:bcid_abs - (bcid_dif - bcid_hit)       " << absTime <<endl;
          cout << "      nbHit for this absolute time:" << evt->GetNbHitPerAbsTime(channel->GetSoftId().GetBoardId(),absTime)<<endl;
      //}
      //else  {
        //if ( i ==0 ) 
        //{
         // cout << "**** Event Information "<< dec << evtNum+1 << "/" << nbEntries<< " ****" <<  endl;
        //  cout << "----     bcId  abs  " << channel->GetBcId_Abs()  <<endl;
        //  cout << "----      absolute time:bcid_abs + (bcid_dif - bcid_hit)       " << absTime <<endl;
        //}
      //}
      
        if ( absTime > evtMax ) evtMax = absTime;
        if ( absTime < evtMin ) evtMin = absTime;
      } // END LOOP ON CHANNELHIT
     cout << "Min Evt absTime: " << evtMin << endl;
     cout << "Max Evt absTime: " << evtMax << endl;
    } // END LOOP ON EVENT
    cout << "End of TTree[" << tree->GetName() << "]" <<endl;
    cout << endl << "** Nombre d Event["<< nbEvt <<"] ** "<< endl;
    cout << endl << "** Nombre de Hits["<< nbHit <<"] ** "<< endl;

  } // end loop on TKEy
return 0 ;
}

---