Inheritance diagram for nnet::PointIdEffTest:

Classes
struct	Config

Public Types
enum	EId { kShower = 0, kTrack = 1, kMichel = 2 }

using	Parameters = art::EDAnalyzer::Table< Config >

Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Public Member Functions
	PointIdEffTest (Parameters const &config)

	PointIdEffTest (PointIdEffTest const &)=delete

	PointIdEffTest (PointIdEffTest &&)=delete

PointIdEffTest &	operator= (PointIdEffTest const &)=delete

PointIdEffTest &	operator= (PointIdEffTest &&)=delete

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 Member Functions
void	beginRun (const art::Run &run) override

void	beginJob () override

void	endJob () override

void	analyze (art::Event const &e) override

void	cleanup ()

void	countTruthDep (const std::vector< sim::SimChannel > &channels, float &emLike, float &trackLike) const

void	countPfpDep (detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, const std::vector< recob::PFParticle > &pfparticles, const art::FindManyP< recob::Cluster > &pfpclus, const art::FindManyP< recob::Hit > &cluhits, float &emLike, float &trackLike) const

bool	isMuonDecaying (const simb::MCParticle &particle, const std::unordered_map< int, const simb::MCParticle * > &particleMap) const

int	testCNN (detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, const std::vector< sim::SimChannel > &channels, const std::vector< art::Ptr< recob::Hit >> &hits, const std::array< float, MVA_LENGTH > &cnn_out, const std::vector< anab::FeatureVector< MVA_LENGTH >> &hit_outs, size_t cidx)

Private Attributes
int	fRun

int	fEvent

float	fMcDepEM

float	fMcDepTrack

float	fMcFractionEM

float	fPfpDepEM

float	fPfpDepTrack

float	fHitEM_0p5

float	fHitTrack_0p5

float	fHitMichel_0p5

float	fHitMcFractionEM

float	fHitEM_mc

float	fpEM

float	fHitMichel_mc

float	fpMichel_hit

float	fpMichel_cl

float	fOutTrk

float	fOutEM

float	fOutNone

float	fHitEM_0p85

float	fHitTrack_0p85

float	fTotHit

float	fCleanHit

float	fHitsEM_OK_0p5 [100]

float	fHitsTrack_OK_0p5 [100]

float	fHitsMichel_OK_0p5 [100]

float	fHitsMichel_False_0p5 [100]

float	fHitsEM_OK_0p85 [100]

float	fHitsTrack_OK_0p85 [100]

float	fHitRecoEM [100]

float	fHitRecoFractionEM [100]

int	fMcPid

int	fClSize

float	fPidValue

int	fTrkOk [100]

int	fTrkBad [100]

int	fShOk [100]

int	fShBad [100]

int	fNone

int	fTotal

double	fElectronsToGeV

int	fTrkLikeIdx

int	fEmLikeIdx

int	fNoneIdx

int	fMichelLikeIdx

TTree *	fEventTree

TTree *	fClusterTree

TTree *	fHitTree

std::ofstream	fHitsOutFile

unsigned int	fView

std::unordered_map< int, const simb::MCParticle * >	fParticleMap

calo::CalorimetryAlg	fCalorimetryAlg

art::InputTag	fSimulationProducerLabel

art::InputTag	fPfpModuleLabel

art::InputTag	fNNetModuleLabel

bool	fSaveHitsFile

Additional Inherited Members
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 55 of file PointIdEffTest_module.cc.

Member Typedef Documentation

using nnet::PointIdEffTest::Parameters = art::EDAnalyzer::Table<Config>

Definition at line 81 of file PointIdEffTest_module.cc.

Member Enumeration Documentation

enum nnet::PointIdEffTest::EId

Enumerator
kShower
kTrack
kMichel

Definition at line 57 of file PointIdEffTest_module.cc.

57 { kShower = 0, kTrack = 1, kMichel = 2 };

nnet::PointIdEffTest::kTrack

Definition: PointIdEffTest_module.cc:57

nnet::PointIdEffTest::kMichel

Definition: PointIdEffTest_module.cc:57

nnet::PointIdEffTest::kShower

Definition: PointIdEffTest_module.cc:57

Constructor & Destructor Documentation

nnet::PointIdEffTest::PointIdEffTest ( Parameters const & config )

explicit

