Hitdumper Class Reference

Inheritance diagram for Hitdumper:

Public Member Functions
	Hitdumper (fhicl::ParameterSet const &p)
virtual	~Hitdumper ()
void	beginJob () override
void	reconfigure (fhicl::ParameterSet const &pset)
void	analyze (const art::Event &evt) override
virtual void	beginSubRun (art::SubRun const &sr) override

Private Member Functions
void	ResetVars ()
	Resets the variables that are saved to the TTree. More...
void	ResetWireHitsVars (int n)
	Resets wire hits tree variables. More...
void	ResetCRTStripsVars (int n)
	Resets crt strips tree variables. More...
void	ResetCRTCustomTracksVars (int n)
	Resets custom crt tracks tree variables. More...
void	ResetCRTTracksVars (int n)
	Resets crt tracks tree variables. More...
void	ResetCRTHitsVars (int n)
	Resets crt hits tree variables. More...
void	ResetOpHitsVars (int n)
	Resets optical hits tree variables. More...
void	ResetPmtTriggerVars (int n)
	Resets pmt hardware trigger variables. More...
void	ResetPmtSoftTriggerVars ()
	Rests pmt software trigger variables. More...
void	ResetCrtSoftTriggerVars ()
	Resets crt software trigger variables. More...
void	ResetMuonTracksVars (int n)
	Resets crossing muon tracks tree variables. More...
void	ResetMuonHitVars (int n)
	Resets crossing muon hit tree variables. More...
void	ResizeMCNeutrino (int nNeutrinos)
	Resize the data structure for MCNeutrino particles. More...
void	ResizeGenie (int nPrimaries)
	Resize the data structure for Genie primaries. More...

Private Attributes
opdet::sbndPDMapAlg	_pd_map
TTree *	fTree
int	_run
	The run number. More...
int	_subrun
	The subrun number. More...
int	_event
	The event number. More...
double	_evttime
	The event time. More...
int	_t0
	The t0. More...
int	_nhits
	Number of reco hits in the event. More...
std::vector< int >	_hit_cryostat
	Cryostat where the hit belongs to. More...
std::vector< int >	_hit_tpc
	TPC where the hit belongs to. More...
std::vector< int >	_hit_plane
	Plane where the hit belongs to. More...
std::vector< int >	_hit_wire
	Wire where the hit belongs to. More...
std::vector< int >	_hit_channel
	Channel where the hit belongs to. More...
std::vector< double >	_hit_peakT
	Hit peak time. More...
std::vector< double >	_hit_charge
	Hit charge. More...
std::vector< double >	_hit_ph
	Hit ph? More...
std::vector< double >	_hit_width
	Hit width. More...
std::vector< double >	_hit_full_integral
	Hit charge integral. More...
std::vector< int >	_waveform_number
	Number for each waveform, to allow for searching. More...
std::vector< double >	_adc_on_wire
	ADC on wire to draw waveform. More...
std::vector< int >	_time_for_waveform
	Time for waveform to plot. More...
int	_adc_count
	Used for plotting waveforms. More...
std::vector< int >	_waveform_integral
	Used to see progression of the waveform integral. More...
std::vector< int >	_adc_count_in_waveform
	Used to view all waveforms on a hitplane together. More...
int	_nstrips
	Number of CRT strips. More...
std::vector< int >	_crt_plane
	CRT plane. More...
std::vector< int >	_crt_module
	CRT module. More...
std::vector< int >	_crt_strip
	CRT strip. More...
std::vector< int >	_crt_orient
	CRT orientation (0 for y (horizontal) and 1 for x (vertical)) More...
std::vector< double >	_crt_time
	CRT time. More...
std::vector< double >	_crt_adc
	CRT adc. More...
std::vector< double >	_crt_pos_x
	CRT position X. More...
std::vector< double >	_crt_pos_y
	CRT position Y. More...
std::vector< double >	_crt_pos_z
	CRT position Z. More...
int	_nctrks
	Number of created CRT tracks. More...
std::vector< double >	_ctrk_x1
	CRT track x1. More...
std::vector< double >	_ctrk_y1
	CRT track y1. More...
std::vector< double >	_ctrk_z1
	CRT track z1. More...
std::vector< double >	_ctrk_t1
	CRT track t1. More...
std::vector< double >	_ctrk_adc1
	CRT track adc1. More...
std::vector< int >	_ctrk_mod1x
	CRT track mod2x. More...
std::vector< double >	_ctrk_x2
	CRT track x2. More...
std::vector< double >	_ctrk_y2
	CRT track y2. More...
std::vector< double >	_ctrk_z2
	CRT track z2. More...
std::vector< double >	_ctrk_t2
	CRT track t2. More...
std::vector< double >	_ctrk_adc2
	CRT track adc2. More...
std::vector< int >	_ctrk_mod2x
	CRT track mod2x. More...
int	_nchits
	Number of CRT hits. More...
std::vector< double >	_chit_x
	CRT hit x. More...
std::vector< double >	_chit_y
	CRT hit y. More...
std::vector< double >	_chit_z
	CRT hit z. More...
std::vector< double >	_chit_time
	CRT hit time. More...
std::vector< int >	_chit_plane
	CRT hit plane. More...
int	_ncts
	Number of CRT tracks. More...
std::vector< double >	_ct_time
	CRT track time. More...
std::vector< double >	_ct_pes
	CRT track PEs. More...
std::vector< double >	_ct_x1
	CRT track x1. More...
std::vector< double >	_ct_y1
	CRT track y1. More...
std::vector< double >	_ct_z1
	CRT track z1. More...
std::vector< double >	_ct_x2
	CRT track x2. More...
std::vector< double >	_ct_y2
	CRT track y2. More...
std::vector< double >	_ct_z2
	CRT track z2. More...
int	_nophits
	Number of Optical Hits. More...
std::vector< int >	_ophit_opch
	OpChannel of the optical hit. More...
std::vector< int >	_ophit_opdet
	OpDet of the optical hit. More...
std::vector< double >	_ophit_peakT
	Peak time of the optical hit. More...
std::vector< double >	_ophit_width
	Width of the optical hit. More...
std::vector< double >	_ophit_area
	Area of the optical hit. More...
std::vector< double >	_ophit_amplitude
	Amplitude of the optical hit. More...
std::vector< double >	_ophit_pe
	PEs of the optical hit. More...
std::vector< double >	_ophit_opdet_x
	OpDet X coordinate of the optical hit. More...
std::vector< double >	_ophit_opdet_y
	OpDet Y coordinate of the optical hit. More...
std::vector< double >	_ophit_opdet_z
	OpDet Z coordinate of the optical hit. More...
std::vector< int >	_ophit_opdet_type
	OpDet tyoe of the optical hit. More...
std::vector< int >	_pmtTrigger_npmtshigh
	number of pmt pairs above threshold, index = time during trigger window (usually beam spill) More...
int	_pmtTrigger_maxpassed
	maximum number of pmt pairs above threshold during trigger window (usually beam spill) More...
bool	_pmtSoftTrigger_foundBeamTrigger
int	_pmtSoftTrigger_tts
	Whether the beam spill was found or not. More...
double	_pmtSoftTrigger_promptPE
	Trigger Time Stamp (TTS), ns (relative to start of beam spill) More...
double	_pmtSoftTrigger_prelimPE
	Total photoelectron count 100 ns after the TTS. More...
int	_pmtSoftTrigger_nAboveThreshold
	Total photoelectron count before the TTS, during the beam spill. More...
int	_crtSoftTrigger_hitsperplane [7]
	number of individual PMTs above ADC threshold (fcl) during the beam spill More...
int	_nmuontrks
	number of muon tracks More...
std::vector< double >	_muontrk_t0
	t0 (time of interaction) More...
std::vector< float >	_muontrk_x1
	x coordinate closer to anode More...
std::vector< float >	_muontrk_y1
	y coordinate closer to anode More...
std::vector< float >	_muontrk_z1
	z coordinate closer to anode More...
std::vector< float >	_muontrk_x2
	x coordinate closer to cathode More...
std::vector< float >	_muontrk_y2
	y coordinate closer to cathode More...
std::vector< float >	_muontrk_z2
	z coordinate closer to cathode More...
std::vector< float >	_muontrk_theta_xz
	theta_xz trajectory angle More...
std::vector< float >	_muontrk_theta_yz
	theta_yz trajectory angle More...
std::vector< int >	_muontrk_tpc
	tpc that muon is located in More...
std::vector< int >	_muontrk_type
	type of muon track More...
int	_nmhits
	Number of muon collection hits per track. More...
std::vector< int >	_mhit_trk
	Track number that the hit belongs to. More...
std::vector< int >	_mhit_tpc
	TPC where the hit belongs to. More...
std::vector< int >	_mhit_wire
	Wire where the hit belongs to. More...
std::vector< int >	_mhit_channel
	Channel where the hit belongs to. More...
std::vector< double >	_mhit_peakT
	Hit peak time. More...
std::vector< double >	_mhit_charge
	Hit charge. More...
size_t	MaxMCNeutrinos
Int_t	mcevts_truth
	! The number of MCNeutrinos there is currently room for More...
std::vector< Int_t >	nuScatterCode_truth
	Scattering code given by Genie for each neutrino. More...
std::vector< Int_t >	nuID_truth
	Unique ID of each true neutrino. More...
std::vector< Int_t >	nuPDG_truth
	neutrino PDG code More...
std::vector< Int_t >	ccnc_truth
	0=CC 1=NC More...
std::vector< Int_t >	mode_truth
	0=QE/El, 1=RES, 2=DIS, 3=Coherent production More...
std::vector< Float_t >	enu_truth
	true neutrino energy More...
std::vector< Float_t >	Q2_truth
	Momentum transfer squared. More...
std::vector< Float_t >	W_truth
	hadronic invariant mass More...
std::vector< Int_t >	hitnuc_truth
	hit nucleon More...
std::vector< Float_t >	nuvtxx_truth
	neutrino vertex x More...
std::vector< Float_t >	nuvtxy_truth
	neutrino vertex y More...
std::vector< Float_t >	nuvtxz_truth
	neutrino vertex z More...
std::vector< Float_t >	nu_dcosx_truth
	neutrino dcos x More...
std::vector< Float_t >	nu_dcosy_truth
	neutrino dcos y More...
std::vector< Float_t >	nu_dcosz_truth
	neutrino dcos z More...
std::vector< Float_t >	lep_mom_truth
	lepton momentum More...
std::vector< Float_t >	lep_dcosx_truth
	lepton dcos x More...
std::vector< Float_t >	lep_dcosy_truth
	lepton dcos y More...
std::vector< Float_t >	lep_dcosz_truth
	lepton dcos z More...
std::vector< Float_t >	tpx_flux
	Px of parent particle leaving BNB target. More...
std::vector< Float_t >	tpy_flux
	Py of parent particle leaving BNB target. More...
std::vector< Float_t >	tpz_flux
	Pz of parent particle leaving BNB target. More...
std::vector< Int_t >	tptype_flux
	Type of parent particle leaving BNB target. More...
size_t	MaxGeniePrimaries = 0
Int_t	genie_no_primaries
std::vector< Int_t >	genie_primaries_pdg
std::vector< Float_t >	genie_Eng
std::vector< Float_t >	genie_Px
std::vector< Float_t >	genie_Py
std::vector< Float_t >	genie_Pz
std::vector< Float_t >	genie_P
std::vector< Int_t >	genie_status_code
std::vector< Float_t >	genie_mass
std::vector< Int_t >	genie_trackID
std::vector< Int_t >	genie_ND
std::vector< Int_t >	genie_mother
TTree *	_sr_tree
	A tree filled per subrun (for POT accounting) More...
