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

Go to the documentation of this file.

/* @version $Revision: 1376 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2011-11-18 16:11:02 +0100 (Fri, 18 Nov 2011) $ */



#include "MicroException.hh"
#include "Mysql.hh"
#include "Toolbox.hh"
#include "CaloRun.hh"
#include "CaloEvent.hh"
#include "CaloHit.hh"

#include "TKey.h"
#include <TKey.h>
#include <TList.h>
#include <TTree.h>
#include <TFile.h>

#include "XMLTool.hh"
#include "geometry/Detector.hh"
#include "geometry/Chamber.hh"
#include "Log.hh"

#include <iostream>
#include <iomanip>

//using namespace std;


int main(int argc, char**argv){
   if ( argc !=2  ) {
   FILE_LOG(logERROR)  << "usage: caloExample rootFile " << endl;
   exit(1);
  }

  // Get detector from SQL and runId
  Detector detector;

  string rootName;
  rootName.assign(argv[1]);
  int nbHit = 0;
  int nbEvt = 0;
  
  TFile f(rootName.c_str());
  TIter nextkey(f.GetListOfKeys());
  TKey *key;

  while (key = (TKey*)nextkey()) 
  {

    TString keyname = key->GetName();

    TTree *tree = (TTree*)key->ReadObj();                
    cout << "---------------------------- NEW TTREE ---------------------------"<< endl;
    nbEvt =0;
    nbHit = 0;

    unsigned int runId = 0;
    try {
      CaloRun *run=dynamic_cast<CaloRun *>(tree->GetUserInfo()->First());
      if ( run != NULL) 
      {
        runId =  run->GetRunId() ;
        cout << "RunId[" << run->GetRunId() << "] " << endl ; 
      }
    }
    catch ( ... ) {}
    Detector detector;
    Toolbox::getDetector(runId,detector);


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

    CaloHit* hit =NULL;
    int nbEntries = tree->GetEntries();
    for ( int evtNum = 0; evtNum < nbEntries ; evtNum++)
    {
      tree->GetEntry(evtNum);
      nbEvt++;
      int nbHits = evt->GetNHits();
      for(int i=0 ; i < nbHits ; i++)
      {
        hit = (CaloHit*)evt->GetHits()->UncheckedAt(i);
        Int_t xpos = hit->GetX();
        Int_t ypos = hit->GetY();
        Int_t layer = hit->GetLayer();
        Chamber& chamber = detector.getChamberById(layer);
        cout << "Hit in chamber id " << layer <<", chamber type " << chamber.getType() << endl;
        
        cout << "\t Nw Hit " << endl;
        hit->GetSolftId().ToString();
        cout << "new hit layer[" << hit->GetLayer() << endl;
        cout << "\thit abs row [" << hit->GetRowInChamber() << setprecision(12) << "], posX[" << xpos << "]"  << endl;
        cout << "\thit abs col [" << hit->GetColInChamber() << setprecision(12) << "], posY[" << ypos << "]"  << endl;
      }


    } // end event
/*
*/
  } // end run ( TTree )

return 0;
}

---