Definition at line 166 of file PointIdEffTest_module.cc.

   : art::EDAnalyzer(config)
   , fMcPid(-1)
   , fClSize(0)
   , fTrkLikeIdx(-1)
   , fEmLikeIdx(-1)
   , fNoneIdx(-1)
   , fMichelLikeIdx(-1)
   , fView(config().View())
   , fCalorimetryAlg(config().CalorimetryAlg())
   , fSimulationProducerLabel(config().SimModuleLabel())
   , fPfpModuleLabel(config().PfpModuleLabel())
   , fNNetModuleLabel(config().NNetModuleLabel())
   , fSaveHitsFile(config().SaveHitsFile())
 {}

nnet::PointIdEffTest::PointIdEffTest ( PointIdEffTest const & )

delete

nnet::PointIdEffTest::PointIdEffTest ( PointIdEffTest && )

delete

Member Function Documentation

void nnet::PointIdEffTest::analyze ( art::Event const & e )

overrideprivatevirtual

Implements art::EDAnalyzer.

Definition at line 311 of file PointIdEffTest_module.cc.

 {
   cleanup(); // remove everything from members
 
   fRun = e.run();
   fEvent = e.id().event();
 
   // access to MC information
 
   // MC particles list
   auto particleHandle = e.getValidHandle<std::vector<simb::MCParticle>>(fSimulationProducerLabel);
   for (auto const& particle : *particleHandle) {
     fParticleMap[particle.TrackId()] = &particle;
   }
 
   // SimChannels
   auto simChannelHandle = e.getValidHandle<std::vector<sim::SimChannel>>(fSimulationProducerLabel);
   countTruthDep(*simChannelHandle, fMcDepEM, fMcDepTrack);
 
   auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
   auto const detProp =
     art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(e, clockData);
 
   // PFParticle selection results
   art::Handle<std::vector<recob::PFParticle>> pfpHandle;
   if (e.getByLabel(fPfpModuleLabel, pfpHandle)) {
     auto cluHandle = e.getValidHandle<std::vector<recob::Cluster>>(fPfpModuleLabel);
     const art::FindManyP<recob::Cluster> clusFromPfps(pfpHandle, e, fPfpModuleLabel);
     const art::FindManyP<recob::Hit> hitsFromClus(cluHandle, e, fPfpModuleLabel);
     countPfpDep(
       clockData, detProp, *pfpHandle, clusFromPfps, hitsFromClus, fPfpDepEM, fPfpDepTrack);
   }
 
   // output from cnn's
 
   anab::MVAReader<recob::Hit, MVA_LENGTH> hitResults(
     e, fNNetModuleLabel); // hit-by-hit outpus just to be dumped to file for debugging
   fTrkLikeIdx = hitResults.getIndex("track");
   fEmLikeIdx = hitResults.getIndex("em");
   fNoneIdx = hitResults.getIndex("none");
   fMichelLikeIdx = hitResults.getIndex("michel");
   if ((fTrkLikeIdx < 0) || (fEmLikeIdx < 0)) {
     throw cet::exception("PointIdEffTest")
       << "No em/track labeled columns in MVA data products." << std::endl;
   }
 
   auto cluResults = anab::MVAReader<recob::Cluster, MVA_LENGTH>::create(
     e, fNNetModuleLabel); // outputs for clusters recovered in not-throwing way
   if (cluResults) {
     const art::FindManyP<recob::Hit> hitsFromClusters(
       cluResults->dataHandle(), e, cluResults->dataTag());
 
     for (size_t c = 0; c < cluResults->size(); ++c) {
       const recob::Cluster& clu = cluResults->item(c);
 
       if (clu.Plane().Plane != fView) continue;
 
       const std::vector<art::Ptr<recob::Hit>>& hits = hitsFromClusters.at(c);
       std::array<float, MVA_LENGTH> cnn_out = cluResults->getOutput(c);
 
       testCNN(clockData,
               detProp,
               *simChannelHandle,
               hits,
               cnn_out,
               hitResults.outputs(),
               c); // test hits in the cluster
     }
 
     if (fTotHit > 0)
       fCleanHit = fCleanHit / fTotHit;
     else
       fCleanHit = 0;
 
     double totMcDep = fMcDepEM + fMcDepTrack;
     if (totMcDep)
       fMcFractionEM = fMcDepEM / totMcDep;
     else
       fMcFractionEM = 0;
 
     double totEmTrk0p5 = fHitEM_0p5 + fHitTrack_0p5;
     if (totEmTrk0p5 > 0)
       fHitMcFractionEM = fHitEM_0p5 / totEmTrk0p5;
     else
       fHitMcFractionEM = 0;
 
     for (size_t i = 0; i < 100; ++i) {
       if (fHitEM_0p5 > 0)
         fHitsEM_OK_0p5[i] /= fHitEM_0p5;
       else
         fHitsEM_OK_0p5[i] = 0;
 
       if (fHitTrack_0p5 > 0)
         fHitsTrack_OK_0p5[i] /= fHitTrack_0p5;
       else
         fHitsTrack_OK_0p5[i] = 0;
 
       if (fHitEM_0p85 > 0)
         fHitsEM_OK_0p85[i] /= fHitEM_0p85;
       else
         fHitsEM_OK_0p85[i] = 0;
 
       if (fHitTrack_0p85 > 0)
         fHitsTrack_OK_0p85[i] /= fHitTrack_0p85;
       else
         fHitsTrack_OK_0p85[i] = 0;
 
       if (totEmTrk0p5 > 0)
         fHitRecoFractionEM[i] = fHitRecoEM[i] / totEmTrk0p5;
       else
         fHitRecoFractionEM[i] = 0;
     }
 
     fEventTree->Fill();
   }
   else {
     mf::LogWarning("TrainingDataAlg") << "MVA FOR CLUSTERS NOT FOUND";
   }
 
   cleanup(); // remove everything from members
 }

