cluster::HoughLineFinderAna Class Reference

Inheritance diagram for cluster::HoughLineFinderAna:

Public Member Functions
	HoughLineFinderAna (fhicl::ParameterSet const &pset)

	~HoughLineFinderAna ()

Private Member Functions
void	analyze (const art::Event &)

void	beginJob ()

Private Attributes
std::string	fHoughModuleLabel

std::string	fDigitModuleLabel

std::string	fHitsModuleLabel

std::string	fDBScanModuleLabel

TTree *	ftree

int	fm_run

unsigned long int	fm_run_timestamp

int	fm_event

int	fm_plane

int	fm_dbsize

int	fm_clusterid

int	fm_wirespan

int	fm_sizeClusterZ

int	fm_sizeHitZ

float	fm_clusterslope

float	fm_clusterintercept

int *	fm_wireZ

int *	fm_hitidZ

float *	fm_mipZ

float *	fm_drifttimeZ

float *	fm_widthZ

float *	fm_upadcZ

Detailed Description

Definition at line 39 of file HoughLineFinderAna_module.cc.

Constructor & Destructor Documentation

cluster::HoughLineFinderAna::HoughLineFinderAna ( fhicl::ParameterSet const & pset )

explicit

Definition at line 84 of file HoughLineFinderAna_module.cc.

     : EDAnalyzer(pset)
     , fHoughModuleLabel (pset.get< std::string >("HoughModuleLabel"))
     , fDigitModuleLabel (pset.get< std::string >("DigitModuleLabel"))
     , fHitsModuleLabel  (pset.get< std::string >("HitsModuleLabel"))
     , fDBScanModuleLabel(pset.get< std::string >("DBScanModuleLabel"))
     , fm_run(0)
     , fm_event(0)
     , fm_plane(0)
     , fm_dbsize(0)
     , fm_clusterid(0)
     , fm_wirespan(0)
     , fm_sizeClusterZ(10000)
     , fm_sizeHitZ(10000)
     , fm_clusterslope(0)
     , fm_clusterintercept(0)
   {
   }

cluster::HoughLineFinderAna::~HoughLineFinderAna ( )

Definition at line 104 of file HoughLineFinderAna_module.cc.

   {
     delete fm_hitidZ;
     delete fm_mipZ;
     delete fm_drifttimeZ;
     delete fm_widthZ;
     delete fm_upadcZ;
     delete fm_wireZ;
   }

Member Function Documentation

void cluster::HoughLineFinderAna::analyze ( const art::Event & evt )

private

Definition at line 147 of file HoughLineFinderAna_module.cc.

   {
 
     art::Handle< std::vector<recob::Cluster> > hlfListHandle;
     evt.getByLabel(fHoughModuleLabel,hlfListHandle);
     art::Handle< std::vector<recob::Hit> > hitListHandle;
     evt.getByLabel(fHitsModuleLabel,hitListHandle);
     art::Handle< std::vector<recob::Cluster> > dbscanListHandle;
     evt.getByLabel(fDBScanModuleLabel,dbscanListHandle);
 
     art::FindManyP<recob::Hit> fmh(dbscanListHandle, evt, fDBScanModuleLabel);
     art::FindManyP<recob::Hit> fmhhl(hlfListHandle, evt, fHoughModuleLabel);
 
     art::PtrVector<recob::Cluster> clusters;
     art::PtrVector<recob::Cluster> dbclusters;
     //   art::PtrVector<recob::Hit> hits;// unused, as yet. EC, 5-Oct-2010.
 
     for (size_t ii = 0; ii <  hlfListHandle->size(); ++ii){
       art::Ptr<recob::Cluster> cluster(hlfListHandle,ii);
       clusters.push_back(cluster);
     }
 
     for (size_t ii = 0; ii <  dbscanListHandle->size(); ++ii){
       art::Ptr<recob::Cluster> dbcluster(dbscanListHandle,ii);
       dbclusters.push_back(dbcluster);
     }
 
     MF_LOG_VERBATIM("HoughLineFinderAna") << "run    : " << evt.id().run();
     //std::cout << "subrun : " << evt.subRun() << std::endl;
     MF_LOG_VERBATIM("HoughLineFinderAna") << "event  : " << evt.id().event();
     fm_run=evt.id().run();
     fm_event=evt.id().event();
     fm_run_timestamp=evt.time().value(); // won't cast, EC, 7-Oct-2010.
     unsigned int firstwire=0;
     unsigned int lastwire=0;
     fm_sizeClusterZ=0;
     fm_sizeHitZ=0;
     fm_dbsize=0;
     art::ServiceHandle<geo::Geometry const> geo;
 
     for(auto view : geo->Views()){
 
       fm_dbsize       = 0;
       fm_sizeClusterZ = clusters.size();
 
       for(size_t j = 0; j < dbclusters.size(); ++j) {
         if(dbclusters[j]->View() == view){
           std::vector< art::Ptr<recob::Hit> > _dbhits = fmh.at(j);
           fm_dbsize += _dbhits.size();
           if(_dbhits.size() > 0) fm_plane   = _dbhits.at(0)->WireID().Plane;
         }
       }
 
       for(size_t j = 0; j < clusters.size(); ++j) {
         if(clusters[j]->View() == view){
           fm_clusterid=clusters[j]->ID();
           std::vector< art::Ptr<recob::Hit> > _hits = fmhhl.at(j);
           fm_clusterslope=(double) std::tan(clusters[j]->StartAngle());
           fm_clusterintercept=(double)clusters[j]->StartTick();
           if(_hits.size()!=0){
             fm_plane   = _hits.at(0)->WireID().Plane;
             firstwire = _hits[0]->WireID().Wire;
             lastwire  = _hits[_hits.size()-1]->WireID().Wire;
             fm_wirespan = lastwire-firstwire;
             fm_sizeHitZ = _hits.size();
 
             for(unsigned int i = 0; i < _hits.size(); ++i){
 
               fm_hitidZ[i]     = i;
               fm_wireZ[i]      = _hits[i]->WireID().Wire;
               fm_mipZ[i]       = (double)_hits[i]->Integral();
               fm_drifttimeZ[i] = (double)_hits[i]->PeakTime();
               fm_widthZ[i]     = (double) (2. * _hits[i]->RMS());
               fm_upadcZ[i]     = (double)_hits[i]->Integral();
             }
 
             ftree->Fill();
           }
         }//end if in the correct view
       }// end loop over clusters
     }// end loop over views
 
   }

void cluster::HoughLineFinderAna::beginJob ( )

private

Definition at line 115 of file HoughLineFinderAna_module.cc.

   {
 
 
     // get access to the TFile service
     art::ServiceHandle<art::TFileService const> tfs;
     ftree= tfs->make<TTree>("HoughTree","HoughTree");
     fm_hitidZ = new int[fm_sizeHitZ];
     fm_mipZ = new float[fm_sizeHitZ];
     fm_drifttimeZ = new float[fm_sizeHitZ];
     fm_widthZ = new float[fm_sizeHitZ];
     fm_upadcZ = new float[fm_sizeHitZ];
     fm_wireZ = new int[fm_sizeHitZ];
     ftree->Branch("run", &fm_run, "run/I");
     ftree->Branch("run_timestamp", &fm_run_timestamp, "run_timestamp/l"); //l is for ULong64_t
     ftree->Branch("event", &fm_event, "event/I");
     ftree->Branch("plane",&fm_plane,"plane/I");
     ftree->Branch("dbsize",&fm_dbsize,"dbsize/I");
     ftree->Branch("clusterid",&fm_clusterid,"clusterid/I");
     ftree->Branch("clusterslope",&fm_clusterslope,"clusterslope/F");
     ftree->Branch("clusterintercept",&fm_clusterintercept,"clusterintecept/F");
     ftree->Branch("wirespan",&fm_wirespan,"wirespan/I");
     ftree->Branch("numberHits",&fm_sizeHitZ,"numberHits/I");
     ftree->Branch("numberClusters",&fm_sizeClusterZ,"numberClusters/I");
     ftree->Branch("hitidZ",fm_hitidZ,"hitidZ[numberHits]/I");
     ftree->Branch("wireZ",fm_wireZ,"wireZ[numberHits]/I");
     ftree->Branch("mipZ",fm_mipZ,"mipZ[numberHits]/F");
     ftree->Branch("drifttimeZ",fm_drifttimeZ,"drifttitmeZ[numberHits]/F");
     ftree->Branch("widthZ",fm_widthZ,"widthZ[numberHits]/F");
   }