Changeset 900 for trunk/MagicSoft/Mars/mhist/MHMcRate.cc
- Timestamp:
- 07/31/01 10:46:18 (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mhist/MHMcRate.cc
r896 r900 51 51 fBackSim=-1.0; // Number of simulated showers for the background 52 52 53 fTrig Rate= -1.0; // Trigger rate in Hz54 fTrig RateError= -1.0; // Estimated error for the trigger rate in Hz53 fTriggerRate= -1.0; // Trigger rate in Hz 54 fTriggerRateError= -1.0; // Estimated error for the trigger rate in Hz 55 55 } 56 56 … … 67 67 fFlux0=-1.0; // dn/dE = fFlux0 * E^{-a} 68 68 69 fShow Rate= -1.0; // Showers rate in Hz70 fShow RateError=0.0; // Estimated error of shower rate in Hz69 fShowerRate= -1.0; // Showers rate in Hz 70 fShowerRateError=0.0; // Estimated error of shower rate in Hz 71 71 } 72 72 … … 85 85 fFlux0=-1.0; // dn/dE = fFlux0 * E^{-a} 86 86 87 fShow Rate= showrate; // Showers rate in Hz88 fShow RateError=sqrt(showrate); // Estimated error of shower rate in Hz87 fShowerRate= showrate; // Showers rate in Hz 88 fShowerRateError=sqrt(showrate); // Estimated error of shower rate in Hz 89 89 } 90 90 … … 103 103 fFlux0=flux0; // dn/dE = fFlux0 * E^{-a} 104 104 105 fShow Rate= -1.0;106 fShow RateError=0.0;105 fShowerRate= -1.0; 106 fShowerRateError=0.0; 107 107 } 108 108 … … 130 130 // update the limits for energy, theta, phi and impact parameter 131 131 // 132 void MHMcRate:: Boundaries(Float_t energy, Float_t theta,133 Float_t phi, Float_t impact)132 void MHMcRate::UpdateBoundaries(Float_t energy, Float_t theta, 133 Float_t phi, Float_t impact) 134 134 { 135 135 // It updates the limit values … … 163 163 const Double_t epowmin = pow(fEnergyMin, specidx); 164 164 165 fShow Rate = fFlux0/specidx*(epowmax-epowmin);165 fShowerRate = fFlux0/specidx*(epowmax-epowmin); 166 166 167 167 if (fPartId!=1) 168 fShow Rate *= (fPhiMax-fPhiMin)*(cos(fThetaMax)-cos(fThetaMin));168 fShowerRate *= (fPhiMax-fPhiMin)*(cos(fThetaMax)-cos(fThetaMin)); 169 169 170 170 const Double_t impactdiff = fImpactMax-fImpactMin; 171 171 172 fShow Rate *= TMath::Pi()*(impactdiff/100.0*impactdiff/100.0);173 174 fShow RateError = sqrt(fShowRate);172 fShowerRate *= TMath::Pi()*(impactdiff/100.0*impactdiff/100.0); 173 174 fShowerRateError = sqrt(fShowerRate); 175 175 176 176 const Double_t anal2 = 1.0-anal*160.0e-9; … … 179 179 // Then the trigger rate and its error is evaluated 180 180 if(fBackTrig<0){ 181 fTrig RateError = sqrt((trig*fShowRate*fShowRate/(simu*simu)) +181 fTriggerRateError = sqrt((trig*fShowerRate*fShowerRate/(simu*simu)) + 182 182 (anal2*anal2*1/(fBackSim*back2*back2))); 183 183 fBackTrig=0; 184 184 } 185 185 else 186 fTrig RateError = sqrt((trig*fShowRate*fShowRate/(simu*simu)) +186 fTriggerRateError = sqrt((trig*fShowerRate*fShowerRate/(simu*simu)) + 187 187 (anal2*anal2*fBackTrig/(back2*back2))); 188 188 189 fTrig Rate = trig*fShowRate/simu + anal2*fBackTrig/back2;189 fTriggerRate = trig*fShowerRate/simu + anal2*fBackTrig/back2; 190 190 } 191 191 … … 194 194 // print the trigger rate 195 195 // 196 void MHMcRate::Print( )197 { 198 *fLog << "Incident rate " << fShow Rate << " Hz " << endl;199 *fLog << "Trigger Rate " << fTrig Rate << " +- " << fTrigRateError << " Hz" << endl;196 void MHMcRate::Print(Option_t *) 197 { 198 *fLog << "Incident rate " << fShowerRate << " Hz " << endl; 199 *fLog << "Trigger Rate " << fTriggerRate << " +- " << fTriggerRateError << " Hz" << endl; 200 200 } 201 201 … … 204 204 // draw the trigger rate 205 205 // 206 void MHMcRate::Draw(Option_t * option)206 void MHMcRate::Draw(Option_t *) 207 207 { 208 208 *fLog << "To be iplemented" << endl;
Note:
See TracChangeset
for help on using the changeset viewer.