void nnet::PointIdEffTest::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 190 of file PointIdEffTest_module.cc.

 {
   art::ServiceHandle<art::TFileService const> tfs;
 
   fEventTree = tfs->make<TTree>("event", "event info");
   fEventTree->Branch("fRun", &fRun, "fRun/I");
   fEventTree->Branch("fEvent", &fEvent, "fEvent/I");
   fEventTree->Branch("fMcDepEM", &fMcDepEM, "fMcDepEM/F");
   fEventTree->Branch("fMcDepTrack", &fMcDepTrack, "fMcDepTrack/F");
   fEventTree->Branch("fMcFractionEM", &fMcFractionEM, "fMcFractionEM/F");
   fEventTree->Branch("fPfpDepEM", &fPfpDepEM, "fPfpDepEM/F");
   fEventTree->Branch("fPfpDepTrack", &fPfpDepTrack, "fPfpDepTrack/F");
   fEventTree->Branch("fHitEM_0p5", &fHitEM_0p5, "fHitEM_0p5/F");
   fEventTree->Branch("fHitMichel_0p5", &fHitMichel_0p5, "fHitMichel_0p5/F");
   fEventTree->Branch("fHitTrack_0p5", &fHitTrack_0p5, "fHitTrack_0p5/F");
   fEventTree->Branch("fHitMcFractionEM", &fHitMcFractionEM, "fHitMcFractionEM/F");
   fEventTree->Branch("fHitEM_0p85", &fHitEM_0p85, "fHitEM_0p85/F");
   fEventTree->Branch("fHitTrack_0p85", &fHitTrack_0p85, "fHitTrack_0p85/F");
   fEventTree->Branch("fCleanHit", &fCleanHit, "fCleanHit/F");
 
   fEventTree->Branch("fHitsEM_OK_0p5", fHitsEM_OK_0p5, "fHitsEM_OK_0p5[100]/F");
   fEventTree->Branch("fHitsTrack_OK_0p5", fHitsTrack_OK_0p5, "fHitsTrack_OK_0p5[100]/F");
   fEventTree->Branch("fHitsMichel_OK_0p5", fHitsMichel_OK_0p5, "fHitsMichel_OK_0p5[100]/F");
   fEventTree->Branch(
     "fHitsMichel_False_0p5", fHitsMichel_False_0p5, "fHitsMichel_False_0p5[100]/F");
 
   fEventTree->Branch("fHitsEM_OK_0p85", fHitsEM_OK_0p85, "fHitsEM_OK_0p85[100]/F");
   fEventTree->Branch("fHitsTrack_OK_0p85", fHitsTrack_OK_0p85, "fHitsTrack_OK_0p85[100]/F");
 
   fEventTree->Branch("fHitRecoEM", fHitRecoEM, "fHitRecoEM[100]/F");
   fEventTree->Branch("fHitRecoFractionEM", fHitRecoFractionEM, "fHitRecoFractionEM[100]/F");
 
   fClusterTree = tfs->make<TTree>("cluster", "clusters info");
   fClusterTree->Branch("fMcPid", &fMcPid, "fMcPid/I");
   fClusterTree->Branch("fClSize", &fClSize, "fClSize/I");
   fClusterTree->Branch("fPidValue", &fPidValue, "fPidValue/F");
   fClusterTree->Branch("fpMichel_cl", &fpMichel_cl, "fpMichel_cl/F");
 
   fHitTree = tfs->make<TTree>("hits", "hits info");
   fHitTree->Branch("fRun", &fRun, "fRun/I");
   fHitTree->Branch("fEvent", &fEvent, "fEvent/I");
   fHitTree->Branch("fHitEM_mc", &fHitEM_mc, "fHitEM_mc/F");
   fHitTree->Branch("fpEM", &fpEM, "fpEM/F");
   fHitTree->Branch("fPidValue", &fPidValue, "fPidValue/F");
   fHitTree->Branch("fHitMichel_mc", &fHitMichel_mc, "fHitMichel_mc/F");
   fHitTree->Branch("fpMichel_hit", &fpMichel_hit, "fpMichel_hit/F");
   fHitTree->Branch("fOutTrk", &fOutTrk, "fOutTrk/F");
   fHitTree->Branch("fOutEM", &fOutEM, "fOutEM/F");
   fHitTree->Branch("fOutNone", &fOutNone, "fOutNone/F");
 
   if (fSaveHitsFile) fHitsOutFile.open("hits_pid.prn");
 
   fNone = 0;
   fTotal = 0;
   for (size_t i = 0; i < 100; ++i) {
     fShOk[i] = 0;
     fShBad[i] = 0;
     fTrkOk[i] = 0;
     fTrkBad[i] = 0;
   }
 }

