CosmicTag.cxx

Go to the documentation of this file.

//
// \brief Definition of CosmicTag analysis object
//
// \author lockwitz@fnal.gov

#include "lardataobj/AnalysisBase/CosmicTag.h"

#include <iomanip>
#include <iostream>

namespace anab {

  //----------------------------------------------------------------------
  CosmicTag::CosmicTag()
    : endPt1(0), endPt2(0), fCosmicScore(0), fCosmicType(CosmicTagID_t::kUnknown)
  {}

  //----------------------------------------------------------------------
  CosmicTag::CosmicTag(float f)
    : endPt1(0), endPt2(0), fCosmicScore(f), fCosmicType(CosmicTagID_t::kUnknown)
  {}

  //----------------------------------------------------------------------
  CosmicTag::CosmicTag(std::vector<float> ePt1,
                       std::vector<float> ePt2,
                       //                      double flashTime, // should go
                       float cScore,
                       CosmicTagID_t cType)
    : endPt1(ePt1), endPt2(ePt2), fCosmicScore(cScore), fCosmicType(cType)
  {}

  //----------------------------------------------------------------------
  // ostream operator.
  //
  std::ostream& operator<<(std::ostream& o, CosmicTag const& a)
  {
    o //<< "Flash Time     : "     <<  a.flashTime     // should go
      << "\n Cosmic Score     : " << a.fCosmicScore << "\n Cosmic Type     : " << a.fCosmicType
      << "\n End Point 1: " << std::endl;
    for (size_t i = 0; i < a.endPt1.size(); i++)
      o << a.endPt1[i] << ", ";
    o << std::endl;
    o << "\n End Point 2: " << std::endl;
    for (size_t i = 0; i < a.endPt2.size(); i++)
      o << a.endPt2[i] << ", ";
    o << std::endl;

    return o;
  }

  float CosmicTag::getXInteraction(float oldX,
                                   float xDrift,
                                   int tSample,
                                   float realTime,
                                   int tick) const
  {
    // fix this
    int t = 3 * tSample;
    if ((tick > 0 * tSample) && (tick < 1 * tSample)) t = 0 * tSample;
    if ((tick > 1 * tSample) && (tick < 2 * tSample)) t = 1 * tSample;
    if ((tick > 2 * tSample) && (tick < 3 * tSample)) t = 2 * tSample;

    float newX = oldX + (xDrift / tSample) * (t - realTime);
    return newX;
  }

}