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

#include <map>
#include "TKey.h"
#include <TROOT.h>
#include <TStyle.h>
#include <TRint.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 "TColor.h"
#include "TFunction.h"
#include "TAttLine.h"
#include "TPaveStats.h"
#include "MicroException.hh"
#include <string>
#include "root/MTChannel.hh"
#include "root/MTEvent.hh"

Include dependency graph for analog.cpp:

Go to the source code of this file.

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

---

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 39 of file analog.cpp.

                              {
        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 ;
}

---