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

Go to the documentation of this file.

/* @version $Revision: 1894 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-09-21 10:57:18 +0200 (Fri, 21 Sep 2012) $ */
#include <map>
#include "TKey.h"
#include <TROOT.h>
#include <TStyle.h>
#include <TRint.h>
#include <TKey.h>
#include <TTree.h>
#include <TCanvas.h>
#include "TH1.h"
#include "TF1.h"
#include <TH1I.h>
#include <TH2I.h>
#include <TFile.h>
#include <stdlib.h>
#include <TSystem.h>
#include <iostream>
#include "TApplication.h"
#include <sstream>
#include "Log.hh"
#include <Rtypes.h>
#include <TCanvas.h>
#include "TColor.h"
#include "TFunction.h"
#include "TAttLine.h"
#include "TPaveStats.h"
#include "MicroException.hh"
#include <string>



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

using namespace std;


int main(int argc, char**argv){
        gROOT->Reset();
        
        gROOT->SetStyle("Plain");
        gStyle->SetTitleBorderSize(0);
        gStyle->SetTitleX(0.35);
        gStyle->SetTitleY(0.995);

        gStyle->SetFrameFillStyle(0);
        gStyle->SetOptStat("");
        gStyle->SetOptFit(0111);
        gStyle->SetStatX(0.9);
        gStyle->SetStatY(0.9);
        gStyle->SetStatH(0.1);
        gStyle->SetStatW(0.1);

        gStyle->SetOptStat("") ;
        gStyle->SetOptFit(0000);        
        
        gStyle->SetOptStat("");
        gStyle->SetOptFit(10001);
        
        
        
        
   if ( argc !=2  ) {
   FILE_LOG(logERROR)  << "usage: analog rootFile " << endl;
   exit(1);
  }

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

  TApplication *theApp  = NULL;

  theApp  = new TRint("App",NULL,NULL );


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


  UInt_t nb_ped=0 ;     // cut on number of pedestal event

  TFile f(rootName.c_str());
  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 < 5 ; evtNum++)
    {
      cout << "New dacPed" << endl;
      TH1I * dacPed = new TH1I("ped","",100000,0,50000);
      tree->GetEntry(evtNum);
      nbEvt++;
      UInt_t channel68 = 0;

                  UInt_t absTime = 0 ;

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

        UInt_t adc = channel->GetAnalogValue();
        UInt_t id = channel->GetSoftId().GetValue()  - 304000000;

        cout << "Id[" << id << "] adc" << adc << endl;

                        dacPed->Fill(id,adc);
      }
      output->cd();
      dacPed->Write("");
    }
  }

        
   theApp->Run();
   delete theApp;
    return 0 ;
}

---