/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/Microroc/Calibration/Calib_chip_microroc.cpp

Go to the documentation of this file.

/* @version $Revision: 1585 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-03-13 11:03:35 +0100 (Tue, 13 Mar 2012) $ */
#include <map>
#include "TKey.h"
#include <TROOT.h>
#include <TRint.h>
#include <TKey.h>
#include <TTree.h>
#include <TFile.h>
#include <TSystem.h>
#include <iostream>
#include <TF1.h>
#include <TGraph.h>
#include <sstream>
#include "Log.hh"

#include "MicroException.hh"

#include "root/MTRun.hh"
#include "root/MTDetector.hh"
#include "root/MTChamber.hh"
#include "root/MTChip.hh"
#include "root/MTEvent.hh"
#include "root/MTChannel.hh"
#include "root/MTHardroc2Chip.hh"

using namespace std;


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

  /*****************************************************************/
  /* ARGUMENTS TO THE PROGRAM **************************************/
  /* 1 - DIRECTORY OF THE ROOT FILE PRODUCED BY THE RECONSTRUCTION */
  /* 2 - NAME OF THE ROOT FILE *************************************/
  /* 3 - THRESHOLD THAT IS MOVED DURING THE CALIBRATION ************/
  /*    0 -> DAC_0 (LOW THRESHOLD) *********************************/
  /*    1 -> DAC_1 (MEDIUM THRESHOLD) ******************************/
  /*    2 -> DAC_2 (HIGH THRESHOLD) ********************************/
  /*****************************************************************/

  TString root_directory = argv[1];
  TString root_filename = argv[2];
  int dac = atoi(argv[3]);

  bool display = 0;

  TString root_file = root_directory + root_filename;
  TString scurve_file = root_directory + "scurve_" + root_filename;

  cout<<"Wish to write "<<scurve_file<<endl;

  int nbChannel = 0;
  float nbEvt = 0;
  int threshold = 0;
  int chipid = 0;
  int hardid = 0;
  int value = 0;
  int content = 0;
 
  TH1I * hvalue = new TH1I("hvalue","",5,0,5);
  TH1I * hnpulse_threshold = new TH1I("hnpulse_threshold","",1024,0,1024);

  TH1I * hscurve[64];
  TString name,title;
 
  for (int j=0;j<64;j++){

    name  = Form("hscurve_%i",j);
    title = Form("SCurve from channel %i",j);
    hscurve[j] = new TH1I(name,title,1024,0,1024);}




  
  
  TFile f(root_file);
  TIter nextkey(f.GetListOfKeys());
  TKey *key;
  

  while (key = (TKey*)nextkey()) {

    TTree *tree = (TTree*)key->ReadObj();                
    MTRun* run = (MTRun*)tree->GetUserInfo()->FindObject("MTRun");

                const MTChip&  chip = run->GetOneChip();
                try{
                        const MTHardroc2Chip &hr2Chip = dynamic_cast<const MTHardroc2Chip &> (chip);
                        cout << "threshold dac 0:[" << hr2Chip.GetThresholdDac_0()  << endl;
                        cout << "threshold dac 1:[" << hr2Chip.GetThresholdDac_1()  << endl;
                        cout << "threshold dac 2:[" << hr2Chip.GetThresholdDac_2()  << endl;
                        if (dac==0){threshold = hr2Chip.GetThresholdDac_0();}
                        if (dac==1){threshold = hr2Chip.GetThresholdDac_1();}
                        if (dac==2){threshold = hr2Chip.GetThresholdDac_2();}
                }
                catch (...) {}


    MTEvent *evt =  new MTEvent();
    TBranch *branch= tree->GetBranch("MTEvent");
    branch->SetAddress(&evt);
    
    MTChannel* channel = NULL;
    
    nbEvt = tree->GetEntries();
    cout<<"Number of events in the tree = "<<nbEvt<<endl;

    hnpulse_threshold->SetBinContent(threshold+1,nbEvt);

    for ( int evtNum = 0; evtNum < nbEvt ; evtNum++){
      
      if (evtNum%1000 == 0){cout<<evtNum/nbEvt*100.<<" %"<<"\r"<<flush;}

      tree->GetEntry(evtNum);
      nbChannel = evt->GetNchannel();

      for(int i=0;i<nbChannel ;i++){    
      
        channel = (MTChannel*)evt->GetChannels()->UncheckedAt(i);
        chipid = channel->GetChipId();
        hardid = channel->GetHardId();

        value = channel->GetDigitalValue();
        hvalue->Fill(value);

        if (value==(dac+1)){

          hscurve[hardid]->Fill(threshold);

            if (display){
              
              cout << "---- digital value:" <<  value <<endl;
              cout << "---- hard id:" << hardid <<endl;
              cout << "---- x:"<< channel->GetX() <<endl;
              cout << "---- y:" << channel->GetY() <<endl;
              cout << "---- chip id:" << chipid <<endl;}}}}}

    

  cout<<endl;
  
 





  TFile * tf = new TFile(scurve_file,"RECREATE");
  hvalue->Write();
  hnpulse_threshold->Write();
  for (int j=0;j<64;j++){
    hscurve[j]->Write();}
  tf->Close();











    return 0 ;}

---