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 "messagefacility/MessageLogger/MessageLogger.h"

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;
  }
  
  
}