Changeset 2015 for trunk/MagicSoft/Mars/mhist/MHSigmaTheta.cc
- Timestamp:
- 04/28/03 09:52:57 (21 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mhist/MHSigmaTheta.cc
r1992 r2015 180 180 { 181 181 Double_t theta = fMcEvt->GetTelescopeTheta()*kRad2Deg; 182 Double_t my Sig = fSigmabar->Calc(*fCam, *fPed, *fEvt);183 184 fSigmaTheta.Fill(theta, my Sig);182 Double_t mysig = fSigmabar->Calc(*fCam, *fPed, *fEvt); 183 184 fSigmaTheta.Fill(theta, mysig); 185 185 186 186 const UInt_t npix = fEvt->GetNumPixels(); … … 188 188 for (UInt_t i=0; i<npix; i++) 189 189 { 190 MCerPhotPix cerpix = (*fEvt)[i]; 191 if (!cerpix.IsPixelUsed()) 192 continue; 193 194 const Int_t id = cerpix.GetPixId(); 195 const MPedestalPix &pix = (*fPed)[id]; 196 197 const Double_t sigma = pix.GetMeanRms(); 198 const Double_t area = fCam->GetPixRatio(id); 199 200 201 fSigmaPixTheta.Fill(theta, (Double_t)id, sigma); 202 203 const Double_t diff = sigma*sigma/area - mySig*mySig; 204 fDiffPixTheta.Fill(theta, (Double_t)id, diff); 190 MCerPhotPix cerpix = (*fEvt)[i]; 191 if (!cerpix.IsPixelUsed()) 192 continue; 193 194 const Int_t id = cerpix.GetPixId(); 195 const MPedestalPix &pix = (*fPed)[id]; 196 197 const Double_t sigma = pix.GetMeanRms(); 198 const Double_t area = fCam->GetPixRatio(id); 199 200 fSigmaPixTheta.Fill(theta, (Double_t)id, sigma); 201 202 const Double_t diff = sigma*sigma/area - mysig*mysig; 203 fDiffPixTheta.Fill(theta, (Double_t)id, diff); 205 204 } 206 205 207 206 return kTRUE; 208 }209 210 // --------------------------------------------------------------------------211 //212 // Draw a copy of the histogram213 //214 TObject *MHSigmaTheta::DrawClone(Option_t *opt)215 {216 return MH::DrawClone(opt, 900, 900);217 207 } 218 208
Note:
See TracChangeset
for help on using the changeset viewer.