Changeset 862 for trunk


Ignore:
Timestamp:
07/10/01 14:19:56 (23 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mmontecarlo/MMcThresholdCalc.cc

    r861 r862  
    8787}
    8888
     89// --------------------------------------------------------------------------
     90//
     91// connect Monte Carlo data with this task
     92//
    8993Bool_t MMcThresholdCalc::PreProcess(MParList* pList)
    9094{
    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
    9496    // condition (MMcTrig) or severl of them (MMcTrig;#.)
    9597
     
    140142}
    141143
     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//
    142150Bool_t MMcThresholdCalc::Process()
    143151{
    144 
    145     // The histograms are filled with log10 of the energy for triggered
    146     // events and weighted with 1/E because it is needed the dN/dE vs. logE
    147     // distribution to get the energy threshold.
    148 
    149152    const Float_t energy   = fMcEvt->GetEnergy();
    150153    const Float_t lg10     = log10(energy);
     
    152155
    153156    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);
    160159
    161160    return kTRUE;
    162161}
    163162
     163// --------------------------------------------------------------------------
     164//
     165// fit the energy distribution to get the threshold
     166// Some iterations are done to be sure the fit parameters converge.
     167//
    164168Bool_t MMcThresholdCalc::PostProcess()
    165169{
    166     // fit the energy distribution to get the threshold
    167     // Some iterations are done to be sure the fit parameters converge.
    168 
    169170    const Float_t sqrt2 = sqrt(2);
    170171
Note: See TracChangeset for help on using the changeset viewer.