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

Go to the documentation of this file.

/* @version $Revision: 1871 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-08-29 09:45:12 +0200 (Wed, 29 Aug 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 "mTypes.h"

#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;
  //TH1I * dacPed = new TH1I("ped","",391643019,820234834,1211877853);
  //TH1I * dacPed = new TH1I("ped","",100000,820234834,1211877853); //1109
  // TH1I * dacPed = new TH1I("ped","",100000,1682148314,2170757373); // 1112
  // TH1I * dacPed = new TH1I("ped","",100000,2437670539,2850616067); // 1114
  // TH1I * dacPed = new TH1I("ped","",100000,3098777110,3548395920); // 1116
   //TH1I * dacPed = new TH1I("ped","",100000,12721389,4293974986); // 1118
   TH1I * dacPed = new TH1I("ped","",100000,0,10000000000); // 1116




  UInt_t min = 0xffffffff;
  UInt_t max = 0;

  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 < nbEntries ; evtNum++)
    {
      tree->GetEntry(evtNum);
      nbEvt++;

      int nbChannel = evt->GetNchannel();
      for(int i=0;i<nbChannel  ;i++)
      {
        nbHit++;
        channel = (MTChannel*)evt->GetChannels()->UncheckedAt(i);
                    UInt_t absTime = (channel->GetBcIdAbs() ); 
        if ( absTime > max ) max = absTime;
        if ( absTime < min ) min = absTime;
        dacPed->Fill(absTime,1);
      }
    }
  }
  cout << "max " << max << endl;
  cout << "min " << min << endl;
 output->cd();
 dacPed->Write("");

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

---