void nnet::PointIdEffTest::beginRun ( const art::Run & run )

overrideprivate

Definition at line 183 of file PointIdEffTest_module.cc.

 {
   art::ServiceHandle<sim::LArG4Parameters const> larParameters;
   fElectronsToGeV = 1. / larParameters->GeVToElectrons();
 }

void nnet::PointIdEffTest::cleanup ( )

private

Definition at line 275 of file PointIdEffTest_module.cc.

 {
   fParticleMap.clear();
 
   fMcDepEM = 0;
   fMcDepTrack = 0;
   fMcFractionEM = 0;
   fPfpDepEM = 0;
   fPfpDepTrack = 0;
   fHitEM_0p5 = 0;
   fHitTrack_0p5 = 0;
   fHitEM_0p85 = 0;
   fHitTrack_0p85 = 0;
   fHitMichel_0p5 = 0;
   fHitMcFractionEM = 0;
   fTotHit = 0;
   fCleanHit = 0;
 
   for (size_t i = 0; i < 100; ++i) {
     fHitsEM_OK_0p5[i] = 0;
     fHitsTrack_OK_0p5[i] = 0;
     fHitsEM_OK_0p85[i] = 0;
     fHitsTrack_OK_0p85[i] = 0;
     fHitsMichel_OK_0p5[i] = 0;
     fHitsMichel_False_0p5[i] = 0;
     fHitRecoEM[i] = 0;
     fHitRecoFractionEM[i] = 0;
   }
 
   fTrkLikeIdx = -1;
   fEmLikeIdx = -1;
   fNoneIdx = -1;
   fMichelLikeIdx = -1;
 }

void nnet::PointIdEffTest::countPfpDep	(	detinfo::DetectorClocksData const &	clockData,
		detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< recob::PFParticle > &	pfparticles,
		const art::FindManyP< recob::Cluster > &	pfpclus,
		const art::FindManyP< recob::Hit > &	cluhits,
		float &	emLike,
		float &	trackLike
	)		const

private

Definition at line 479 of file PointIdEffTest_module.cc.

 {
   emLike = 0;
   trackLike = 0;
   for (size_t i = 0; i < pfparticles.size(); ++i) {
     const auto& pfp = pfparticles[i];
     const auto& clus = pfpclus.at(i);
 
     float hitdep = 0;
     for (const auto& c : clus) {
       const auto& hits = cluhits.at(c.key());
       for (const auto& h : hits) {
         if (h->View() == fView) {
           hitdep +=
             h->SummedADC() * fCalorimetryAlg.LifetimeCorrection(clockData, detProp, h->PeakTime());
         }
       }
     }
 
     if ((pfp.PdgCode() == 11) || pfp.PdgCode() == -11) { emLike += hitdep; }
     else {
       trackLike += hitdep;
     }
   }
 }

