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

Go to the documentation of this file.

/* @version $Revision: 1566 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-03-08 14:00:10 +0100 (Thu, 08 Mar 2012) $ */
/* 
 * File:   beamProfile.cpp
 * Author: jacquem
 *
 * Created on May 6, 2010, 10:03 AM
 */

// System includes
#include <iostream>
#include <stdlib.h>
#include <map>
#include <string>

// Root includes
#include "TROOT.h"
#include "TGraph.h"
#include "TStyle.h"
#include <TFile.h>
#include <TTree.h>
#include <TLegend.h>
#include <TKey.h>
#include <TH1I.h>
#include <TH2I.h>
#include <TText.h>
#include <TCanvas.h>
#include <TPaveLabel.h>

// Framework inlcudes
#include "Run.hh"
#include "Detector.hh"
#include "XMLTool.hh"
#include "Toolbox.hh"
#include "MicroException.hh"
#include "Log.hh"

#include "root/MTRun.hh"
#include "root/MTEvent.hh"
#include "root/MTChannel.hh"
#include "root/CaloEvent.hh"
#include "root/CaloHit.hh"
#include "tools/Arguments.hh"

#include "Chamber.hh"
#include "Detector.hh"
using namespace std;



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


  FILELog::ReportingLevel() = FILELog::FromString(INFO);

  Arguments arg(argc,argv,"monitoring inputFile outputFile -l evtNum");
  if (arg.getNbOfArgs() != 2 ) { arg.usage() ; exit(-1);}
  string rootName = arg.getArgument(1);
  string outputFileName = arg.getArgument(2);

  TFile *output= new TFile(outputFileName.c_str(),"RECREATE");


  Char_t flag[20];
  TTree *flagTree = new TTree("FLAG", "evant flag");
  flagTree->Branch("Flag", flag,"flag/C");
   

  UInt_t eventLimit = 0;
  if ( arg.getOption("-l").size() != 0 )
  {
    eventLimit = atoi(arg.getOption("-l").c_str());
  }

  bool construct = false;

  TFile f(rootName.c_str());


  TH2I* xy = NULL;
  TIter nextkey(f.GetListOfKeys());
  TKey *key;

  TH2I* hists[100];
  while (key = (TKey*)nextkey()) 
  {
    TTree *tree = (TTree*)key->ReadObj();
    FILE_LOG(logINFO  ) << "New TTree [" << tree->GetName() << endl;
    Int_t max = 10400*96; // 10400 micro * 96 pad

    CaloEvent *evt =  new CaloEvent();
    TBranch *branch= tree->GetBranch("CaloEvent");
    branch->SetAddress(&evt);
    CaloHit* channel =NULL;
    int nbEntries = tree->GetEntries();
    if ( eventLimit == 0 ) { eventLimit = 100; }
    unsigned int interval = (unsigned int) eventLimit / 10;
    if (interval == 0) interval = 1;
    for ( int evtNum = 0; evtNum <eventLimit ; evtNum++)
    {
      tree->GetEntry(evtNum);
      string name =  "event_" ; 
      Toolbox::addIntToString(name,evt->GetEventId());
      TH2I* xy = new TH2I("name", "title",96,0,max,96,0,max);
      hists[evtNum] = xy;
      int nbHits = evt->GetNHits();

      for(int i=0 ; i < nbHits ; i++)
      {
        channel = (CaloHit*)evt->GetHits()->UncheckedAt(i);
        xy->Fill(channel->GetY(),channel->GetX(),1);
      }
        FILE_LOG(logINFO  ) << "Write Event[" << evtNum << "] \t\t\t\t\t\t\t\t\t\t\t\r" << flush;
    }
 }



  output->cd();


  for ( unsigned int i =  0 ; i < eventLimit ; i++)
  {
    hists[i]->Write();
  }

  output->Close();
  f.Close(); 

  cout <<endl << endl << "Done" << endl;
  
return 0;
}


  

---