/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/root/analog.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) $ */
#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 );



  TH1I * dacAll = new TH1I("all","",512,0,4096);
  TH1I * dacPed = new TH1I("ped","",512,0,4096);
  TH1I * dac0 = new TH1I("c0","",2048,0,4096);
  TH1I * dac1 = new TH1I("c1","",2048,0,4096);
  TH1I * dac2 = new TH1I("c2","",2048,0,4096);
  TH2I * histo2 = new TH2I("analog2","",3,1,4,4096,0,4096);

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

      UInt_t channel66 = 0;
      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);

//        if (channel->GetChipId() == 67 && channel->GetHardId() == 53 ) { channel66 =  channel->GetDigitalValue(); ;} 
//        if (channel->GetChipId() ==67 && channel->GetHardId() == 55 ) { channel68 =  channel->GetDigitalValue();}


                absTime = channel->GetBcId_Abs() - (channel->GetBcId_Dif() - channel->GetBcId_Hit() ); // see src/analyse/root/example.cpp
                UInt_t nbHitPerAbsTime = evt->GetNbHitPerAbsTime(channel->GetBoardId(), absTime) ; // see src
                

        if ( channel->GetAnalogValue() != 0 )//&& channel->GetDigitalValue() > 0 )// && channel->GetDigitalValue() < 3)
       {
          if ( channel->GetChipId() ==  67 )
          {
            if ( channel->GetHardId() == 54 )
            {   UInt_t adc = 4096-channel->GetAnalogValue();
            
                if (  channel->GetDigitalValue() == 0  && nb_ped < 300) 
                                {       dacPed->Fill(adc);
                                        //dacAll->Fill(adc);
                                        nb_ped++;
                                }

            
            if (nbHitPerAbsTime < 4)
                {
            
 //             if (  channel->GetDigitalValue() == 1  && adc < 1400) 
              if (  channel->GetDigitalValue() == 1  ) 
              {     
                dac0->Fill(adc);
                histo2->Fill(channel->GetDigitalValue(),adc);
                dacAll->Fill(adc);
              }
              //if (  channel->GetDigitalValue() == 2  && adc >= 1200 && adc <  1700) 
              if (  channel->GetDigitalValue() == 2 ) 
              {
                dac1->Fill(adc);
                histo2->Fill(channel->GetDigitalValue(),adc);
                dacAll->Fill(adc);
              }
              //if (  channel->GetDigitalValue() == 3  && adc  >=1600)
              if (  channel->GetDigitalValue() == 3 )
              {
                /* dac2->Fill(adc);
                 histo2->Fill(channel->GetDigitalValue(),adc);*/
                // dacAll->Fill(adc);
              }
            }
          }
        }
      }} //break;
    }
  }
   TCanvas * c0 = new TCanvas("c0","Prototype m2");//,0,0,800,600);
//   c0->Divide(3,2);
  
   //c0->cd(1);
   //dac0->Draw();
   //c0->cd(2);
   //dac1->Draw();
   //c0->cd(3);
   //dac2->Draw();
   //c0->cd(4);
   dacPed->Draw("");
   dacPed->Fit("gaus","R", "SAME", 1140, 1200);
//   dacPed->Fit("gaus");
   dacPed->GetFunction("gaus")->SetLineColor(kBlue+1);
   dacPed->GetYaxis()->SetRangeUser(0,210);
   dacPed->GetXaxis()->SetTitle("ADC count");
   dacPed->GetYaxis()->SetTitle("N");
   dacAll->Draw("sames");
   dacAll->Fit("landau","R", "SAME", 1100, 1800);
   dacAll->GetFunction("landau")->SetLineColor(kOrange+5);
        


        //dacAll->GetYaxis()->SetRangeUser(0,300);
        
        c0->Update();   // Create the TGRAPH. IF not present, following lines

        TPaveStats *p1 = (TPaveStats*) dacPed->GetListOfFunctions()->FindObject("stats");
//      dacPed->GetListOfFunctions()->Print();
//      p1->Print();
        dacPed->GetListOfFunctions()->Remove(p1);
        dacPed->SetStats(0);
        p1->SetTextColor(kBlue+1);
        p1->SetX1NDC(0.75);
        p1->SetX2NDC(0.9);
        p1->SetY1NDC(0.8);
        p1->SetY2NDC(0.9);
        p1->Draw();


        // NEED SAMES (!=SAME) in Drawin histo !!!!!
//      dacAll->GetListOfFunctions()->Print();
        TPaveStats *p2 = (TPaveStats*) dacAll->GetListOfFunctions()->FindObject("stats");
//      p2->Print();
        dacAll->GetListOfFunctions()->Remove(p2);
        dacAll->SetStats(0);
        p2->SetTextColor(kOrange+5);
        p2->SetX1NDC(0.75);
        p2->SetX2NDC(0.9);
        p2->SetY1NDC(0.7);
        p2->SetY2NDC(0.8);
        p2->Draw();
        

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

---