Changeset 2255 for trunk/MagicSoft/Mars/manalysis
- Timestamp:
- 07/01/03 10:46:07 (21 years ago)
- Location:
- trunk/MagicSoft/Mars/manalysis
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/manalysis/MCT1PadONOFF.cc
r2206 r2255 1401 1401 Int_t j = pix.GetPixId(); 1402 1402 1403 Double_t ratioArea = fCam->GetPixRatio(j); 1403 // GetPixRatio returns (area of pixel 0 / area of current pixel) 1404 Double_t ratioArea = 1.0 / fCam->GetPixRatio(j); 1404 1405 1405 1406 MPedestalPix &ppix = (*fPed)[j]; -
trunk/MagicSoft/Mars/manalysis/MCT1PadSchweizer.cc
r2237 r2255 540 540 Int_t j = pix.GetPixId(); 541 541 542 Double_t ratioArea = fCam->GetPixRatio(j); 542 // GetPixRatio returns (area of pixel 0 / area of current pixel) 543 Double_t ratioArea = 1.0 / fCam->GetPixRatio(j); 543 544 544 545 MPedestalPix &ppix = (*fPed)[j]; -
trunk/MagicSoft/Mars/manalysis/MSigmabar.cc
r2237 r2255 136 136 137 137 const Int_t idx = cerpix.GetPixId(); 138 const Double_t area = geom.GetPixRatio(idx);139 140 138 if (idx == 0) 141 139 { … … 145 143 } 146 144 145 // ratio is the area of pixel 0 146 // divided by the area of the current pixel 147 const Double_t ratio = geom.GetPixRatio(idx); 148 147 149 const MGeomPix &gpix = geom[idx]; 148 150 … … 157 159 continue; 158 160 159 if (area < 1.5) 161 //if (area < 1.5) 162 if (ratio > 0.5) 160 163 { 161 164 innerPixels[sector]++; 162 innerSquaredSum[sector]+= sigma*sigma / area;165 innerSquaredSum[sector]+= sigma*sigma * ratio; 163 166 } 164 167 else 165 168 { 166 169 outerPixels[sector]++; 167 outerSquaredSum[sector]+= sigma*sigma / area;170 outerSquaredSum[sector]+= sigma*sigma * ratio; 168 171 } 169 172 } … … 243 246 244 247 } 248 249 250 251 252 253 254 255 256
Note:
See TracChangeset
for help on using the changeset viewer.