JRECONSTRUCTION::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 710 of file JReconstruction/JEvtToolkit.hh.

Constructor & Destructor Documentation

JRECONSTRUCTION::JAtmosphericMuon::JAtmosphericMuon ( )

inline

Default constructor.

Definition at line 715 of file JReconstruction/JEvtToolkit.hh.

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

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

inline

Constructor.

Parameters

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

Definition at line 727 of file JReconstruction/JEvtToolkit.hh.

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

Member Function Documentation

double JRECONSTRUCTION::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 741 of file JReconstruction/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 774 of file JReconstruction/JEvtToolkit.hh.

    {
      double theta1, theta2;

      in >> theta1 >> theta2;

      object.dot1 = cos(theta1 * JMATH::PI/180.0);
      object.dot2 = cos(theta2 * JMATH::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 794 of file JReconstruction/JEvtToolkit.hh.

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

      return out;
    }

Member Data Documentation

double JRECONSTRUCTION::JAtmosphericMuon::dot1

Definition at line 803 of file JReconstruction/JEvtToolkit.hh.

double JRECONSTRUCTION::JAtmosphericMuon::dot2

Definition at line 804 of file JReconstruction/JEvtToolkit.hh.

---

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

• JReconstruction/JEvtToolkit.hh