/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/analyse/Damien/evt_cadre.C

Go to the documentation of this file.

{
  gROOT->Reset();
  gROOT->SetStyle("Bold");
  gStyle->SetPalette(1);
  gSystem->Load("~/Micromegas/Analyse/trunk/lib/libMicro.so");
  
  
  //////////////////////
  //User's parameters
  
  //which event do you want to plot (0 = no plot, a number = the event you want)
  //Int_t plot = 0;
  Int_t plot = 4555;
  
  //Do you want some info ? (0 = no, 1 = yes a little, 2 = yes a lot, 3 = yes all for debug purpose)
  Int_t blabla = 2;
  
  //what are the id's of our asu's ? (assuming they go from first_asu to first_asu+5)
  Int_t asu[6];
  for(int it=0;it<6;it++) asu[it] = it+13;
  
  //square events
  Int_t nb_square = 20;
  Int_t list_square[nb_square] = {105,375,558,1090,1356,1600,2212,2563,2970,2985,2985,3481,3481,3633,3736,4481,4562,4606,5329,5662};
  
  //threshold
  Int_t threshold = 30;

  //input data
  string file = "05082011_1055_MGT.root";  
  //////////////////////
  
  
  //our data file
  TFile *f = new TFile(file.c_str());
  if(blabla==1||blabla==2||blabla==3){
    cout<<endl;
    cout<<"reading File : "<<file<<endl;
  }
  
  //key iterator
  TIter nextkey(f.GetListOfKeys());
  TKey *key;
  
  //counters
  UInt_t counter = 0;
  UInt_t channels_asu[6][128];
  UInt_t nb_border[6][128];
  UInt_t nb_tot;
  
  Int_t xpos;
  Int_t ypos;

//Là y'a un truc bizarre:
//si la déclaration de hxy est dans le if -> ca rame, si elle est hors du if -> ca tourne bien
//je comprends pas du tout pourquoi ??? 
  //plotting
  //TH2I *hxy = new TH2I("hxy","Hit position distribution;y (cm);x (cm)",96,0,96,96,0,96);
  if(plot != 0){
    TH2I *hxy = new TH2I("hxy","Hit position distribution;y (cm);x (cm)",96,0,96,96,0,96);
    TCanvas *c0 = new TCanvas("c0", "Prototype m2", 100, 50, 500, 500);
  }
  
  if(blabla==1||blabla==2||blabla==3) cout<<endl;
  
  
  //loop on all the events of the key
  while (key = (TKey*)nextkey()) {
    
    counter++;
    if(blabla==1||blabla==2||blabla==3) cout<<"key number = "<<counter<<endl;
    
    TTree *tree = (TTree*)key->ReadObj();                
    MTRun* run = (MTRun*)tree->GetUserInfo()->FindObject("MTRun");
    
    Int_t nEvent = tree.GetEntries();
    if(blabla==1||blabla==2||blabla==3) cout <<"Nevent = "<<nEvent<<endl;
    
    //Int_t azu = run->GetDetector()->GetChambers()->GetBoards(); 
    
    MTEvent *evt = new MTEvent();
    
    TBranch *branch= tree->GetBranch("MTEvent");
    branch->SetAddress(&evt);
    
    UInt_t list_event[nEvent];
    
    //at each new key, reset the event border counters
    for (int bd=0;bd<=5;bd++) for (int ti=0;ti<128;ti++) nb_border[bd][ti] = 0; 
    nb_tot=0;
    
    
    //loop on all the events of the key
    for (int i=0;i<nEvent-1;i++){

      //test if this event is square (if it is square, my algorythm would see it as a false border event because there are hits on the borders...)
      bool I_am_a_square_event = 0;
      for (int it=0;it<nb_square;it++) if (i+1==list_square[it]) I_am_a_square_event = 1;
      if (I_am_a_square_event==1) continue;
      
      //at each new event, reset the channels counter of each board_id
      for (int bd=0;bd<=5;bd++) for (int ti=0;ti<128;ti++) channels_asu[bd][ti] = 0; 
      
      tree.GetEntry(i);
      
      UInt_t nchannel = evt->GetNchannel();
      UInt_t board_id;
      UInt_t mem_order;
      
      if(blabla==3) cout<<"Event "<<i+1<<" / "<<nEvent<<"; nb of channels=  "<<nchannel<<endl;

      
      //loop on all the channels of the event
      for (int j=0;j<nchannel;j++){
        
        //Reset board_id and memory order
        board_id = 0;
        mem_order = 0;
        
        //create this channel
        MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
        
        //get the board_id of this channel
        board_id = channel->GetBoardId();

        //get the memory order of this hit
        mem_order = channel->GetMemoryOrder();
        
        //if the event is on the borders of the asu upgrade the hit counter of the good board at the good time
        if (board_id==asu[0]) if (xpos == 000000 || xpos == 310000 || ypos == 000000 || ypos == 470000) channels_asu[0][mem_order]++;
        if (board_id==asu[1]) if (xpos == 000000 || xpos == 310000 || ypos == 480000 || ypos == 960000) channels_asu[1][mem_order]++;
        if (board_id==asu[2]) if (xpos == 320000 || xpos == 640000 || ypos == 000000 || ypos == 470000) channels_asu[2][mem_order]++;
        if (board_id==asu[3]) if (xpos == 320000 || xpos == 640000 || ypos == 480000 || ypos == 960000) channels_asu[3][mem_order]++;
        if (board_id==asu[4]) if (xpos == 650000 || xpos == 960000 || ypos == 000000 || ypos == 470000) channels_asu[4][mem_order]++;
        if (board_id==asu[5]) if (xpos == 650000 || xpos == 960000 || ypos == 480000 || ypos == 960000) channels_asu[5][mem_order]++;
        
        //get the position of this hit/channel
        xpos = channel->GetX()/10000;
        ypos = channel->GetY()/10000;
        
        if(blabla==3) cout<<"Channel : "<<channel<<"; hit on asu "<<board_id<<"; at memory order "<<mem_order<<", x= "<<xpos<<", y= "<<ypos<<endl;
        
        //if we want to plot this event, fill the histogramm
        if (plot==i+1) hxy->Fill(ypos,xpos); 
      }
      
      //On each board, at each time, identify the events which are on the border
      for (int bd=0;bd<=5;bd++) {
         for (int ti=0;ti<128;ti++) {
           if (channels_asu[bd][ti]>= threshold) {

             if(blabla==3) cout <<bd<<" "<<ti;
             if(blabla==2||blabla==3) cout <<" -> event: "<<i+1<<" = border in asu no "<<asu[bd]<<" at memory order "<<ti<<" with Nhits = "<<channels_asu[bd][ti]<<endl;
             
             //update the event border counter of the good board at the good time
             nb_border[bd][ti]++;
             //update the total event border counter of this key
             nb_tot++;

            //update the event list
            list_event[nb_tot-1] = i+1;
           }
         }
      }
      //if we want to plot this event, let's do it now
      if (plot==i+1){
        hxy->SetTitle("Hits position");
        hxy.GetXaxis()->SetTitle("X (cm)");
        hxy.GetYaxis()->SetTitle("Y (cm)");
        hxy->Draw("zcol");
      }  
    }

    //Now the Prints: at each key print nb of evts border
    if(blabla==2||blabla==3) cout<<endl;
    if(blabla==1||blabla==2||blabla==3) {
      //at each key print nb of evts border
      cout<<"In key "<<counter<<" Nb tot evt border = "<<nb_tot<<endl;
      //at each key print print the liste of border events
      cout<<"List of border events: "<<endl;
      for(int it=0;it<nb_tot;it++) cout<<" "<<list_event[it];
      cout<<endl<<endl;
    }
    if(blabla==2||blabla==3) {
      cout<<"Hereafter I show the nb of event border in each board/asu at each time:"<<endl<<endl<<"Board/asu i"<<endl<<"memory order:";
      for (int ti=0;ti<128;ti++) cout<<" "<<ti;
      cout<<endl;
      for (int bd=0;bd<=5;bd++) {
        cout<<endl;
        cout<<"Board/asu "<<asu[bd]<<endl;
        for (int ti=0;ti<128;ti++) {
          cout<<" "<<nb_border[bd][ti];
          nb_border[bd][ti] = 0; 
          
        }
        cout<<endl;
      }
      cout<<endl;
    }
  }
}

---