/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/TB2010/Testbeam_june_10/Hardroc2/Beam_profile/Beam_profile_hardroc2.C

Go to the documentation of this file.

{
  gROOT->Reset();
   gStyle->SetPalette(1);
  gSystem->Load("/home/micromegas/Framework/libMicro.so");






  TString filename = "/lapp_data/LC/Detecteurs/MicroMegas/data/TB2010/SPS_H4_june_2010/Root_files/acq_HR2_14062010_1844_SLAB2_1.root"; 


  Long64_t t,t1,t2,t3,dt;
  int nchannel = 0;
  int hardid = 0;
  int softid = 0;
  int chipid = 0;
  int nhit = 0;
  int xpos = 0;  
  int ypos = 0;  
  int zpos = 0;  
  int adc = 0;  

  TH1I * hdt = new TH1I("hdt","",100,-10,90);
  TH2I * hxy_m2 = new TH2I("hxy_m2","",96,0,96,96,0,96);
  TH2I * hxy_m2_cut = new TH2I("hxy_m2_cut","",96,0,96,96,0,96);
  
  


  TFile *f = new TFile(filename);
  cout<<"reading File : "<<filename<<endl;

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

  while (key = (TKey*)nextkey()) {
    
    TTree *tree = (TTree*)key->ReadObj();                
    MTRun* run = (MTRun*)tree->GetUserInfo()->FindObject("MTRun");
    //run->Info();

  int nEvent = tree.GetEntries();
  cout <<"nEvent="<<nEvent<<endl;

  MTEvent *evt = new MTEvent();

  TBranch *branch= tree->GetBranch("MTEvent");
  branch->SetAddress(&evt);







  for (int i=0;i<nEvent;i++){  
  //for (int i=113;i<114;i++){  

    tree.GetEntry(i);
    nchannel = evt->GetNchannel();

    //cout<<"event number = "<<i<<"  Number of hits = "<<nchannel<<endl;

    MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(0);
    
    for (int j=0;j<nchannel;j++){
      
      MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
      
      xpos = channel->GetX();
      ypos = channel->GetY();

      hardid = channel->GetHardId();
      chipid = channel->GetChipId();
      
      hxy_m2->Fill(xpos,ypos);
      
      t1 = channel->GetBcId_Abs();
      t2 = channel->GetBcId_Dif();
      t3 = channel->GetBcId_Hit();
      dt = (t2-t3);
      t = (t1 - (t2-t3));
      
      hdt->Fill(dt);
      
      if (dt>=4 && dt<=10){hxy_m2_cut->Fill(xpos,ypos);}}}}


hdt->Draw();


}

---