Inheritance diagram for lightana::FlashT0SelectedChannels:

Classes
struct	Config

Public Member Functions
	FlashT0SelectedChannels (art::ToolConfigTable< Config > const &config)

double	GetFlashT0 (double flash_peaktime, LiteOpHitArray_t ophit_list) override

Private Attributes
double	fPDFraction

double	fPreWindow

double	fPostWindow

double	fMinHitPE

Additional Inherited Members
Private Member Functions inherited from lightana::FlashT0Base
virtual	~FlashT0Base () noexcept=default

Detailed Description

Definition at line 24 of file FlashT0SelectedChannels_tool.cc.

Constructor & Destructor Documentation

lightana::FlashT0SelectedChannels::FlashT0SelectedChannels ( art::ToolConfigTable< Config > const & config )

explicit

Definition at line 69 of file FlashT0SelectedChannels_tool.cc.

     : fPDFraction { config().PDFraction() },
     fPreWindow  { config().PreWindow()  },
     fPostWindow { config().PostWindow() },
     fMinHitPE   { config().MinHitPE()   }
   {
   }

Member Function Documentation

double lightana::FlashT0SelectedChannels::GetFlashT0	(	double	flash_peaktime,
		LiteOpHitArray_t	ophit_list
	)

overridevirtual

Implements lightana::FlashT0Base.

Definition at line 77 of file FlashT0SelectedChannels_tool.cc.

                                                                                           {
 
     std::vector< std::pair<double, double> > selected_hits;
     double pe_sum = 0;
 
     // fill vector with selected hits in the specified window
     for(auto const& hit : ophit_list) {
       if( hit.peak_time<flash_time+fPostWindow && hit.peak_time>flash_time-fPreWindow && hit.pe>fMinHitPE){
         selected_hits.push_back( std::make_pair(hit.pe, hit.peak_time));
         pe_sum += hit.pe ;
       }
     }
 
     if(pe_sum>0){
       // sort vector by number of #PE (ascending order)
       std::sort( selected_hits.begin(), selected_hits.end(), std::greater< std::pair<double, double> >() );
 
       double flasht0_mean=0, pe_count=0;
       int nophits=0;
 
       // loop over selected ophits
       for (size_t ix=0; ix<selected_hits.size(); ix++) {
         pe_count += selected_hits[ix].first;
         flasht0_mean += selected_hits[ix].second;
         nophits++;
         if( pe_count/pe_sum>fPDFraction ) break;
       }
 
       return flasht0_mean/nophits;
     }
     else
       return flash_time;
   }