/Users/jacquem/NetBeansProjects/MicromegasFramework/trunk/src/analyse/Gassiplex/Shower_profile/Event_viewer.C

Go to the documentation of this file.

void View_event(int i){
  
  gSystem->Load("/home1/chefdevi/micromegasFrameWork/trunk/libMicro.so");
  



  
  //TString path = "/lapp_data/LC/Detecteurs/MicroMegas/data/TB2009/May/Rootfiles/dev/";
  TString path = "/home1/chefdevi/micromegasFrameWork/trunk/src/analyse/Gassiplex/Shower_profile/";
  TString fileName= "electron_2gev_6absorber.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;
  int NC = 672;


 

  /*Distributions of the number of hits in the xy-plane*/

  TH2I * hxy1 = new TH2I("hxy1","",32,0,32,32,0,32);
  TH2I * hxy2 = new TH2I("hxy2","",32,0,32,32,0,32);
  TH2I * hxy3 = new TH2I("hxy3","",32,0,32,32,0,32);
  TH2I * hxy4 = new TH2I("hxy4","",32,0,32,32,0,32);




  /*Define variables*/

  int xpos1,ypos1,adc1,sum_adc1;
  int xpos2,ypos2,adc2,sum_adc2;
  int xpos3,ypos3,adc3,sum_adc3;
  int xpos4,ypos4,adc4,sum_adc4;

  vector <int> vxpos1,vypos1,vadc1;
  vector <int> vxpos2,vypos2,vadc2;
  vector <int> vxpos3,vypos3,vadc3;
  vector <int> vxpos4,vypos4,vadc4;

  float avgx[3],avgy[3];
  float posz1 = 6;
  float posz2 = 13;
  float posz3 = 19;
  float posz4 = 98;
  float avgz[3] = {posz1,posz2,posz3};
  float avgx1=0,avgx2=0,avgx3=0,avgx4=0;
  float avgy1=0,avgy2=0,avgy3=0,avgy4=0;

  int threshold = 35;
  float nhit1=0,nhit2=0,nhit3=0,nhit4=0;
  int nclean = 0;
  bool clean = 0;

  int event = 0;
  int crap_event = 0;

  float x0 = 6.89;
  float y0 = 16.48;
  float deltax = 0,deltay = 0,radial = 0;





  
  cout<<"Event # "<<i<<endl;
  
  tree->GetEntry(i);
  NC = evt->GetNchannel();
  if (NC != 672){cout<<"Crap event"<<endl;}
  
  clean=0;
  nhit1=0,nhit2=0,nhit3=0,nhit4=0;
  sum_adc1=0;sum_adc2=0;sum_adc3=0;sum_adc4=0;
  
  

  /*Chamber beta 2.1 D*/
  
  for(int j=0;j<96;j++){
    
    MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
    xpos1 = channel->GetX();
    ypos1 = channel->GetY();
    adc1  = channel->GetValue();
    
    if (adc1>=threshold){
      
      cout<<" x1="<<xpos1<<"  y1="<<ypos1<<"  adc1="<<adc1<<endl;
      nhit1++;
      hxy1->SetBinContent(xpos1+1,ypos1+1,adc1);}}
  

    /*Chamber beta 2.1 E*/

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

      MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
      xpos2 = channel->GetX();
      ypos2 = channel->GetY();
      adc2  = channel->GetValue();

      if (adc2>=threshold){

        cout<<" x2="<<xpos2<<"  y2="<<ypos2<<"  adc2="<<adc2<<endl;
        nhit2++;
        hxy2->SetBinContent(xpos2+1,ypos2+1,adc2);}}


    /*Chamber beta 2.1 A*/

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

      MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
      xpos3 = channel->GetX();
      ypos3 = channel->GetY();
      adc3  = channel->GetValue();

      if (adc3>=threshold){

        cout<<" x3="<<xpos3<<"  y3="<<ypos3<<"  adc3="<<adc3<<endl;
        nhit3++;
        hxy3->SetBinContent(xpos3+1,ypos3+1,adc3);}}


  
    /*chamber beta 2.4*/

    for(int j=287;j<672;j++){
      MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
      xpos4 = channel->GetX();
      ypos4 = channel->GetY();
      adc4  = channel->GetValue();

      if (adc4>=threshold){
        
        cout<<" x4="<<xpos4<<"  y4="<<ypos4<<"  adc4="<<adc4<<endl;
        nhit4++;
        hxy4->SetBinContent(xpos4+1,ypos4+1,adc4);}}


    cout<<nhit1<<" "<<nhit2<<" "<<nhit3<<" "<<nhit4<<endl;
    

    
    /*Plot distributions*/
    
    
    TCanvas * myc3 = new TCanvas("myc3");
    myc3->Divide(2,2);
    myc3->cd(1);
    hxy1->Draw("zcol");
    myc3->cd(2);
    hxy2->Draw("zcol");
    myc3->cd(3);
    hxy3->Draw("zcol");
    myc3->cd(4);
    hxy4->Draw("zcol");

  }

---