Changeset 5127 for trunk/MagicSoft/Mars/mtemp/mmpi/macros/calibration_shape.C

Timestamp:

09/24/04 18:15:19 (20 years ago)

Author:

hbartko

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mmpi/macros/calibration_shape.C (modified) (8 diffs)

trunk/MagicSoft/Mars/mtemp/mmpi/macros/calibration_shape.C

@@ -30 +30 @@
 
 Int_t pedr = 36038;
-Int_t calr = 36040;
+Int_t calr = 36042;
 //Int_t datar = 12517;
 
 
 
-void calibration_shape(const TString inpath = "/.magic/magicserv01/scratch/hbartko/rootdata/2004_09_06/", Int_t pedruns = pedr, Int_t calruns = calr, Int_t ipix=297)
+void calibration_shape(const TString inpath = "/.magic/magicserv01/scratch/hbartko/rootdata/2004_09_06/", Int_t pedruns = pedr, Int_t calruns = calr, Int_t ipix=316)
 
 {
@@ -112 +112 @@
   plist2.AddToList(&badcam);
   
-  gLog << endl;;
-  gLog << "Calculate MCalibrationCam from Runs " << cruns.GetRunsAsString() << endl;
-  gLog << endl;
-  
+   
   MReadMarsFile read2("Events");
   read2.DisableAutoScheme();
@@ -149 +146 @@
   MExtractTimeHighestIntegral   timeext;
   timeext.SetRange(0, 14, 0, 14);
+  //timeext.SetWindowSize(10,10);
   timeext.SetWindowSize(WindowSize,WindowSize);
   
@@ -195 +193 @@
   int count = 1; 
   
-  TH2F * shape = new TH2F("shape","shape",30,0,30,1000,-20,250);
+  TH2F * shape = new TH2F("shape","shape",300,0,30,1000,-20,250);
   TH2F * shape_corr = new TH2F("shape_corr","shape_corr",300,0,30,1000,-20,250);
   TH2F * shape_corr_all = new TH2F("shape_corr_all","shape_corr_all",300,0,30,1000,-20,250);
@@ -225 +223 @@
     sigcam.GetPixelContent(sig, ipix, geomcam, 0);
     
+    Byte_t sat = (sigcam)[ipix].GetNumHiGainSaturated();
+
+
     htime->Fill(time);
     hsig->Fill(sig);
@@ -243 +244 @@
     const Float_t ABoffs = pedcam[ipix].GetPedestalABoffset();
     
+   
     Float_t PedMean[2];
     PedMean[0] = ped_mean + ABoffs;
@@ -254 +256 @@
       shape->Fill(sample+15.5,logains[sample]-PedMean[(sample+ABFlag)&0x1]);
       shape_corr->Fill(sample+0.5+5.-time,higains[sample]-PedMean[(sample+ABFlag)&0x1]);
-      shape_corr->Fill(sample+15.5+5.-time,logains[sample]-PedMean[(sample+ABFlag)&0x1]);
-      shape_corr_all->Fill(sample+0.5+5.-time,(higains[sample]-PedMean[(sample+ABFlag)&0x1])/sig*250);
-      shape_corr_all->Fill(sample+15.5+5.-time,(logains[sample]-PedMean[(sample+ABFlag)&0x1])/sig*250);
+      shape_corr->Fill(sample+15.5+5.-time,logains[sample]-PedMean[(sample+nh+ABFlag)&0x1]);
+      if (sat==0){
+        shape_corr_all->Fill(sample+0.5+5.-time,(higains[sample]-PedMean[(sample+ABFlag)&0x1])/sig*250);
+        shape_corr_all->Fill(sample+15.5+5.-time,(logains[sample]-PedMean[(sample+nh+ABFlag)&0x1])/sig*250);
+      }
     }
         
@@ -350 +354 @@
   htime->Draw();
 
+
+  TCanvas * c7 = new TCanvas("c7","c7",600,400);
+  c7->cd();
+  c7->SetGridx();
+  c7->SetGridy();
+  c7->SetFillColor(0);
+  //  c7->SetBordermode(0);
+    
+  TProfile * shape_corr_all_pfx = (TProfile*)shape_corr_all->ProfileX();
+  
+  shape_corr_all_pfx->SetStats(0);
+  TString title6 = "Average Calibration Signal Shape, Amplitude + Arrival Time Corrected, Pixel";
+  title6+=ipix;
+  shape_corr_all_pfx->SetTitle(title6);
+  shape_corr_all_pfx->SetXTitle("FADC sample no. + (<t_{arrival}> - t_{arrival}) / T_{ADC}");
+  shape_corr_all_pfx->SetYTitle("signal [a.u.]");
+  shape_corr_all_pfx->SetMarkerStyle(20);
+  shape_corr_all_pfx->SetMarkerColor(4);
+  shape_corr_all_pfx->SetMarkerSize(0.5);
+  shape_corr_all_pfx->Draw();
 }