Ignore:
Timestamp:
08/04/08 11:44:11 (17 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mjobs/MJCalibrateSignal.cc

    r9038 r9077  
    540540    // an error is raised if more than 50% of the events were skipped.
    541541    //
     542    // Since unsuitable pixels are not counted in the ideal case
     543    // there is no empty pixel at all.
     544    //
    542545    MFCosmics fcosmicscal;
    543546    fcosmicscal.SetNamePedestalCam("MPedestalFundamental"); //CORRECT?
     
    548551
    549552    //------------------------------
    550     /*
    551      MFCosmics fcosmicsped; // Def=20% empty pixels
    552      fcosmicsped.SetNamePedestalCam("MPedestalFundamental"); //CORRECT?
    553      MContinue contcosped(&fcosmicsped, "ContCosmicsPed");
    554      contcosped.SetInverted();
    555     */
     553    //
     554    // Remove events which are too bright. These events could either
     555    // be calibration events or events which produce many many
     556    // saturating low-gains and thus spoil the dead pixel plot.
     557    // These events don't seem to be physical events so we can fairly
     558    // remove them (they are no analysable gammas anyway)
     559    // So we remove all events which have less than 5% empty pixels.
     560    //
     561    MFCosmics fcosmicsbright;
     562    fcosmicsbright.SetNamePedestalCam("MPedestalFundamental");
     563    fcosmicsbright.SetMaxEmptyPixels(0.10);
     564    fcosmicsbright.SetMinAcceptedFraction(0.9);
     565
     566    MContinue contbright(&fcosmicsbright, "ContBrightEvts");
     567    contbright.SetInverted();
    556568
    557569    //------------------------------
     
    845857        tlist2.AddToList(&taskenv2);
    846858
     859    // Check if we have an extremely bright event (remove them,
     860    // they are presumably errornous events or calibration events)
     861    tlist2.AddToList(&contbright);
     862
    847863    /*
    848864    if (fIsHiLoCalibration)
Note: See TracChangeset for help on using the changeset viewer.