RockEvent.cxx

Go to the documentation of this file.

// RockEvent.cxx
// A reduced version of NuEvent for storing the rock events

#include <iostream>
#include <cmath>
#include "MessageService/MsgService.h"
#include "RockMuons/RockEvent.h"
#include "NtupleUtils/NuEvent.h"
#include "NtupleUtils/NuCutImps.h"
#include "NtupleUtils/NuReco.h"

#include "TMath.h"

CVSID("");

//......................................................................
RockEvent::RockEvent()
{
  this->Reset();
}

//......................................................................

RockEvent::RockEvent(const NuEvent& nu) : neuEnMCBin(-1), binNumber(-1)
{
  this->Fill(nu);
}

//......................................................................

void RockEvent::Reset(void)
{
  // Create a new, blank NuEvent so that the reset values will always be the same
  this->Fill(NuEvent());

  this->neuEnMCBin = -1;
  this->binNumber    = -1;
}

//......................................................................

void RockEvent::Fill(NuEvent nu, const bool dohoveshift)
{
  // Apply geometry corrections.  This is not the ideal place to do it, but
  // Justin and I (Matt) decided that this was the best way that is feasible
  // at the moment.  For the next analysis round, this should be handled
  // in NuReco::ApplyReweights in the same way as energy shifts are.
  // Make sure this is done only for monte carlo and fake data (which is,
  // of course, monte carlo and so still gets SimFlag::kMC)
  if(dohoveshift && nu.simFlag == SimFlag::kMC)
    ShiftHoveStrip(nu);

  // Event variables
  this->rw = nu.rw;

  // Binning/histogram selecting variables

  this->regionTrkVtx = nu.regionTrkVtx;

  this->edgeRegionTrkVtx = nu.edgeRegionTrkVtx;
  this->stripHoveNumTrkVtx = nu.stripHoveNumTrkVtx;
  this->stripTrkBegPerpFlag = nu.stripTrkBegPerpFlag;
  this->parallelStripTrkVtx = nu.parallelStripTrkVtx;
  this->planeTrkBegu = nu.planeTrkBegu;
  this->planeTrkBegv = nu.planeTrkBegv;

  // edgemid = Edge Medium Identification
  // Possible media are:
  // 0 = rock-like
  // 1 = detector-like
  //
  // this has a value whether or not a track vertex is near an edge
  // Vertices far from an edge should all be labeled detector-like,
  // although near the ends of each supermodule this is not reliable
  this->edgemid = NuCutImps::EdgeMID(nu.edgeRegionTrkVtx, nu.stripHoveNumTrkVtx,
                                     nu.parallelStripTrkVtx);

  // "Best" track energy of the "primary" track, where "best" and "primary"
  // are defined upstream.  Force the correct definition of containment
  // for choosing range or curvature.
  const float trkEn = (nu.containmentFlagCC0250Std == 1) ? nu.trkEnRange : nu.trkEnCurv;
  if(fabs(trkEn-nu.trkEn) > 1e-3){
    MAXMSG("RockEvent", Msg::kWarning, 10)
      << "Adjusting trkEn because it doesn't agree with "
      << "containmentFlagCC0250Std. Are you using old DSTs?" << endl;
  }
  this->trkEn = trkEn;

  this->dirCosNu = nu.dirCosNu;

  // Angle of the track with respect to the beam direction, in degrees.
  // It's easier to bin sensibly than dirCosNu
  this->dang = TMath::ACos(nu.dirCosNu)*180/TMath::Pi();

  // Truth variables
  this->neuEnMC = nu.neuEnMC;
  this->inu = nu.inu;
  this->iaction = nu.iaction;
  this->trkStartEnMC = nu.trkStartEnMC;
}
//......................................................................

std::ostream& operator<<(std::ostream& out,const RockEvent& val)
{
  out << "Rock Event  Energy Trk: " << val.trkEn
      << " MC: " << val.neuEnMC << "\n"
         "            Region: " << val.regionTrkVtx << " rw: " << val.rw
      << "            inu:  " << val.inu;
  return out;
}

---