void nnet::PointIdEffTest::countTruthDep	(	const std::vector< sim::SimChannel > &	channels,
		float &	emLike,
		float &	trackLike
	)		const

private

Definition at line 435 of file PointIdEffTest_module.cc.

 {
   emLike = 0;
   trackLike = 0;
   for (auto const& channel : channels) {
     // for every time slice in this channel:
     auto const& timeSlices = channel.TDCIDEMap();
     for (auto const& timeSlice : timeSlices) {
       // loop over the energy deposits.
       auto const& energyDeposits = timeSlice.second;
       for (auto const& energyDeposit : energyDeposits) {
         int trackID = energyDeposit.trackID;
 
         double energy = energyDeposit.numElectrons * fElectronsToGeV * 1000;
 
         if (trackID < 0) { emLike += energy; }
         else if (trackID > 0) {
           auto search = fParticleMap.find(trackID);
           bool found = true;
           if (search == fParticleMap.end()) {
             mf::LogWarning("TrainingDataAlg") << "PARTICLE NOT FOUND";
             found = false;
           }
 
           int pdg = 0;
           if (found) {
             const simb::MCParticle& particle = *((*search).second);
             if (!pdg) pdg = particle.PdgCode(); // not EM activity so read what PDG it is
           }
 
           if ((pdg == 11) || (pdg == -11) || (pdg == 22))
             emLike += energy;
           else
             trackLike += energy;
         }
       }
     }
   }
 }

void nnet::PointIdEffTest::endJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 253 of file PointIdEffTest_module.cc.

 {
   if (fSaveHitsFile) fHitsOutFile.close();
 
   art::ServiceHandle<art::TFileService const> tfs;
 
   TTree* thrTree = tfs->make<TTree>("threshold", "error rate vs threshold");
 
   float thr, shErr, trkErr;
   thrTree->Branch("thr", &thr, "thr/F");
   thrTree->Branch("shErr", &shErr, "shErr/F");
   thrTree->Branch("trkErr", &trkErr, "trkErr/F");
 
   for (size_t i = 0; i < 100; ++i) {
     thr = 0.01 * i;
     shErr = fShBad[i] / float(fShBad[i] + fShOk[i]);
     trkErr = fTrkBad[i] / float(fTrkBad[i] + fTrkOk[i]);
     thrTree->Fill();
   }
 }

bool nnet::PointIdEffTest::isMuonDecaying	(	const simb::MCParticle &	particle,
		const std::unordered_map< int, const simb::MCParticle * > &	particleMap
	)		const

private

Definition at line 513 of file PointIdEffTest_module.cc.

 {
   bool hasElectron = false, hasNuMu = false, hasNuE = false;
 
   int pdg = abs(particle.PdgCode());
   if ((pdg == 13) && (particle.EndProcess() == "FastScintillation")) // potential muon decay at rest
   {
     unsigned int nSec = particle.NumberDaughters();
     for (size_t d = 0; d < nSec; ++d) {
       auto d_search = particleMap.find(particle.Daughter(d));
       if (d_search != particleMap.end()) {
         auto const& daughter = *((*d_search).second);
         int d_pdg = abs(daughter.PdgCode());
         if (d_pdg == 11)
           hasElectron = true;
         else if (d_pdg == 14)
           hasNuMu = true;
         else if (d_pdg == 12)
           hasNuE = true;
       }
     }
   }
 
   return (hasElectron && hasNuMu && hasNuE);
 }

PointIdEffTest& nnet::PointIdEffTest::operator= ( PointIdEffTest const & )

delete

PointIdEffTest& nnet::PointIdEffTest::operator= ( PointIdEffTest && )

delete

int nnet::PointIdEffTest::testCNN	(	detinfo::DetectorClocksData const &	clockData,
		detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< sim::SimChannel > &	channels,
		const std::vector< art::Ptr< recob::Hit >> &	hits,
		const std::array< float, MVA_LENGTH > &	cnn_out,
		const std::vector< anab::FeatureVector< MVA_LENGTH >> &	hit_outs,
		size_t	cidx
	)

private