int	_sr_run
int	_sr_subrun
double	_sr_begintime
double	_sr_endtime
double	_sr_pot
	POTs in each subrun. More...
int	_max_hits
	maximum number of hits (to be set via fcl) More...
int	_max_ophits
	maximum number of hits (to be set via fcl) More...
int	_max_samples
	maximum number of samples (to be set via fcl) More...
int	_max_chits
	maximum number of CRT hits (to be set via fcl) More...
int	_max_nctrks
	maximum number of CRT tracks (to be set via fcl) More...
std::string	fHitsModuleLabel
	Label for Hit dataproduct (to be set via fcl) More...
std::string	fLArG4ModuleLabel
	Label for LArG4 dataproduct (to be set via fcl) More...
std::string	fCRTStripModuleLabel
	Label for CRTStrip dataproduct (to be set via fcl) More...
std::string	fCRTHitModuleLabel
	Label for CRTHit dataproduct (to be set via fcl) More...
std::string	fCRTTrackModuleLabel
	Label for CRTTrack dataproduct (to be set via fcl) More...
std::string	fpmtTriggerModuleLabel
	Label for pmtTrigger dataproduct (to be set vis fcl) More...
std::string	fpmtSoftTriggerModuleLabel
	Label for pmt software trigger data product (to be set via fcl) More...
std::string	fcrtSoftTriggerModuleLabel
	Label for crt software trigger data product (to be set via fcl) More...
std::string	fMuonTrackModuleLabel
	Label for MuonTrack dataproduct (to be set via fcl) More...
std::string	fDigitModuleLabel
	Label for digitizer (to be set via fcl) More...
std::string	fGenieGenModuleLabel
	Label for Genie dataproduct (to be set via fcl) More...
std::vector< std::string >	fOpHitsModuleLabels
	Labels for OpHit dataproducts (to be set via fcl) More...
bool	fkeepCRThits
	Keep the CRT hits (to be set via fcl) More...
bool	fkeepCRTstrips
	Keep the CRT strips (to be set via fcl) More...
bool	fmakeCRTtracks
	Make the CRT tracks (to be set via fcl) More...
bool	freadCRTtracks
	Keep the CRT tracks (to be set via fcl) More...
bool	freadOpHits
	Add OpHits to output (to be set via fcl) More...
bool	freadMuonTracks
	Add MuonTracks to output (to be set via fcl) More...
bool	freadMuonHits
	Add MuonTrack hits to output(to be set via fcl) More...
bool	freadTruth
	Add Truth info to output (to be set via fcl) More...
bool	freadpmtTrigger
	Add pmt hardware trigger info to output (to be set via fcl) More...
bool	freadpmtSoftTrigger
	Add pmt software trigger info to output (to be set via fcl) More...
bool	freadcrtSoftTrigger
	Add crt software trigger info to output (to be set via fcl) More...
bool	fsavePOTInfo
	Add POT info to output (to be set via fcl) More...
bool	fcheckTransparency
	Checks for wire transprency (to be set via fcl) More...
bool	fUncompressWithPed
	Uncompresses the waveforms if true (to be set via fcl) More...
int	fWindow
bool	fSkipInd
	If true, induction planes are not saved (to be set via fcl) More...
std::vector< int >	fKeepTaggerTypes = {0, 1, 2, 3, 4, 5, 6}
	Taggers to keep (to be set via fcl) More...
sbnd::CRTHitRecoAlg	hitAlg
geo::GeometryCore const *	fGeometryService
art::ServiceHandle < geo::AuxDetGeometry >	fAuxDetGeoService
const geo::AuxDetGeometry *	fAuxDetGeo
const geo::AuxDetGeometryCore *	fAuxDetGeoCore

Detailed Description

Definition at line 99 of file HitDumper_module.cc.

Constructor & Destructor Documentation

Hitdumper::Hitdumper ( fhicl::ParameterSet const &  p )

explicit

Definition at line 385 of file HitDumper_module.cc.

: EDAnalyzer(pset)
{

  fGeometryService = lar::providerFrom<geo::Geometry>();
  // fDetectorClocks = lar::providerFrom<detinfo::DetectorClocksService>();
  // fDetectorProperties = lar::providerFrom<detinfo::DetectorPropertiesService>();
  // fTrigClock = fDetectorClocks->TriggerClock();
  fAuxDetGeo = &(*fAuxDetGeoService);
  fAuxDetGeoCore = fAuxDetGeo->GetProviderPtr();

  // Read in the parameters from the .fcl file.
  this->reconfigure(pset);
}

Hitdumper::~Hitdumper ( )

virtual

Definition at line 444 of file HitDumper_module.cc.

{
  // Clean up dynamic memory and other resources here.
}

Member Function Documentation

void Hitdumper::analyze ( const art::Event &  evt )

override

Definition at line 450 of file HitDumper_module.cc.

