disambigcheck::disambigcheck Class Reference

Inheritance diagram for disambigcheck::disambigcheck:

Public Member Functions
	disambigcheck (fhicl::ParameterSet const &pset)
void	analyze (const art::Event &evt)
void	beginJob ()
void	beginRun (const art::Run &run)
void	endJob ()
void	reconfigure (fhicl::ParameterSet const &pset)
virtual	~disambigcheck ()
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
art::ServiceHandle< geo::Geometry >	fGeom
art::ServiceHandle< cheat::BackTrackerService >	bt_serv
art::ServiceHandle< cheat::ParticleInventoryService >	pi_serv
std::string	fChanHitDisambig
std::string	fChanHitCheater
TH1D *	fCorrect
TH1D *	fMissed
TH1D *	fIncorrect
TH2D *	fXCorrect
TH2D *	fXMissed
TH2D *	fXIncorrect
TH2D *	fYCorrect
TH2D *	fYMissed
TH2D *	fZIncorrect
TH2D *	fZCorrect
TH2D *	fZMissed
TH2D *	fYIncorrect
TH2D *	fThetaCorrect
TH2D *	fThetaMissed
TH2D *	fThetaIncorrect
TH2D *	fPhiCorrect
TH2D *	fPhiMissed
TH2D *	fPhiIncorrect

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 60 of file disambigcheck_module.cc.

Constructor & Destructor Documentation

disambigcheck::disambigcheck::disambigcheck ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 99 of file disambigcheck_module.cc.

    : EDAnalyzer (pset)  
  {
    this->reconfigure(pset);
    //  produces< std::vector<recob::Hit> >();
  }

disambigcheck::disambigcheck::~disambigcheck ( )

virtual

Definition at line 106 of file disambigcheck_module.cc.

{}

Member Function Documentation

void disambigcheck::disambigcheck::analyze

(

const art::Event &

evt

)

