WebSVN – f9daq – Diff – //cvi/instr/AitSipmDAQ/data/aitana.cxx

 #include <TStyle.h>
 #include <TCanvas.h>
 #include <TLegend.h>
 // Program za analizo podatkov zajetih z AITSiPMDAQ
-int aitana( const char *fname="a1.dat", int range=10, float cut=500, int *offset = NULL)
+int aitana( const char *fname="a1.dat", int range=5, float cut=300, int *offset = NULL)
+{
         int debug=0;
        //int off[4]={507,521,497,459};
+                //int off[4]={198,203,190,150}; /// run 1
-                int off[4]={87,125,95,80};
+                int off[4]={120,120,120,120}; /// run 6
                 if (offset) for (int i=0;i<4;i++){off[i]=offset[i];};
                 gStyle->SetOptStat(0);
         char rootname[0xFF];
                 sprintf(rootname, "%s.root", fname);
         TFile *froot = new  TFile(rootname,"RECREATE");
             TLegend * leg = new TLegend(0.8, 0.8, 1, 1);
                 TH1F *hadc[4];
                 TH1F *hadcsum;
                 TH1F *hadcsumx;
                 TH1F *hadcsumy;
-                TH2F *hxy;
+                TH2F *hxy,*hxy1,*hxy2;
                 TH2F *h2d=NULL;
                 TH2F *h2dx=NULL;
                 TH2F *h2dy=NULL;
                 TH2F *hxcorr=NULL;
                 TH2F *hycorr=NULL;
          hadc[i]=new TH1F(histname,histname,200,0.5,200*range+0.5);
+                }
                 hadcsum=new TH1F("adcsum","adcsum",200,0.5,8*200*range+0.5);
                 hadcsumx=new TH1F("adcsumx","adcsumx",200,0.5,200*range+0.5);
                 hadcsumy=new TH1F("adcsumy","adcsumy",200,0.5,200*range+0.5);
-        hxy=new TH2F("hxy","Center of gravity;ycog(a.u.);xcog(a.u.)",200,0,1,200,0,1);
+        hxy=new TH2F("hxy","Center of gravity;ycog(a.u.);xcog(a.u.)",200,-0.1,1.1,200,-0.1,1.1);
+        hxy1=new TH2F("hxy1","Center of gravity for position 1;ycog(a.u.);xcog(a.u.)",200,-0.1,1.1,200,-0.1,1.1);
+        hxy2=new TH2F("hxy2","Center of gravity for position 2;ycog(a.u.);xcog(a.u.)",200,-0.1,1.1,200,-0.1,1.1);
                 FILE *fin=fopen(fname,"rb");
         if (fin==NULL) {
                 printf("Error opening file %s\n",fname);
                 return 0;
+        }
                                                            sum=0;
                                                            for (int j=0;j<4;j++){
                                                                  a[j]=adc[j]&0xFFF;
                                                                  //sum+=a[j];
                                                                  hadc[j]->Fill(a[j]);
-                                 b[j]=a[j]-off[j]+10;
+                                 b[j]=a[j]-off[j];
+                                                                 if (b[j]< 0) b[j]=0;
                                  sum+=b[j];
+                                                           }
                                                            hadcsum->Fill(sum);
                                                            hadcsumx->Fill(a[0]+a[1]);
                                                            hadcsumy->Fill(a[2]+a[3]);
                                                            if (h2dx) h2dx->Fill(posx,posy,b[0]+b[1]);
                                                            if (h2dy) h2dy->Fill(posx,posy,b[2]+b[3]);
                                                            if (sum>cut) {
                                                                    hxy->Fill(y,x);
+                                   if (pos[3]==poshdr[6]/4 && pos[4]==poshdr[7]/4) hxy1->Fill(y,x);
+                                   if (pos[3]==poshdr[6]/4*3 && pos[4]==poshdr[7]/4*3) hxy2->Fill(y,x);
                                                                    if (hxcorr) hxcorr->Fill(pos[0],y);
                                                                    if (hycorr) hycorr->Fill(pos[1],x);
+                                                           }
+                                                   }
                 */
                 c->Print(rootname,"pdf");
                 c= new TCanvas("c1","ait reconstruction",750,50,700,700);
-                c->Divide(2,3);
+                c->Divide(2,4);
                 if (!h2d){
               c->Clear();
-                  c->Divide(1,4);
+                  c->Divide(2,3);
+                }
                 hxy->SetMinimum(-1);
                 hxy->SetMaximum(1000);
                 c->cd(1)->SetLogz(); hxy->DrawCopy("colz");
+        c->cd(3)->SetLogz(); hxy1->DrawCopy("colz");
+        c->cd(4)->SetLogz(); hxy2->DrawCopy("colz");
-                c->cd(4)->SetLogy(); hadcsum->DrawCopy();
+                c->cd(2)->SetLogy(); hadcsum->DrawCopy();
                 if (h2d){
                   hxcorr->SetMinimum(-1);
                   hycorr->SetMinimum(-1);
-                  c->cd(5); hxcorr->DrawCopy("colz");
+                  c->cd(7); hxcorr->DrawCopy("colz");
-                  c->cd(6); hycorr->DrawCopy("colz");
+                  c->cd(8); hycorr->DrawCopy("colz");
+                }
                 TH1D *px = hxy->ProjectionX("_px",0,-1);
-                c->cd(2); px->DrawCopy();
+                c->cd(5); px->DrawCopy();
             TH1D *py = hxy->ProjectionY("_py",0,-1);
-                c->cd(3); py->DrawCopy();
+                c->cd(6); py->DrawCopy();
                 c->Modified();
                 c->Update();
                 c->Print(rootname,"pdf");

Subversion Repositories f9daq

(root)//cvi/instr/AitSipmDAQ/data/aitana.cxx – Rev 218 → 315