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...

typedef JClonable< JClonable_t >::clone_type	clone_type

Public Member Functions
	JEvtWeightFactorMupage ()
	Default constructor.

bool	is_valid () const override final
	Check whether this formula is valid.

double	getFactor (const Evt &evt) const override final
	Get weighting factor for given event.

virtual clone_type	clone () const override
	Get clone of this object.

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 40 of file JEvtWeightFactorMupage.hh.

Member Typedef Documentation

◆ clone_type

template<class JClonable_t , class JDerived_t >

JClonable<JClonable_t>::clone_type JLANG::JClonable< JClonable_t, JDerived_t >::clone_type

inherited

Definition at line 61 of file JClonable.hh.

Member Enumeration Documentation

◆ variables

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 46 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

◆ JEvtWeightFactorMupage()

JAANET::JEvtWeightFactorMupage::JEvtWeightFactorMupage ( )

inline

Default constructor.

Definition at line 60 of file JEvtWeightFactorMupage.hh.

61 {}

Member Function Documentation

◆ is_valid()

bool JAANET::JEvtWeightFactorMupage::is_valid ( ) const

inlinefinaloverride

Check whether this formula is valid.

Returns: true if valid; else false

Definition at line 69 of file JEvtWeightFactorMupage.hh.

    {
      const int N = this->GetNdim();
 
      return N >= 0 && N <= (int) NUMBER_OF_VARIABLES;
    }

◆ getFactor()

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

inlinefinaloverride

Get weighting factor for given event.

Parameters

evt event

Returns: weighting factor

Definition at line 83 of file JEvtWeightFactorMupage.hh.

    {
      using namespace std;
      using namespace JPP;
      
      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)) {
          muons.push_back(getTrack(*i));
        }
      }
 
      if (muons.empty()) {
        THROW(JValueOutOfRange, "JEvtWeightFactorMupage::operator(): No muon for event " << evt.id << '.' << endl);
      }
 
      Double_t  vars[NUMBER_OF_VARIABLES] = { 0.0 };
 
      double    Etot = 0.0;
      JVector3D dir;
 
      for (vector<JTrack3E>::const_iterator i = muons.begin(); i != muons.end(); ++i) {
        Etot += i->getE();
        dir  += i->getDirection();
      }
        
      const JDirection3D D(dir);
      const JRotation3D  R(D);
 
      for (vector<JTrack3E>::iterator i = muons.begin(); i != muons.end(); ++i) {
        i->rotate(R);
      }
 
      vector<Trk>::const_iterator iBundle = find_if(evt.mc_trks.cbegin(), evt.mc_trks.cend(),
                                                    make_predicate(&Trk::status, TRK_ST_MUONBUNDLE));
 
      vars[MUON_MULTIPLICITY] = (double) muons.size();
      vars[TOTAL_MUON_ENERGY] = (iBundle != evt.mc_trks.cend() ? iBundle->E : Etot);
      vars[MEAN_ZENITH_ANGLE] = D.getDZ();
 
      if (this->GetNdim() > (int) LATERAL_SPREAD) {
 
        const JCircle2D circle = JCircle2D(muons.cbegin(), muons.cend()); // smallest enclosing circle
 
        vars[LATERAL_SPREAD] = circle.getRadius();
      }
          
      return this->EvalPar(&vars[0]);
    }

◆ clone()

template<class JClonable_t , class JDerived_t >

virtual clone_type JLANG::JClonable< JClonable_t, JDerived_t >::clone ( ) const

inlineoverridevirtualinherited

Get clone of this object.

Returns: pointer to newly created object

Reimplemented in JPHYSICS::JPD0Transformer_t< JArgument_t >, JPHYSICS::JPDFTransformer< 2, JArgument_t >, JPHYSICS::JPDFTransformer< 3, JArgument_t >, JPHYSICS::JPDFTransformer< 4, JArgument_t >, JPHYSICS::JPDFTransformer< 5, JArgument_t >, JPHYSICS::JPDFTransformer_t< JArgument_t >, JPHYSICS::JPDGTransformer_t< JArgument_t >, and JTOOLS::JMultiMapTransformer< N, JArgument_t >::JMultiMapDefaultTransformer.

Definition at line 69 of file JClonable.hh.

    {
      return new JDerived_t(static_cast<const JDerived_t&>(*this));
    }

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

JEvtWeightFactorMupage.hh

Public Types

Public Member Functions