Definition at line 543 of file PointIdEffTest_module.cc.

 {
   fClSize = hits.size();
 
   std::unordered_map<int, int> mcHitPid;
 
   fPidValue = 0;
   double p_trk_or_sh = cnn_out[fTrkLikeIdx] + cnn_out[fEmLikeIdx];
   if (p_trk_or_sh > 0) { fPidValue = cnn_out[fTrkLikeIdx] / p_trk_or_sh; }
 
   double p_michel = 0;
   if (fMichelLikeIdx >= 0) { fpMichel_cl = cnn_out[fMichelLikeIdx]; }
 
   double totEnSh = 0, totEnTrk = 0, totEnMichel = 0;
   for (auto const& hit : hits) {
     // the channel associated with this hit.
     auto hitChannelNumber = hit->Channel();
 
     double hitEn = 0, hitEnSh = 0, hitEnTrk = 0, hitEnMichel = 0;
 
     auto const& vout = hit_outs[hit.key()];
     fOutTrk = vout[fTrkLikeIdx];
     fOutEM = vout[fEmLikeIdx];
     if (fNoneIdx >= 0) { fOutNone = vout[fNoneIdx]; }
     if (fMichelLikeIdx >= 0) { p_michel = vout[fMichelLikeIdx]; }
 
     for (auto const& channel : channels) {
       if (channel.Channel() != hitChannelNumber) continue;
 
       // for every time slice in this channel:
       auto const& timeSlices = channel.TDCIDEMap();
       for (auto const& timeSlice : timeSlices) {
         int time = timeSlice.first;
         if (std::abs(hit->TimeDistanceAsRMS(time)) < 1.0) {
           // loop over the energy deposits.
           auto const& energyDeposits = timeSlice.second;
 
           for (auto const& energyDeposit : energyDeposits) {
             int trackID = energyDeposit.trackID;
 
             double energy = energyDeposit.numElectrons * fElectronsToGeV * 1000;
             hitEn += energy;
 
             if (trackID < 0) { hitEnSh += energy; } // EM activity
             else if (trackID > 0) {
               auto search = fParticleMap.find(trackID);
               if (search != fParticleMap.end()) {
                 const simb::MCParticle& particle = *((*search).second);
                 int pdg = particle.PdgCode(); // not EM activity so read what PDG it is
 
                 if ((pdg == 11) || (pdg == -11) || (pdg == 22))
                   hitEnSh += energy;
                 else
                   hitEnTrk += energy;
 
                 if (pdg == 11) // electron, check if it is Michel
                 {
                   auto msearch = fParticleMap.find(particle.Mother());
                   if (msearch != fParticleMap.end()) {
                     auto const& mother = *((*msearch).second);
                     if (isMuonDecaying(mother, fParticleMap)) { hitEnMichel += energy; }
                   }
                 }
               }
               else {
                 mf::LogWarning("TrainingDataAlg") << "PARTICLE NOT FOUND";
               }
             }
           }
         }
       }
     }
     totEnSh += hitEnSh;
     totEnTrk += hitEnTrk;
     totEnMichel += hitEnMichel;
 
     double hitAdc =
       hit->SummedADC() * fCalorimetryAlg.LifetimeCorrection(clockData, detProp, hit->PeakTime());
     fTotHit += hitAdc;
 
     int hitPidMc_0p5 = -1;
     if (hitEnSh > hitEnTrk) {
       fHitEM_0p5 += hitAdc;
       hitPidMc_0p5 = nnet::PointIdEffTest::kShower;
     }
     else {
       fHitTrack_0p5 += hitAdc;
       hitPidMc_0p5 = nnet::PointIdEffTest::kTrack;
     }
     mcHitPid[hit.key()] = hitPidMc_0p5;
     auto const& hout = hit_outs[hit.key()];
     fpEM = 0;
     float hit_trk_or_sh = hout[fTrkLikeIdx] + hout[fEmLikeIdx];
     if (hit_trk_or_sh > 0) fpEM = hout[fEmLikeIdx] / hit_trk_or_sh;
     fHitEM_mc = hitEnSh / (hitEnSh + hitEnTrk);
 
     int hitPidMc_0p85 = -1;
     double hitDep = hitEnSh + hitEnTrk;
     if (hitEnSh > 0.85 * hitDep) {
       fHitEM_0p85 += hitAdc;
       fCleanHit += hitAdc;
       hitPidMc_0p85 = nnet::PointIdEffTest::kShower;
     }
     else if (hitEnTrk > 0.85 * hitDep) {
       fHitTrack_0p85 += hitAdc;
       fCleanHit += hitAdc;
       hitPidMc_0p85 = nnet::PointIdEffTest::kTrack;
     }
 
     bool mcMichel = false;
     fpMichel_hit = p_michel;
     fHitMichel_mc = hitEnMichel / hitEn;
     if (fHitMichel_mc > 0.5) {
       fHitMichel_0p5 += hitAdc;
       mcMichel = true;
     }
 
     fHitTree->Fill();
 
     for (size_t i = 0; i < 100; ++i) {
       double thr = 0.01 * i;
 
       if (p_michel > thr) {
         if (mcMichel) { fHitsMichel_OK_0p5[i] += hitAdc; }
         else {
           fHitsMichel_False_0p5[i] += hitAdc;
         }
       }
 
       int recoPid = -1;
       if (fPidValue < thr) {
         recoPid = nnet::PointIdEffTest::kShower;
         fHitRecoEM[i] += hitAdc;
       }
       else
         recoPid = nnet::PointIdEffTest::kTrack;
 
       if ((recoPid == nnet::PointIdEffTest::kShower) &&
           (hitPidMc_0p5 == nnet::PointIdEffTest::kShower)) {
         fHitsEM_OK_0p5[i] += hitAdc;
       }
       else if ((recoPid == nnet::PointIdEffTest::kTrack) &&
                (hitPidMc_0p5 == nnet::PointIdEffTest::kTrack)) {
         fHitsTrack_OK_0p5[i] += hitAdc;
       }
 
       if ((recoPid == nnet::PointIdEffTest::kShower) &&
           (hitPidMc_0p85 == nnet::PointIdEffTest::kShower)) {
         fHitsEM_OK_0p85[i] += hitAdc;
       }
       else if ((recoPid == nnet::PointIdEffTest::kTrack) &&
                (hitPidMc_0p85 == nnet::PointIdEffTest::kTrack)) {
         fHitsTrack_OK_0p85[i] += hitAdc;
       }
     }
   }
 
   // ************ count clusters *************
 
   fMcPid = -1;
   if (totEnSh > 1.5 * totEnTrk) // major energy deposit from EM activity
   {
     fMcPid = nnet::PointIdEffTest::kShower;
   }
   else if (totEnTrk > 1.5 * totEnSh) {
     fMcPid = nnet::PointIdEffTest::kTrack;
   }
 
   for (size_t i = 0; i < 100; ++i) {
     double thr = 0.01 * i;
 
     int recoPid = -1;
     if (fPidValue < thr)
       recoPid = nnet::PointIdEffTest::kShower;
     else
       recoPid = nnet::PointIdEffTest::kTrack;
 
     if ((recoPid == nnet::PointIdEffTest::kShower) && (fMcPid == nnet::PointIdEffTest::kShower)) {
       fShOk[i] += fClSize;
     }
     else if ((recoPid == nnet::PointIdEffTest::kTrack) &&
              (fMcPid == nnet::PointIdEffTest::kTrack)) {
       fTrkOk[i] += fClSize;
     }
     else if ((recoPid == nnet::PointIdEffTest::kShower) &&
              (fMcPid == nnet::PointIdEffTest::kTrack)) {
       fTrkBad[i] += fClSize;
     }
     else if ((recoPid == nnet::PointIdEffTest::kTrack) &&
              (fMcPid == nnet::PointIdEffTest::kShower)) {
       fShBad[i] += fClSize;
     }
     else {
       fNone++;
     }
   }
   fTotal++;
 
   if (fSaveHitsFile) {
     for (auto const& h : hits) {
       auto const& vout = hit_outs[h.key()];
       double hitPidValue = 0;
       double h_trk_or_sh = vout[fTrkLikeIdx] + vout[fEmLikeIdx];
       if (h_trk_or_sh > 0) hitPidValue = vout[fTrkLikeIdx] / h_trk_or_sh;
 
       fHitsOutFile << fRun << " " << fEvent << " " << h->WireID().TPC << " " << h->WireID().Wire
                    << " " << h->PeakTime() << " "
                    << h->SummedADC() *
                         fCalorimetryAlg.LifetimeCorrection(clockData, detProp, h->PeakTime())
                    << " " << mcHitPid[h.key()] << " " << fPidValue << " " << hitPidValue;
 
       if (fMichelLikeIdx >= 0) {
         fHitsOutFile << " " << vout[fMichelLikeIdx]; // is michel?
       }
 
       fHitsOutFile << " " << cidx << std::endl;
     }
   }
 
   fClusterTree->Fill();
   return fMcPid;
 }