{
  // Reset the TTree variables
  ResetVars();

  _run = evt.run();
  _subrun = evt.subRun();
  _event = evt.id().event();

  _t0 = 0.;
  // t0 = detprop->TriggerOffset();  // units of TPC ticks

  //
  // Hits
  //
  art::Handle<std::vector<recob::Hit>> hitListHandle;
  std::vector<art::Ptr<recob::Hit>> hitlist;
  if (evt.getByLabel(fHitsModuleLabel,hitListHandle)) {
    art::fill_ptr_vector(hitlist, hitListHandle);
    _nhits = hitlist.size();

    // Calculate how many hits we will save if skipping the induction planes
    if (fSkipInd) {
      _nhits = 0;
      for (auto h : hitlist) {
        if (h->WireID().Plane == 2) {
          _nhits++;
        }
      }
    }
  }
  else {
    std::cout << "Failed to get recob::Hit data product." << std::endl;
    _nhits = 0;
  }

  if (_nhits > _max_hits) {
    std::cout << "Available hits are " << _nhits
              << ", which is above the maximum number allowed to store." << std::endl;
    std::cout << "Will only store " << _max_hits << "hits." << std::endl;
    _nhits = _max_hits;
  }

  ResetWireHitsVars(_nhits);

  size_t counter = 0;
  for (size_t i = 0; i < hitlist.size(); ++i) {
    geo::WireID wireid = hitlist[i]->WireID();
    if (fSkipInd && wireid.Plane != 2) {
      continue;
    }

    _hit_cryostat[counter] = wireid.Cryostat;
    _hit_tpc[counter] = wireid.TPC;
    _hit_plane[counter] = wireid.Plane;
    _hit_wire[counter] = wireid.Wire;
    _hit_channel[counter] = hitlist[i]->Channel();
    // peak time needs plane dependent offset correction applied.
    _hit_peakT[counter] = hitlist[i]->PeakTime();
    _hit_charge[counter] = hitlist[i]->Integral();
    _hit_ph[counter] = hitlist[i]->PeakAmplitude();
    _hit_width[counter] = hitlist[i]->RMS();
    counter ++;
  }

  //
  // CRT strips
  //
  int _nstr = 0;
  art::Handle<std::vector<sbnd::crt::CRTData> > crtStripListHandle;
  std::vector<art::Ptr<sbnd::crt::CRTData> > striplist;
  // art::Handle< std::vector<sbnd::crt::CRTData> > crtStripListHandle;
  // std::vector< art::Ptr<sbnd::crt::CRTData> > striplist;
  if (evt.getByLabel(fCRTStripModuleLabel, crtStripListHandle))  {
    art::fill_ptr_vector(striplist, crtStripListHandle);
    _nstr = striplist.size();
  } else {
    std::cout << "Failed to get sbnd::crt::CRTData data product." << std::endl;
  }

  int ns = 0;
  if (_nstr > _max_chits) _nstr = _max_chits;
  // strips are always in pairs, one entry for each sipm (2 sipms per strip)

  ResetCRTStripsVars(_nstr);

  for (int i = 0; i < _nstr; i += 2){
    uint32_t chan = striplist[i]->Channel();

    //    std::pair<std::string,unsigned> tagger = CRTHitRecoAlg::ChannelToTagger(chan);
    std::pair<std::string,unsigned> tagger = hitAlg.ChannelToTagger(chan);
    sbnd::CRTPlane ip = sbnd::CRTCommonUtils::GetPlaneIndex(tagger.first);

    bool keep_tagger = false;
    for (auto t : fKeepTaggerTypes) {
      if (ip == t) {
        keep_tagger = true;
      }
    }
    // std::cout << "Tagger name " << tagger.first << ", ip " << ip << ", kept? " << (keep_tagger ? "yes" : "no") << std::endl;

    if (ip != sbnd::kCRTNotDefined && keep_tagger) {

      uint32_t ttime = striplist[i]->T0();
      float ctime = (int)ttime * 0.001; // convert form ns to us
      // recover simulation bug where neg times are sotred as unsigned integers
      // if (ttime > MAX_INT) ctime = 0.001 * (ttime - MAX_INT * 2);
      if (ctime < 1600. && ctime > -1400.) {
        uint32_t adc1 = striplist[i]->ADC();
        uint32_t adc2 = striplist[i+1]->ADC();
        if (adc1 > 4095) adc1 = 4095;
        if (adc2 > 4095) adc2 = 4095;
        //    std::cout << tagger.first << " " << tagger.second << std::endl;
        //    int sipm = chan & 1;  // 0 or 1
        int strip = (chan >> 1) & 15;
        int module = (chan>> 5);
        //
        std::string name = fGeometryService->AuxDet(module).TotalVolume()->GetName();
        TVector3 center = fAuxDetGeoCore->AuxDetChannelToPosition(2*strip, name);
        _crt_plane.push_back(ip);
        _crt_module.push_back(module);
        _crt_strip.push_back(strip);
        _crt_orient.push_back(tagger.second);
        _crt_time.push_back(ctime);
        _crt_adc.push_back(adc1 + adc2 - 127.2); // -127.2/131.9 correct for gain and 2*ped to get pe
        _crt_pos_x.push_back(center.X());
        _crt_pos_y.push_back(center.Y());
        _crt_pos_z.push_back(center.Z());
        ns++;
      }
    }
  }
  _nstrips = ns;



  //
  // CRT Custom Tracks
  //
  ResetCRTCustomTracksVars(_nstrips);
  _nctrks = 0;
  if (fmakeCRTtracks) {
    int ntr = 0;
    int iflag[1000] = {0};
    for (int i = 0; i < (ns - 1); ++i) {
      if (iflag[i] == 0) {
        iflag[i] = 1;
        float plane1x = 0, plane2x = 0;
        float plane1y = 0, plane2y = 0;
        float plane1tx = 0, plane2tx = 0;
        float plane1ty = 0, plane2ty = 0;
        float plane1xm = -1, plane2xm = -1;
        float plane1ym = -1, plane2ym = -1;
        int  nh1x = 0, nh2x = 0;
        int  nh1y = 0, nh2y = 0;
        float adc1x = 0, adc2x = 0;
        float adc1y = 0, adc2y = 0;
        if (_crt_plane[i] == sbnd::kCRTFaceSouth) { // 1
          if (_crt_orient[i] == kCRTVertical && _crt_adc[i] > 500) { // < 500 hardcoded
            if (nh1x == 0 || (_crt_module[i] == plane1xm)) {
              nh1x++;
              if (_crt_adc[i] > adc1x) {
                plane1tx = _crt_time[i];
                adc1x += _crt_adc[i];
                plane1x = _crt_pos_x[i];
                plane1xm = _crt_module[i];
              }
            }
          }
          else if (_crt_orient[i]==kCRTHorizontal && _crt_adc[i]>500) { // < 500 hardcoded
            if (nh1y==0 || (_crt_module[i]==plane1ym)) {
              nh1y++;
              if (_crt_adc[i]>adc1y) {
                plane1ty=_crt_time[i];
                adc1y+=_crt_adc[i];
                plane1y=_crt_pos_y[i];
                plane1ym=_crt_module[i];
              }
            }
          }
        }
        else {
          if (_crt_orient[i]==kCRTVertical && _crt_adc[i]>500) { // < 500 hardcoded
            if (nh2x==0 ||  (_crt_module[i]==plane2xm)) {
              nh2x++;
              if (_crt_adc[i]>adc2x) {
                plane2tx=_crt_time[i];
                adc2x+=_crt_adc[i];
                plane2x=_crt_pos_x[i];
                plane2xm=_crt_module[i];
              }
            }
          }
          else if (_crt_orient[i]==kCRTHorizontal && _crt_adc[i]>500) { // < 500 hardcoded
            if (nh2y==0 ||  (_crt_module[i]==plane2ym)) {
              nh2y++;
              if (_crt_adc[i]>adc2y) {
                plane2ty=_crt_time[i];
                adc2y+=_crt_adc[i];
                plane2y=_crt_pos_y[i];
                plane2ym=_crt_module[i];
              }
            }
          }
        }
        for (int j=i+1;j<ns;++j) {
          float tdiff = fabs(_crt_time[i]-_crt_time[j]);
          if (tdiff<0.1) {
            iflag[j]=1;
            if (_crt_plane[j]==sbnd::kCRTFaceSouth) {
              if (_crt_orient[j]==kCRTVertical && _crt_adc[j]>1000) {
                if (nh1x==0 ||  (_crt_module[j]==plane1xm)) {
                  nh1x++;
                  if (_crt_adc[j]>adc1x) {
                    plane1tx=_crt_time[j];
                    adc1x+=_crt_adc[j];
                    plane1x=_crt_pos_x[j];
                    plane1xm=_crt_module[j];
                  }
                }
              }
              else if (_crt_orient[j]==kCRTHorizontal && _crt_adc[j]>1000) {
                if (nh1y==0 ||  (_crt_module[j]==plane1ym)) {
                  nh1y++;
                  if (_crt_adc[j]>adc1y) {
                    plane1ty=_crt_time[j];
                    adc1y+=_crt_adc[j];
                    plane1y=_crt_pos_y[j];
                    plane1ym=_crt_module[j];
                  }
                }
              }
            }
            else {
              if (_crt_orient[j]==kCRTVertical && _crt_adc[j]>1000) {
                if (nh2x==0 ||  (_crt_module[j]==plane2xm)) {
                  nh2x++;
                  if (_crt_adc[j]>adc2x) {
                    plane2tx=_crt_time[j];
                    adc2x+=_crt_adc[j];
                    plane2x=_crt_pos_x[j];
                    plane2xm=_crt_module[j];
                  }
                }
              }
              else if (_crt_orient[j]==kCRTHorizontal && _crt_adc[j]>1000) {
                if (nh2y==0 ||  (_crt_module[j]==plane2ym)) {
                  nh2y++;
                  if (_crt_adc[j]>adc2y) {
                    plane2ty=_crt_time[j];
                    adc2y+=_crt_adc[j];
                    plane2y=_crt_pos_y[j];
                    plane2ym=_crt_module[j];
                  }
                }
              }
            }
          }
              } // look for hits at the same time as hit i
              if (nh1x>0 && nh1y>0 && nh2x>0 && nh2y>0 && adc1x<9000 && adc1y<9000 && adc2x<9000 && adc2y<9000) {
              // make a track!
          _ctrk_x1.push_back(plane1x);
          _ctrk_y1.push_back(plane1y);
          _ctrk_z1.push_back(-239.95);
          _ctrk_t1.push_back(0.5*(plane1tx+plane1ty));
          _ctrk_adc1.push_back(adc1x+adc1y);
          _ctrk_mod1x.push_back(plane1xm);
          _ctrk_x2.push_back(plane2x);
          _ctrk_y2.push_back(plane2y);
          _ctrk_z2.push_back(656.25);
          _ctrk_t2.push_back(0.5*(plane2tx+plane2ty));
          _ctrk_adc2.push_back(adc2x+adc2y);
          _ctrk_mod2x.push_back(plane2xm);
          ntr++;
                // std::cout << "track " << ntr << std::endl;
                // std::cout <<  "x y t adc: plane 1 " << plane1x << " " << plane1y << " " <<
                //   0.5*(plane1tx+plane1ty) << " " << adc1x << " " << adc1y << std::endl;
                // std::cout <<  "         : plane 2 " << plane2x << " " << plane2y << " " <<
                //   0.5*(plane2tx+plane2ty) << " " << adc2x << " " << adc2y << std::endl;
        }
      } // i is the first hit with this time
    } // loop over hits
    _nctrks = ntr;
  }  // end if make tracks

  //
  // CRT hits
  //
  if (fkeepCRThits) {
    art::Handle<std::vector<sbn::crt::CRTHit> > crtHitListHandle;
    std::vector<art::Ptr<sbn::crt::CRTHit> > chitlist;
    // art::Handle< std::vector<sbnd::crt::CRTData> > crtStripListHandle;
    // std::vector< art::Ptr<sbnd::crt::CRTData> > striplist;
    if (evt.getByLabel(fCRTHitModuleLabel, crtHitListHandle))  {
      art::fill_ptr_vector(chitlist, crtHitListHandle);
      _nchits = chitlist.size();
    }
    else {
      std::cout << "Failed to get sbn::crt::CRTHit data product." << std::endl;
      _nchits = 0;
    }

    if (_nchits > _max_chits) {
      std::cout << "Available CRT hits are " << _nchits
                << ", which is above the maximum number allowed to store." << std::endl;
      std::cout << "Will only store " << _max_chits << "CRT hits." << std::endl;
      _nchits = _max_chits;
    }

    ResetCRTHitsVars(_nchits);

    for (int i = 0; i < _nchits; ++i){
      // int ip = kNotDefined;
      sbnd::CRTPlane ip = sbnd::CRTCommonUtils::GetPlaneIndex(chitlist[i]->tagger);

      _chit_time[i]=chitlist[i]->ts1_ns*0.001;
      if (chitlist[i]->ts1_ns > MAX_INT) {
        _chit_time[i] = 0.001 * (chitlist[i]->ts1_ns - TIME_CORRECTION);
      }

      _chit_x[i] = chitlist[i]->x_pos;
      _chit_y[i] = chitlist[i]->y_pos;
      _chit_z[i] = chitlist[i]->z_pos;
      _chit_plane[i] = ip;
    }
  }

  //
  // CRT tracks
  //
  _ncts = 0;
  if (freadCRTtracks) {
    art::Handle<std::vector<sbn::crt::CRTTrack> > crtTrackListHandle;
    std::vector<art::Ptr<sbn::crt::CRTTrack> > ctrklist;
    if (evt.getByLabel(fCRTTrackModuleLabel, crtTrackListHandle))  {
      art::fill_ptr_vector(ctrklist, crtTrackListHandle);
      _ncts = ctrklist.size();
      if (_ncts > _max_nctrks) _ncts = _max_nctrks;

      ResetCRTTracksVars(_ncts);

      for (int i = 0; i < _ncts; ++i){
        _ct_pes[i] = ctrklist[i]->peshit;
        _ct_time[i] = ctrklist[i]->ts1_ns*0.001;
        if (ctrklist[i]->ts1_ns > MAX_INT) {
          _ct_time[i] = 0.001 * (ctrklist[i]->ts1_ns - TIME_CORRECTION);
        }
        _ct_x1[i] = ctrklist[i]->x1_pos;
        _ct_y1[i] = ctrklist[i]->y1_pos;
        _ct_z1[i] = ctrklist[i]->z1_pos;
        _ct_x2[i] = ctrklist[i]->x2_pos;
        _ct_y2[i] = ctrklist[i]->y2_pos;
        _ct_z2[i] = ctrklist[i]->z2_pos;
      }
    } else {
      std::cout << "Failed to get sbn::crt::CRTTrack data product." << std::endl;
    }
  }


  //
  // Optical Hits
  //
  if (freadOpHits) {
    _nophits = 0;
    size_t previous_nophits = 0;

    // Loop over all the ophits labels
    for (auto ophit_label : fOpHitsModuleLabels) {

      art::Handle<std::vector<recob::OpHit>> ophitListHandle;
      std::vector<art::Ptr<recob::OpHit>> ophitlist;
      if (evt.getByLabel(ophit_label, ophitListHandle)) {
        art::fill_ptr_vector(ophitlist, ophitListHandle);
        _nophits += ophitlist.size();
      }
      else {
        std::cout << "Failed to get recob::OpHit data product." << std::endl;
      }

      if (_nophits > _max_ophits) {
        std::cout << "Available optical hits are " << _nophits << ", which is above the maximum number allowed to store." << std::endl;
        std::cout << "Will only store " << _max_ophits << " optical hits." << std::endl;
        _nophits = _max_ophits;
      }

      ResetOpHitsVars(_nophits);

      for (size_t i = 0; i < ophitlist.size(); ++i) {
        size_t index = previous_nophits + i;
        _ophit_opch[index] = ophitlist.at(i)->OpChannel();
        _ophit_opdet[index] = fGeometryService->OpDetFromOpChannel(ophitlist.at(i)->OpChannel());
        _ophit_peakT[index] = ophitlist.at(i)->PeakTime();
        _ophit_width[index] = ophitlist.at(i)->Width();
        _ophit_area[index] = ophitlist.at(i)->Area();
        _ophit_amplitude[index] = ophitlist.at(i)->Amplitude();
        _ophit_pe[index] = ophitlist.at(i)->PE();
        auto opdet_center = fGeometryService->OpDetGeoFromOpChannel(ophitlist.at(i)->OpChannel()).GetCenter();
        _ophit_opdet_x[index] = opdet_center.X();
        _ophit_opdet_y[index] = opdet_center.Y();
        _ophit_opdet_z[index] = opdet_center.Z();
        auto pd_type = _pd_map.pdType(ophitlist.at(i)->OpChannel());
        if (pd_type == "pmt_coated") {_ophit_opdet_type[index] = kPMTCoated;}
        else if (pd_type == "pmt_uncoated") {_ophit_opdet_type[index] = kPMTUnCoated;}
        else if (pd_type == "xarapuca_vis") {_ophit_opdet_type[index] = kXArapucaVis;}
        else if (pd_type == "xarapuca_vuv") {_ophit_opdet_type[index] = kXArapucaVuv;}
        else {_ophit_opdet_type[index] = kPDNotDefined;}
      }
      previous_nophits = _nophits;
    }
  }

  //
  // pmt hardware trigger
  //

  if (freadpmtTrigger){
    art::Handle<std::vector<sbnd::comm::pmtTrigger> > pmtTriggerListHandle;
    std::vector<art::Ptr<sbnd::comm::pmtTrigger> > pmttriggerlist;

    if (evt.getByLabel(fpmtTriggerModuleLabel, pmtTriggerListHandle)){
      art::fill_ptr_vector(pmttriggerlist, pmtTriggerListHandle);
      ResetPmtTriggerVars( (int)pmttriggerlist[0]->numPassed.size());

      for (int i=0; i < (int)pmttriggerlist[0]->numPassed.size(); i++){
        _pmtTrigger_npmtshigh[i] = pmttriggerlist[0]->numPassed[i];
      }
      _pmtTrigger_maxpassed = pmttriggerlist[0]->maxPMTs;

    }
    else{
      std::cout << "Failed to get sbnd::comm::pmtTrigger data product" << std::endl;
    }
  }

  //
  // PMT Software Trigger
  //
  if (freadpmtSoftTrigger){
    art::Handle<sbnd::trigger::pmtSoftwareTrigger> pmtSoftTriggerHandle;
    if (evt.getByLabel(fpmtSoftTriggerModuleLabel, pmtSoftTriggerHandle)){
      const sbnd::trigger::pmtSoftwareTrigger &pmtSoftTriggerMetrics = (*pmtSoftTriggerHandle);
      ResetPmtSoftTriggerVars();
      _pmtSoftTrigger_foundBeamTrigger = pmtSoftTriggerMetrics.foundBeamTrigger;
      _pmtSoftTrigger_tts = pmtSoftTriggerMetrics.triggerTimestamp;
      _pmtSoftTrigger_promptPE = pmtSoftTriggerMetrics.promptPE;
      _pmtSoftTrigger_prelimPE = pmtSoftTriggerMetrics.prelimPE;
      _pmtSoftTrigger_nAboveThreshold = pmtSoftTriggerMetrics.nAboveThreshold;
      // _pmtSoftTrigger_pmtInfoVec = pmtsofttriggerlist[0]->pmtInfoVec;
    }
    else{
      std::cout << "Failed to get sbnd::trigger::pmtSoftwareTrigger data product" << std::endl;
    }
  }

  //
  // CRT Software Trigger
  //
  if (freadcrtSoftTrigger){
    art::Handle<sbndaq::CRTmetric> crtSoftTriggerHandle;
    if (evt.getByLabel(fcrtSoftTriggerModuleLabel, crtSoftTriggerHandle)){
      const sbndaq::CRTmetric &crtSoftTriggerMetrics = (*crtSoftTriggerHandle);
      ResetCrtSoftTriggerVars();
      for (int i=0; i<7; i++){
        _crtSoftTrigger_hitsperplane[i] = crtSoftTriggerMetrics.hitsperplane[i];
      }
    }
    else{
      std::cout << "Failed to get sbndaq::crtMetric data product" << std::endl;
    }

  }

  //
  // Muon tracks 
  //
  _nmuontrks = 0; 
  if (freadMuonTracks){
    art::Handle<std::vector<sbnd::comm::MuonTrack> > muonTrackListHandle;
    std::vector<art::Ptr<sbnd::comm::MuonTrack> > muontrklist;

    if (evt.getByLabel(fMuonTrackModuleLabel, muonTrackListHandle)){
      art::fill_ptr_vector(muontrklist, muonTrackListHandle); 
      _nmuontrks = muontrklist.size();
      ResetMuonTracksVars(_nmuontrks);

      for (int i=0; i < _nmuontrks; i++){ 
        _muontrk_t0[i] = muontrklist[i]->t0_us; 
        _muontrk_x1[i] = muontrklist[i]->x1_pos;
        _muontrk_y1[i] = muontrklist[i]->y1_pos;
        _muontrk_z1[i] = muontrklist[i]->z1_pos;
        _muontrk_x2[i] = muontrklist[i]->x2_pos;
        _muontrk_y2[i] = muontrklist[i]->y2_pos;
        _muontrk_z2[i] = muontrklist[i]->z2_pos; 
        _muontrk_theta_xz[i] = muontrklist[i]->theta_xz; 
        _muontrk_theta_yz[i] = muontrklist[i]->theta_yz;
        _muontrk_tpc[i] = muontrklist[i]->tpc; 
        _muontrk_type[i] = muontrklist[i]->type;
      }
      if (freadMuonHits){
        art::FindMany<recob::Hit> muontrkassn(muonTrackListHandle, evt, fMuonTrackModuleLabel);
        ResetMuonHitVars(3000); //estimate of maximum collection hits
        _nmhits = 0;
        for (int i=0; i < _nmuontrks; i++){ 
        std::vector< const recob::Hit*> muonhitsVec = muontrkassn.at(i);
          _nmhits += (muonhitsVec.size()); 
          for (size_t j=0; j<muonhitsVec.size(); j++){
            auto muonhit = muonhitsVec.at(j);
            geo::WireID wireid = muonhit->WireID();
            _mhit_trk.push_back(i);
            _mhit_tpc.push_back(wireid.TPC);
            _mhit_wire.push_back(wireid.Wire);
            _mhit_channel.push_back(muonhit->Channel());
            _mhit_peakT.push_back(muonhit->PeakTime());
            _mhit_charge.push_back(muonhit->Integral());
          }
        }
      }
    }
    else{
      std::cout << "Failed to get sbnd::comm::MuonTrack data product" << std::endl;
    }
  }

  if (fcheckTransparency) {

    _waveform_number.resize(_max_hits*_max_samples, -9999.);
    _adc_on_wire.resize(_max_hits*_max_samples, -9999.);
    _time_for_waveform.resize(_max_hits*_max_samples, -9999.);
    _waveform_integral.resize(_max_hits*_max_samples, -9999.);
    _adc_count_in_waveform.resize(_max_hits*_max_samples, -9999.);

    art::Handle<std::vector<raw::RawDigit>> digitVecHandle;

    bool retVal = evt.getByLabel(fDigitModuleLabel, digitVecHandle);
    if(retVal == true) {
      mf::LogInfo("HitDumper")    << "I got fDigitModuleLabel: "         << fDigitModuleLabel << std::endl;
    } else {
      mf::LogWarning("HitDumper") << "Could not get fDigitModuleLabel: " << fDigitModuleLabel << std::endl;
    }

    int waveform_number_tracker = 0;
    int adc_counter = 1;
    _adc_count = _nhits * (fWindow * 2 + 1);

    // loop over waveforms
    for(size_t rdIter = 0; rdIter < digitVecHandle->size(); ++rdIter) {

      //GET THE REFERENCE TO THE CURRENT raw::RawDigit.
      art::Ptr<raw::RawDigit> digitVec(digitVecHandle, rdIter);
      int channel   = digitVec->Channel();
      auto fDataSize = digitVec->Samples();
      std::vector<short> rawadc;      //UNCOMPRESSED ADC VALUES.
      rawadc.resize(fDataSize);

      // see if there is a hit on this channel
      for (int ihit = 0; ihit < _nhits; ++ihit) {
        if (_hit_channel[ihit] == channel) {

          int pedestal = (int)digitVec->GetPedestal();
          //UNCOMPRESS THE DATA.
          if (fUncompressWithPed) {
            raw::Uncompress(digitVec->ADCs(), rawadc, pedestal, digitVec->Compression());
          }
          else {
            raw::Uncompress(digitVec->ADCs(), rawadc, digitVec->Compression());
          }

          unsigned int bin = _hit_peakT[ihit];
          unsigned int low_edge,high_edge;
          if((int)bin > fWindow and _hit_plane[ihit] == 0) {
            low_edge = bin - (2*fWindow);
          }
          else if ((int)bin>fWindow) {
            low_edge = bin-fWindow;
          }
          else {
            low_edge = 0;
          }
          high_edge = bin + fWindow;
          if (high_edge > (fDataSize-1)) {
            high_edge = fDataSize - 1;
          }
          double integral = 0.0;
          waveform_number_tracker++;
          int counter_for_adc_in_waveform = 0;
          for (size_t ibin = low_edge; ibin <= high_edge; ++ibin) {
            _adc_count_in_waveform[adc_counter] = counter_for_adc_in_waveform;
            counter_for_adc_in_waveform++;
            _waveform_number[adc_counter] = waveform_number_tracker;
            _adc_on_wire[adc_counter] = rawadc[ibin]-pedestal;
            _time_for_waveform[adc_counter] = ibin;
            //std::cout << "DUMP: " << _waveform_number[adc_counter] << " " << _adc_count << " " << _hit_plane[ihit] << " " << _hit_wire[ihit] << " " <<ibin << " " << (rawadc[ibin]-pedestal) << " " << _time_for_waveform[adc_counter] << " " << _adc_on_wire[adc_counter] << std::endl;
            integral+=_adc_on_wire[adc_counter];
            _waveform_integral[adc_counter] = integral;
            adc_counter++;
          }
          std::cout << "DUMP SUM: " << _hit_tpc[ihit] << " " << _hit_plane[ihit] << " " << _hit_wire[ihit] << " " <<  integral << " " << waveform_number_tracker << std::endl;
          _hit_full_integral[ihit] = integral;
        } // if hit channel matches waveform channel
      } //end loop over hits
    }// end loop over waveforms
  }// end if fCheckTrasparency


  if (freadTruth){

    int nGeniePrimaries = 0, nMCNeutrinos = 0;
    art::Handle< std::vector<simb::MCTruth> > mctruthListHandle;
    std::vector<art::Ptr<simb::MCTruth> > mclist;
    if (evt.getByLabel(fGenieGenModuleLabel,mctruthListHandle)){
      art::fill_ptr_vector(mclist, mctruthListHandle);
    }else {
      std::cout << "Failed to get Genie data product." << std::endl;
    }

    art::Ptr<simb::MCTruth> mctruth;

      if (!mclist.empty()) {//at least one mc record

        mctruth = mclist[0];

        if (mctruth->NeutrinoSet()) nGeniePrimaries = mctruth->NParticles();

    } // if have MC truth
      MF_LOG_DEBUG("HitDumper") << "Expected " << nGeniePrimaries << " GENIE particles";

    //Initially call the number of neutrinos to be stored the number of MCTruth objects.  This is not strictly true i.e. BNB + cosmic overlay but we will count the number of neutrinos later
    nMCNeutrinos = mclist.size();

    ResizeGenie(nGeniePrimaries);
    ResizeMCNeutrino(nMCNeutrinos);

    mcevts_truth = mclist.size();
    //Brailsford 2017/10/16
    //Issue 17917
    //To keep a 1:1 between neutrinos and 'flux' we need the assns
    art::FindManyP<simb::MCFlux> fmFluxNeutrino(mctruthListHandle, evt, fGenieGenModuleLabel);
    // Get GTruth information for scattering code
    art::FindManyP< simb::GTruth > fmgt( mctruthListHandle, evt, fGenieGenModuleLabel );

    if (mcevts_truth > 0){//at least one mc record

      //Brailsford 2017/10/16
      //Issue 17917
      //Loop over every truth in the spill rather than just looking at the first one.
      //Because MCTruth could be a neutrino OR something else (e.g. cosmics) we are going to have to count up how many neutrinos there are
      mcevts_truth = 0;
      for (unsigned int i_mctruth = 0; i_mctruth < mclist.size(); i_mctruth++){
        //fetch an mctruth
        art::Ptr<simb::MCTruth> curr_mctruth = mclist[i_mctruth];
        //Check if it's a neutrino
        if (!curr_mctruth->NeutrinoSet()) continue;

        // Genie Truth association only for the neutrino
        if (fmgt.size()>i_mctruth) {
          std::vector< art::Ptr<simb::GTruth> > mcgtAssn = fmgt.at(i_mctruth);

          nuScatterCode_truth[i_mctruth] = mcgtAssn[0]->fGscatter;
        } else {
          nuScatterCode_truth[i_mctruth] = -1.;
        }

        nuPDG_truth[i_mctruth] = curr_mctruth->GetNeutrino().Nu().PdgCode();
        ccnc_truth[i_mctruth] = curr_mctruth->GetNeutrino().CCNC();
        mode_truth[i_mctruth] = curr_mctruth->GetNeutrino().Mode();
        Q2_truth[i_mctruth] = curr_mctruth->GetNeutrino().QSqr();
        W_truth[i_mctruth] = curr_mctruth->GetNeutrino().W();
        hitnuc_truth[i_mctruth] = curr_mctruth->GetNeutrino().HitNuc();
        enu_truth[i_mctruth] = curr_mctruth->GetNeutrino().Nu().E();
        nuvtxx_truth[i_mctruth] = curr_mctruth->GetNeutrino().Nu().Vx();
        nuvtxy_truth[i_mctruth] = curr_mctruth->GetNeutrino().Nu().Vy();
        nuvtxz_truth[i_mctruth] = curr_mctruth->GetNeutrino().Nu().Vz();
        if (curr_mctruth->GetNeutrino().Nu().P()){
          nu_dcosx_truth[i_mctruth] = curr_mctruth->GetNeutrino().Nu().Px()/curr_mctruth->GetNeutrino().Nu().P();
          nu_dcosy_truth[i_mctruth] = curr_mctruth->GetNeutrino().Nu().Py()/curr_mctruth->GetNeutrino().Nu().P();
          nu_dcosz_truth[i_mctruth] = curr_mctruth->GetNeutrino().Nu().Pz()/curr_mctruth->GetNeutrino().Nu().P();
        }
        lep_mom_truth[i_mctruth] = curr_mctruth->GetNeutrino().Lepton().P();
        if (curr_mctruth->GetNeutrino().Lepton().P()){
          lep_dcosx_truth[i_mctruth] = curr_mctruth->GetNeutrino().Lepton().Px()/curr_mctruth->GetNeutrino().Lepton().P();
          lep_dcosy_truth[i_mctruth] = curr_mctruth->GetNeutrino().Lepton().Py()/curr_mctruth->GetNeutrino().Lepton().P();
          lep_dcosz_truth[i_mctruth] = curr_mctruth->GetNeutrino().Lepton().Pz()/curr_mctruth->GetNeutrino().Lepton().P();
        }
        //Brailsford
        //2017/10/17
        //Issue 12918
        //Use the art::Ptr key as the neutrino's unique ID
        nuID_truth[i_mctruth] = curr_mctruth.key();
        //We need to also store N 'flux' neutrinos per event so now check that the FindOneP is valid and, if so, use it!
        if (fmFluxNeutrino.isValid()){
          if (fmFluxNeutrino.at(0).size()>i_mctruth){
          art::Ptr<simb::MCFlux> curr_mcflux = fmFluxNeutrino.at(0).at(i_mctruth);
          tpx_flux[i_mctruth] = curr_mcflux->ftpx;
          tpy_flux[i_mctruth] = curr_mcflux->ftpy;
          tpz_flux[i_mctruth] = curr_mcflux->ftpz;
          tptype_flux[i_mctruth] = curr_mcflux->ftptype;
          }
        }

        //Let's increase the neutrino count
        mcevts_truth++;
      }

      if (mctruth->NeutrinoSet()){
        //genie particles information
        genie_no_primaries = mctruth->NParticles();

        size_t StoreParticles = std::min((size_t) genie_no_primaries, MaxGeniePrimaries);
        if (genie_no_primaries > (int) StoreParticles) {
          // got this error? it might be a bug,
          // since the structure should have enough room for everything
          mf::LogError("HitDumper") << "event has "
            << genie_no_primaries << " MC particles, only "
            << StoreParticles << " stored in tree";
        }
        for(size_t iPart = 0; iPart < StoreParticles; ++iPart){
          const simb::MCParticle& part(mctruth->GetParticle(iPart));
          genie_primaries_pdg[iPart]=part.PdgCode();
          genie_Eng[iPart]=part.E();
          genie_Px[iPart]=part.Px();
          genie_Py[iPart]=part.Py();
          genie_Pz[iPart]=part.Pz();
          genie_P[iPart]=part.P();
          genie_status_code[iPart]=part.StatusCode();
          genie_mass[iPart]=part.Mass();
          genie_trackID[iPart]=part.TrackId();
          genie_ND[iPart]=part.NumberDaughters();
          genie_mother[iPart]=part.Mother();
        } // for particle
      } //if neutrino set
    }//if (mcevts_truth)
  }//if (fReadTruth){



  fTree->Fill();

}

