GammaFilter.h

Go to the documentation of this file.

#ifndef GW_GAMMAFILTER_H
#define GW_GAMMAFILTER_H

#ifndef ROOT_TH1F
#include <TH1F.h>
#endif

#ifndef GW_GAMMALINK_H
#include <GammaLink.h>
#endif

#ifndef GW_BASEFILTER_H
#include <BaseFilter.h>
#endif

namespace Gw {


class GammaFilter : public BaseFilter
{
private:

protected:
        TH1F *fComptonShape;    // What is the energy after a compton scattering
        TH1F *fCompton;         // histogram that is used to do compton scattering
        TH1F *fDoppler;         // doppler broadring
        TH1F *fEfficiency;      // efficiency curve
        
public:
        GammaFilter();

        GammaFilter(Axis_t xmin, Axis_t xmax, Float_t eff = 1, Float_t PT = 0.5, Float_t dop = 1.25);
        virtual ~GammaFilter();


        Bool_t IsEffective();
        

        // void SetEfficiency(TH1F &);

        
        Bool_t  InitFilter(const char *);
        

        Bool_t Compton ( GammaLink *, Float_t &);

        Float_t Doppler ( GammaLink * );

        bool IsRejected ( GammaLink * );


        virtual Int_t ApplyE(const TSeqCollection &, TBuffer &);
        virtual Int_t ApplyE(const TSeqCollection &, Float_t *);
                        
        ClassDef(GammaFilter,1); // An experimental filter to real gamma-ray cascades
};
// GammaFilter inline members //////////////////////////////////////////////////////
inline bool GammaFilter::IsRejected ( GammaLink *gam )
{ return fEfficiency->GetBinContent(fEfficiency->FindBin(gam->GetEnergy().GetValue())) < fRand->Current()->Rndm();}
inline Float_t GammaFilter::Doppler ( GammaLink *gam )
{
        Float_t e = gam->GetEnergy().GetValue();
        return fRand->Current()->Gaus(e,fDoppler->GetBinContent(fDoppler->FindBin(e)));
}
// End GammaFilter inline members //////////////////////////////////////////////////

}
#endif