JAANET::JEvtWeightFactorMupage Struct Reference

Implementation of reweighting factor for mupage events according to a specifiable ROOT TFormula. More...

#include <JEvtWeightFactorMupage.hh>

Inheritance diagram for JAANET::JEvtWeightFactorMupage:

Public Types
enum	variables {   MUON_MULTIPLICITY, MEAN_ZENITH_ANGLE, TOTAL_MUON_ENERGY, LATERAL_SPREAD,   NUMBER_OF_VARIABLES }
	Indices of reweighting variables for MUPAGE. More...

Public Member Functions
	JEvtWeightFactorMupage ()
	Default constructor. More...
	JEvtWeightFactorMupage (const char *name, const char *formula)
	Constructor. More...
double	operator() (const Evt &evt) const
	Get weighting factor for given event. More...

Detailed Description

Implementation of reweighting factor for mupage events according to a specifiable ROOT TFormula.

Note: The ROOT TFormula may assume any number of parameters, but should be restricted to
the physical variables listed among JEvtWeightFactorMupage::variables.
These variables may be specified within the ROOT TFormula as 'x[<index>]',
where <index> corresponds to the index of the desired physical variable within JEvtWeightFactorMupage::variables.

Definition at line 33 of file JEvtWeightFactorMupage.hh.

Member Enumeration Documentation

enum JAANET::JEvtWeightFactorMupage::variables

Indices of reweighting variables for MUPAGE.

Enumerator
MUON_MULTIPLICITY	Muon multiplicity.
MEAN_ZENITH_ANGLE	Average cosine of zenith angle.
TOTAL_MUON_ENERGY	Muon bundle total energy [GeV].
LATERAL_SPREAD	Muon bundle lateral spread [m].
NUMBER_OF_VARIABLES	Number of reweighting variables; N.B. This enum value needs to be specified last!

Definition at line 39 of file JEvtWeightFactorMupage.hh.

                   {
      MUON_MULTIPLICITY,  //!< Muon multiplicity      
      MEAN_ZENITH_ANGLE,  //!< Average cosine of zenith angle
      TOTAL_MUON_ENERGY,  //!< Muon bundle total energy         [GeV]
      LATERAL_SPREAD,     //!< Muon bundle lateral spread       [m]
      
      NUMBER_OF_VARIABLES //!< Number of reweighting variables; \n
                          //!< N.B.\ This enum value needs to be specified last!
    }; 

Constructor & Destructor Documentation

JAANET::JEvtWeightFactorMupage::JEvtWeightFactorMupage ( )

inline

Default constructor.

Definition at line 53 of file JEvtWeightFactorMupage.hh.

                             :
      TFormula()
    {}

JAANET::JEvtWeightFactorMupage::JEvtWeightFactorMupage ( const char *  name, const char *  formula  )

inline

Constructor.

Parameters

name	name
formula	formula

Definition at line 64 of file JEvtWeightFactorMupage.hh.

                                                :
      TFormula(name, formula)
    {}

Member Function Documentation

double JAANET::JEvtWeightFactorMupage::operator() ( const Evt &  evt ) const

inline

Get weighting factor for given event.

Parameters

evt

event

Returns

weighting factor

Definition at line 76 of file JEvtWeightFactorMupage.hh.

    {
      using namespace std;
      using namespace JPP;
      
      vector<Double_t> vars(NUMBER_OF_VARIABLES, 0.0);

      // Compute average zenith-angle of final state muons and convert to JTrack3E

      double Etot = 0.0;
      
      vector<JTrack3E> muons;
      
      for (vector<Trk>::const_iterator i = evt.mc_trks.cbegin(); i != evt.mc_trks.cend(); ++i) {
        
        if (is_muon(*i) && is_finalstate(*i)) {

          Etot += i->E;
          
          vars[MEAN_ZENITH_ANGLE] += i->dir.z / i->dir.len();

          muons.push_back(getTrack(*i));
        }
      }

      vector<Trk>::const_iterator iBundle = find_if(evt.mc_trks.cbegin(), evt.mc_trks.cend(),
                                                    JLANG::make_predicate(&Trk::status, TRK_ST_MUONBUNDLE));

      vars[MUON_MULTIPLICITY] = (iBundle != evt.mc_trks.cend() && iBundle->len > 0 ?
                                 iBundle->len : (Double_t) muons.size());

      vars[TOTAL_MUON_ENERGY] = (iBundle != evt.mc_trks.cend() ?
                                 iBundle->E : Etot);

      vars[MEAN_ZENITH_ANGLE] /= (Double_t) muons.size();

      if (muons.size() > 0) {

        if (this->GetNdim() > (int) LATERAL_SPREAD) {

          JCircle2D circle     = JCircle2D(muons.cbegin(), muons.cend()); // smallest enclosing circle
          vars[LATERAL_SPREAD] = (Double_t) circle.getRadius();
        }
          
        return this->DoEval(&vars[0]);

      } else {
        
        THROW(JNullPointerException, "JEvtWeightFactorMupage::operator(): No muon for event " << evt.id << '.' << endl);        
      }
    }

---

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

• JEvtWeightFactorMupage.hh