void Hitdumper::beginJob ( )

override

Definition at line 1190 of file HitDumper_module.cc.

 {
  // Implementation of optional member function here.
  art::ServiceHandle<art::TFileService> tfs;
  fTree = tfs->make<TTree>("hitdumpertree","analysis tree");
  fTree->Branch("run",&_run,"run/I");
  fTree->Branch("subrun",&_subrun,"subrun/I");
  fTree->Branch("event",&_event,"event/I");
  fTree->Branch("evttime",&_evttime,"evttime/D");
  fTree->Branch("t0",&_t0,"t0/I");

  fTree->Branch("nhits", &_nhits, "nhits/I");
  fTree->Branch("hit_cryostat", &_hit_cryostat);
  fTree->Branch("hit_tpc", &_hit_tpc);
  fTree->Branch("hit_plane", &_hit_plane);
  fTree->Branch("hit_wire", &_hit_wire);
  fTree->Branch("hit_channel", &_hit_channel);
  fTree->Branch("hit_peakT", &_hit_peakT);
  fTree->Branch("hit_charge", &_hit_charge);
  fTree->Branch("hit_ph", &_hit_ph);
  fTree->Branch("hit_width", &_hit_width);
  fTree->Branch("hit_full_integral", &_hit_full_integral);

  if (fcheckTransparency) {
    fTree->Branch("adc_count", &_adc_count,"adc_count/I");
    fTree->Branch("waveform_number", &_waveform_number);
    fTree->Branch("time_for_waveform",&_time_for_waveform);
    fTree->Branch("adc_on_wire", &_adc_on_wire);
    fTree->Branch("waveform_integral", &_waveform_integral);
    fTree->Branch("adc_count_in_waveform", &_adc_count_in_waveform);
  }

  if (fkeepCRTstrips) {
    fTree->Branch("nstrips", &_nstrips, "nstrips/I");
    fTree->Branch("crt_plane", &_crt_plane);
    fTree->Branch("crt_module", &_crt_module);
    fTree->Branch("crt_strip", &_crt_strip);
    fTree->Branch("crt_orient", &_crt_orient);
    fTree->Branch("crt_time", &_crt_time);
    fTree->Branch("crt_adc", &_crt_adc);
    fTree->Branch("crt_pos_x", &_crt_pos_x);
    fTree->Branch("crt_pos_y", &_crt_pos_y);
    fTree->Branch("crt_pos_z", &_crt_pos_z);
  }

  if (fmakeCRTtracks) {
    fTree->Branch("nctrks", &_nctrks, "nctrks/I");
    fTree->Branch("ctrk_x1", &_ctrk_x1);
    fTree->Branch("ctrk_y1", &_ctrk_y1);
    fTree->Branch("ctrk_z1", &_ctrk_z1);
    fTree->Branch("ctrk_t1", &_ctrk_t1);
    fTree->Branch("ctrk_adc1", &_ctrk_adc1);
    fTree->Branch("ctrk_mod1x", &_ctrk_mod1x);
    fTree->Branch("ctrk_x2", &_ctrk_x2);
    fTree->Branch("ctrk_y2", &_ctrk_y2);
    fTree->Branch("ctrk_z2", &_ctrk_z2);
    fTree->Branch("ctrk_t2", &_ctrk_t2);
    fTree->Branch("ctrk_adc2", &_ctrk_adc2);
    fTree->Branch("ctrk_mod2x", &_ctrk_mod2x);
  }
  if (fkeepCRThits) {
    fTree->Branch("nchits", &_nchits, "nchits/I");
    fTree->Branch("chit_x", &_chit_x);
    fTree->Branch("chit_y", &_chit_y);
    fTree->Branch("chit_z", &_chit_z);
    fTree->Branch("chit_time", &_chit_time);
    fTree->Branch("chit_plane", &_chit_plane);
  }
  if (freadCRTtracks) {
    fTree->Branch("ncts", &_ncts, "ncts/I");
    fTree->Branch("ct_x1", &_ct_x1);
    fTree->Branch("ct_y1", &_ct_y1);
    fTree->Branch("ct_z1", &_ct_z1);
    fTree->Branch("ct_x2", &_ct_x2);
    fTree->Branch("ct_y2", &_ct_y2);
    fTree->Branch("ct_z2", &_ct_z2);
    fTree->Branch("ct_time", &_ct_time);
    fTree->Branch("ct_pes", &_ct_pes);
  }

  if (freadOpHits) {
    fTree->Branch("nophits", &_nophits, "nophits/I");
    fTree->Branch("ophit_opch", &_ophit_opch);
    fTree->Branch("ophit_opdet", &_ophit_opdet);
    fTree->Branch("ophit_peakT", &_ophit_peakT);
    fTree->Branch("ophit_width", &_ophit_width);
    fTree->Branch("ophit_area", &_ophit_area);
    fTree->Branch("ophit_amplitude", &_ophit_amplitude);
    fTree->Branch("ophit_pe", &_ophit_pe);
    fTree->Branch("ophit_opdet_x", &_ophit_opdet_x);
    fTree->Branch("ophit_opdet_y", &_ophit_opdet_y);
    fTree->Branch("ophit_opdet_z", &_ophit_opdet_z);
    fTree->Branch("ophit_opdet_type", &_ophit_opdet_type);
  }

  if (freadpmtTrigger){
    fTree->Branch("pmtTrigger_npmtshigh", &_pmtTrigger_npmtshigh);
    fTree->Branch("pmtTrigger_maxpassed", &_pmtTrigger_maxpassed, "pmtTrigger_maxpassed/I");
  }

  if (freadpmtSoftTrigger){
    fTree->Branch("pmtSoftTrigger_foundBeamTrigger", &_pmtSoftTrigger_foundBeamTrigger);
    fTree->Branch("pmtSoftTrigger_tts", &_pmtSoftTrigger_tts);
    fTree->Branch("pmtSoftTrigger_promptPE", &_pmtSoftTrigger_promptPE);
    fTree->Branch("pmtSoftTrigger_prelimPE", &_pmtSoftTrigger_prelimPE);
    fTree->Branch("pmtSoftTrigger_nAboveThreshold", &_pmtSoftTrigger_nAboveThreshold);
    // fTree->Branch("pmtSoftTrigger_", &_pmtSoftTigger_)
  }

  if (freadcrtSoftTrigger){
    fTree->Branch("crtSoftTrigger_hitsperplane", &_crtSoftTrigger_hitsperplane,"crtSoftTrigger_hitsperplane[7]/I");
  }

  if (freadMuonTracks) {
    fTree->Branch("nmuontrks", &_nmuontrks, "nmuontrks/I");
    fTree->Branch("muontrk_t0", &_muontrk_t0);
    fTree->Branch("muontrk_x1", &_muontrk_x1);
    fTree->Branch("muontrk_y1", &_muontrk_y1);
    fTree->Branch("muontrk_z1", &_muontrk_z1);
    fTree->Branch("muontrk_x2", &_muontrk_x2);
    fTree->Branch("muontrk_y2", &_muontrk_y2);
    fTree->Branch("muontrk_z2", &_muontrk_z2);
    fTree->Branch("muontrk_theta_xz", &_muontrk_theta_xz);
    fTree->Branch("muontrk_theta_yz", &_muontrk_theta_yz);
    fTree->Branch("muontrk_tpc", &_muontrk_tpc); 
    fTree->Branch("muontrk_type", &_muontrk_type); 
  }

    if (freadMuonHits) {
    fTree->Branch("nmhits", &_nmhits, "nmhits/I");
    fTree->Branch("mhit_trk", &_mhit_trk);
    fTree->Branch("mhit_tpc", &_mhit_tpc);
    fTree->Branch("mhit_wire", &_mhit_wire); 
    fTree->Branch("mhit_channel", &_mhit_channel);
    fTree->Branch("mhit_peakT", &_mhit_peakT);
    fTree->Branch("mhit_charge", &_mhit_charge); 
  }

  if (freadTruth) {
    fTree->Branch("mcevts_truth",&mcevts_truth,"mcevts_truth/I");
    fTree->Branch("nuScatterCode_truth",&nuScatterCode_truth);
    fTree->Branch("nuID_truth",&nuID_truth);
    fTree->Branch("nuPDG_truth",&nuPDG_truth);
    fTree->Branch("ccnc_truth",&ccnc_truth);
    fTree->Branch("mode_truth",&mode_truth);
    fTree->Branch("enu_truth",&enu_truth);
    fTree->Branch("Q2_truth",&Q2_truth);
    fTree->Branch("W_truth",&W_truth);
    fTree->Branch("hitnuc_truth",&hitnuc_truth);
    fTree->Branch("nuvtxx_truth",&nuvtxx_truth);
    fTree->Branch("nuvtxy_truth",&nuvtxy_truth);
    fTree->Branch("nuvtxz_truth",&nuvtxz_truth);
    fTree->Branch("nu_dcosx_truth",&nu_dcosx_truth);
    fTree->Branch("nu_dcosy_truth",&nu_dcosy_truth);
    fTree->Branch("nu_dcosz_truth",&nu_dcosz_truth);
    fTree->Branch("lep_mom_truth",&lep_mom_truth);
    fTree->Branch("lep_dcosx_truth",&lep_dcosx_truth);
    fTree->Branch("lep_dcosy_truth",&lep_dcosy_truth);
    fTree->Branch("lep_dcosz_truth",&lep_dcosz_truth);

    fTree->Branch("tpx_flux",&tpx_flux);
    fTree->Branch("tpy_flux",&tpy_flux);
    fTree->Branch("tpz_flux",&tpz_flux);
    fTree->Branch("tptype_flux",&tptype_flux);

    fTree->Branch("genie_no_primaries",&genie_no_primaries);
    fTree->Branch("genie_primaries_pdg",&genie_primaries_pdg);
    fTree->Branch("genie_Eng",&genie_Eng);
    fTree->Branch("genie_Px",&genie_Px);
    fTree->Branch("genie_Py",&genie_Py);
    fTree->Branch("genie_Pz",&genie_Pz);
    fTree->Branch("genie_P",&genie_P);
    fTree->Branch("genie_status_code",&genie_status_code);
    fTree->Branch("genie_mass",&genie_mass);
    fTree->Branch("genie_trackID",&genie_trackID);
    fTree->Branch("genie_ND",&genie_ND);
    fTree->Branch("genie_mother",&genie_mother);
  }

  if (fsavePOTInfo) {
    _sr_tree = tfs->make<TTree>("pottree","");
    _sr_tree->Branch("run", &_sr_run, "run/I");
    _sr_tree->Branch("subrun", &_sr_subrun, "subrun/I");
    _sr_tree->Branch("begintime", &_sr_begintime, "begintime/D");
    _sr_tree->Branch("endtime", &_sr_endtime, "endtime/D");
    _sr_tree->Branch("pot", &_sr_pot, "pot/D");
  }

}

