- Timestamp:
- 07/10/01 14:19:56 (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc
r861 r862 87 87 } 88 88 89 // -------------------------------------------------------------------------- 90 // 91 // connect Monte Carlo data with this task 92 // 89 93 Bool_t MMcThresholdCalc::PreProcess(MParList* pList) 90 94 { 91 // connect Monte Carlo data with this task 92 93 // This task has into accout if the root file has one trigger 95 // This task takes into accout if the root file has one trigger 94 96 // condition (MMcTrig) or severl of them (MMcTrig;#.) 95 97 … … 140 142 } 141 143 144 // -------------------------------------------------------------------------- 145 // 146 // The histograms are filled with log10 of the energy for triggered 147 // events and weighted with 1/E because it is needed the dN/dE vs. logE 148 // distribution to get the energy threshold. 149 // 142 150 Bool_t MMcThresholdCalc::Process() 143 151 { 144 145 // The histograms are filled with log10 of the energy for triggered146 // events and weighted with 1/E because it is needed the dN/dE vs. logE147 // distribution to get the energy threshold.148 149 152 const Float_t energy = fMcEvt->GetEnergy(); 150 153 const Float_t lg10 = log10(energy); … … 152 155 153 156 for (unsigned int i=0; i<fDimension; i++) 154 { 155 if (fMcTrig[i]->GetFirstLevel()<=0) 156 continue; 157 158 fHMcEnergy[i]->Fill(lg10, reciproc); 159 } 157 if (fMcTrig[i]->GetFirstLevel()>0) 158 fHMcEnergy[i]->Fill(lg10, reciproc); 160 159 161 160 return kTRUE; 162 161 } 163 162 163 // -------------------------------------------------------------------------- 164 // 165 // fit the energy distribution to get the threshold 166 // Some iterations are done to be sure the fit parameters converge. 167 // 164 168 Bool_t MMcThresholdCalc::PostProcess() 165 169 { 166 // fit the energy distribution to get the threshold167 // Some iterations are done to be sure the fit parameters converge.168 169 170 const Float_t sqrt2 = sqrt(2); 170 171
Note:
See TracChangeset
for help on using the changeset viewer.