/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/root/beamProfile.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 <vector>

// Root includes
#include "TString.h"
#include "TApplication.h"
#include <TStyle.h>
#include <TFile.h>
#include <TTree.h>
#include <TKey.h>

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

#include "Board.hh"
#include "Chip.hh"
#include "Dif.hh"

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

using namespace std;



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


FILELog::ReportingLevel() = FILELog::FromString(INFO);
string filename ;
ui32 nbEvents = 0;
ui32 nbChannels = 0;

// Manage command options
//if ( argc == 1)
//{
    //cout << endl  << "usage: beamProfile   [root file Name]" << endl ;
    //cout << "help:  beamProfile -h " << endl;
    //exit(0);
//}
BeamProfile manager(argc,argv);

filename = manager.getFileName();

if ( filename.empty())
{
    cout << "a root file name is mandatory. Please use -f option to set it. " << endl;
    exit(0);
}

// Build detector from XML file
        Detector detector;
        // Get detector information
        if ( ! manager.getXMLFileName().empty())
        {
            SteerDesc steerDesc;
            XMLTool xml;
            xml.parse(steerDesc, manager.getXMLFileName());
            detector.build(steerDesc);
            manager.SetDetector(&detector);
        }


// Create ROOT Application
TApplication theApp("App", 0, NULL);
gStyle->SetPalette(1);

// Open root file
int nchannel = 0;

TFile f(filename.c_str());
cout<<"reading File : "<<filename<<endl;

// Get TTree list form root file
TIter nextkey(f.GetListOfKeys());
TKey *key;

// Loop over TTree ( one TTree by rawfile )
while (key = (TKey*)nextkey())
{
    TTree *tree = (TTree*)key->ReadObj();
    if ( tree == NULL) { exit(0);} ;

          if  (manager.getShowRunInfo())
                {
                    //Get Run information ( name, date, dif and chip confiration )
        MTRun* run = (MTRun*)tree->GetUserInfo()->FindObject("MTRun");
        if ( run != NULL) run->Info();
          }
                 FILE_LOG(logINFO  ) << endl << "TTree[" <<tree->GetName() << "]" << endl;

                // Get number of event for current  TTree
    int nEvent = tree->GetEntries();

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


                // Loop over events
    for (int i=0;i<nEvent;i++)
    {
        nbEvents++;
                         if ( nbEvents % 100 == 0) FILE_LOG(logINFO  ) << "Event " << i+1 <<" / " << nEvent << "\r" << flush;
                // Print computate status
                //if (i%10 == 0){cout<<(float)(i*100)/(float)(nEvent)<<" %"<<"\r"<<flush;}

                // Get ith event
        tree->GetEntry(i);
                                // Get number of MTChannel for this event
        nchannel = evt->GetNchannel();
                    manager.FillHotChamber(*evt);

                                // loop over MTChannel
        for (int j=0;j<nchannel;j++)
        {
            nbChannels++;
            MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
                                                // Get x,y position

                             manager.FillXyDist(*channel);

                             manager.FillHitPerChipDist(*channel);

                       manager.FillHitPerChannel(*channel);

                             manager.FillDeltaT(*channel);

                            manager.FillHotChip(*channel);

            //manager.FillTimestamp(*channel);

        }
    } // End of event loop.
} // end while over TTree

manager.DrawXyDist();
manager.DrawHitPerChipDist();
manager.DrawHitPerChannelDist();
manager.DrawDeltaT();
manager.DrawHotChamber();


manager.PrintHotChip();
cout << endl <<endl << "----------- SUCCESS ---------------" << endl;
cout << endl << "----------- " << nbEvents <<" events, " << nbChannels << " channels. ---------------" << endl;
cout << " Close Root  ( or  press Ctrl^C ) to finish."<< endl;

theApp.Run();
return (EXIT_SUCCESS);
}

---