void Hitdumper::beginSubRun ( art::SubRun const &  sr )

overridevirtual

Definition at line 1591 of file HitDumper_module.cc.

                                               {

  if (!fsavePOTInfo) {
    return;
  }

  _sr_run       = sr.run();
  _sr_subrun    = sr.subRun();
  _sr_begintime = sr.beginTime().value();
  _sr_endtime   = sr.endTime().value();

  art::Handle<sumdata::POTSummary> pot_handle;
  sr.getByLabel(fGenieGenModuleLabel, pot_handle);

  if (pot_handle.isValid()) {
    _sr_pot = pot_handle->totpot;
  } else {
    _sr_pot = 0.;
  }

  _sr_tree->Fill();

}

void Hitdumper::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 400 of file HitDumper_module.cc.

{

  _max_hits = p.get<int>("MaxHits", 50000);
  _max_ophits = p.get<int>("MaxOpHits", 50000);
  _max_samples = p.get<int>("MaxSamples", 5001);
  _max_chits = p.get<int>("MaxCRTHits", 5000);
  _max_nctrks = p.get<int>("MaxCRTTracks", 10);

  fHitsModuleLabel     = p.get<std::string>("HitsModuleLabel");
  fDigitModuleLabel    = p.get<std::string>("DigitModuleLabel", "daq");
  fLArG4ModuleLabel    = p.get<std::string>("LArG4ModuleLabel", "largeant");
  fCRTStripModuleLabel = p.get<std::string>("CRTStripModuleLabel", "crt");
  fCRTHitModuleLabel   = p.get<std::string>("CRTHitModuleLabel", "crthit");
  fCRTTrackModuleLabel = p.get<std::string>("CRTTrackModuleLabel", "crttrack");
  fOpHitsModuleLabels  = p.get<std::vector<std::string>>("OpHitsModuleLabel");
  fpmtTriggerModuleLabel = p.get<std::string>("pmtTriggerModuleLabel", "pmttriggerproducer");
  fpmtSoftTriggerModuleLabel = p.get<std::string>("pmtSoftTriggerModuleLabel", "pmtSoftwareTrigger");
  fcrtSoftTriggerModuleLabel = p.get<std::string>("crtSoftTriggerModuleLabel", "MetricProducer");
  fMuonTrackModuleLabel  = p.get<std::string>("MuonTrackModuleLabel", "MuonTrackProducer");
  fGenieGenModuleLabel = p.get<std::string>("GenieGenModuleLabel", "generator");

  fkeepCRThits       = p.get<bool>("keepCRThits",true);
  fkeepCRTstrips     = p.get<bool>("keepCRTstrips",false);
  fmakeCRTtracks     = p.get<bool>("makeCRTtracks",true);
  freadCRTtracks     = p.get<bool>("readCRTtracks",true);
  freadOpHits        = p.get<bool>("readOpHits",true);
  freadpmtTrigger    = p.get<bool>("readpmtTrigger",true);
  freadpmtSoftTrigger= p.get<bool>("readpmtSoftTrigger",true);
  freadcrtSoftTrigger= p.get<bool>("readcrtSoftTrigger",false);
  freadMuonTracks    = p.get<bool>("readMuonTracks",true);
  freadMuonHits      = p.get<bool>("readMuonHits",false);
  fcheckTransparency = p.get<bool>("checkTransparency",false);
  freadTruth         = p.get<bool>("readTruth",true);
  fsavePOTInfo       = p.get<bool>("savePOTinfo",true);
  fUncompressWithPed = p.get<bool>("UncompressWithPed",false);

  fWindow            = p.get<int>("window",100);
  fKeepTaggerTypes   = p.get<std::vector<int>>("KeepTaggerTypes");

  fSkipInd           = p.get<bool>("SkipInduction",false);
}

