/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/Pi0RecAlg/Pi0RecAlg-00-00-06/Pi0RecAlg/Criteria.h

Go to the documentation of this file.

#ifndef CRITERIA__H
#define CRITERAI__H

#include "GaudiKernel/SmartDataPtr.h"

#include "EventModel/EventHeader.h"
#include "EvtRecEvent/EvtRecEvent.h"
#include "EvtRecEvent/EvtRecTrack.h"

HepLorentzVector getP4(RecEmcShower* gTrk)
{

   double eraw = gTrk->energy();
   double phi =  gTrk->phi();
   double the =  gTrk->theta();

   return HepLorentzVector( eraw * sin(the) * cos(phi),
                            eraw * sin(the) * sin(phi),
                            eraw * cos(the),
                            eraw );
}
bool isGoodGamma(EvtRecTrack* gtrk)
{

}
struct GammaPair : public std::pair<EvtRecTrack*,EvtRecTrack*>
{
        GammaPair(EvtRecTrack* g1, EvtRecTrack* g2, double m/*, double s, double c*/)
                :First(g1), Second(g2), inv_m(m)/*, chisq(s), cos_theta(c)*/{;}
        EvtRecTrack* First;
        EvtRecTrack* Second;
        double inv_m;
//      double chisq;
//      double cos_theta;
};
/*struct less_chisq: public binary_function <GammaPair, GammaPair, bool>
{
        bool operator()(const GammaPair& lgp, const GammaPair& rgp)
        {
                return lgp.chisq < rgp.chisq;
        }
};*/
struct high_momentum: public binary_function <GammaPair, GammaPair, bool>
{
        bool operator()(const GammaPair& lgp, const GammaPair& rgp)
        {
                double p1 = (getP4(lgp.First->emcShower())+ getP4(lgp.Second->emcShower())).rho();
                double p2 = (getP4(rgp.First->emcShower())+ getP4(rgp.Second->emcShower())).rho();
                return p1 > p2;
        }
};
struct Criteria
{
        Criteria(double inv_min = 0.10, double inv_max = 0.18):_inv_min(inv_min), _inv_max(inv_max){;}
        double _inv_min;
        double _inv_max;
        virtual bool check(const GammaPair& gp) const 
        {
                //                std::cout<<"gp.inv_m"<<gp.inv_m<<" inv_min"<<_inv_min<<std::endl;
                if( gp.inv_m> _inv_min && gp.inv_m < _inv_max) return true;
                else return false;
        }
};
/*struct KFitCriteria : public Criteria
{
        KFitCriteria(double chisq_cut = 20., double costheta_cut = 0.98 ): _chisq_cut(chisq_cut),_costheta_cut(costheta_cut){;}
        double _chisq_cut;
        double _costheta_cut;
        virtual bool check(const GammaPair& gp) const 
        {
                //                std::cout<<"gp.chisq"<<gp.chisq<<" chisq_cut"<<_chisq_cut<<std::endl;
                if(gp.chisq > _chisq_cut) return false;
                if(fabs(gp.cos_theta) > _costheta_cut) return false;
                return true;
        }
};*/
/*struct helper : public Criteria
{
        helper(){;}
        double mass(const GammaPair& gp) const{
                EvtRecTrack* g1Trk = gp.First;
                EvtRecTrack* g2Trk = gp.Second;

                RecEmcShower* g1Shower = g1Trk->emcShower();
                RecEmcShower* g2Shower = g2Trk->emcShower();
                return (getP4(g1Shower)+getP4(g2Shower)).m();
        }
        double angle(const GammaPair& gp) const{
                EvtRecTrack* g1Trk =  gp.First;
                EvtRecTrack* g2Trk =  gp.Second;

                RecEmcShower* g1Shower = g1Trk->emcShower();
                RecEmcShower* g2Shower = g2Trk->emcShower();
                return getP4(g1Shower).angle(getP4(g2Shower).vect());;
        }
};*/
#endif

---