Member Data Documentation

calo::CalorimetryAlg nnet::PointIdEffTest::fCalorimetryAlg

private

Definition at line 157 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fCleanHit

private

Definition at line 129 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fClSize

private

Definition at line 138 of file PointIdEffTest_module.cc.

TTree * nnet::PointIdEffTest::fClusterTree

private

Definition at line 149 of file PointIdEffTest_module.cc.

double nnet::PointIdEffTest::fElectronsToGeV

private

Definition at line 145 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fEmLikeIdx

private

Definition at line 147 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fEvent

private

Definition at line 121 of file PointIdEffTest_module.cc.

TTree* nnet::PointIdEffTest::fEventTree

private

Definition at line 149 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitEM_0p5

private

Definition at line 124 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitEM_0p85

private

Definition at line 128 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitEM_mc

private

Definition at line 125 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitMcFractionEM

private

Definition at line 124 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitMichel_0p5

private

Definition at line 124 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitMichel_mc

private

Definition at line 126 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitRecoEM[100]

private

Definition at line 135 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitRecoFractionEM[100]

private

Definition at line 135 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitsEM_OK_0p5[100]

private

Definition at line 132 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitsEM_OK_0p85[100]

private

Definition at line 134 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitsMichel_False_0p5[100]

private

Definition at line 133 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitsMichel_OK_0p5[100]