void Hitdumper::ResetCRTCustomTracksVars ( int  n )

private

Resets custom crt tracks tree variables.

Definition at line 1416 of file HitDumper_module.cc.

                                              {
  _ctrk_x1.clear();
  _ctrk_y1.clear();
  _ctrk_z1.clear();
  _ctrk_t1.clear();
  _ctrk_adc1.clear();
  _ctrk_mod1x.clear();
  _ctrk_x2.clear();
  _ctrk_y2.clear();
  _ctrk_z2.clear();
  _ctrk_t2.clear();
  _ctrk_adc2.clear();
  _ctrk_mod2x.clear();

  _ctrk_x1.reserve(n);
  _ctrk_y1.reserve(n);
  _ctrk_z1.reserve(n);
  _ctrk_t1.reserve(n);
  _ctrk_adc1.reserve(n);
  _ctrk_mod1x.reserve(n);
  _ctrk_x2.reserve(n);
  _ctrk_y2.reserve(n);
  _ctrk_z2.reserve(n);
  _ctrk_t2.reserve(n);
  _ctrk_adc2.reserve(n);
  _ctrk_mod2x.reserve(n);
}

void Hitdumper::ResetCRTHitsVars ( int  n )

private

Resets crt hits tree variables.

Definition at line 1455 of file HitDumper_module.cc.

                                      {
  _chit_plane.assign(n, DEFAULT_VALUE);
  _chit_time.assign(n, DEFAULT_VALUE);
  _chit_x.assign(n, DEFAULT_VALUE);
  _chit_y.assign(n, DEFAULT_VALUE);
  _chit_z.assign(n, DEFAULT_VALUE);
}

void Hitdumper::ResetCrtSoftTriggerVars ( )

private

Resets crt software trigger variables.

Definition at line 1490 of file HitDumper_module.cc.

                                       {
  for (int i=0; i<7; i++){
    _crtSoftTrigger_hitsperplane[i] = 0;
  }
}

void Hitdumper::ResetCRTStripsVars ( int  n )

private

Resets crt strips tree variables.

Definition at line 1394 of file HitDumper_module.cc.

                                        {
  _crt_plane.clear();
  _crt_module.clear();
  _crt_strip.clear();
  _crt_orient.clear();
  _crt_time.clear();
  _crt_adc.clear();
  _crt_pos_x.clear();
  _crt_pos_y.clear();
  _crt_pos_z.clear();

  _crt_plane.reserve(n);
  _crt_module.reserve(n);
  _crt_strip.reserve(n);
  _crt_orient.reserve(n);
  _crt_time.reserve(n);
  _crt_adc.reserve(n);
  _crt_pos_x.reserve(n);
  _crt_pos_y.reserve(n);
  _crt_pos_z.reserve(n);
}

void Hitdumper::ResetCRTTracksVars ( int  n )

private

Resets crt tracks tree variables.

Definition at line 1444 of file HitDumper_module.cc.

                                        {
  _ct_x1.assign(n, DEFAULT_VALUE);
  _ct_y1.assign(n, DEFAULT_VALUE);
  _ct_z1.assign(n, DEFAULT_VALUE);
  _ct_time.assign(n, DEFAULT_VALUE);
  _ct_pes.assign(n, DEFAULT_VALUE);
  _ct_x2.assign(n, DEFAULT_VALUE);
  _ct_y2.assign(n, DEFAULT_VALUE);
  _ct_z2.assign(n, DEFAULT_VALUE);
}

void Hitdumper::ResetMuonHitVars ( int  n )

private

Resets crossing muon hit tree variables.

Definition at line 1510 of file HitDumper_module.cc.

                                     {
  _mhit_trk.clear(); 
  _mhit_tpc.clear();
  _mhit_wire.clear();
  _mhit_channel.clear();
  _mhit_peakT.clear();
  _mhit_charge.clear();

  _mhit_trk.reserve(n);
  _mhit_tpc.reserve(n);
  _mhit_wire.reserve(n);
  _mhit_channel.reserve(n);
  _mhit_peakT.reserve(n);
  _mhit_charge.reserve(n);
}

void Hitdumper::ResetMuonTracksVars ( int  n )

private

Resets crossing muon tracks tree variables.

Definition at line 1496 of file HitDumper_module.cc.

                                        {
  _muontrk_t0.assign(n, DEFAULT_VALUE);
  _muontrk_x1.assign(n, DEFAULT_VALUE);
  _muontrk_y1.assign(n, DEFAULT_VALUE);
  _muontrk_z1.assign(n, DEFAULT_VALUE);
  _muontrk_x2.assign(n, DEFAULT_VALUE);
  _muontrk_y2.assign(n, DEFAULT_VALUE);
  _muontrk_z2.assign(n, DEFAULT_VALUE); 
  _muontrk_theta_xz.assign(n, DEFAULT_VALUE); 
  _muontrk_theta_yz.assign(n, DEFAULT_VALUE);
  _muontrk_tpc.assign(n, DEFAULT_VALUE);
  _muontrk_type.assign(n, DEFAULT_VALUE);
}

void Hitdumper::ResetOpHitsVars ( int  n )

private

Resets optical hits tree variables.

