lightana::FlashGeoBarycenter Class Reference

Inheritance diagram for lightana::FlashGeoBarycenter:

Classes
struct	Config

Public Member Functions
	FlashGeoBarycenter (art::ToolConfigTable< Config > const &config)
void	GetFlashLocation (std::vector< double > pePerOpChannel, double &Ycenter, double &Zcenter, double &Ywidth, double &Zwidth) override

Private Attributes
unsigned int	fWeightExp

Additional Inherited Members
Private Member Functions inherited from lightana::FlashGeoBase
virtual	~FlashGeoBase () noexcept=default

Detailed Description

Definition at line 20 of file FlashGeoBarycenter_tool.cc.

Constructor & Destructor Documentation

lightana::FlashGeoBarycenter::FlashGeoBarycenter ( art::ToolConfigTable< Config > const &  config )

explicit

Definition at line 51 of file FlashGeoBarycenter_tool.cc.

    : fWeightExp{ config().WeightExp() }
  {
  }

Member Function Documentation

void lightana::FlashGeoBarycenter::GetFlashLocation ( std::vector< double >  pePerOpChannel, double &  Ycenter, double &  Zcenter, double &  Ywidth, double &  Zwidth  )

overridevirtual

Implements lightana::FlashGeoBase.

Definition at line 57 of file FlashGeoBarycenter_tool.cc.

  {

    // Reset variables
    Ycenter = Zcenter = 0.;
    Ywidth  = Zwidth  = -999.;
    double totalPE = 0.;
    double sumy = 0., sumz = 0., sumy2 = 0., sumz2 = 0.;
    double weight =1.;
    for (unsigned int opch = 0; opch < pePerOpChannel.size(); opch++) {
      // Get physical detector location for this opChannel
      double PMTxyz[3];
      ::lightana::OpDetCenterFromOpChannel(opch, PMTxyz);

      // Get weight for this channel
      if(fWeightExp==1) weight = pePerOpChannel[opch];
      else if(fWeightExp==2) weight = pePerOpChannel[opch]*pePerOpChannel[opch];
      else weight = std::pow(pePerOpChannel[opch], fWeightExp);

      // Add up the position, weighting with PEs
      sumy    += weight*PMTxyz[1];
      sumy2   += weight*PMTxyz[1]*PMTxyz[1];
      sumz    += weight*PMTxyz[2];
      sumz2   += weight*PMTxyz[2]*PMTxyz[2];
      totalPE += weight;
    }

    Ycenter = sumy/totalPE;
    Zcenter = sumz/totalPE;

    // This is just sqrt(<x^2> - <x>^2)
    if ( (sumy2*totalPE - sumy*sumy) > 0. )
      Ywidth = std::sqrt(sumy2*totalPE - sumy*sumy)/totalPE;

    if ( (sumz2*totalPE - sumz*sumz) > 0. )
      Zwidth = std::sqrt(sumz2*totalPE - sumz*sumz)/totalPE;
  }

Member Data Documentation

unsigned int lightana::FlashGeoBarycenter::fWeightExp

private

Definition at line 46 of file FlashGeoBarycenter_tool.cc.

---

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

• FlashGeoBarycenter_tool.cc