pdsp::EMCNNCheck Class Reference

Inheritance diagram for pdsp::EMCNNCheck:

Public Member Functions
	EMCNNCheck (fhicl::ParameterSet const &p)
	EMCNNCheck (EMCNNCheck const &)=delete
	EMCNNCheck (EMCNNCheck &&)=delete
EMCNNCheck &	operator= (EMCNNCheck const &)=delete
EMCNNCheck &	operator= (EMCNNCheck &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
std::string	fGeneratorTag
std::string	fCNNTag
protoana::ProtoDUNEBeamCuts	beam_cuts
protoana::ProtoDUNEBeamlineUtils	fBeamlineUtils
TTree *	ftree
int	run
int	subrun
int	event
int	beampdg
double	average_score_em
double	average_score_trk
double	average_score_mic
double	track_endz
int	ndaughterhits
double	average_daughter_score_mic
double	vtxx
double	vtxy
double	vtxz
double	endx
double	endy
double	endz
double	dirx
double	diry
double	dirz
std::vector< short >	channel
std::vector< short >	tpc
std::vector< short >	plane
std::vector< short >	wire
std::vector< double >	charge
std::vector< double >	peakt
std::vector< double >	score_em
std::vector< double >	score_trk
std::vector< double >	score_mic
std::vector< short >	daughter_channel
std::vector< short >	daughter_tpc
std::vector< short >	daughter_plane
std::vector< short >	daughter_wire
std::vector< double >	daughter_charge
std::vector< double >	daughter_peakt
std::vector< double >	daughter_score_em
std::vector< double >	daughter_score_trk
std::vector< double >	daughter_score_mic
std::vector< int >	pdg
std::vector< int >	origin
std::vector< std::string >	process

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 46 of file EMCNNCheck_module.cc.

Constructor & Destructor Documentation

pdsp::EMCNNCheck::EMCNNCheck ( fhicl::ParameterSet const &  p )

explicit

Definition at line 113 of file EMCNNCheck_module.cc.

  : EDAnalyzer{p},
//fselectpdg(p.get<int>("selectpdg")),
  fGeneratorTag(p.get<std::string>("GeneratorTag")),
  fCNNTag(p.get<std::string>("CNNTag","emtrkmichelid:emtrkmichel")),
  beam_cuts(p.get<fhicl::ParameterSet>("BeamCuts")),
  fBeamlineUtils(p.get<fhicl::ParameterSet>("BeamlineUtils"))
  //BeamCuts(p.get<fhicl::ParameterSet>("BeamCuts"))
  {
    //beam_cuts = protoana::ProtoDUNEBeamCuts(BeamCuts);
  }

pdsp::EMCNNCheck::EMCNNCheck ( EMCNNCheck const &  )

delete

pdsp::EMCNNCheck::EMCNNCheck ( EMCNNCheck &&  )

delete

Member Function Documentation

void pdsp::EMCNNCheck::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 125 of file EMCNNCheck_module.cc.

{

  run = e.run();
  subrun = e.subRun();
  event = e.id().event();
  beampdg = 0;
  average_score_em  = 0.;
  average_score_trk = 0.;
  average_score_mic = 0.;
  track_endz = -1;
  ndaughterhits = 0;
  average_daughter_score_mic = 0.;
  vtxx = -9999;
  vtxy = -9999;
  vtxz = -9999;
  endx = -9999;
  endy = -9999;
  endz = -9999;
  dirx = -9999;
  diry = -9999;
  dirz = -9999;
  channel.clear();
  tpc.clear();
  plane.clear();
  wire.clear();
  charge.clear();
  peakt.clear();
  score_em.clear();
  score_trk.clear();
  score_mic.clear();

  daughter_channel.clear();
  daughter_tpc.clear();
  daughter_plane.clear();
  daughter_wire.clear();
  daughter_charge.clear();
  daughter_peakt.clear();
  daughter_score_em.clear();
  daughter_score_trk.clear();
  daughter_score_mic.clear();

  pdg.clear();
  origin.clear();
  process.clear();

  //Services
  art::ServiceHandle<cheat::BackTrackerService> bt_serv;
  art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
  auto const* SCE = lar::providerFrom<spacecharge::SpaceChargeService>();
  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService>()->DataFor(e);
  auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService>()->DataFor(e, clockData);
  art::ServiceHandle<geo::Geometry const> geom;

  std::vector < art::Ptr < recob::Slice > > sliceList;
  auto sliceListHandle = e.getHandle < std::vector < recob::Slice > > ("pandora");
  if (sliceListHandle) {
    art::fill_ptr_vector(sliceList, sliceListHandle);
  }
  else return;

  // Get all pfparticles
  std::vector < art::Ptr < recob::PFParticle > > pfpList;
  auto pfpListHandle = e.getHandle < std::vector < recob::PFParticle > >("pandora");
  if (pfpListHandle) {
    art::fill_ptr_vector(pfpList, pfpListHandle);
  }

  // Get all clusters
  std::vector < art::Ptr < recob::Cluster > > cluList;
  auto cluListHandle = e.getHandle < std::vector < recob::Cluster > >("pandora");
  if (cluListHandle) {
    art::fill_ptr_vector(cluList, cluListHandle);
  }

  // Get all tracks
  std::vector < art::Ptr < recob::Track > > trkList;
  auto trkListHandle = e.getHandle < std::vector < recob::Track > >("pandoraTrack");
  if (trkListHandle) {
    art::fill_ptr_vector(trkList, trkListHandle);
  }

  // Get all hits
  std::vector < art::Ptr < recob::Hit > > hitList;
  auto hitListHandle = e.getHandle < std::vector < recob::Hit > >("hitpdune");
  if (hitListHandle) {
    art::fill_ptr_vector(hitList, hitListHandle);
  }

  // Get cluster-PFParticle association
  art::FindManyP<recob::Cluster> fmcpfp(pfpListHandle, e, "pandora");

  // Get vertex-PFParticle association
  art::FindManyP<recob::Vertex> fmvpfp(pfpListHandle, e, "pandora");

  // Get hit-cluster association
  art::FindManyP<recob::Hit> fmhc(cluListHandle, e, "pandora");

  art::FindManyP <recob::Hit> hitsFromSlice(sliceListHandle, e, "pandora");

  // Get track-hit association
  art::FindManyP<recob::Hit, recob::TrackHitMeta> fmthm(trkListHandle, e,"pandoraTrack"); // to associate tracks and hits

  // Get hit-track association
  art::FindManyP<recob::Track> thass(hitListHandle, e, "pandoraTrack"); //to associate hit just trying

  anab::MVAReader<recob::Hit,4> hitResults(e, fCNNTag);

  if (!e.isRealData()){
    // Get the truth utility to help us out
    protoana::ProtoDUNETruthUtils truthUtil;
    // Firstly we need to get the list of MCTruth objects from the generator. The standard protoDUNE
    // simulation has fGeneratorTag = "generator"
    auto mcTruths = e.getValidHandle<std::vector<simb::MCTruth>>(fGeneratorTag);
    // mcTruths is basically a pointer to an std::vector of simb::MCTruth objects. There should only be one
    // of these, so we pass the first element into the function to get the good particle
    const simb::MCParticle* geantGoodParticle = truthUtil.GetGeantGoodParticle((*mcTruths)[0],e);
    if(geantGoodParticle != 0x0){
      //std::cout << "Found GEANT particle corresponding to the good particle with pdg = " << geantGoodParticle->PdgCode() << std::endl;
//      if (fselectpdg==211){
//        if (geantGoodParticle->PdgCode()!=211 &&
//            geantGoodParticle->PdgCode()!=13){
//          return;
//        }
//      }
//      else{
//        if (geantGoodParticle->PdgCode()!=fselectpdg) return;
//      }
      beampdg = geantGoodParticle->PdgCode();
    }
  }
  else{
    //Access the Beam Event
    auto beamHandle = e.getValidHandle<std::vector<beam::ProtoDUNEBeamEvent>>("beamevent");
  
    std::vector<art::Ptr<beam::ProtoDUNEBeamEvent>> beamVec;
    if( beamHandle.isValid()){
      art::fill_ptr_vector(beamVec, beamHandle);
    }
    
    if (beamVec.empty()) return;

    const beam::ProtoDUNEBeamEvent & beamEvent = *(beamVec.at(0)); //Should just have one
    /////////////////////////////////////////////////////////////
  
  
    //Check the quality of the event
    std::cout << "Timing Trigger: " << beamEvent.GetTimingTrigger() << std::endl; 
    std::cout << "Is Matched: "     << beamEvent.CheckIsMatched() << std::endl << std::endl;
    
    if( !fBeamlineUtils.IsGoodBeamlineTrigger( e ) ){
      std::cout << "Failed quality check" << std::endl;
      return;
    }
    //Access PID
    std::vector< int > pids = fBeamlineUtils.GetPID( beamEvent, 1. );
//    bool foundparticle = false;
//    for( size_t i = 0; i < pids.size(); ++i ){
//      //std::cout << pids[i] << std::endl;
//      if (pids[i] == fselectpdg) foundparticle = true;
//    }
//    if (!foundparticle) return;
    if (pids.empty()) return;
    beampdg = pids[0];
  }

  if (!beampdg) return;

  //std::cout<<"Found pion"<<std::endl;

  // Get the PFParticle utility
  protoana::ProtoDUNEPFParticleUtils pfpUtil;

  std::vector<const recob::PFParticle*> beamParticles = pfpUtil.GetPFParticlesFromBeamSlice(e,"pandora");

  if(beamParticles.size() == 0){
    std::cerr << "We found no beam particles for this event... moving on" << std::endl;
    return;
  }

  // We can now look at these particles
  int trackid = -1;
  int endwire = -1;
  int endtpc = -1;
  double endpeakt = -1;
  std::vector<int> wirekeys;
  for(const recob::PFParticle* particle : beamParticles){

    const recob::Track* thisTrack = pfpUtil.GetPFParticleTrack(*particle,e,"pandora","pandoraTrack");
    const recob::Shower* thisShower = pfpUtil.GetPFParticleShower(*particle,e,"pandora","pandoraShower");
    if (thisTrack){
      //if (!beam_cuts.IsBeamlike(*thisTrack, e, "1")) return;
      // Track ID
      trackid = thisTrack->ID();
      // Track end point z
      track_endz = thisTrack->End().Z();
      vtxx = thisTrack->Vertex().X();
      vtxy = thisTrack->Vertex().Y();
      vtxz = thisTrack->Vertex().Z();
      if (!geom->FindTPCAtPosition(geo::Point_t(vtxx, vtxy, vtxz)).isValid) return;
      geo::Vector_t offset = {0., 0., 0.};
      if (SCE->EnableCalSpatialSCE()){
        offset = SCE->GetCalPosOffsets(geo::Point_t(vtxx, vtxy, vtxz), (geom->FindTPCAtPosition(geo::Point_t(vtxx, vtxy, vtxz))).TPC);
      }
      //std::cout<<"track "<<offset.X()<<" "<<offset.Y()<<" "<<offset.Z()<<std::endl;
      vtxx -= offset.X();
      vtxy += offset.Y();
      vtxz += offset.Z();
      endx = thisTrack->End().X();
      endy = thisTrack->End().Y();
      endz = thisTrack->End().Z();
      if (!geom->FindTPCAtPosition(geo::Point_t(endx, endy, endz)).isValid) return;
      offset = {0., 0., 0.};
      if (SCE->EnableCalSpatialSCE()){
        offset = SCE->GetCalPosOffsets(geo::Point_t(endx, endy, endz), (geom->FindTPCAtPosition(geo::Point_t(endx, endy, endz))).TPC);
      }
      endx -= offset.X();
      endy += offset.Y();
      endz += offset.Z();
      TVector3 dir(endx-vtxx, endy-vtxy, endz-vtxz);
      dir = dir.Unit();
      dirx = dir.X();
      diry = dir.Y();
      dirz = dir.Z();
      // Find the last wire number and peak time on the track
      if (fmthm.isValid()){
        float zlast0=-99999;
        auto vhit=fmthm.at(trackid);
        auto vmeta=fmthm.data(trackid);
        for (size_t ii = 0; ii<vhit.size(); ++ii){ //loop over all meta data hit
          bool fBadhit = false;
          if (vmeta[ii]->Index() == static_cast<unsigned int>(std::numeric_limits<int>::max())){
            fBadhit = true;
            continue;
          }
          if (vmeta[ii]->Index()>=thisTrack->NumberTrajectoryPoints()){
            throw cet::exception("Calorimetry_module.cc") << "Requested track trajectory index "<<vmeta[ii]->Index()<<" exceeds the total number of trajectory points "<<thisTrack->NumberTrajectoryPoints()<<" for track index "<<trackid<<". Something is wrong with the track reconstruction. Please contact tjyang@fnal.gov!!";
          }
          if (!thisTrack->HasValidPoint(vmeta[ii]->Index())){
            fBadhit = true;
            continue;
          }
          auto loc = thisTrack->LocationAtPoint(vmeta[ii]->Index());
          if (fBadhit) continue; //HY::If BAD hit, skip this hit and go next
          if (loc.Z()<-100) continue; //hit not on track
          if(vhit[ii]->WireID().Plane==2){
            wirekeys.push_back(vhit[ii].key());
            float zlast=loc.Z();
            if(zlast>zlast0){
              zlast0=zlast;
              endwire=vhit[ii]->WireID().Wire;
              endpeakt=vhit[ii]->PeakTime();
              endtpc=vhit[ii]->WireID().TPC;
            }
          }
        }
      }
    }
    if (thisShower){
      //if (!beam_cuts.IsBeamlike(*thisShower, e, "1")) return;
      vtxx = thisShower->ShowerStart().X();
      vtxy = thisShower->ShowerStart().Y();
      vtxz = thisShower->ShowerStart().Z();
      endx = vtxx + thisShower->Direction().X();
      endy = vtxy + thisShower->Direction().Y();
      endz = vtxz + thisShower->Direction().Z();
      if (!geom->FindTPCAtPosition(geo::Point_t(vtxx, vtxy, vtxz)).isValid) return;
      geo::Vector_t offset = {0., 0., 0.};
      if (SCE->EnableCalSpatialSCE()){
        offset = SCE->GetCalPosOffsets(geo::Point_t(vtxx, vtxy, vtxz), (geom->FindTPCAtPosition(geo::Point_t(vtxx, vtxy, vtxz))).TPC);
      }
      //std::cout<<"shower "<<offset.X()<<" "<<offset.Y()<<" "<<offset.Z()<<std::endl;
      vtxx -= offset.X();
      vtxy += offset.Y();
      vtxz += offset.Z();
      if (!geom->FindTPCAtPosition(geo::Point_t(endx, endy, endz)).isValid) return;
      offset = {0., 0., 0.};
      if (SCE->EnableCalSpatialSCE()){
        offset = SCE->GetCalPosOffsets(geo::Point_t(endx, endy, endz), (geom->FindTPCAtPosition(geo::Point_t(endx, endy, endz))).TPC);
      }
      endx -= offset.X();
      endy += offset.Y();
      endz += offset.Z();
      TVector3 dir(endx-vtxx, endy-vtxy, endz-vtxz);
      dir = dir.Unit();
      dirx = dir.X();
      diry = dir.Y();
      dirz = dir.Z();
    }
  }

  //int sliceid = pfpUtil.GetBeamSlice(e, "pandora");
  
  //if (sliceid!=9999){
  if (fmcpfp.isValid()){
    // Get clusters associated with pfparticle
    auto const& clusters = fmcpfp.at(beamParticles[0]->Self());
    for (auto const & cluster : clusters){
      if (fmhc.isValid()){
        // Get hits associated with cluster
        auto const& hits = fmhc.at(cluster.key());
        //auto const& hits = hitsFromSlice.at(sliceid);
        //std::cout<<hits.size()<<std::endl;
        for (auto & hit : hits){
          std::array<float,4> cnn_out = hitResults.getOutput(hit);
          if (hit->WireID().Plane == 2){
            channel.push_back(hit->Channel());
            tpc.push_back(hit->WireID().TPC);
            plane.push_back(hit->WireID().Plane);
            wire.push_back(hit->WireID().Wire);
            charge.push_back(hit->Integral());
            peakt.push_back(hit->PeakTime());     
            score_em.push_back(cnn_out[hitResults.getIndex("em")]);
            score_trk.push_back(cnn_out[hitResults.getIndex("track")]);
            score_mic.push_back(cnn_out[hitResults.getIndex("michel")]);
            int this_pdg = 0;
            int this_origin = -1;
            std::string this_process = "null";
            if (!e.isRealData()){
              int TrackID = 0;
              std::map<int,double> trkide;
              std::vector<sim::TrackIDE> TrackIDs = bt_serv->HitToEveTrackIDEs(clockData, hit);
              for(size_t e = 0; e < TrackIDs.size(); ++e){
                trkide[TrackIDs[e].trackID] += TrackIDs[e].energy;
              }     
              
              // Work out which IDE despoited the most charge in the hit if there was more than one.
              double maxe = -1;
              double tote = 0;
              for (std::map<int,double>::iterator ii = trkide.begin(); ii!=trkide.end(); ++ii){
                tote += ii->second;
                if ((ii->second)>maxe){
                  maxe = ii->second;
                  TrackID = ii->first;
                }
              }
              // Now have trackID, so get PdG code and T0 etc.
              const simb::MCParticle *particle = pi_serv->TrackIdToParticle_P(TrackID);
              if (particle){
                this_pdg = particle->PdgCode();
                this_process = particle->Process();
                this_origin = pi_serv->ParticleToMCTruth_P(particle)->Origin();
              }
            }
            pdg.push_back(this_pdg);
            origin.push_back(this_origin);
            process.push_back(this_process);
          }
        }
      }
    }
  }

  // Get the average of the collection plane scores
  unsigned int nCollectionHits = 0;
  for(unsigned int h = 0; h < plane.size(); ++h){
    if(plane.at(h) == 2){
      ++nCollectionHits;
      average_score_em += score_em.at(h);
      average_score_trk += score_trk.at(h);
      average_score_mic += score_mic.at(h);
    }
  }
  if(nCollectionHits > 0){
    average_score_em /= static_cast<double>(nCollectionHits);
    average_score_trk /= static_cast<double>(nCollectionHits);
    average_score_mic /= static_cast<double>(nCollectionHits);
  }

  // Get the hits near the track end (Michel candidate)
  //std::cout<<trackid<<std::endl;
  if (trackid!=-1){
    for(size_t hitl=0;hitl<hitList.size();hitl++){
      std::array<float,4> cnn_out=hitResults.getOutput(hitList[hitl]);
      auto & tracks = thass.at(hitList[hitl].key());
      // hit not on the track
      if (std::find(wirekeys.begin(), wirekeys.end(), hitl) != wirekeys.end()) continue;
      // hit not on a long track
      if (!tracks.empty() && int(tracks[0].key()) != trackid && trkList[tracks[0].key()]->Length()>25) continue;
      int planeid=hitList[hitl]->WireID().Plane;
      if (planeid!=2) continue;
      int tpcid=hitList[hitl]->WireID().TPC;
      if (tpcid!=endtpc) continue;
      float peakth1=hitList[hitl]->PeakTime();
      int wireh1=hitList[hitl]->WireID().Wire;
      if(std::abs(wireh1-endwire)<15 && std::abs(peakth1-endpeakt)<100 && tpcid==endtpc){
        daughter_channel.push_back(hitList[hitl]->Channel());
        daughter_tpc.push_back(hitList[hitl]->WireID().TPC);
        daughter_plane.push_back(hitList[hitl]->WireID().Plane);
        daughter_wire.push_back(hitList[hitl]->WireID().Wire);
        daughter_charge.push_back(hitList[hitl]->Integral());
        daughter_peakt.push_back(hitList[hitl]->PeakTime());     
        daughter_score_em.push_back(cnn_out[hitResults.getIndex("em")]);
        daughter_score_trk.push_back(cnn_out[hitResults.getIndex("track")]);
        daughter_score_mic.push_back(cnn_out[hitResults.getIndex("michel")]);
        
        ++ndaughterhits;
        average_daughter_score_mic += cnn_out[hitResults.getIndex("michel")];
        //std::cout<<hitList[hitl]->WireID().Wire<<" "<<hitList[hitl]->PeakTime()<<" "<<hitList[hitl]->Integral()<<" "<<cnn_out[hitResults.getIndex("michel")]<<std::endl;
      }
    }
  }
  if (ndaughterhits) average_daughter_score_mic /= ndaughterhits;
  
  if (!channel.empty()) ftree->Fill();
}

void pdsp::EMCNNCheck::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 532 of file EMCNNCheck_module.cc.

                             {

  art::ServiceHandle<art::TFileService> fileServiceHandle;
  ftree = fileServiceHandle->make<TTree>("ftree", "hit info");
  ftree->Branch("run", &run, "run/I");
  ftree->Branch("event", &event, "event/I");
  ftree->Branch("beampdg", &beampdg, "beampdg/I");
  ftree->Branch("average_score_em" , &average_score_em , "average_score_em/D");
  ftree->Branch("average_score_trk", &average_score_trk, "average_score_trk/D");
  ftree->Branch("average_score_mic", &average_score_mic, "average_score_mic/D");
  ftree->Branch("track_endz", &track_endz, "track_endz/D");
  ftree->Branch("ndaughterhits", &ndaughterhits, "ndaughterhits/I");
  ftree->Branch("average_daughter_score_mic", &average_daughter_score_mic, "average_daughter_score_mic/D");
  ftree->Branch("vtxx", &vtxx, "vtxx/D");
  ftree->Branch("vtxy", &vtxy, "vtxy/D");
  ftree->Branch("vtxz", &vtxz, "vtxz/D");
  ftree->Branch("endx", &endx, "endx/D");
  ftree->Branch("endy", &endy, "endy/D");
  ftree->Branch("endz", &endz, "endz/D");
  ftree->Branch("dirx", &dirx, "dirx/D");
  ftree->Branch("diry", &diry, "diry/D");
  ftree->Branch("dirz", &dirz, "dirz/D");
  ftree->Branch("channel", &channel);
  ftree->Branch("tpc", &tpc);
  ftree->Branch("plane", &plane);
  ftree->Branch("wire", &wire);
  ftree->Branch("charge", &charge);
  ftree->Branch("peakt", &peakt);
  ftree->Branch("score_em", &score_em);
  ftree->Branch("score_trk", &score_trk);
  ftree->Branch("score_mic", &score_mic);

  ftree->Branch("daughter_channel", &daughter_channel);
  ftree->Branch("daughter_tpc", &daughter_tpc);
  ftree->Branch("daughter_plane", &daughter_plane);
  ftree->Branch("daughter_wire", &daughter_wire);
  ftree->Branch("daughter_charge", &daughter_charge);
  ftree->Branch("daughter_peakt", &daughter_peakt);
  ftree->Branch("daughter_score_em", &daughter_score_em);
  ftree->Branch("daughter_score_trk", &daughter_score_trk);
  ftree->Branch("daughter_score_mic", &daughter_score_mic);

  ftree->Branch("pdg", &pdg);
  ftree->Branch("origin", &origin);
  ftree->Branch("process", &process);

}

EMCNNCheck& pdsp::EMCNNCheck::operator= ( EMCNNCheck const &  )

delete

EMCNNCheck& pdsp::EMCNNCheck::operator= ( EMCNNCheck &&  )

delete

Member Data Documentation

double pdsp::EMCNNCheck::average_daughter_score_mic

private

Definition at line 82 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::average_score_em

private

Definition at line 77 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::average_score_mic

private

Definition at line 79 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::average_score_trk

private

Definition at line 78 of file EMCNNCheck_module.cc.

protoana::ProtoDUNEBeamCuts pdsp::EMCNNCheck::beam_cuts

private

Definition at line 69 of file EMCNNCheck_module.cc.

int pdsp::EMCNNCheck::beampdg

private

Definition at line 76 of file EMCNNCheck_module.cc.

std::vector<short> pdsp::EMCNNCheck::channel

private

Definition at line 86 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::charge

private

Definition at line 90 of file EMCNNCheck_module.cc.

std::vector<short> pdsp::EMCNNCheck::daughter_channel

private

Definition at line 96 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::daughter_charge

private

Definition at line 100 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::daughter_peakt

private

Definition at line 101 of file EMCNNCheck_module.cc.

std::vector<short> pdsp::EMCNNCheck::daughter_plane

private

Definition at line 98 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::daughter_score_em

private

Definition at line 102 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::daughter_score_mic

private

Definition at line 104 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::daughter_score_trk

private

Definition at line 103 of file EMCNNCheck_module.cc.

std::vector<short> pdsp::EMCNNCheck::daughter_tpc

private

Definition at line 97 of file EMCNNCheck_module.cc.

std::vector<short> pdsp::EMCNNCheck::daughter_wire

private

Definition at line 99 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::dirx

private

Definition at line 85 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::diry

private

Definition at line 85 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::dirz

private

Definition at line 85 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::endx

private

Definition at line 84 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::endy

private

Definition at line 84 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::endz

private

Definition at line 84 of file EMCNNCheck_module.cc.

int pdsp::EMCNNCheck::event

private

Definition at line 75 of file EMCNNCheck_module.cc.

protoana::ProtoDUNEBeamlineUtils pdsp::EMCNNCheck::fBeamlineUtils

private

Definition at line 70 of file EMCNNCheck_module.cc.

std::string pdsp::EMCNNCheck::fCNNTag

private

Definition at line 67 of file EMCNNCheck_module.cc.

std::string pdsp::EMCNNCheck::fGeneratorTag

private

Definition at line 66 of file EMCNNCheck_module.cc.

TTree* pdsp::EMCNNCheck::ftree

private

Definition at line 72 of file EMCNNCheck_module.cc.

int pdsp::EMCNNCheck::ndaughterhits

private

Definition at line 81 of file EMCNNCheck_module.cc.

std::vector<int> pdsp::EMCNNCheck::origin

private

Definition at line 107 of file EMCNNCheck_module.cc.

std::vector<int> pdsp::EMCNNCheck::pdg

private

Definition at line 106 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::peakt

private

Definition at line 91 of file EMCNNCheck_module.cc.

std::vector<short> pdsp::EMCNNCheck::plane

private

Definition at line 88 of file EMCNNCheck_module.cc.

std::vector<std::string> pdsp::EMCNNCheck::process

private

Definition at line 108 of file EMCNNCheck_module.cc.

int pdsp::EMCNNCheck::run

private

Definition at line 73 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::score_em

private

Definition at line 92 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::score_mic

private

Definition at line 94 of file EMCNNCheck_module.cc.

std::vector<double> pdsp::EMCNNCheck::score_trk

private

Definition at line 93 of file EMCNNCheck_module.cc.

int pdsp::EMCNNCheck::subrun

private

Definition at line 74 of file EMCNNCheck_module.cc.

std::vector<short> pdsp::EMCNNCheck::tpc

private

Definition at line 87 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::track_endz

private

Definition at line 80 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::vtxx

private

Definition at line 83 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::vtxy

private

Definition at line 83 of file EMCNNCheck_module.cc.

double pdsp::EMCNNCheck::vtxz

private

Definition at line 83 of file EMCNNCheck_module.cc.

std::vector<short> pdsp::EMCNNCheck::wire

private

Definition at line 89 of file EMCNNCheck_module.cc.

---

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

• protoduneana/protoduneana/singlephase/EMTaskForce/EMCNNCheck_module.cc