Definition at line 1463 of file HitDumper_module.cc.

                                     {
  _ophit_opch.resize(n, DEFAULT_VALUE);
  _ophit_opdet.resize(n, DEFAULT_VALUE);
  _ophit_peakT.resize(n, DEFAULT_VALUE);
  _ophit_width.resize(n, DEFAULT_VALUE);
  _ophit_area.resize(n, DEFAULT_VALUE);
  _ophit_amplitude.resize(n, DEFAULT_VALUE);
  _ophit_pe.resize(n, DEFAULT_VALUE);
  _ophit_opdet_x.resize(n, DEFAULT_VALUE);
  _ophit_opdet_y.resize(n, DEFAULT_VALUE);
  _ophit_opdet_z.resize(n, DEFAULT_VALUE);
  _ophit_opdet_type.resize(n, DEFAULT_VALUE);
}

void Hitdumper::ResetPmtSoftTriggerVars ( )

private

Rests pmt software trigger variables.

Definition at line 1482 of file HitDumper_module.cc.

                                       {
  _pmtSoftTrigger_foundBeamTrigger = false;
  _pmtSoftTrigger_tts = 0;
  _pmtSoftTrigger_promptPE = 0;
  _pmtSoftTrigger_prelimPE = 0;
  _pmtSoftTrigger_nAboveThreshold = 0;
}

void Hitdumper::ResetPmtTriggerVars ( int  n )

private

Resets pmt hardware trigger variables.

Definition at line 1477 of file HitDumper_module.cc.

                                        {
  _pmtTrigger_npmtshigh.assign(n, DEFAULT_VALUE);
  _pmtTrigger_maxpassed = 0;
}

void Hitdumper::ResetVars ( )

private

Resets the variables that are saved to the TTree.

Definition at line 1526 of file HitDumper_module.cc.

                          {


  _run = -99999;
  _subrun = -99999;
  _event = -99999;
  _evttime = -99999;
  _t0 = -99999;

  mcevts_truth = 0;
  genie_no_primaries = 0;
}

void Hitdumper::ResetWireHitsVars ( int  n )

private

Resets wire hits tree variables.

Definition at line 1381 of file HitDumper_module.cc.

                                       {
  _hit_cryostat.assign(n, DEFAULT_VALUE);
  _hit_tpc.assign(n, DEFAULT_VALUE);
  _hit_plane.assign(n, DEFAULT_VALUE);
  _hit_wire.assign(n, DEFAULT_VALUE);
  _hit_channel.assign(n, DEFAULT_VALUE);
  _hit_peakT.assign(n, DEFAULT_VALUE);
  _hit_charge.assign(n, DEFAULT_VALUE);
  _hit_ph.assign(n, DEFAULT_VALUE);
  _hit_width.assign(n, DEFAULT_VALUE);
  _hit_full_integral.assign(n, DEFAULT_VALUE);
}

void Hitdumper::ResizeGenie ( int  nPrimaries )

private

Resize the data structure for Genie primaries.

Definition at line 1571 of file HitDumper_module.cc.

                                          {

  // minimum size is 1, so that we always have an address
  MaxGeniePrimaries = (size_t) std::max(nPrimaries, 1);

  genie_primaries_pdg.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_Eng.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_Px.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_Py.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_Pz.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_P.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_status_code.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_mass.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_trackID.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_ND.assign(MaxGeniePrimaries, DEFAULT_VALUE);
  genie_mother.assign(MaxGeniePrimaries, DEFAULT_VALUE);

}

void Hitdumper::ResizeMCNeutrino ( int  nNeutrinos )

private

Resize the data structure for MCNeutrino particles.

Definition at line 1539 of file HitDumper_module.cc.

                                               {

  //min size is 1, to guarantee an address
  MaxMCNeutrinos = (size_t) std::max(nNeutrinos, 1);

  nuScatterCode_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  nuID_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  nuPDG_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  ccnc_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  mode_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  enu_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  Q2_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  W_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  hitnuc_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  nuvtxx_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  nuvtxy_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  nuvtxz_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  nu_dcosx_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  nu_dcosy_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  nu_dcosz_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  lep_mom_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  lep_dcosx_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  lep_dcosy_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  lep_dcosz_truth.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  //Also resize the flux information here as it's a 1:1 with the MCNeutrino
  tpx_flux.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  tpy_flux.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  tpz_flux.assign(MaxMCNeutrinos, DEFAULT_VALUE);
  tptype_flux.assign(MaxMCNeutrinos, DEFAULT_VALUE);

}

Member Data Documentation

int Hitdumper::_adc_count

private

Used for plotting waveforms.

Definition at line 180 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_adc_count_in_waveform

private

Used to view all waveforms on a hitplane together.

Definition at line 182 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_adc_on_wire

private

ADC on wire to draw waveform.

Definition at line 178 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_chit_plane

private

CRT hit plane.

Definition at line 219 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_chit_time

private

CRT hit time.

Definition at line 217 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_chit_x

private

CRT hit x.

Definition at line 214 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_chit_y

private

CRT hit y.

Definition at line 215 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_chit_z

private

CRT hit z.

Definition at line 216 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_crt_adc

private

CRT adc.

Definition at line 192 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_crt_module

private

CRT module.

Definition at line 188 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_crt_orient

private

CRT orientation (0 for y (horizontal) and 1 for x (vertical))

Definition at line 190 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_crt_plane

private

CRT plane.

Definition at line 187 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_crt_pos_x

private

CRT position X.

Definition at line 193 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_crt_pos_y

private

CRT position Y.

Definition at line 194 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_crt_pos_z

private

CRT position Z.

Definition at line 195 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_crt_strip

private

CRT strip.

Definition at line 189 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_crt_time

private

CRT time.

Definition at line 191 of file HitDumper_module.cc.

int Hitdumper::_crtSoftTrigger_hitsperplane[7]

private

number of individual PMTs above ADC threshold (fcl) during the beam spill

Number of (very low level) CRT hits per plane

Definition at line 258 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ct_pes

private

CRT track PEs.

Definition at line 224 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ct_time

private

CRT track time.

Definition at line 223 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ct_x1

private

CRT track x1.

Definition at line 225 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ct_x2

private

CRT track x2.

Definition at line 228 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ct_y1

private

CRT track y1.

Definition at line 226 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ct_y2

private

CRT track y2.

Definition at line 229 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ct_z1

private

CRT track z1.

Definition at line 227 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ct_z2

private

CRT track z2.

Definition at line 230 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_adc1

private

CRT track adc1.

Definition at line 203 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_adc2

private

CRT track adc2.

Definition at line 209 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_ctrk_mod1x

private

CRT track mod2x.

Definition at line 204 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_ctrk_mod2x

private

CRT track mod2x.

Definition at line 210 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_t1

private

CRT track t1.

Definition at line 202 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_t2

private

CRT track t2.

Definition at line 208 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_x1

private

CRT track x1.

Definition at line 199 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_x2

private

CRT track x2.

Definition at line 205 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_y1

private

CRT track y1.

Definition at line 200 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_y2

private

CRT track y2.

Definition at line 206 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_z1

private

CRT track z1.

Definition at line 201 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ctrk_z2

private

CRT track z2.

Definition at line 207 of file HitDumper_module.cc.

int Hitdumper::_event

private

The event number.

Definition at line 161 of file HitDumper_module.cc.

double Hitdumper::_evttime

private

The event time.

Definition at line 162 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_hit_channel

private

Channel where the hit belongs to.

Definition at line 171 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_hit_charge

private

Hit charge.

Definition at line 173 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_hit_cryostat

private

Cryostat where the hit belongs to.

Definition at line 167 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_hit_full_integral

private

Hit charge integral.

Definition at line 176 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_hit_peakT

private

Hit peak time.

Definition at line 172 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_hit_ph

private

Hit ph?

Definition at line 174 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_hit_plane

private

Plane where the hit belongs to.

Definition at line 169 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_hit_tpc

private

TPC where the hit belongs to.

Definition at line 168 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_hit_width

private

Hit width.

Definition at line 175 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_hit_wire

private

Wire where the hit belongs to.

Definition at line 170 of file HitDumper_module.cc.

int Hitdumper::_max_chits

private

maximum number of CRT hits (to be set via fcl)

Definition at line 337 of file HitDumper_module.cc.

int Hitdumper::_max_hits

private

maximum number of hits (to be set via fcl)

Definition at line 334 of file HitDumper_module.cc.

int Hitdumper::_max_nctrks

private

maximum number of CRT tracks (to be set via fcl)

Definition at line 338 of file HitDumper_module.cc.

int Hitdumper::_max_ophits

private

maximum number of hits (to be set via fcl)

Definition at line 335 of file HitDumper_module.cc.

int Hitdumper::_max_samples

private

maximum number of samples (to be set via fcl)

Definition at line 336 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_mhit_channel

private

Channel where the hit belongs to.

Definition at line 279 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_mhit_charge

private

Hit charge.

Definition at line 281 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_mhit_peakT

private

Hit peak time.

Definition at line 280 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_mhit_tpc

private

TPC where the hit belongs to.

Definition at line 277 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_mhit_trk

private

Track number that the hit belongs to.

Definition at line 276 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_mhit_wire

private

Wire where the hit belongs to.

Definition at line 278 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_muontrk_t0

private

t0 (time of interaction)

Definition at line 262 of file HitDumper_module.cc.

std::vector<float> Hitdumper::_muontrk_theta_xz

private

theta_xz trajectory angle

Definition at line 269 of file HitDumper_module.cc.

std::vector<float> Hitdumper::_muontrk_theta_yz

private

theta_yz trajectory angle

Definition at line 270 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_muontrk_tpc

private

tpc that muon is located in

Definition at line 271 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_muontrk_type

private

type of muon track

Definition at line 272 of file HitDumper_module.cc.

std::vector<float> Hitdumper::_muontrk_x1

private

x coordinate closer to anode

Definition at line 263 of file HitDumper_module.cc.

std::vector<float> Hitdumper::_muontrk_x2

private

x coordinate closer to cathode

Definition at line 266 of file HitDumper_module.cc.

std::vector<float> Hitdumper::_muontrk_y1

private

y coordinate closer to anode

Definition at line 264 of file HitDumper_module.cc.

std::vector<float> Hitdumper::_muontrk_y2

private

y coordinate closer to cathode

Definition at line 267 of file HitDumper_module.cc.

std::vector<float> Hitdumper::_muontrk_z1

private

z coordinate closer to anode

Definition at line 265 of file HitDumper_module.cc.

std::vector<float> Hitdumper::_muontrk_z2

private

z coordinate closer to cathode

Definition at line 268 of file HitDumper_module.cc.

int Hitdumper::_nchits

private

Number of CRT hits.

Definition at line 213 of file HitDumper_module.cc.

int Hitdumper::_nctrks

private

Number of created CRT tracks.

Definition at line 198 of file HitDumper_module.cc.

int Hitdumper::_ncts

private

Number of CRT tracks.

Definition at line 222 of file HitDumper_module.cc.

int Hitdumper::_nhits

private

Number of reco hits in the event.

Definition at line 166 of file HitDumper_module.cc.

int Hitdumper::_nmhits

private

Number of muon collection hits per track.

Definition at line 275 of file HitDumper_module.cc.

int Hitdumper::_nmuontrks

private

number of muon tracks

Definition at line 261 of file HitDumper_module.cc.

int Hitdumper::_nophits

private

Number of Optical Hits.

Definition at line 232 of file HitDumper_module.cc.

int Hitdumper::_nstrips

private

Number of CRT strips.

Definition at line 186 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ophit_amplitude

private

Amplitude of the optical hit.

Definition at line 238 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ophit_area

private

Area of the optical hit.

Definition at line 237 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_ophit_opch

private

OpChannel of the optical hit.

Definition at line 233 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_ophit_opdet

private

OpDet of the optical hit.

Definition at line 234 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_ophit_opdet_type

private

OpDet tyoe of the optical hit.

Definition at line 243 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ophit_opdet_x

