/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/root/Test_time.C

Go to the documentation of this file.

{
  gROOT->Reset();
  
  gSystem->Load("libMicro.so");
  


  /*Total number of events is 62046*/
  /*Number of event with Nchannel != 672 is 17888*/
  /*Number of good event (Nch=672) with one single hit (ADC>35) per chamber is 3570*/



  bool display = 0;

  
  TString path = "/lapp_data/LC/Detecteurs/MicroMegas/data/TB2009/May/Rootfiles/dev/temp/";
  //TString path = "/home1/chefdevi/micromegasFrameWork/trunk/src/analyse/Gassiplex/Beam_profile/";
  //TString fileName= "beam_profile_alice_pm.root";
  TString fileName= "bck_060609_064948_1.cfg.root";
  
  TFile *f = new TFile(path+fileName);
  cout<<"reading File : "<<path<<fileName<<endl;
  TTree *tree;
  tree  = (TTree*) f.Get("t1");cout<<"tree built\n";
  int nEvent = tree.GetEntries();
  MTEvent *evt = new MTEvent();
  TBranch *branch= tree->GetBranch("MTEvent");
  branch->SetAddress(&evt);
  cout <<"nEvent="<<nEvent<<endl;


  TH1F * histo = new TH1F("histo","",1e7,0,1e7);

  //for( int i=0;i<nEvent;i++){
    
    if (display){cout<<"Event # "<<i<<endl;}

    tree.GetEntry(0);
    NC = evt->GetNchannel();

//    Int_t nchannel = evt->GetChannels()->GetEntries();
    Int_t nchannel = evt->GetNchannel();

    for (int j=0;j<nchannel;j++){

      MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);

       evt->SortChannel();
       UInt_t t1 = channel->GetBcId_Abs();
       UInt_t t2 = channel->GetBcId_Dif();
       UInt_t t3 = channel->GetBcId_Hit();
       UInt_t value = channel->GetValue();
        
       cout << "t3:[" << t3 << "]" ;
       cout << "value:[" << value << "]" ;
       cout << "softId:[" << channel->GetSoftId() << "]" << endl;
//       histo->Fill(t3);
    }



    //histo->Draw();


}

---