JFIT::JAtmosphericMuon Class Reference

Auxiliary class to evaluate atmospheric muon hypothesis. More...

#include <JEvtToolkit.hh>

Public Member Functions
	JAtmosphericMuon ()
	Default constructor. More...
	JAtmosphericMuon (const double theta1, const double theta2)
	Constructor. More...
double	operator() (JEvt::const_iterator __begin, JEvt::const_iterator __end) const
	Test is event is atmospheric muon. More...

Public Attributes
double	dot1
double	dot2

Friends
std::istream &	operator>> (std::istream &in, JAtmosphericMuon &object)
	Read atmospheric muon analyser from input. More...
std::ostream &	operator<< (std::ostream &out, const JAtmosphericMuon &object)
	Write atmospheric muon analyser to output. More...

Detailed Description

Auxiliary class to evaluate atmospheric muon hypothesis.

The hypothesis is tested by means of the difference in quality between the best upward and best downward track.

Definition at line 624 of file JEvtToolkit.hh.

Constructor & Destructor Documentation

JFIT::JAtmosphericMuon::JAtmosphericMuon ( )

inline

Default constructor.

Definition at line 629 of file JEvtToolkit.hh.

                       :
      dot1(0.0),
      dot2(0.0)
    {}

JFIT::JAtmosphericMuon::JAtmosphericMuon ( const double  theta1, const double  theta2  )

inline

Constructor.

Parameters

theta1	upper hemisphere angle [deg]
theta2	lower hemisphere angle [deg]

Definition at line 641 of file JEvtToolkit.hh.

                                          :
      dot1(cos(theta1 * JTOOLS::PI/180.0)),
      dot2(cos(theta2 * JTOOLS::PI/180.0))
    {}

Member Function Documentation

double JFIT::JAtmosphericMuon::operator() ( JEvt::const_iterator  __begin, JEvt::const_iterator  __end  ) const

inline

Test is event is atmospheric muon.

Parameters

__begin	begin of data
__end	end of data

Returns

negative if preferably down / positive if preferably up

Definition at line 655 of file JEvtToolkit.hh.

    {
      double Qup   = 0.0;
      double Qdown = 0.0;

      for (JEvt::const_iterator i = __begin; i != __end; ++i) {

        if        (i->getDZ() >= dot1) {
            
          if (i->getQ() > Qup) {
            Qup   = i->getQ();
          }
          
        } else if (i->getDZ() <= dot2) {
          
          if (i->getQ() > Qdown) {
            Qdown = i->getQ();
          }
        }
      }

      return Qup - Qdown;
    }

Friends And Related Function Documentation

std::istream& operator>> ( std::istream &  in, JAtmosphericMuon &  object  )

friend

Read atmospheric muon analyser from input.

Parameters

in	input stream
object	atmospheric muon analyser

Returns

input stream

Definition at line 688 of file JEvtToolkit.hh.

    {
      double theta1, theta2;

      in >> theta1 >> theta2;

      object.dot1 = cos(theta1 * JTOOLS::PI/180.0);
      object.dot2 = cos(theta2 * JTOOLS::PI/180.0);

      return in;
    }

std::ostream& operator<< ( std::ostream &  out, const JAtmosphericMuon &  object  )

friend

Write atmospheric muon analyser to output.

Parameters

out	output stream
object	atmospheric muon analyser

Returns

output stream

Definition at line 708 of file JEvtToolkit.hh.

    {
      out << acos(object.dot1) * 180.0 / JTOOLS::PI << ' '
          << acos(object.dot2) * 180.0 / JTOOLS::PI;

      return out;
    }

Member Data Documentation

double JFIT::JAtmosphericMuon::dot1

Definition at line 717 of file JEvtToolkit.hh.

double JFIT::JAtmosphericMuon::dot2

Definition at line 718 of file JEvtToolkit.hh.

---

The documentation for this class was generated from the following file:

• JEvtToolkit.hh