private

OpDet X coordinate of the optical hit.

Definition at line 240 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ophit_opdet_y

private

OpDet Y coordinate of the optical hit.

Definition at line 241 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ophit_opdet_z

private

OpDet Z coordinate of the optical hit.

Definition at line 242 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ophit_pe

private

PEs of the optical hit.

Definition at line 239 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ophit_peakT

private

Peak time of the optical hit.

Definition at line 235 of file HitDumper_module.cc.

std::vector<double> Hitdumper::_ophit_width

private

Width of the optical hit.

Definition at line 236 of file HitDumper_module.cc.

opdet::sbndPDMapAlg Hitdumper::_pd_map

private

Definition at line 155 of file HitDumper_module.cc.

bool Hitdumper::_pmtSoftTrigger_foundBeamTrigger

private

Definition at line 250 of file HitDumper_module.cc.

int Hitdumper::_pmtSoftTrigger_nAboveThreshold

private

Total photoelectron count before the TTS, during the beam spill.

Definition at line 254 of file HitDumper_module.cc.

double Hitdumper::_pmtSoftTrigger_prelimPE

private

Total photoelectron count 100 ns after the TTS.

Definition at line 253 of file HitDumper_module.cc.

double Hitdumper::_pmtSoftTrigger_promptPE

private

Trigger Time Stamp (TTS), ns (relative to start of beam spill)

Definition at line 252 of file HitDumper_module.cc.

int Hitdumper::_pmtSoftTrigger_tts

private

Whether the beam spill was found or not.

Definition at line 251 of file HitDumper_module.cc.

int Hitdumper::_pmtTrigger_maxpassed

private

maximum number of pmt pairs above threshold during trigger window (usually beam spill)

Definition at line 247 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_pmtTrigger_npmtshigh

private

number of pmt pairs above threshold, index = time during trigger window (usually beam spill)

Definition at line 246 of file HitDumper_module.cc.

int Hitdumper::_run

private

The run number.

Definition at line 159 of file HitDumper_module.cc.

double Hitdumper::_sr_begintime

private

Definition at line 330 of file HitDumper_module.cc.

double Hitdumper::_sr_endtime

private

Definition at line 330 of file HitDumper_module.cc.

double Hitdumper::_sr_pot

private

POTs in each subrun.

Definition at line 331 of file HitDumper_module.cc.

int Hitdumper::_sr_run

private

Definition at line 329 of file HitDumper_module.cc.

int Hitdumper::_sr_subrun

private

Definition at line 329 of file HitDumper_module.cc.

TTree* Hitdumper::_sr_tree

private

A tree filled per subrun (for POT accounting)

Definition at line 328 of file HitDumper_module.cc.

int Hitdumper::_subrun

private

The subrun number.

Definition at line 160 of file HitDumper_module.cc.

int Hitdumper::_t0

private

The t0.

Definition at line 163 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_time_for_waveform

private

Time for waveform to plot.

Definition at line 179 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_waveform_integral

private

Used to see progression of the waveform integral.

Definition at line 181 of file HitDumper_module.cc.

std::vector<int> Hitdumper::_waveform_number

private

Number for each waveform, to allow for searching.

Definition at line 177 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::ccnc_truth

private

0=CC 1=NC

Definition at line 289 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::enu_truth

private

true neutrino energy

Definition at line 291 of file HitDumper_module.cc.

const geo::AuxDetGeometry* Hitdumper::fAuxDetGeo

private

Definition at line 380 of file HitDumper_module.cc.

const geo::AuxDetGeometryCore* Hitdumper::fAuxDetGeoCore

private

Definition at line 381 of file HitDumper_module.cc.

art::ServiceHandle<geo::AuxDetGeometry> Hitdumper::fAuxDetGeoService

private

Definition at line 379 of file HitDumper_module.cc.

bool Hitdumper::fcheckTransparency

private

Checks for wire transprency (to be set via fcl)

Definition at line 367 of file HitDumper_module.cc.

std::string Hitdumper::fCRTHitModuleLabel

private

Label for CRTHit dataproduct (to be set via fcl)

Definition at line 343 of file HitDumper_module.cc.

std::string Hitdumper::fcrtSoftTriggerModuleLabel

private

Label for crt software trigger data product (to be set via fcl)

Definition at line 347 of file HitDumper_module.cc.

std::string Hitdumper::fCRTStripModuleLabel

private

Label for CRTStrip dataproduct (to be set via fcl)

Definition at line 342 of file HitDumper_module.cc.

std::string Hitdumper::fCRTTrackModuleLabel

private

Label for CRTTrack dataproduct (to be set via fcl)

Definition at line 344 of file HitDumper_module.cc.

std::string Hitdumper::fDigitModuleLabel

private

Label for digitizer (to be set via fcl)

Definition at line 349 of file HitDumper_module.cc.

std::string Hitdumper::fGenieGenModuleLabel

private

Label for Genie dataproduct (to be set via fcl)

Definition at line 350 of file HitDumper_module.cc.

geo::GeometryCore const* Hitdumper::fGeometryService

private

Definition at line 377 of file HitDumper_module.cc.

std::string Hitdumper::fHitsModuleLabel

private

Label for Hit dataproduct (to be set via fcl)

Definition at line 340 of file HitDumper_module.cc.

bool Hitdumper::fkeepCRThits

private

Keep the CRT hits (to be set via fcl)

Definition at line 355 of file HitDumper_module.cc.

bool Hitdumper::fkeepCRTstrips

private

Keep the CRT strips (to be set via fcl)

Definition at line 356 of file HitDumper_module.cc.

std::vector<int> Hitdumper::fKeepTaggerTypes = {0, 1, 2, 3, 4, 5, 6}

private

Taggers to keep (to be set via fcl)

Definition at line 373 of file HitDumper_module.cc.

std::string Hitdumper::fLArG4ModuleLabel

private

Label for LArG4 dataproduct (to be set via fcl)

Definition at line 341 of file HitDumper_module.cc.

bool Hitdumper::fmakeCRTtracks

private

Make the CRT tracks (to be set via fcl)

Definition at line 357 of file HitDumper_module.cc.

std::string Hitdumper::fMuonTrackModuleLabel

private

Label for MuonTrack dataproduct (to be set via fcl)

Definition at line 348 of file HitDumper_module.cc.

std::vector<std::string> Hitdumper::fOpHitsModuleLabels

private

Labels for OpHit dataproducts (to be set via fcl)

Definition at line 351 of file HitDumper_module.cc.

std::string Hitdumper::fpmtSoftTriggerModuleLabel

private

Label for pmt software trigger data product (to be set via fcl)

Definition at line 346 of file HitDumper_module.cc.

std::string Hitdumper::fpmtTriggerModuleLabel

private

Label for pmtTrigger dataproduct (to be set vis fcl)

Definition at line 345 of file HitDumper_module.cc.

bool Hitdumper::freadcrtSoftTrigger

private

Add crt software trigger info to output (to be set via fcl)

Definition at line 365 of file HitDumper_module.cc.

bool Hitdumper::freadCRTtracks

private

Keep the CRT tracks (to be set via fcl)

Definition at line 358 of file HitDumper_module.cc.

bool Hitdumper::freadMuonHits

private

Add MuonTrack hits to output(to be set via fcl)

Definition at line 361 of file HitDumper_module.cc.

bool Hitdumper::freadMuonTracks

private

Add MuonTracks to output (to be set via fcl)

Definition at line 360 of file HitDumper_module.cc.

bool Hitdumper::freadOpHits

private

Add OpHits to output (to be set via fcl)

Definition at line 359 of file HitDumper_module.cc.

bool Hitdumper::freadpmtSoftTrigger

private

Add pmt software trigger info to output (to be set via fcl)

Definition at line 364 of file HitDumper_module.cc.

bool Hitdumper::freadpmtTrigger

private

Add pmt hardware trigger info to output (to be set via fcl)

Definition at line 363 of file HitDumper_module.cc.

bool Hitdumper::freadTruth

private

Add Truth info to output (to be set via fcl)

Definition at line 362 of file HitDumper_module.cc.

bool Hitdumper::fsavePOTInfo

private

Add POT info to output (to be set via fcl)

Definition at line 366 of file HitDumper_module.cc.

bool Hitdumper::fSkipInd

private

If true, induction planes are not saved (to be set via fcl)

Definition at line 370 of file HitDumper_module.cc.

TTree* Hitdumper::fTree

private

Definition at line 157 of file HitDumper_module.cc.

bool Hitdumper::fUncompressWithPed

private

Uncompresses the waveforms if true (to be set via fcl)

Definition at line 368 of file HitDumper_module.cc.

int Hitdumper::fWindow

private

Definition at line 369 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::genie_Eng

private

Definition at line 316 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::genie_mass

private

Definition at line 322 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::genie_mother

private

Definition at line 325 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::genie_ND

private

Definition at line 324 of file HitDumper_module.cc.

Int_t Hitdumper::genie_no_primaries

private

Definition at line 314 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::genie_P

private

Definition at line 320 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::genie_primaries_pdg

private

Definition at line 315 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::genie_Px

private

Definition at line 317 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::genie_Py

private

Definition at line 318 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::genie_Pz

private

Definition at line 319 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::genie_status_code

private

Definition at line 321 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::genie_trackID

private

Definition at line 323 of file HitDumper_module.cc.

sbnd::CRTHitRecoAlg Hitdumper::hitAlg

private

Definition at line 375 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::hitnuc_truth

private

hit nucleon

Definition at line 294 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::lep_dcosx_truth

private

lepton dcos x

Definition at line 302 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::lep_dcosy_truth

private

lepton dcos y

Definition at line 303 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::lep_dcosz_truth

private

lepton dcos z

Definition at line 304 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::lep_mom_truth

private

lepton momentum

Definition at line 301 of file HitDumper_module.cc.

size_t Hitdumper::MaxGeniePrimaries = 0

private

Definition at line 313 of file HitDumper_module.cc.

size_t Hitdumper::MaxMCNeutrinos

private

Definition at line 284 of file HitDumper_module.cc.

Int_t Hitdumper::mcevts_truth

private

! The number of MCNeutrinos there is currently room for

number of neutrino Int_teractions in the spill

Definition at line 285 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::mode_truth

private

0=QE/El, 1=RES, 2=DIS, 3=Coherent production

Definition at line 290 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::nu_dcosx_truth

private

neutrino dcos x

Definition at line 298 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::nu_dcosy_truth

private

neutrino dcos y

Definition at line 299 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::nu_dcosz_truth

private

neutrino dcos z

Definition at line 300 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::nuID_truth

private

Unique ID of each true neutrino.

Definition at line 287 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::nuPDG_truth

private

neutrino PDG code

Definition at line 288 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::nuScatterCode_truth

private

Scattering code given by Genie for each neutrino.

Definition at line 286 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::nuvtxx_truth

private

neutrino vertex x

Definition at line 295 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::nuvtxy_truth

private

neutrino vertex y

Definition at line 296 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::nuvtxz_truth

private

neutrino vertex z

Definition at line 297 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::Q2_truth

private

Momentum transfer squared.

Definition at line 292 of file HitDumper_module.cc.

std::vector<Int_t> Hitdumper::tptype_flux

private

Type of parent particle leaving BNB target.

Definition at line 310 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::tpx_flux

private

Px of parent particle leaving BNB target.

Definition at line 307 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::tpy_flux

private

Py of parent particle leaving BNB target.

Definition at line 308 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::tpz_flux

private

Pz of parent particle leaving BNB target.

Definition at line 309 of file HitDumper_module.cc.

std::vector<Float_t> Hitdumper::W_truth

private

hadronic invariant mass

Definition at line 293 of file HitDumper_module.cc.

---

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

• HitDumper_module.cc