Changeset 1180 for trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc
- Timestamp:
- 01/15/02 11:48:08 (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/manalysis/MCerPhotCalc.cc
r1145 r1180 124 124 while (pixel.Next()) 125 125 { 126 127 UInt_t pixid = pixel.GetPixelId(); 126 const UInt_t pixid = pixel.GetPixelId(); 128 127 129 128 const MPedestalPix &ped = (*fPedestals)[pixid]; 130 129 131 130 // 132 // sanity check 131 // sanity check (old MC files sometimes have pixids>577) 133 132 // 134 if (!fPedestals->CheckBounds(pixid)){ 135 *fLog<<inf<<"Pixel ID larger than camera ... Skip the event"<<endl; 133 if (!fPedestals->CheckBounds(pixid)) 134 { 135 *fLog << inf << "Pixel ID larger than camera... skipping event." << endl; 136 136 return kCONTINUE; 137 137 } … … 139 139 Float_t nphot = (Float_t)pixel.GetSumHiGainSamples(); 140 140 141 // 141 142 // We check that the pixel is not empty, if it is empty 142 143 // we won't substract the pedestal. Empty means that it has 143 144 // 0 signal in all the slices. 144 145 if(nphot!=0) nphot -= fRunHeader->GetNumSamplesHiGain()*ped.GetMean(); 145 // 146 if (nphot!=0) 147 nphot -= fRunHeader->GetNumSamplesHiGain()*ped.GetMean(); 146 148 147 149 fCerPhotEvt->AddPixel(pixid, nphot, sqrt(fRunHeader->GetNumSamplesHiGain())*ped.GetSigma());
Note:
See TracChangeset
for help on using the changeset viewer.