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

Go to the documentation of this file.

/* @version $Revision: 2001 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-12-18 13:39:25 +0100 (Tue, 18 Dec 2012) $ */
#include <map>
#include "TKey.h"
#include "TList.h"
#include <TROOT.h>
#include <TRint.h>
#include <TKey.h>
#include <TList.h>
#include <TTree.h>
#include <TFile.h>
#include <stdlib.h>
#include <TSystem.h>
#include <iostream>
#include <sstream>
#include "Log.hh"

#include "MicroException.hh"
#include <string>

#include "mysql/Mysql.hh"


#include "root/MTRun.hh"
#include "root/MTTrack.hh"
#include "root/MTChamber.hh"
#include "root/MTDetector.hh"
#include "root/MTBoard.hh"
#include "root/MTChip.hh"
#include "root/MTChannel.hh"
#include "root/MTEvent.hh"
#include "root/MTDif.hh"
#include "root/MTHardroc1Chip.hh"
#include "root/MTHardroc2Chip.hh"
#include "root/MTDiracChip.hh"
#include "root/MTMicrorocChip.hh"


using namespace std;


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

/*
  Mysql mysql;
  MYSQL_ROW row;

  string sql = "SELECT * FROM RUN";

  MYSQL_RES* result = mysql.executeSelect(sql);

  int num_fields = mysql_num_fields(result);

  while ((row = mysql_fetch_row(result)))
  {
      for( int i = 0; i < num_fields; i++)
      {
          printf("%s ", row[i] ? row[i] : "NULL");
      }
      printf("\n");
  }

  int id = -1;
  try {
    id  = mysql.getReconstructionSoftwareId("114");
  }
  catch (MicroException e) {}

  cout << "Id[" << id << "]" <<endl ;
  id = mysql.addReconstructionSoftware("$Revision: 2001 $", "comments");
  cout << "Id[" << id << "]" <<endl ;
  try {
    id  = mysql.getReconstructionSoftwareId("114");
  }
  catch (MicroException e) {}
  cout << "Id[" << id << "]" <<endl ;

*/

   if ( argc !=2  ) {
   FILE_LOG(logERROR)  << "usage: example rootFile " << endl;
   exit(1);
  }

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

  TFile f(rootName.c_str());
  f.ls();



  TIter nextkey(f.GetListOfKeys());
  TKey *key;
  while (key = (TKey*)nextkey()) 
  {
    TTree *tree = (TTree*)key->ReadObj();                
    MTRun *run=NULL;
    try 
    {
      run=dynamic_cast<MTRun *>(tree->GetUserInfo()->First());
      cout << "---------------------------- NEW TTREE ---------------------------"<< endl;
      nbEvt =0;
      nbHit = 0;
      //cout << "RunId[" << run->GetRunId() << "]" << endl;

      bool showConfig = false;
      if ( run != NULL && showConfig) 
      {
        run->Info(); 
        MTDetector* det = run->GetDetector();
        // loop over chamber config
          
        const std::map<UInt_t,MTChamber*>& chambers = det->GetChambers();
        for( std::map<UInt_t,MTChamber*>::const_iterator ii=chambers.begin(); ii!=chambers.end(); ++ii)
        {
          const MTChamber& chamber= *((*ii).second);
          cout << "chamber id: " << chamber.GetId() << endl;
          cout << "chamber type: " << chamber.GetType() << endl;

          // loop over dif config
          //const std::map<UInt_t,MTDif*>& difs = chamber.GetDifs();
          //for( std::map<UInt_t,MTDif*>::const_iterator ii=difs.begin(); ii!=difs.end(); ++ii)
          //{
            //const MTDif& dif= *((*ii).second);
            //dif.Info();
          //}
          // loop over board config
          const std::map<UInt_t,MTBoard*>& boards = chamber.GetBoards();
          for( std::map<UInt_t,MTBoard*>::const_iterator ii=boards.begin(); ii!=boards.end(); ++ii)
          {
            const MTBoard& board= *((*ii).second);
            board.Info();

            cout << endl << endl << "\t\t\t**** Chips configuration ****" << endl ;
cout << "DEBUG board.GetChips().size()[" << board.GetChips().size() << "]" << endl; 
            map<UInt_t,MTChip*> chips = board.GetChips();
            for( map<UInt_t,MTChip*>::const_iterator iiChip=chips.begin(); iiChip!=chips.end(); ++iiChip)
            {
              const MTChip& chip  = *((*iiChip).second);
              cout << "chip[" << chip.GetId() << "] bypass[" << chip.GetBypassed() << "] " << endl;
              if ( chip.IsConfigured() )
              {
              cout << endl << "   -- New Chip configuration:" <<  chip.GetSoftId().toString() << " softId[" << chip.GetSoftId().GetValue() <<"], serial Number[" << chip.GetSerialNumber() << "]"<< endl <<endl  <<endl ;
              try {
                const MTDiracChip &diracChip = dynamic_cast<const MTDiracChip &> (chip);
                cout << "   ** Dirac Chip config **" << endl;
                cout << "      threshold B High: " << diracChip.GetBHiThreshold()  << endl;
                cout << "      threshold B Mid : " << diracChip.GetBMidThreshold()  << endl;
                cout << "      threshold B Low : " << diracChip.GetBLoThreshold()  << endl;

                cout << "      threshold A High :" << diracChip.GetAHiThreshold()  << endl;
                cout << "      threshold A Mid : " << diracChip.GetAMidThreshold()  << endl;
                cout << "      threshold A Low : " << diracChip.GetALoThreshold()  << endl;
                cout << "   ** END Dirac Chip config **" << endl;
              }
              catch (...) {}

              try{
                const MTHardroc1Chip &hr1Chip = dynamic_cast<const MTHardroc1Chip &> (chip);
                cout << "   ** Hardroc1 Chip config **" << endl;
                cout << "      threshold dac 0: " << hr1Chip.GetThresholdDac_0()  << endl;
                cout << "      threshold dac 1: " << hr1Chip.GetThresholdDac_1()  << endl;
                cout <<    "   gain   :         " << hr1Chip.GetGain() << endl;
                cout << "   ** END Hardroc1 Chip config **" << endl;
              }
              catch (...) {}

              try{
                const MTHardroc2Chip &hr2Chip = dynamic_cast<const MTHardroc2Chip &> (chip);
                cout << "   ** Hardroc2 Chip config **" << endl;
                cout << "      threshold dac 0: " << hr2Chip.GetThresholdDac_0()  << endl;
                cout << "      threshold dac 1: " << hr2Chip.GetThresholdDac_1()  << endl;
                cout << "      threshold dac 2: " << hr2Chip.GetThresholdDac_2()  << endl;
                cout <<    "   gain   :         " << hr2Chip.GetGain() << endl;
                cout << "   ** END Hardroc2 Chip config **" << endl;
              }
              catch (...) {}

              try{
                const MTMicrorocChip &micro2Chip = dynamic_cast<const MTMicrorocChip &> (chip);
                cout << "   ** MicrorocChip Chip config **" << endl;
                cout << dec << "      Serial number : " << micro2Chip.GetSerialNumber()  << endl;
                cout << dec << "      low shaper : " << micro2Chip.GetShaper_lg()  << endl;
                cout << dec << "      high shaper : " << micro2Chip.GetShaper_hg()  << endl;
                cout << dec << "      threshold dac 0: " << micro2Chip.GetThresholdDac_0()  << endl;
                cout << dec << "      threshold dac 1: " << micro2Chip.GetThresholdDac_1()  << endl;
                cout << dec << "      threshold dac 2: " << micro2Chip.GetThresholdDac_2()  << endl;
                cout << "   ** END MicrorocChip Chip config **" << endl;
              }
              catch (...) {}

              cout << endl << endl << "* channels config **" << endl << endl;
              cout <<"      Enable channel:" ;
              for ( UInt_t num = 0 ; num <=63 ; num++)
              {
                cout << "[" << dec <<num << ":"<< chip.GetEnableChannel(num) << "]";
              }
              cout << endl <<endl << endl <<  "    Stimulate channel:";
              for ( UInt_t num = 0 ; num <=63 ; num++)
              {
                cout << "[" << dec <<num << ":"<< chip.GetStimulatedChannel(num) << "]";
              }

              cout << endl <<endl << endl <<  "    Mask ( Microroc Only  ) channel:";
              for ( UInt_t num = 0 ; num <=63 ; num++)
              {
                for ( UInt_t index = 0 ; index < 3 ; index++ )
                {
                  cout << "[" << dec <<num << ":MASK" << index << "-> "<< chip.GetMaskChannel(num,index) << "]";
                }
              }


              cout << endl <<endl << endl  << "    pedastal correction channel:";
              for ( UInt_t num = 0 ; num <=63 ; num++)
              {
                cout << "[" <<dec << num << ":"<< chip.GetPedestalChannel(num) << "]";
              }
              cout << endl << "   -- End Chip configuration:" << endl;
              } // end if chip is configured

              }  // end loop for chip
                cout  <<endl << "\t\t\t**** END Chips configuration **** "<< endl  << endl << endl;
          }
        }
      }
    }
    catch (MicroException e) {}


  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 nbPad = 1300;
    //UInt_t my_asu;
    //UInt_t my_abstime;
    //cout<<evt->IsSquare(nbPad,my_asu,my_abstime)<<endl;

/*
  //    cout << "    event num in TTree:" << evtNum <<endl;
    cout << "Validity_flag: [0x" << hex << evt->GetValidFlag() << "]" << dec << evt->GetValidFlagString() << endl;

    bool spark = false;
    std::map< UInt_t, std::map<UInt_t, UInt_t>* >& nbHitPerAbsTime = evt->GetNbHitPerAbsTime();
    for (std::map< UInt_t, std::map<UInt_t, UInt_t>* >::iterator iter = nbHitPerAbsTime.begin(); iter != nbHitPerAbsTime.end(); iter++)
    {
    map<UInt_t, UInt_t> *byBoard = iter->second;
    for ( map<UInt_t,UInt_t>::iterator it = byBoard->begin(); it != byBoard->end(); it++)
    {
    UInt_t absTime = it->first;
    UInt_t nbHit   = it->second;
        cout << "Board[" << iter->first << "] Absolute time [" << absTime  << "]nbHit[" << nbHit <<"]"    << endl;
    if (nbHit > 1) { spark = true; }
    }
    }

*/

    //if (!evt->GetCrcIsCorrect())
    {
    cout << "    event id:" << evt->GetEventId() <<endl;

    cout << "    event Dif suynchro:" << evt->GetDifSynchro() <<endl;
    cout << "    event total Trigger counter:" << evt->GetGlobalTriggerCounter() <<endl;
    cout << "    event Dif Trigger counter:" << evt->GetDifTriggerCounter() <<endl;
    const map<unsigned int , float>& tempAsu1 = evt->GetTemperatureAsu1();
    for ( map<unsigned int , float>::const_iterator it = tempAsu1.begin() ; it != tempAsu1.end(); it++)
    {
      cout << "    temperature difId: " << it->first << " Asu1:" << it->second <<endl;
    }
    const map<unsigned int , float>& tempAsu2 = evt->GetTemperatureAsu2();
    for ( map<unsigned int , float>::const_iterator it = tempAsu2.begin() ; it != tempAsu2.end(); it++)
    {
      cout << "    temperature difId: " << it->first << " Asu2:" << it->second <<endl;
    }

    const map<unsigned int , float>& tempDif = evt->GetTemperatureDif();
    cout << "evt->GetTemperatureDif().size[" << evt->GetTemperatureDif().size() << "]" <<endl;
    for ( map<unsigned int , float>::const_iterator it = tempDif.begin() ; it != tempDif.end(); it++)
    {
      cout << "    temperature difId: " << it->first << " dif:" << it->second <<endl;
      cout << "DEBUG 1" << endl;
    }
cout << "evt->GetTemperatureAsu1().size[" << evt->GetTemperatureAsu1().size() << "]" <<endl;
cout << "evt->GetTemperatureAsu2().size[" << evt->GetTemperatureAsu2().size() << "]" <<endl;
    cout << "    temperature:" << evt->GetTemperature() <<endl;
    cout << "    pressure:"<< evt->GetPressure() <<endl;
    cout << "    timestamp:"<< evt->GetTimestamp() <<endl;
    cout << "    dif_synchro:"<< evt->GetDifSynchro() <<endl;
    cout << "    CRC is correct:"<< evt->GetCrcIsCorrect() <<endl;

    std::vector<UInt_t>& dtCands = evt->GetDtCandidate();
    for ( std::vector<UInt_t>::const_iterator it =dtCands.begin(); it!=dtCands.end(); it++)
    {
      cout << "dt canditate: [" << *it << "] " << endl; 
    }


    const std::map<UInt_t, std::map<UInt_t, std::map<UInt_t, Float_t>* >* >& meshs = evt->GetMeshVoltage();
    for ( std::map<UInt_t, std::map<UInt_t, std::map<UInt_t, Float_t>* >* >::const_iterator iterChamber = meshs.begin(); iterChamber != meshs.end(); iterChamber++)
    {
      std::map<UInt_t, std::map<UInt_t, Float_t>* >&difContainer = *(iterChamber)->second;
      for ( std::map<UInt_t, std::map<UInt_t, Float_t>* >::const_iterator iterDif = difContainer.begin(); iterDif != difContainer.end(); iterDif++)
      {
        std::map<UInt_t, Float_t> &boardContainer = *(iterDif)->second;
        for ( std::map<UInt_t, Float_t>::const_iterator iterVoltage = boardContainer.begin(); iterVoltage != boardContainer.end(); iterVoltage++)
        {
          cout << "    MeshVoltage for chamber [" << iterChamber->first << "]Dif = ["<< iterDif->first << "] board[" << iterVoltage->first << "] :" <<iterVoltage->second << "]" << endl;
        }
      }

    }


    const std::map<UInt_t, std::map<UInt_t, std::map<UInt_t, Float_t>* >* >& drifts = evt->GetDriftVoltage();
    for ( std::map<UInt_t, std::map<UInt_t, std::map<UInt_t, Float_t>* >* >::const_iterator iterChamber = drifts.begin(); iterChamber != drifts.end(); iterChamber++)
    {
    std::map<UInt_t, std::map<UInt_t, Float_t>* >&difContainer = *(iterChamber)->second;
    for ( std::map<UInt_t, std::map<UInt_t, Float_t>* >::const_iterator iterDif = difContainer.begin(); iterDif != difContainer.end(); iterDif++)
    {
    std::map<UInt_t, Float_t> &boardContainer = *(iterDif)->second;
    for ( std::map<UInt_t, Float_t>::const_iterator iterVoltage = boardContainer.begin(); iterVoltage != boardContainer.end(); iterVoltage++)
    {
    // mergedEvt->SetDriftMeshVoltage(iterChamber->first,iterDif->first,iterVoltage->first,iterVoltage->second);
    cout << "    DriftVoltage for chamber [" << iterChamber->first << "]Dif = ["<< iterDif->first << "] board[" << iterVoltage->first << "] :" <<iterVoltage->second << "]" << endl;

    }
    }

    }
    }
/*
      if ( run->GetProcessed() == true )
*/
      {
        cout << "    is square[" << evt->GetSquareEvt() << "]"  << endl;
        cout << "    is time info correct[" << evt->GetTimeInfo() << "]"  << endl;
        MTTrack *track = evt->GetTrack();
        if ( track != NULL )
        {
          cout << "    straight: " << track->GetStraight()
          << ", from Strip: " << track->GetFromStrip()
          << ", from Pad: " << track->GetFromPad()
          << ", single hit: " << track->GetSingleHit()
          << ", ax: " << track->GetAx()
          << ", bx: " << track->GetBx()
          << ", ay: " << track->GetAy()
          << ", by: " << track->GetBy()  << endl;
        }
      }
      cout << "    number of channel hit:" << evt->GetNchannel() <<endl;
      int nbChannel = evt->GetNchannel();
/*
      for(int i=0;i<nbChannel  ;i++)
      {
        nbHit++;
        channel = (MTChannel*)evt->GetChannels()->UncheckedAt(i);


        cout << "****** Channels Information "<< i+1 << "/" << nbChannel << " ******" <<  endl;
        cout << "      Hit num:            " << i << endl;
        cout << "      chamber id:         " << channel->GetChamberId() <<endl;
        cout << "      dif id:             " << channel->GetDifId() <<endl;
        cout << "      boardId id:         " << channel->GetBoardId() <<endl;
        cout << "      chip id(format dec):" << dec << channel->GetChipId() <<endl;
        cout << "      hard id    :        " << dec << channel->GetHardId() <<endl;
        cout << "      row inside chip :   " << channel->GetRowInChip() <<endl;
        cout << "      col inside chip:    " << channel->GetColInChip() <<endl ;
        cout << "      soft id:            " << channel->GetSoftId().GetValue() << endl<< endl;
        cout << "      soft id:            " << channel->GetSoftId().GetValue() << endl<< endl;
        cout << "      asu row:            " << channel->GetSoftId().GetBoardRow() << endl<< endl;
        cout << "      asu col:            " << channel->GetSoftId().GetBoardCol() << endl<< endl;
        cout << "      channel row inside chamber: " << channel->GetRowInChamber() <<endl;
        cout << "      channel col inside chamber: " << channel->GetColInChamber() <<endl ;
        cout << "      analog value:       " << channel->GetAnalogValue() <<endl;
        cout << "      digital value:      " <<  channel->GetDigitalValue() <<endl;
        cout << "      x:                  " <<  channel->GetX() <<endl;
        cout << "      y:                  " << channel->GetY() <<endl;
        cout << "      z:                  " << channel->GetZ() <<endl;
        cout << "      order:              " << channel->GetOrder() <<endl;
        cout << "      timestamp:          " <<  channel->GetTimestamp() <<endl;
        cout << "      bcId abs:           " << channel->GetBcIdAbs() <<endl;
        cout << "      bcId dif:           " << channel->GetBcIdDif() <<endl;
        cout << "      bcId hit:           " << channel->GetBcIdHit() <<endl;
        UInt_t absTime =  channel->GetBcIdAbs() + (channel->GetBcIdDif() - channel->GetBcIdHit() );
        cout << "DEBUG      absolute time:bcid_abs - (bcid_dif - bcid_hit)       " << absTime <<endl;
        cout << "DEBUG      dt (bcid_dif - bcid_hit)       " << channel->GetBcIdDif() - channel->GetBcIdHit() << endl;
        cout << "      nbHit for this absolute time:" << evt->GetNbHitPerAbsTime(channel->GetSoftId().GetBoardId(),absTime)<<endl;
        cout << "      bcId hit Fine:      " << channel->GetBcIdHitFine() <<endl;
        cout << "      memory order :      " << channel->GetMemoryOrder() <<endl;
      } // END LOOP ON CHANNELHIT
*/

    } // END LOOP ON EVENT
    cout << "End of TTree[" << tree->GetName() << "]" <<endl;
    cout << endl << "** Nombre d Event["<< nbEvt <<"] ** "<< endl;
    cout << endl << "** Nombre de Hits["<< nbHit <<"] ** "<< endl;

  } // end loop on TKEy
return 0 ;
}

---