AlgDeMuxBeam.h

Go to the documentation of this file.

// $Id: AlgDeMuxBeam.h,v 1.22 2004/11/01 20:12:03 brebel Exp $
//
// AlgDeMuxBeam
//
// An Algorithm class that DeMultiplexes an event and adjusts the
//weights in the PlexSEIdAltLists of the CandDigits from that event
//
// Author:  B. Rebel 11/2000

#ifndef ALGDEMUXBEAM_H
#define ALGDEMUXBEAM_H

#include "Algorithm/AlgBase.h"
#include "Validity/VldContext.h"

class CandDeMuxDigitHandle;
class DmxPlaneItr;
class DmxStatus;
class CandDeMuxDigitHandleItr;

class AlgDeMuxBeam : public AlgBase
{

        public:

        // Constructors and destructors
        AlgDeMuxBeam();
        virtual ~AlgDeMuxBeam();

        // State testing methods
        virtual void Trace(const char *c) const;

        // State changing methods
        virtual void RunAlg(AlgConfig &acd, CandHandle &ch, CandContext &cx);

        private:
                
                Float_t fFinalShowerRegionSlope;
                Float_t fFinalShowerRegionIntercept;
                Float_t fInitialMuonRegionSlope;
                Float_t fInitialMuonRegionIntercept;
                Float_t fDigitResetCut;
                Int_t fMaxBadResidual;
                Float_t fMaxResidual;
                Float_t fMaxMuonRemovalFraction;
                Int_t fMaxStripAdjustment;
                Int_t fMinMuonPlanesRequired;
                Float_t fMinMatedChargeFrac;
                Int_t fMuonPlanesInSet;
                Float_t fMuonSlopeLimit;
                Int_t fMuonStartPlaneNumber; //starting plane of a muon track
                Float_t fMuonStartPlaneZPos;
                Int_t fPlanesInSet; //number of planes in sliding window
                Int_t fStrayCut;
                Int_t fStrayPlanes;
                Int_t fNumberOfHypotheses;
                Float_t fShowerMuonTransverseDifference;
                VldContext fVldContext; //current context

                void DeMuxFirstNPlanesTest(DmxPlaneItr &planeItr, DmxPlaneItr &validPlaneItr);
                void FindShowerWindowSolution(DmxPlaneItr &windowPlaneItr, Float_t cog, 
                                                                          Int_t firstPlane, Int_t lastPlane);
                void SetFirstNPlanes(DmxPlaneItr &planeItr, Float_t slope, Float_t intercept);
                void SetMuonRegionWindow(DmxPlaneItr &planeItr, Float_t slope, Float_t intercept, bool setAll);
                void SetShowerRegionWindow(DmxPlaneItr &planeItr, Float_t slope, Float_t intercept);
                void UseShowerSlidingWindow(DmxPlaneItr &planeItr, DmxPlaneItr &validPlaneItr, 
                                                                        DmxStatus *status);
                void UseMuonSlidingWindow(DmxPlaneItr &planeItr, DmxPlaneItr &validPlaneItr, 
                                                                  DmxStatus *status);
                void UseReverseMuonSlidingWindow(DmxPlaneItr &planeItr, DmxPlaneItr &validPlaneItr, 
                                                                                 DmxStatus *status);
                void RegionLinearFit(Float_t &slope, Float_t &intercept, Float_t &chiSqr,
                                                         Float_t *x, Float_t *y, Float_t *sigma,
                                                         Float_t *residual, Int_t numPoints);
                void FindFitOverNPlanes(Float_t *x, Float_t *y, Float_t *sigma, Float_t *residual, 
                                                                Float_t &slope, Float_t &intercept, Int_t numPoints,
                                                                Int_t *planeNum);
                Float_t FindRegionMatedChargeFraction(DmxPlaneItr &validPlaneItr);
                void ReconcileShowerAndMuonRegions(DmxPlaneItr &planeItr,
                                                                                   CandDeMuxDigitHandleItr &cdhit);
   
                //depricated methods left in for now
                void DeMuxFirstNPlanes(DmxPlaneItr &planeItr, DmxPlaneItr &validPlaneItr, 
                                                           CandDeMuxDigitHandleItr &cdhit);
                Float_t FindDigitCoG(DmxPlaneItr &validPlaneItr, CandDeMuxDigitHandleItr &cdhit);
                void FindFit(DmxPlaneItr &planeItr, DmxPlaneItr &validPlaneItr, 
                                         Double_t &a1, Double_t &a2, Double_t &a3, 
                                         Double_t &a4, Double_t &bestChi2, 
                                         Int_t leverArm, Int_t *hypotheses);
                void FindFitForHypothesisSet(DmxPlaneItr &planeItr, Double_t &prevChiSq, 
                                                                         Double_t &a1, Double_t &a2, Double_t &a3, 
                                                                         Double_t &a4, Int_t* hypotheses, 
                                                                         Int_t leverArm);       
                Bool_t FindPlanesToDropFromFit(DmxPlaneItr &planeItr, Double_t a1, Double_t a2,
                                                                           Double_t a3, Double_t a4,
                                                                           Int_t leverArm);
                void FindRMSOfPlanes(DmxPlaneItr &planeItr, Int_t leverArm, Int_t *hypsToUse, 
                                                         Float_t &cog, Float_t &rms);

                ClassDef(AlgDeMuxBeam,8)                      // Beam DeMux Algorithm Class

};

#endif                                                  // ALGDEMUXBEAM_H

---