private

Definition at line 133 of file PointIdEffTest_module.cc.

std::ofstream nnet::PointIdEffTest::fHitsOutFile

private

Definition at line 151 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitsTrack_OK_0p5[100]

private

Definition at line 132 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitsTrack_OK_0p85[100]

private

Definition at line 134 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitTrack_0p5

private

Definition at line 124 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fHitTrack_0p85

private

Definition at line 128 of file PointIdEffTest_module.cc.

TTree * nnet::PointIdEffTest::fHitTree

private

Definition at line 149 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fMcDepEM

private

Definition at line 122 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fMcDepTrack

private

Definition at line 122 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fMcFractionEM

private

Definition at line 122 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fMcPid

private

Definition at line 137 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fMichelLikeIdx

private

Definition at line 147 of file PointIdEffTest_module.cc.

art::InputTag nnet::PointIdEffTest::fNNetModuleLabel

private

Definition at line 160 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fNone

private

Definition at line 143 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fNoneIdx

private

Definition at line 147 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fOutEM

private

Definition at line 127 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fOutNone

private

Definition at line 127 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fOutTrk

private

Definition at line 127 of file PointIdEffTest_module.cc.

std::unordered_map<int, const simb::MCParticle*> nnet::PointIdEffTest::fParticleMap

private

Definition at line 155 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fpEM

private

Definition at line 125 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fPfpDepEM

private

Definition at line 123 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fPfpDepTrack

private

Definition at line 123 of file PointIdEffTest_module.cc.

art::InputTag nnet::PointIdEffTest::fPfpModuleLabel

private

Definition at line 159 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fPidValue

private

Definition at line 139 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fpMichel_cl

private

Definition at line 126 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fpMichel_hit

private

Definition at line 126 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fRun

private

Definition at line 121 of file PointIdEffTest_module.cc.

bool nnet::PointIdEffTest::fSaveHitsFile

private

Definition at line 161 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fShBad[100]

private

Definition at line 142 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fShOk[100]

private

Definition at line 142 of file PointIdEffTest_module.cc.

art::InputTag nnet::PointIdEffTest::fSimulationProducerLabel

private

Definition at line 158 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fTotal

private

Definition at line 143 of file PointIdEffTest_module.cc.

float nnet::PointIdEffTest::fTotHit

private

Definition at line 129 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fTrkBad[100]

private

Definition at line 141 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fTrkLikeIdx

private

Definition at line 147 of file PointIdEffTest_module.cc.

int nnet::PointIdEffTest::fTrkOk[100]

private

Definition at line 141 of file PointIdEffTest_module.cc.

unsigned int nnet::PointIdEffTest::fView

private

Definition at line 153 of file PointIdEffTest_module.cc.

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

larrecodnn/larrecodnn/ImagePatternAlgs/Modules/PointIdEffTest_module.cc

Classes

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation