Changeset 9521 for trunk/MagicSoft
- Timestamp:
- 10/26/09 15:31:08 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/msimcamera/MSimRandomPhotons.cc
r9425 r9521 34 34 // FIXME: We should add the wavelength distribution. 35 35 // 36 // The artificial night sky background rate is calculated as follows: 37 // 38 // * The photon detection efficiency vs. wavelength of the detector is obtained 39 // from "PhotonDetectionEfficiency" of type "MParSpline" 40 // 41 // * The angular acceptance of the light collectors is obtained 42 // from "ConesAngularAcceptance" of type "MParSpline" 43 // 44 // * The spectral acceptance of the light collectors is obtained 45 // from "ConesTransmission" of type "MParSpline" 46 // 47 // * The reflectivity of the mirrors vs wavelength is obtained 48 // from "MirrorReflectivity" of type "MParSpline" 49 // 50 // The rate is then calculated as 51 // 52 // R = R0 * Ai * f 53 // 54 // R0 is the night sky background rate as given in Eckart's paper. Ai the 55 // area of the cones acceptance window, f is given as: 56 // 57 // f = nm * Min(Ar, sr*d^2) 58 // 59 // with 60 // 61 // nm being the integral of the product of the mirror reflectivity, the cone 62 // transmission and the photon detection efficiency. 63 // 64 // d the distance of the focal plane to the mirror 65 // 66 // Ar is the total reflective area of the reflector 67 // 68 // sr is the effective solid angle corresponding to the integral of the 69 // cones angular acceptance 70 // 71 // 36 72 // Input Containers: 37 73 // fNameGeomCam [MGeomCam] … … 151 187 } 152 188 189 *fLog << warn << "FIXME: A check is needed that the PDE is 0 at both ends!" << endl; 190 153 191 const MParSpline *s1 = (MParSpline*)pList->FindObject("PhotonDetectionEfficiency", "MParSpline"); 154 192 const MParSpline *s2 = (MParSpline*)pList->FindObject("ConesAngularAcceptance", "MParSpline"); 155 193 const MParSpline *s3 = (MParSpline*)pList->FindObject("MirrorReflectivity", "MParSpline"); 194 const MParSpline *s5 = (MParSpline*)pList->FindObject("ConesTransmission", "MParSpline"); 156 195 157 196 const Double_t d2 = fGeom->GetCameraDist()*fGeom->GetCameraDist(); 158 const Double_t pde = s1 && s1->GetSpline() ? s1->GetSpline()->Integral() : 1;197 // const Double_t pde = s1 && s1->GetSpline() ? s1->GetSpline()->Integral() : 1; 159 198 const Double_t sr = s2 && s2->GetSpline() ? s2->GetSpline()->IntegralSolidAngle() : 1; 160 const Double_t mir = s3 && s3->GetSpline() ? s3->GetSpline()->Integral() : 1;199 // const Double_t mir = s3 && s3->GetSpline() ? s3->GetSpline()->Integral() : 1; 161 200 const Double_t Ar = r->GetA()/1e4; 162 201 … … 168 207 MParSpline *s4 = (MParSpline*)s1->Clone(); 169 208 s4->Multiply(*s3->GetSpline()); 209 s4->Multiply(*s5->GetSpline()); 170 210 171 211 const Double_t nm = s4 && s4->GetSpline() ? s4->GetSpline()->Integral() : 1; … … 180 220 *fLog << "Reflector area: " << Form("%.2f", Ar) << "m^2" << endl; 181 221 *fLog << "Resulting eff. collection area: " << Form("%.1f", TMath::Min(Ar, sr*d2)) << "m^2" << endl; 182 *fLog << "Eff. wavelength band (PDE): " << Form("%.1f", pde) << "nm" << endl;183 *fLog << "Eff. wavelength band (Mirror): " << Form("%.1f", mir) << "nm" << endl;184 *fLog << "Eff. wavelength band ( PDE+MIR): " << Form("%.1f", nm) << "nm" << endl;222 // *fLog << "Eff. wavelength band (PDE): " << Form("%.1f", pde) << "nm" << endl; 223 // *fLog << "Eff. wavelength band (Mirror): " << Form("%.1f", mir) << "nm" << endl; 224 *fLog << "Eff. wavelength band (MIR+Cone+PDE): " << Form("%.1f", nm) << "nm" << endl; 185 225 *fLog << "Pixel area of " << fNameGeomCam << "[0]: " << Form("%.2e", (*fGeom)[0].GetA()*conv*conv) << "sr" << endl; 186 226 //*fLog << "Effective angular acceptance: " << sr << "sr" << endl; … … 290 330 // 291 331 // FrequencyFixed: 0.040 292 // FrequencyNSB: 0.040332 // FrequencyNSB: 5.8 293 333 // 294 334 // The fixed frequency is given in units fitting the units of the time.
Note:
See TracChangeset
for help on using the changeset viewer.