Definition at line 159 of file disambigcheck_module.cc.

  {
    
    std::unique_ptr<std::vector<recob::Hit> > hcol(new std::vector<recob::Hit>);
    std::vector< art::Ptr<recob::Hit> >  ChHitsDisambig;
    auto ChannelHitsDisambig = evt.getHandle< std::vector<recob::Hit> >(fChanHitDisambig);
    if (ChannelHitsDisambig)
      art::fill_ptr_vector(ChHitsDisambig, ChannelHitsDisambig);

    // Make unambiguous collection hits
    std::vector< art::Ptr<recob::Hit> >  ChHitsCheater;
    auto ChannelHitsCheater = evt.getHandle< std::vector<recob::Hit> >(fChanHitCheater);
    if (ChannelHitsCheater)
      art::fill_ptr_vector(ChHitsCheater, ChannelHitsCheater);
    
    int correcthits = 0;
    int incorrecthits = 0; 
    int missinghits = 0;
    // Map of MCParticle Track IDs, Vector 1 ( matched, incorrect, missed ),  Pair 1 ( Pair 2 (Theta, Phi) , Vector 2 direction cosines ( x, y, z ) ) 
    std::map< int, std::pair< TVector3, std::pair < std::pair< double,double >, TVector3 > > > ParticleMap;     
    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
    if (ChannelHitsDisambig.isValid()) {
      for( size_t h = 0; h < ChHitsCheater.size(); h++ ){
        if(ChHitsCheater[h]->View() == geo::kZ ) continue;
        uint32_t cheatchannel = ChHitsCheater[h]->Channel();
        double cheatpeaktime = ChHitsCheater[h]->PeakTime();
        uint32_t cheatwire = ChHitsCheater[h]->WireID().Wire;

        std::vector<sim::IDE> ides;
        std::vector<sim::TrackIDE> TrackIDs = bt_serv->HitToTrackIDEs(clockData, ChHitsCheater[h]);
        double MaxE = 0;
        int TrackID;
        for(size_t e = 0; e < TrackIDs.size(); ++e){
          if ( TrackIDs[e].energy > MaxE ) {
            TrackID = TrackIDs[e].trackID;
            MaxE    = TrackIDs[e].energy;
          }
        }
        // Now have trackID, so get PdG code and T0 etc.
        const simb::MCParticle *particle = pi_serv->TrackIdToParticle_P(TrackID);
        if ( ParticleMap.count(TrackID) == 0) {
          //std::cout << "Looking at key " << TrackID << " for the first time so setting the angles." << std::endl;
          TVector3 MC_NormMom;
          MC_NormMom[0] = ( particle->Px() / particle->P() );
          MC_NormMom[1] = ( particle->Py() / particle->P() );
          MC_NormMom[2] = ( particle->Pz() / particle->P() );
          double MC_Theta = MC_NormMom.Theta();
          double MC_Phi   = MC_NormMom.Phi();
          ParticleMap[TrackID].second.second[0] = MC_NormMom[0];
          ParticleMap[TrackID].second.second[1] = MC_NormMom[1];
          ParticleMap[TrackID].second.second[2] = MC_NormMom[2];
          ParticleMap[TrackID].second.first.first  = MC_Theta;
          ParticleMap[TrackID].second.first.second = MC_Phi;
        } // If first instance of this trackID fill quantities
        bool found = false;
        for( size_t w = 0; w < ChHitsDisambig.size(); w++ ) {
          uint32_t disambigchannel = ChHitsDisambig[w]->Channel();
          double disambigpeaktime = ChHitsDisambig[w]->PeakTime();
          uint32_t disambigwire = ChHitsDisambig[w]->WireID().Wire;
        
          if((cheatchannel == disambigchannel)&&(TMath::Abs(cheatpeaktime - disambigpeaktime)< 0.1)) {
          
            found=true;
            if(cheatwire==disambigwire){
              correcthits++;   
              ++ParticleMap[TrackID].first[0];
            }
            else{
              incorrecthits++;
              ++ParticleMap[TrackID].first[1];
            }   
          }     
        } // Loop through disambiguated hits
        if(!found){
          missinghits++;
          ++ParticleMap[TrackID].first[2];
        } 
        //std::cout << "Particle Map Angles...." << ParticleMap[TrackID].first[0] << " " << ParticleMap[TrackID].first[1] << " " << ParticleMap[TrackID].first[2] <<"\n"
        //      << "And it now has " << ParticleMap[TrackID].second[0] << " correct, " << ParticleMap[TrackID].second[1] << " incorrect and " << ParticleMap[TrackID].second[2] << " missed hits." << std::endl;
      } // size of cheated hits

      int totalhits = (correcthits + incorrecthits + missinghits); 
           
      //double disambighitsFraction = ddisambighits / dtotalhits;
      if ( totalhits != 0 ) {
        fCorrect  ->Fill( (double)correcthits   / (double)totalhits );
        fMissed   ->Fill( (double)missinghits / (double)totalhits ); 
        fIncorrect->Fill( (double)incorrecthits   / (double)totalhits ); 
        if ( (double)correcthits / (double)totalhits == 0 ) std::cout << "WHY IS THIS EVENT NOT MATCHING ANY HITS!!!????" << std::endl;
      } // totalhits != 0
      std::cout << "\nTotal hits " << totalhits << ", correct hits " << correcthits << ", incorrect hits " << incorrecthits << ", missing hits " << missinghits << std::endl;

      for (std::map< int, std::pair< TVector3, std::pair < std::pair< double,double >, TVector3 > > >::iterator ii = ParticleMap.begin(); ii!=ParticleMap.end(); ++ii){
        int MapTotalHits = ii->second.first[0] + ii->second.first[1] + ii->second.first[2];
        std::cout << "Total hits for this TrackID " << ii->first << " is; " << MapTotalHits << ", Correct " << ii->second.first[0] << ", Incorrect " << ii->second.first[1] << ", Missed " << ii->second.first[2] << std::endl;
        if ( MapTotalHits != 0 ) {
          fXCorrect   -> Fill ( ii->second.second.second[0], (double)ii->second.first[0] / (double)MapTotalHits );
          fXMissed    -> Fill ( ii->second.second.second[0], (double)ii->second.first[1] / (double)MapTotalHits );
          fXIncorrect -> Fill ( ii->second.second.second[0], (double)ii->second.first[2] / (double)MapTotalHits );
          fYCorrect   -> Fill ( ii->second.second.second[1], (double)ii->second.first[0] / (double)MapTotalHits );
          fYMissed    -> Fill ( ii->second.second.second[1], (double)ii->second.first[1] / (double)MapTotalHits );
          fYIncorrect -> Fill ( ii->second.second.second[1], (double)ii->second.first[2] / (double)MapTotalHits );
          fZCorrect   -> Fill ( ii->second.second.second[2], (double)ii->second.first[0] / (double)MapTotalHits );
          fZMissed    -> Fill ( ii->second.second.second[2], (double)ii->second.first[1] / (double)MapTotalHits );
          fZIncorrect -> Fill ( ii->second.second.second[2], (double)ii->second.first[2] / (double)MapTotalHits );

          fThetaCorrect   -> Fill ( ii->second.second.first.first , (double)ii->second.first[0] / (double)MapTotalHits );
          fThetaMissed    -> Fill ( ii->second.second.first.first , (double)ii->second.first[1] / (double)MapTotalHits );
          fThetaIncorrect -> Fill ( ii->second.second.first.first , (double)ii->second.first[2] / (double)MapTotalHits );
          fPhiCorrect     -> Fill ( ii->second.second.first.second, (double)ii->second.first[0] / (double)MapTotalHits );
          fPhiMissed      -> Fill ( ii->second.second.first.second, (double)ii->second.first[1] / (double)MapTotalHits );
          fPhiIncorrect   -> Fill ( ii->second.second.first.second, (double)ii->second.first[2] / (double)MapTotalHits );       

          //Some cout statements to have a look at......
          /*
            if ( (double)ii->second.first[0] / (double)MapTotalHits < 0.8 ) {
            std::cout << "\n\nThis track has efficiency " << (double)ii->second.first[0] / (double)MapTotalHits*100 
            << "\nTheta is " << ii->second.second.first.first << ", and phi is " << ii->second.second.first.second 
            << std::endl;
            if ( ii->second.second.first.first < 1.7 && ii->second.second.first.first > 1.45 ) {        
            std::cout << "THIS TRACK HAS LOWER EFFICIENCY AT THETA IS ROUGHLY PI/2." << std::endl;
            }
            if ( ii->second.second.first.second < -1.45 && ii->second.second.first.second > -1.7 ) {
            std::cout << "THIS TRACK HAS LOWER EFFICIENCY AT PHI IS ROUGHLY -PI/2." << std::endl;
            }
            std::cout <<"\n\n" << std::endl;
            }//*/
        } // MapTotalHits != 0
      } // Loop through Map
    } // If hits valid   
    ParticleMap.clear();
  }

void disambigcheck::disambigcheck::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 119 of file disambigcheck_module.cc.

  {
    art::ServiceHandle<art::TFileService> tfs;
    
    fCorrect   = tfs->make<TH1D>("correct"  ,"Correct disambiguated"  , 100, 0, 1.01);
    fMissed    = tfs->make<TH1D>("missed"   ,"Missed hits"            , 100, 0, 1.01);
    fIncorrect = tfs->make<TH1D>("incorrect","Incorrect disambiguated", 100, 0, 1.01);

    fXCorrect   = tfs->make<TH2D>("Xcorrect"  ,"Correct disambiguated"  , 40, -1.01, 1.01, 40, 0, 1.01);
    fXMissed    = tfs->make<TH2D>("Xmissed"   ,"Missed hits"            , 40, -1.01, 1.01, 40, 0, 1.01);
    fXIncorrect = tfs->make<TH2D>("Xincorrect","Incorrect disambiguated", 40, -1.01, 1.01, 40, 0, 1.01);
    fYCorrect   = tfs->make<TH2D>("Ycorrect"  ,"Correct disambiguated"  , 40, -1.01, 1.01, 40, 0, 1.01);
    fYMissed    = tfs->make<TH2D>("Ymissed"   ,"Missed hits"            , 40, -1.01, 1.01, 40, 0, 1.01);
    fYIncorrect = tfs->make<TH2D>("Yincorrect","Incorrect disambiguated", 40, -1.01, 1.01, 40, 0, 1.01);
    fZCorrect   = tfs->make<TH2D>("Zcorrect"  ,"Correct disambiguated"  , 40, -1.01, 1.01, 40, 0, 1.01);
    fZMissed    = tfs->make<TH2D>("Zmissed"   ,"Missed hits"            , 40, -1.01, 1.01, 40, 0, 1.01);
    fZIncorrect = tfs->make<TH2D>("Zincorrect","Incorrect disambiguated", 40, -1.01, 1.01, 40, 0, 1.01); 

    fThetaCorrect   = tfs->make<TH2D>("Thetacorrect"  ,"Correct disambiguated"  , 20, 0    , 3.15, 40, 0, 1.01);
    fThetaMissed    = tfs->make<TH2D>("Thetamissed"   ,"Missed hits"            , 20, 0    , 3.15, 40, 0, 1.01);
    fThetaIncorrect = tfs->make<TH2D>("Thetaincorrect","Incorrect disambiguated", 20, 0    , 3.15, 40, 0, 1.01);
    fPhiCorrect     = tfs->make<TH2D>("Phicorrect"    ,"Correct disambiguated"  , 40, -3.15, 3.15, 40, 0, 1.01);
    fPhiMissed      = tfs->make<TH2D>("Phimissed"     ,"Missed hits"            , 40, -3.15, 3.15, 40, 0, 1.01);
    fPhiIncorrect   = tfs->make<TH2D>("Phiincorrect"  ,"Incorrect disambiguated", 40, -3.15, 3.15, 40, 0, 1.01);
  
    //cXCorrect = tfs->make<TCanvas>("XCorrect", "", 600, 500);
}

void disambigcheck::disambigcheck::beginRun

(

const art::Run &

run

)

Definition at line 149 of file disambigcheck_module.cc.

{}

void disambigcheck::disambigcheck::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 150 of file disambigcheck_module.cc.

  {
    //cXCorrect -> cd(0);
    //fXCorrect -> Draw("colz");
    //cXCorrect -> Write();
  }

void disambigcheck::disambigcheck::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 109 of file disambigcheck_module.cc.

  {
    
    fChanHitDisambig =  pset.get< std::string >("ChanHitDisambig");
    fChanHitCheater  =  pset.get< std::string >("ChanHitCheater");
    
  }  

Member Data Documentation

art::ServiceHandle<cheat::BackTrackerService> disambigcheck::disambigcheck::bt_serv

private

Definition at line 78 of file disambigcheck_module.cc.

std::string disambigcheck::disambigcheck::fChanHitCheater

private

Definition at line 82 of file disambigcheck_module.cc.

std::string disambigcheck::disambigcheck::fChanHitDisambig

private

Definition at line 81 of file disambigcheck_module.cc.

TH1D* disambigcheck::disambigcheck::fCorrect

private

Definition at line 84 of file disambigcheck_module.cc.

art::ServiceHandle<geo::Geometry> disambigcheck::disambigcheck::fGeom

private

Definition at line 77 of file disambigcheck_module.cc.

TH1D* disambigcheck::disambigcheck::fIncorrect

private

Definition at line 84 of file disambigcheck_module.cc.

TH1D* disambigcheck::disambigcheck::fMissed

private

Definition at line 84 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fPhiCorrect

private

Definition at line 89 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fPhiIncorrect

private

Definition at line 89 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fPhiMissed

private

Definition at line 89 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fThetaCorrect

private

Definition at line 88 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fThetaIncorrect

private

Definition at line 88 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fThetaMissed

private

Definition at line 88 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fXCorrect

private

Definition at line 85 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fXIncorrect

private

Definition at line 85 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fXMissed

private

Definition at line 85 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fYCorrect

private

Definition at line 86 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fYIncorrect

private

Definition at line 87 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fYMissed

private

Definition at line 86 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fZCorrect

private

Definition at line 87 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fZIncorrect

private

Definition at line 86 of file disambigcheck_module.cc.

TH2D* disambigcheck::disambigcheck::fZMissed

private

Definition at line 87 of file disambigcheck_module.cc.

art::ServiceHandle<cheat::ParticleInventoryService> disambigcheck::disambigcheck::pi_serv

private

Definition at line 79 of file disambigcheck_module.cc.

---

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

• dunereco/dunereco/HitFinderDUNE/disambigcheck_module.cc