Micromegas LAPP reconstruction

00001 {gROOT->Reset();
00002   
00003   gSystem->Load("/home1/chefdevi/micromegasFrameWork/trunk/libMicro.so");
00004   
00005 
00006   bool display = 0;
00007 
00008   TString fileName= "electron_2gev_4absorber.root";
00009 
00010   
00011   TString path = "/lapp_data/LC/Detecteurs/MicroMegas/data/TB2009/May/Rootfiles/dev/";
00012   //TString path = "/home1/chefdevi/micromegasFrameWork/trunk/src/analyse/Gassiplex/Shower_profile/";
00013   
00014   TFile *f = new TFile(path+fileName);
00015   cout<<"reading File : "<<path<<fileName<<endl;
00016   TTree *tree;
00017   tree  = (TTree*) f.Get("t1");cout<<"tree built\n";
00018   int nEvent = tree.GetEntries();
00019   MTEvent *evt = new MTEvent();
00020   TBranch *branch= tree->GetBranch("MTEvent");
00021   branch->SetAddress(&evt);
00022   cout <<"nEvent="<<nEvent<<endl;
00023   int NC = 672;
00024 
00025 
00026 
00027   /*Rootfile to save the histos*/
00028 
00029   TString fileName2= "histo_" + fileName;
00030   TFile * f2 = new TFile(fileName2,"recreate");
00031 
00032 
00033 
00034   /*Distributions of the number of hits*/
00035 
00036   TH1F * hnhit1 = new TH1F("hnhit1","",96,0,96);
00037   TH1F * hnhit2 = new TH1F("hnhit2","",96,0,96);
00038   TH1F * hnhit3 = new TH1F("hnhit3","",96,0,96);
00039   TH1F * hnhit4 = new TH1F("hnhit4","",96*4,0,96*4);
00040 
00041   TH1F * hnhit1_clean = new TH1F("hnhit1_clean","",96*4,0,96*4);
00042   TH1F * hnhit2_clean = new TH1F("hnhit2_clean","",96*4,0,96*4);
00043   TH1F * hnhit3_clean = new TH1F("hnhit3_clean","",96*4,0,96*4);
00044   TH1F * hnhit4_clean = new TH1F("hnhit4_clean","",96*4,0,96*4);
00045 
00046 
00047   /*Distributions of the energy (or ADC)*/
00048 
00049   TH1F * hadc1 = new TH1F("hadc1","",100,0,1024*5);
00050   TH1F * hadc2 = new TH1F("hadc2","",100,0,1024*5);
00051   TH1F * hadc3 = new TH1F("hadc3","",100,0,1024*5);
00052   TH1F * hadc4 = new TH1F("hadc4","",100,0,1024*10);
00053 
00054   TH1F * hadc1_clean = new TH1F("hadc1_clean","",100,0,1024*5);
00055   TH1F * hadc2_clean = new TH1F("hadc2_clean","",100,0,1024*5);
00056   TH1F * hadc3_clean = new TH1F("hadc3_clean","",100,0,1024*5);
00057   TH1F * hadc4_clean = new TH1F("hadc4_clean","",100,0,1024*10);
00058 
00059   TH1F * hadc1_no_sat = new TH1F("hadc1_no_sat","",100,0,1024*5);
00060   TH1F * hadc2_no_sat = new TH1F("hadc2_no_sat","",100,0,1024*5);
00061   TH1F * hadc3_no_sat = new TH1F("hadc3_no_sat","",100,0,1024*5);
00062   TH1F * hadc4_no_sat = new TH1F("hadc4_no_sat","",100,0,1024*10);
00063   /*
00064   TH1F * hadc1_no_sat_clean = new TH1F("hadc1_no_sat_clean","",100,0,1024*5);
00065   TH1F * hadc2_no_sat_clean = new TH1F("hadc2_no_sat_clean","",100,0,1024*5);
00066   TH1F * hadc3_no_sat_clean = new TH1F("hadc3_no_sat_clean","",100,0,1024*5);
00067   TH1F * hadc4_no_sat_clean = new TH1F("hadc4_no_sat_clean","",100,0,1024*10);
00068   */
00069 
00070   /*Distributions of the hit positions in x and y-planes*/
00071 
00072   TH1F * hxpos1 = new TH1F("hxpos1","",40,0,40);
00073   TH1F * hxpos2 = new TH1F("hxpos2","",40,0,40);
00074   TH1F * hxpos3 = new TH1F("hxpos3","",40,0,40);
00075   TH1F * hxpos4 = new TH1F("hxpos4","",40,0,40);
00076 
00077   TH1F * hxpos1_clean = new TH1F("hxpos1_clean","",40,0,40);
00078   TH1F * hxpos2_clean = new TH1F("hxpos2_clean","",40,0,40);
00079   TH1F * hxpos3_clean = new TH1F("hxpos3_clean","",40,0,40);
00080   TH1F * hxpos4_clean = new TH1F("hxpos4_clean","",40,0,40);
00081 
00082   TH1F * hypos1 = new TH1F("hypos1","",40,0,40);
00083   TH1F * hypos2 = new TH1F("hypos2","",40,0,40);
00084   TH1F * hypos3 = new TH1F("hypos3","",40,0,40);
00085   TH1F * hypos4 = new TH1F("hxpos4","",40,0,40);
00086 
00087   TH1F * hypos1_clean = new TH1F("hxpos1_clean","",40,0,40);
00088   TH1F * hypos2_clean = new TH1F("hxpos2_clean","",40,0,40);
00089   TH1F * hypos3_clean = new TH1F("hxpos3_clean","",40,0,40);
00090   TH1F * hypos4_clean = new TH1F("hxpos4_clean","",40,0,40);
00091 
00092   TH1F * hxpos4_all = new TH1F("hxpos4_all","",40,0,40);
00093   TH1F * hypos4_all = new TH1F("hypos4_all","",40,0,40);
00094 
00095 
00096   /*Distributions of the number of hits in the xy-plane*/
00097 
00098   TH2I * hxy1 = new TH2I("hxy1","",32,0,32,32,0,32);
00099   TH2I * hxy2 = new TH2I("hxy2","",32,0,32,32,0,32);
00100   TH2I * hxy3 = new TH2I("hxy3","",32,0,32,32,0,32);
00101   TH2I * hxy4 = new TH2I("hxy4","",32,0,32,32,0,32);
00102 
00103   TH2I * hxy1_all = new TH2I("hxy1_all","",32,0,32,32,0,32);
00104   TH2I * hxy2_all = new TH2I("hxy2_all","",32,0,32,32,0,32);
00105   TH2I * hxy3_all = new TH2I("hxy3_all","",32,0,32,32,0,32);
00106   TH2I * hxy4_all = new TH2I("hxy4_all","",32,0,32,32,0,32);
00107 
00108 
00109   /*Radial distribution of hits w.r.t. beam axis*/
00110 
00111   TH1F * hradial_profile = new TH1F("hradial_profile","",400,-50,50);
00112 
00113 
00114 
00115   /*Set Line Width*/
00116   hnhit1->SetLineWidth(2);
00117   hnhit2->SetLineWidth(2);
00118   hnhit3->SetLineWidth(2);
00119   hnhit4->SetLineWidth(2);
00120 
00121   hnhit1_clean->SetLineWidth(2);
00122   hnhit2_clean->SetLineWidth(2);
00123   hnhit3_clean->SetLineWidth(2);
00124   hnhit4_clean->SetLineWidth(2);
00125 
00126   hadc1->SetLineWidth(2);
00127   hadc2->SetLineWidth(2);
00128   hadc3->SetLineWidth(2);
00129   hadc4->SetLineWidth(2);
00130 
00131   hadc1_clean->SetLineWidth(2);
00132   hadc2_clean->SetLineWidth(2);
00133   hadc3_clean->SetLineWidth(2);
00134   hadc4_clean->SetLineWidth(2);
00135 
00136   hadc1_no_sat->SetLineWidth(2);
00137   hadc2_no_sat->SetLineWidth(2);
00138   hadc3_no_sat->SetLineWidth(2);
00139   hadc4_no_sat->SetLineWidth(2);
00140   /*
00141   hadc1_no_sat_clean->SetLineWidth(2);
00142   hadc2_no_sat_clean->SetLineWidth(2);
00143   hadc3_no_sat_clean->SetLineWidth(2);
00144   hadc4_no_sat_clean->SetLineWidth(2);
00145   */
00146   hxpos1->SetLineWidth(2);
00147   hypos1->SetLineWidth(2);
00148   hxpos2->SetLineWidth(2);
00149   hypos2->SetLineWidth(2);
00150   hxpos3->SetLineWidth(2);
00151   hypos3->SetLineWidth(2);
00152   hxpos4->SetLineWidth(2);
00153   hypos4->SetLineWidth(2);
00154 
00155   hxpos1_clean->SetLineWidth(2);
00156   hypos1_clean->SetLineWidth(2);
00157   hxpos2_clean->SetLineWidth(2);
00158   hypos2_clean->SetLineWidth(2);
00159   hxpos3_clean->SetLineWidth(2);
00160   hypos3_clean->SetLineWidth(2);
00161   hxpos4_clean->SetLineWidth(2);
00162   hypos4_clean->SetLineWidth(2);
00163 
00164 
00165   /*Set Line Color*/
00166   hnhit2->SetLineColor(2);
00167   hnhit3->SetLineColor(4);
00168   hnhit4->SetLineColor(8);
00169 
00170   hnhit2_clean->SetLineColor(2);
00171   hnhit3_clean->SetLineColor(4);
00172   hnhit4_clean->SetLineColor(8);
00173 
00174   hadc2->SetLineColor(2);
00175   hadc3->SetLineColor(4);
00176   hadc4->SetLineColor(8);
00177 
00178   hadc2_clean->SetLineColor(2);
00179   hadc3_clean->SetLineColor(4);
00180   hadc4_clean->SetLineColor(8);
00181 
00182   hxpos2->SetLineColor(2);
00183   hypos2->SetLineColor(2);
00184   hxpos3->SetLineColor(4);
00185   hypos3->SetLineColor(4);
00186   hxpos4->SetLineColor(8);
00187   hypos4->SetLineColor(8);
00188 
00189   hxpos2_clean->SetLineColor(2);
00190   hypos2_clean->SetLineColor(2);
00191   hxpos3_clean->SetLineColor(4);
00192   hypos3_clean->SetLineColor(4);
00193   hxpos4_clean->SetLineColor(8);
00194   hypos4_clean->SetLineColor(8);
00195 
00196 
00197   /*Set Line Style*/
00198   hxpos1_clean->SetLineStyle(2);
00199   hypos1_clean->SetLineStyle(2);
00200   hxpos2_clean->SetLineStyle(2);
00201   hypos2_clean->SetLineStyle(2);
00202 
00203   hxpos3_clean->SetLineStyle(2);
00204   hypos3_clean->SetLineStyle(2);
00205   hxpos4_clean->SetLineStyle(2);
00206   hypos4_clean->SetLineStyle(2);
00207 
00208   hnhit1_clean->SetLineStyle(2);
00209   hnhit2_clean->SetLineStyle(2);
00210   hnhit3_clean->SetLineStyle(2);
00211   hnhit4_clean->SetLineStyle(2);
00212 
00213   hadc1_clean->SetLineStyle(2);
00214   hadc2_clean->SetLineStyle(2);
00215   hadc3_clean->SetLineStyle(2);
00216   hadc4_clean->SetLineStyle(2);
00217 
00218   hadc1_no_sat->SetLineStyle(2);
00219   hadc2_no_sat->SetLineStyle(2);
00220   hadc3_no_sat->SetLineStyle(2);
00221   hadc4_no_sat->SetLineStyle(2);
00222   /*
00223   hadc1_no_sat_clean->SetLineStyle(2);
00224   hadc2_no_sat_clean->SetLineStyle(2);
00225   hadc3_no_sat_clean->SetLineStyle(2);
00226   hadc4_no_sat_clean->SetLineStyle(2);
00227   */
00228 
00229 
00230 
00231   /*Define variables*/
00232 
00233   int xpos1,ypos1,adc1,sum_adc1;
00234   int xpos2,ypos2,adc2,sum_adc2;
00235   int xpos3,ypos3,adc3,sum_adc3;
00236   int xpos4,ypos4,adc4,sum_adc4;
00237 
00238   vector <int> vxpos1,vypos1,vadc1;
00239   vector <int> vxpos2,vypos2,vadc2;
00240   vector <int> vxpos3,vypos3,vadc3;
00241   vector <int> vxpos4,vypos4,vadc4;
00242 
00243   float avgx[3],avgy[3];
00244   float posz1 = 6;
00245   float posz2 = 13;
00246   float posz3 = 19;
00247   float posz4 = 98;
00248   float avgz[3] = {posz1,posz2,posz3};
00249   float avgx1=0,avgx2=0,avgx3=0,avgx4=0;
00250   float avgy1=0,avgy2=0,avgy3=0,avgy4=0;
00251 
00252   int threshold = 35;
00253   float nhit1=0,nhit2=0,nhit3=0,nhit4=0;
00254   int nclean = 0;
00255   bool clean = 0;
00256   int adc_sat = 1000;
00257   bool no_sat1 = 1,no_sat2 = 1;
00258   bool no_sat3 = 1,no_sat4 = 1;
00259 
00260   int event = 0;
00261   int crap_event = 0;
00262 
00263   float x0 = 6.89;
00264   float y0 = 16.48;
00265   float deltax = 0,deltay = 0,radial = 0;
00266 
00267 
00268 
00269 
00270 
00271   /*Loop over all events*/
00272 
00273   //for( int i=0;i<nEvent;i++){
00274     for( int i=0;i<10;i++){
00275   //for( int i=event;i<event+1;i++){
00276     
00277     if (display){cout<<"Event # "<<i<<endl;}
00278 
00279     tree.GetEntry(i);
00280     NC = evt->GetNchannel();
00281 
00282     if (NC == 672){
00283     
00284 
00285     clean=0;
00286     no_sat1 = 1;no_sat2 = 1;no_sat3 = 1;no_sat4 = 1;
00287     nhit1=0,nhit2=0,nhit3=0,nhit4=0;
00288     sum_adc1=0;sum_adc2=0;sum_adc3=0;sum_adc4=0;
00289 
00290     vxpos1.clear();vypos1.clear();vadc1.clear();
00291     vxpos2.clear();vypos2.clear();vadc2.clear();
00292     vxpos3.clear();vypos3.clear();vadc3.clear();
00293     vxpos4.clear();vypos4.clear();vadc4.clear();
00294 
00295     hxy1->Reset();
00296     hxy2->Reset();
00297     hxy3->Reset();
00298     hxy4->Reset();
00299 
00300 
00301     /*Chamber beta 2.1 D*/
00302 
00303     for(int j=0;j<96;j++){
00304 
00305       MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
00306       xpos1 = channel.GetX();
00307       ypos1 = channel.GetY();
00308       adc1  = channel.GetValue();
00309       
00310       if (adc1>=threshold){
00311         
00312         vxpos1.push_back(xpos1);
00313         vypos1.push_back(ypos1);
00314         vadc1.push_back(adc1);
00315         sum_adc1+=adc1;
00316 
00317         /*no saturation if all adcs in chamber are below 1000*/
00318         if (adc1>adc_sat){no_sat1 = 0;}
00319 
00320         if (display){cout<<" x1="<<xpos1<<"  y1="<<ypos1<<"  adc1="<<adc1<<endl;}
00321 
00322         nhit1++;
00323         hxpos1->Fill(xpos1);
00324         hypos1->Fill(ypos1);
00325         hxy1->SetBinContent(xpos1+1,ypos1+1,hxy1->GetBinContent(xpos1+1,ypos1+1)+1);
00326         hxy1_all->SetBinContent(xpos1+1,ypos1+1,hxy1_all->GetBinContent(xpos1+1,ypos1+1)+1);}}
00327 
00328 
00329     /*Chamber beta 2.1 E*/
00330 
00331     for(int j=0+96;j<96+96;j++){
00332 
00333       MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
00334       xpos2 = channel.GetX();
00335       ypos2 = channel.GetY();
00336       adc2  = channel.GetValue();
00337 
00338       if (adc2>=threshold){
00339 
00340         vxpos2.push_back(xpos2);
00341         vypos2.push_back(ypos2);
00342         vadc2.push_back(adc2);
00343         sum_adc2+=adc2;
00344 
00345         /*no saturation if all adcs in chamber are below 1000*/
00346         if (adc2>adc_sat){no_sat2 = 0;}
00347 
00348         if (display){cout<<" x2="<<xpos2<<"  y2="<<ypos2<<"  adc2="<<adc2<<endl;}
00349 
00350         nhit2++;
00351         hxpos2->Fill(xpos2);
00352         hypos2->Fill(ypos2);
00353         hxy2->SetBinContent(xpos2+1,ypos2+1,hxy2->GetBinContent(xpos2+1,ypos2+1)+1);
00354         hxy2_all->SetBinContent(xpos2+1,ypos2+1,hxy2_all->GetBinContent(xpos2+1,ypos2+1)+1);}}
00355 
00356 
00357     /*Chamber beta 2.1 A*/
00358 
00359     for(int j=0+96+96;j<96+96+96;j++){
00360 
00361       MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
00362       xpos3 = channel.GetX();
00363       ypos3 = channel.GetY();
00364       adc3  = channel.GetValue();
00365 
00366       if (adc3>=threshold){
00367 
00368         vxpos3.push_back(xpos3);
00369         vypos3.push_back(ypos3);
00370         vadc3.push_back(adc3);
00371         sum_adc3+=adc3;
00372 
00373         /*no saturation if all adcs in chamber are below 1000*/
00374         if (adc3>adc_sat){no_sat3 = 0;}
00375 
00376         if (display){cout<<" x3="<<xpos3<<"  y3="<<ypos3<<"  adc3="<<adc3<<endl;}
00377 
00378         nhit3++;
00379         hxpos3->Fill(xpos3);
00380         hypos3->Fill(ypos3);
00381         hxy3->SetBinContent(xpos3+1,ypos3+1,hxy3->GetBinContent(xpos3+1,ypos3+1)+1);
00382         hxy3_all->SetBinContent(xpos3+1,ypos3+1,hxy3_all->GetBinContent(xpos3+1,ypos3+1)+1);}}
00383 
00384 
00385     
00386 
00387     /*clean event if one and only one hit in each beta 2.1 chamber*/
00388     clean = (nhit1==1) && (nhit2==1) && (nhit3==1);
00389     
00390     /*clean event if at least one hit in each beta 2.1 chamber*/
00391     //clean = (nhit1>=1) && (nhit2>=1) && (nhit3>=1);
00392 
00393     /*clean event if at least one hit in each beta 2.1 chamber but less than 10*/
00394     //clean = (nhit1>=1) && (nhit2>=1) && (nhit3>=1) && (nhit1<10) && (nhit2<10) && (nhit3<10);        
00395   
00396 
00397 
00398 
00399     /*chamber beta 2.4*/
00400 
00401     for(int j=287;j<672;j++){
00402     //for(int j=287;j<287+1;j++){
00403       MTChannel* channel = (MTChannel*)evt->GetChannels()->UncheckedAt(j);
00404       xpos4 = channel.GetX();
00405       ypos4 = channel.GetY();
00406       adc4  = channel.GetValue();
00407 
00408       if (adc4>=threshold){
00409         
00410         vxpos4.push_back(xpos4);
00411         vypos4.push_back(ypos4);
00412         vadc4.push_back(adc4);
00413         sum_adc4+=adc4;
00414 
00415         /*no saturation if all adcs in chamber are below 1000*/
00416         if (adc4>adc_sat){no_sat4 = 0;}
00417 
00418         if (display){cout<<" x4="<<xpos4<<"  y4="<<ypos4<<"  adc4="<<adc4<<endl;}
00419 
00420         nhit4++;
00421         hxpos4->Fill(xpos4);
00422         hypos4->Fill(ypos4);
00423         hxy4->SetBinContent(xpos4+1,ypos4+1,hxy4->GetBinContent(xpos4+1,ypos4+1)+1);
00424         hxy4_all->SetBinContent(xpos4+1,ypos4+1,hxy4_all->GetBinContent(xpos4+1,ypos4+1)+1);
00425         hxpos4_all->Fill(xpos4);
00426         hypos4_all->Fill(ypos4);}}
00427 
00428 
00429 
00430 
00431     
00432     hnhit1->Fill(nhit1);
00433     hnhit2->Fill(nhit2);
00434     hnhit3->Fill(nhit3);
00435     hnhit4->Fill(nhit4);
00436 
00437     hadc1->Fill(sum_adc1);
00438     hadc2->Fill(sum_adc2);
00439     hadc3->Fill(sum_adc3);
00440     hadc4->Fill(sum_adc4);
00441 
00442 
00443     if (no_sat1){hadc1_no_sat->Fill(sum_adc1);}
00444     if (no_sat2){hadc2_no_sat->Fill(sum_adc2);}
00445     if (no_sat3){hadc3_no_sat->Fill(sum_adc3);}
00446     if (no_sat4){hadc4_no_sat->Fill(sum_adc4);}
00447     
00448 
00449     
00450     if (clean){
00451       /*
00452       if (no_sat1){hadc1_no_sat_clean->Fill(sum_adc1);}
00453       if (no_sat2){hadc2_no_sat_clean->Fill(sum_adc2);}
00454       if (no_sat3){hadc3_no_sat_clean->Fill(sum_adc3);}
00455       if (no_sat4){hadc4_no_sat_clean->Fill(sum_adc4);}
00456       */
00457       if (display){cout<<"clean event!"<<endl;}
00458 
00459       nclean++;
00460       
00461       avgx1=0;avgx2=0;avgx3=0;avgx4=0;
00462       avgy1=0;avgy2=0;avgy3=0;avgy4=0;
00463 
00464       hnhit1_clean->Fill(nhit1);
00465       hnhit2_clean->Fill(nhit2);
00466       hnhit3_clean->Fill(nhit3);
00467       hnhit4_clean->Fill(nhit4);
00468 
00469       hadc1_clean->Fill(sum_adc1);
00470       hadc2_clean->Fill(sum_adc2);
00471       hadc3_clean->Fill(sum_adc3);
00472       hadc4_clean->Fill(sum_adc4);
00473 
00474       for (int k=0;k<nhit1;k++){
00475         hxpos1_clean->Fill(vxpos1[k]);
00476         hypos1_clean->Fill(vypos1[k]);}
00477 
00478       for (int k=0;k<nhit2;k++){
00479         hxpos2_clean->Fill(vxpos2[k]);
00480         hypos2_clean->Fill(vypos2[k]);}
00481 
00482       for (int k=0;k<nhit3;k++){
00483         hxpos3_clean->Fill(vxpos3[k]);
00484         hypos3_clean->Fill(vypos3[k]);}
00485 
00486       for (int k=0;k<nhit4;k++){
00487         
00488         deltax = (vxpos4[k]+0.5) - x0;
00489         deltay = (vypos4[k]+0.5) - y0;
00490 
00491         radial = deltay/fabs(deltay) * sqrt(pow(deltax,2) + pow(deltay,2));
00492 
00493         hradial_profile->Fill(radial);
00494         hxpos4_clean->Fill(vxpos4[k]);
00495         hypos4_clean->Fill(vypos4[k]);}}
00496 
00497 
00498     if (display){cout<<nhit1<<" "<<nhit2<<" "<<nhit3<<" "<<nhit4<<endl;}}
00499 
00500 
00501     else{crap_event++;}}
00502 
00503 
00504       
00505 
00506 
00507 
00508   cout<<"ntotal = "<<nEvent<<"  nclean = "<<nclean<<"  ncrap = "<<crap_event<<endl;
00509 
00510 
00511 
00512   /*save histo in a rootfile*/
00513   f2->Write();
00514   f2->Close();
00515 
00516 
00517 
00518   //hxy4_all->Draw("zcol");
00519   //hradial_profile->Rebin(4);
00520   //hradial_profile->Draw();
00521   //hypos4->Draw();
00522 
00523 
00524 
00525 
00526   /*Plot distributions*/
00527 
00528 
00529   bool nhit_distribution = 0;
00530   bool adc_distribution = 0;
00531   bool chamber_xy_display = 0;
00532   bool xyz_plot = 0;
00533 
00534   bool number_hit_distrib_4_chamber = 0;
00535   bool number_hit_distrib_beta24 = 0;
00536   bool x_y_hit_distrib_beta21 = 0;
00537   bool x_y_hit_distrib_beta22 = 0;
00538   bool x_y_hit_distrib_beta23 = 0;
00539   bool x_y_hit_distrib_beta24 = 0;
00540 
00541 
00542   if (nhit_distribution){
00543 
00544     TCanvas * myc = new TCanvas("myc");
00545     myc->Divide(2,2);
00546     myc->cd(1);
00547     hnhit1->Draw();
00548     hnhit1_clean->Draw("same");
00549     myc->cd(2);
00550     hnhit2->Draw();
00551     hnhit2_clean->Draw("same");
00552     myc->cd(3);
00553     hnhit3->Draw();
00554     hnhit3_clean->Draw("same");
00555     myc->cd(4);
00556     hnhit4->Draw();
00557     hnhit4_clean->Draw("same");}
00558 
00559   if (adc_distribution){
00560 
00561     TCanvas * myc = new TCanvas("myc");
00562     myc->Divide(2,2);
00563     myc->cd(1);
00564     hadc1->Draw();
00565     hadc1_clean->Draw("same");
00566     myc->cd(2);
00567     hadc2->Draw();
00568     hadc2_clean->Draw("same");
00569     myc->cd(3);
00570     hadc3->Draw();
00571     hadc3_clean->Draw("same");
00572     myc->cd(4);
00573     hadc4->Draw();
00574     hadc4_clean->Draw("same");}
00575 
00576 
00577   if (chamber_xy_display){
00578 
00579     TCanvas * myc3 = new TCanvas("myc3");
00580     myc3->Divide(2,2);
00581     myc3->cd(1);
00582     hxy1->Draw("zcol");
00583     myc3->cd(2);
00584     hxy2->Draw("zcol");
00585     myc3->cd(3);
00586     hxy3->Draw("zcol");
00587     myc3->cd(4);
00588     hxy4->Draw("zcol");}
00589 
00590 
00591   if (xyz_plot){
00592     
00593     TCanvas * myc4 = new TCanvas("myc4");
00594 
00595     tgxy->SetMarkerStyle(21);
00596     tgzx->SetMarkerStyle(21);
00597     tgzy->SetMarkerStyle(21);
00598     
00599     myc4->Divide(2,2);
00600     myc4->cd(1);
00601     tgxy->Draw("ap");
00602     myc4->cd(2);
00603     tgzx->Draw("ap");
00604     myc4->cd(3);
00605     tgzy->Draw("ap");}
00606 
00607 
00608   if (number_hit_distrib_4_chamber){
00609     
00610     hnhit1->DrawNormalized();
00611     hnhit2->DrawNormalized("same");
00612     hnhit3->DrawNormalized("same");
00613     hnhit4->DrawNormalized("same");}
00614 
00615   if (number_hit_distrib_beta24){
00616 
00617     hnhit4->DrawNormalized();
00618     hnhit4_clean->DrawNormalized("same");}
00619 
00620   if (x_y_hit_distrib_beta21){
00621 
00622     hxpos1->DrawNormalized();
00623     hxpos1_clean->DrawNormalized("same");
00624     hypos1->DrawNormalized("same");
00625     hypos1_clean->DrawNormalized("same");}
00626 
00627   if (x_y_hit_distrib_beta22){
00628 
00629     hxpos2->DrawNormalized();
00630     hxpos2_clean->DrawNormalized("same");
00631     hypos2->DrawNormalized("same");
00632     hypos2_clean->DrawNormalized("same");}
00633 
00634   if (x_y_hit_distrib_beta23){
00635 
00636     hxpos3->DrawNormalized();
00637     hxpos3_clean->DrawNormalized("same");
00638     hypos3->DrawNormalized("same");
00639     hypos3_clean->DrawNormalized("same");}
00640 
00641   if (x_y_hit_distrib_beta24){
00642 
00643     hxpos4->DrawNormalized();
00644     hxpos4_clean->DrawNormalized("same");
00645     hypos4->DrawNormalized("same");
00646     hypos4_clean->DrawNormalized("same");}
00647 
00648 
00649 }
/Users/jacquem/NetBeansProjects/MicromegasFramework/trunk/src/analyse/Gassiplex/Shower_profile/Profile_no_tracking.C
