sn::SupernovaAna Class Reference

Inheritance diagram for sn::SupernovaAna:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer

Public Member Functions
	SupernovaAna (fhicl::ParameterSet const &p)
	SupernovaAna (SupernovaAna const &)=delete
	SupernovaAna (SupernovaAna &&)=delete
SupernovaAna &	operator= (SupernovaAna const &)=delete
SupernovaAna &	operator= (SupernovaAna &&)=delete
void	analyze (art::Event const &e) override
void	beginJob ()
void	endJob ()
void	FillPlaneSeperations (art::Ptr< rb::Cluster > &aSlice, int &min_cell_sep, double &min_plane_sep, double &min_sep)
std::string	workerType () const
bool	modifiesEvent () const
void	registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)
std::string const &	processName () const
bool	wantAllEvents () const
bool	wantEvent (Event const &e)
fhicl::ParameterSetID	selectorConfig () const
art::Handle< art::TriggerResults >	getTriggerResults (Event const &e) const
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)
base_engine_t &	createEngine (seed_t seed)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)
base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const
detail::CachedProducts &	cachedProducts ()
void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
void	prepareForJob (fhicl::ParameterSet const &pset)
void	showMissingConsumes () const

Private Attributes
std::string	fCellHitsLabel
	Label of module making the cell hits. More...
std::string	fSlicerLabel
	Instance label making the slices. More...
bool	fVerbose
	Verbose output. More...
int	nEvent = 0
TH1F *	hNSlices
TH1F *	hSlicesNHits
TH1F *	hSlicesSumADC
TH1F *	hSlicesMeanADC
TH1F *	hSlicesMinInPlaneSeperation
TH1F *	hSlicesMinCrossPlaneSeperation
TH1F *	hSlicesMinSeperation
TH1F *	hNNoiseSlices
TH1F *	hNoiseSlicesNHits
TH1F *	hNoiseSlicesSumADC
TH1F *	hNoiseSlicesMeanADC
TH1F *	hNoiseSlicesMinInPlaneSeperation
TH1F *	hNoiseSlicesMinCrossPlaneSeperation
TH1F *	hNoiseSlicesMinSeperation
TH1F *	h_t_SumADC
TH1F *	h_t_Nhits
TH1F *	h_t_MeanADC
TH1F *	h_t_NSlices
TH1F *	h_t_SliceNHits
TH1F *	h_t_SliceMeanHits
TH1F *	h_t_NNoiseSlices
TH1F *	h_t_NoiseNHits
TH1F *	h_t_NoiseMeanHits
TH2F *	hADC_vs_mb
TH2F *	hSlice_hits_vs_mb
TH2F *	hMeanADC_vs_Slice_hits
TH2F *	hNoiseSlice_hits_vs_mb
TH2F *	hNoiseMeanADC_vs_Slice_hits

Detailed Description

Definition at line 38 of file SupernovaAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 39 of file EDAnalyzer.h.

using art::EDAnalyzer::WorkerType = WorkerT<EDAnalyzer>

inherited

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

sn::SupernovaAna::SupernovaAna ( fhicl::ParameterSet const &  p )

explicit

Definition at line 99 of file SupernovaAna_module.cc.

References om::cout, allTimeWatchdog::endl, fCellHitsLabel, fSlicerLabel, fVerbose, h_t_MeanADC, h_t_Nhits, h_t_NNoiseSlices, h_t_NoiseMeanHits, h_t_NoiseNHits, h_t_NSlices, h_t_SliceMeanHits, h_t_SliceNHits, h_t_SumADC, hADC_vs_mb, hMeanADC_vs_Slice_hits, hNNoiseSlices, hNoiseMeanADC_vs_Slice_hits, hNoiseSlice_hits_vs_mb, hNoiseSlicesMeanADC, hNoiseSlicesMinCrossPlaneSeperation, hNoiseSlicesMinInPlaneSeperation, hNoiseSlicesMinSeperation, hNoiseSlicesNHits, hNoiseSlicesSumADC, hNSlices, hSlice_hits_vs_mb, hSlicesMeanADC, hSlicesMinCrossPlaneSeperation, hSlicesMinInPlaneSeperation, hSlicesMinSeperation, hSlicesNHits, hSlicesSumADC, and art::TFileDirectory::make().

                                                        : EDAnalyzer(p),
  fCellHitsLabel (p.get<std::string >("CellHitsLabel")),
  fSlicerLabel  (p.get<std::string >("SlicerLabel")),
  fVerbose      (p.get<bool >("verbose"))
{
    std::cout << "--- sn::SupernovaAna instantiate" << std::endl;
    std::cout << "\t CellHitModuleLabel:     " << fCellHitsLabel    << std::endl;
    std::cout << "\t SlicerLabel:            " << fSlicerLabel            << std::endl;
    std::cout << "\t verbose:                " << fVerbose                << std::endl;
    
    art::ServiceHandle<art::TFileService> tfs;
    hNSlices= tfs->make<TH1F>("NSlices",  ";NSlices;Entries",200,0,7000);
    hSlicesNHits= tfs->make<TH1F>("SlicesNHits",  ";SlicesNHits;Entries",500,0,500);
    hSlicesSumADC= tfs->make<TH1F>("SlicesSumADC",  ";SlicesSumADC;Entries",400,0,20000);
    hSlicesMeanADC= tfs->make<TH1F>("SlicesMeanADC",  ";SlicesMeanADC;Entries",410,0.5,4100.5);
    hSlicesMinInPlaneSeperation= tfs->make<TH1F>("SlicesMinInPlaneSeperation",  ";hSlicesMinInPlaneSeperation;;Entries",22,-0.5,20.5);
    hSlicesMinCrossPlaneSeperation= tfs->make<TH1F>("SlicesMinCrossPlaneSeperation",  ";hSlicesMinCrossPlaneSeperation;Entries",200,-0.5,20.5);
    hSlicesMinSeperation= tfs->make<TH1F>("SlicesMinSeperation",  ";hSlicesMinSeperation;Entries",200,-0.5,20.5);
    
    hNNoiseSlices= tfs->make<TH1F>("NNoiseSlices",  ";NSlices;Entries",200,0,7000);
    hNoiseSlicesNHits= tfs->make<TH1F>("NoiseSlicesNHits",  ";SlicesNHits;Entries",200,0,200);
    hNoiseSlicesSumADC= tfs->make<TH1F>("NoiseSlicesSumADC",  ";SlicesSumADC;Entries",200,0,10000);
    hNoiseSlicesMeanADC= tfs->make<TH1F>("NoiseSlicesMeanADC",  ";SlicesMeanADC;Entries",410,0.5,4100.5);
    hNoiseSlicesMinInPlaneSeperation= tfs->make<TH1F>("NoiseSlicesMinInPlaneSeperation",  ";hSlicesMinInPlaneSeperation;;Entries",11,-0.5,10.5);
    hNoiseSlicesMinCrossPlaneSeperation= tfs->make<TH1F>("NoiseSlicesMinCrossPlaneSeperation",  ";hSlicesMinCrossPlaneSeperation;Entries",100,-0.5,10.5);
    hNoiseSlicesMinSeperation= tfs->make<TH1F>("NoiseSlicesMinSeperation",  ";hSlicesMinSeperation;Entries",100,-0.5,10.5);
    
    // time dependent
    h_t_SumADC = tfs->make<TH1F>("t_SumADC",   "SumADC;Event;SumADC",  400,-.5,399.5);
    h_t_Nhits  = tfs->make<TH1F>("t_Nhits",    "Nhits;Event;Nhits",    400,-.5,399.5);
    h_t_MeanADC= tfs->make<TH1F>("t_MeanADC",  "MeanADC;Event;MeanADC",400,-.5,399.5);
    h_t_NSlices= tfs->make<TH1F>("t_NSlices",  "NSlices;Event;NSlices",400,-.5,399.5);
    h_t_SliceNHits= tfs->make<TH1F>("t_SliceNHits",  "SliceNHits;Event;NSliceHits",400,-.5,399.5);
    h_t_SliceMeanHits= tfs->make<TH1F>("t_SliceMeanHits",  "SliceMeanHits;Event;SliceMeanHits",400,-.5,399.5);
    h_t_NNoiseSlices= tfs->make<TH1F>("t_NNoiseSlices",  "NNoiseSlices;Event;NNoiseSlices",400,-.5,399.5);
    h_t_NoiseNHits= tfs->make<TH1F>("t_NoiseNHits",  "NoiseNHits;Event;NoiseNHits",400,-.5,399.5);
    h_t_NoiseMeanHits= tfs->make<TH1F>("t_NoiseMeanHits",  "NoiseMeanHits;Event;NoiseMeanHits",400,-.5,399.5);
    
    hADC_vs_mb = tfs->make<TH2F>("ADC_vs_mb", "ADC_vs_mb;mb;ADC;entries",400,-.5,399.5,410,0.5,4100.5);
    
    hSlice_hits_vs_mb = tfs->make<TH2F>("Slice_hits_vs_mb", "Slice_hits_vs_mb;mb;Slice Hits;entries",400,-.5,399.5,101,-0.5,100.5);
    hMeanADC_vs_Slice_hits = tfs->make<TH2F>("MeanADC_vs_Slice_hits", "MeanADC_vs_Slice_hits;MeanADC;Slice Hits;entries",410,0.5,4100.5,101,-0.5,100.5);
    
    hNoiseSlice_hits_vs_mb = tfs->make<TH2F>("NoiseSlice_hits_vs_mb", "NoiseSlice_hits_vs_mb;mb;Noise Slice Hits;entries",200,-.5,199.5,101,-0.5,100.5);
    hNoiseMeanADC_vs_Slice_hits = tfs->make<TH2F>("NoiseMeanADC_vs_Slice_hits", "NoiseMeanADC_vs_Slice_hits;MeanADC;Slice Hits;entries",410,0.5,4100.5,101,-0.5,100.5);
}

sn::SupernovaAna::SupernovaAna ( SupernovaAna const &  )

delete

sn::SupernovaAna::SupernovaAna ( SupernovaAna &&  )

delete

Member Function Documentation

void sn::SupernovaAna::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 160 of file SupernovaAna_module.cc.

References rawdata::RawDigit::ADC(), getBrightness::cell, rawdata::RawDigit::Channel(), visualisationForPaperMasterPlot::cmap, om::cout, allTimeWatchdog::endl, art::EventID::event(), fCellHitsLabel, FillPlaneSeperations(), fSlicerLabel, fVerbose, art::DataViewImpl::getByLabel(), cmap::dataprov::CMap::GetCell(), cmap::dataprov::CMap::GetPlane(), daqchannelmap::DAQChannelMap::getView(), h_t_MeanADC, h_t_Nhits, h_t_NNoiseSlices, h_t_NoiseMeanHits, h_t_NoiseNHits, h_t_NSlices, h_t_SliceMeanHits, h_t_SliceNHits, h_t_SumADC, hADC_vs_mb, hits(), hMeanADC_vs_Slice_hits, hNNoiseSlices, hNoiseMeanADC_vs_Slice_hits, hNoiseSlice_hits_vs_mb, hNoiseSlicesMeanADC, hNoiseSlicesMinCrossPlaneSeperation, hNoiseSlicesMinInPlaneSeperation, hNoiseSlicesMinSeperation, hNoiseSlicesNHits, hNoiseSlicesSumADC, hNSlices, hSlice_hits_vs_mb, hSlicesMeanADC, hSlicesMinCrossPlaneSeperation, hSlicesMinInPlaneSeperation, hSlicesMinSeperation, hSlicesNHits, hSlicesSumADC, art::Event::id(), makeTrainCVSamples::int, rb::Cluster::IsNoise(), cmap::dataprov::CMap::Map(), rb::Cluster::NCell(), nEvent, NDAPDHVSetting::plane, string, rawdata::RawDigit::TDC(), rb::Cluster::TotalADC(), POTSpillRate::view, and daqchannelmap::X_VIEW.

{
    std::cout << "--- sn::SupernovaAna analyzer. Event: " 
              << e.id().event() 
              << ", event counter: " << nEvent
              << std::endl;
              
    ++nEvent;
    
    art::ServiceHandle<cmap::CMap> cmap;
    
    // Get Cell Hits
    art::Handle< std::vector<rb::CellHit> > hits;
    e.getByLabel(fCellHitsLabel, hits);
    
    const unsigned int nHits = hits->size();
    
    if (fVerbose) std::cout << "\tgot " << nHits << " cell hits" << std::endl;
    
    int sum_adc = 0;
    int min_TDC = 0;
    int max_TDC = 0;
    
    if (nHits > 0){
        min_TDC = hits->at(0).TDC();
        max_TDC = hits->at(0).TDC();
    }
    
    int min_x = 999;
    int max_x = 0;
    int min_y = 999;
    int max_y = 0;
    int min_z = 999;
    int max_z = 0;
    
    // Loop over Cell Hits
    for(unsigned int hitIdx = 0; hitIdx < nHits; ++hitIdx){
        const rb::CellHit* aHit = &(*hits)[hitIdx];
        
        if (aHit->TDC() < min_TDC) min_TDC = aHit->TDC();
        if (aHit->TDC() > max_TDC) max_TDC = aHit->TDC();
        
        // Get plane/cell from cmap
        const int plane = cmap->GetPlane(aHit);
        const int cell  = cmap->GetCell(aHit);
        
        // Set max/min plane
        if (plane < min_z) min_z = plane;
        if (plane > max_z) max_z = plane;
        
        // Set max/min cell
        std::string view = "y";
        if (cmap->Map()->getView(aHit->Channel()) == daqchannelmap::X_VIEW){
            view = "x";
            if (cell < min_x) min_x = cell;
            if (cell > max_x) max_x = cell;
        } else {
            if (cell < min_y) min_y = cell;
            if (cell > max_y) max_y = cell;
        }
        
        sum_adc = sum_adc + int(aHit->ADC());
        hADC_vs_mb->Fill(nEvent,aHit->ADC());
    } // End of loop on cell hits
    
    if (fVerbose){
        std::cout << "\tTDC: " << min_TDC
                  << " -> "    << max_TDC
                  << ", dTDC: "<< max_TDC - min_TDC
                  << std::endl;
        std::cout << "\tX: " << min_x
                  << " -> "  << max_x
                  << ", Y: " << min_y
                  << " -> "  << max_y
                  << ", Z: " << min_z
                  << " -> "  << max_z
                  <<std::endl;
        std::cout << "\tADC sum: " << sum_adc << std::endl;
    } // End of if on verbose
    
    h_t_SumADC->SetBinContent(nEvent,sum_adc);
    h_t_Nhits->SetBinContent(nEvent,nHits);
    
    if (nHits > 0){ 
        h_t_MeanADC->SetBinContent(nEvent, (float(sum_adc) / nHits));
    } else {
        h_t_MeanADC->SetBinContent(nEvent, 0);
    }
    
    
    // Slices
    art::Handle< std::vector< rb::Cluster > > slices;
    e.getByLabel(fSlicerLabel, slices);
    
    double nNoiseSlices  = 0.;
    double nNoiseHits    = 0.;
    double nHitsInSlices = 0.;
    
    const unsigned int nSlices = slices->size();
    
    for(unsigned int sliceIdx = 0; sliceIdx < nSlices; ++sliceIdx){
        art::Ptr<rb::Cluster> aSlice(slices, sliceIdx);
        
        int    min_cell_sep  = 999;
        double min_plane_sep = 999;
        double min_sep       = 999;
        
        if(aSlice->IsNoise()){
            ++nNoiseSlices;
            nNoiseHits += aSlice->NCell();
            
            hNoiseSlice_hits_vs_mb->Fill(nEvent, aSlice->NCell());
            hNoiseMeanADC_vs_Slice_hits->Fill((aSlice->TotalADC() / (double)aSlice->NCell()), aSlice->NCell());
            hNoiseSlicesNHits->Fill(aSlice->NCell());
            hNoiseSlicesSumADC->Fill(aSlice->TotalADC());
            hNoiseSlicesMeanADC->Fill((aSlice->TotalADC() / (double)aSlice->NCell()));
            //FillPlaneSeperations(aSlice, min_cell_sep, min_plane_sep, min_sep);
            hNoiseSlicesMinInPlaneSeperation->Fill(min_cell_sep);
            hNoiseSlicesMinCrossPlaneSeperation->Fill(min_plane_sep);
            hNoiseSlicesMinSeperation->Fill(min_sep);
        } else {
            nHitsInSlices += aSlice->NCell();
            
            hSlice_hits_vs_mb->Fill(nEvent, aSlice->NCell());
            hMeanADC_vs_Slice_hits->Fill((aSlice->TotalADC() / (double)aSlice->NCell()), aSlice->NCell());
            hSlicesNHits->Fill(aSlice->NCell());
            hSlicesSumADC->Fill(aSlice->TotalADC());
            hSlicesMeanADC->Fill((aSlice->TotalADC() / (double)aSlice->NCell()));
            FillPlaneSeperations(aSlice, min_cell_sep, min_plane_sep, min_sep);
            hSlicesMinInPlaneSeperation->Fill(min_cell_sep);
            hSlicesMinCrossPlaneSeperation->Fill(min_plane_sep);
            hSlicesMinSeperation->Fill(min_sep);
        }
    } // End loop on slices
    
    if (fVerbose){
        std::cout << "\tSlices:             " << nSlices << std::endl;
        std::cout << "\tNoise slices:       " << nNoiseSlices << std::endl;
        std::cout << "\tMean hits:          " << nNoiseHits / nNoiseSlices << std::endl;
        std::cout << "\tNon-noise slices:   " << nSlices - nNoiseSlices << std::endl;
        std::cout << "\tMean hits:          " << nHitsInSlices / (nSlices - nNoiseSlices) << std::endl;
    } // End of if on verbose
    
    // Slices
    hNSlices->Fill(nSlices - nNoiseSlices);
    h_t_NSlices->SetBinContent(nEvent, nSlices - nNoiseSlices);
    h_t_SliceNHits->SetBinContent(nEvent, nHitsInSlices);
    if( (nSlices - nNoiseSlices) > 0 ){
        h_t_SliceMeanHits->SetBinContent(nEvent, nHitsInSlices/(nSlices - nNoiseSlices));
    } else {
        h_t_SliceMeanHits->SetBinContent(nEvent, 0);
    }
    
    // Noise slice
    hNNoiseSlices->Fill(nNoiseSlices);
    h_t_NNoiseSlices->SetBinContent(nEvent, nNoiseSlices);
    h_t_NoiseNHits->SetBinContent(nEvent, nNoiseHits);
    if (nNoiseSlices > 0){
        h_t_NoiseMeanHits->SetBinContent(nEvent, nNoiseHits / nNoiseSlices);
    } else {
        h_t_NoiseMeanHits->SetBinContent(nEvent, 0);
    }
    
}

void sn::SupernovaAna::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 147 of file SupernovaAna_module.cc.

References om::cout, and allTimeWatchdog::endl.

{
    std::cout << "--- sn::SupernovaAna begin job" << std::endl;
}

detail::CachedProducts& art::EventObserverBase::cachedProducts ( )

inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

  {
    return selectors_;
  }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it )

inherited

Definition at line 146 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 161 of file Consumer.h.

References T.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it )

inherited

Definition at line 171 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed )

inherited

Referenced by evgen::CosmicPionGen::CosmicPionGen(), evgen::CosmicsGen::CosmicsGen(), fuzz::FuzzyKVertex::FuzzyKVertex(), g4n::G4Gen::G4Gen(), photrans::ImprovedTransport::ImprovedTransport(), art::MixFilter< T >::initEngine_(), evgen::IsoGen::IsoGen(), rsim::MakeNoiseSpectrumFile::MakeNoiseSpectrumFile(), photrans::PhotonSystRewriter::PhotonSystRewriter(), photrans::PhotonTransport::PhotonTransport(), rsim::ReadoutSim::ReadoutSim(), evgen::SingleGen::SingleGen(), and evgen::SupernovaGen::SupernovaGen().

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make  )

inherited

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed, std::string const &  kind_of_engine_to_make, label_t const &  engine_label  )

inherited

CurrentProcessingContext const* art::EDAnalyzer::currentContext ( ) const

protectedinherited

Referenced by art::EDAnalyzer::workerType().

void sn::SupernovaAna::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 153 of file SupernovaAna_module.cc.

References om::cout, allTimeWatchdog::endl, and nEvent.

{
    std::cout << "--- sn::SupernovaAna end job" << std::endl;
    std::cout << "\tevents:                       " << nEvent << std::endl;
}

void sn::SupernovaAna::FillPlaneSeperations	(	art::Ptr< rb::Cluster > &	aSlice,
		int &	min_cell_sep,
		double &	min_plane_sep,
		double &	min_sep
	)

Definition at line 326 of file SupernovaAna_module.cc.

References rb::CellHit::Cell(), rb::Cluster::Cell(), DEFINE_ART_MODULE(), stan::math::fabs(), MECModelEnuComparisons::i, calib::j, rb::Cluster::NCell(), rb::CellHit::Plane(), util::pythag(), sep, and rb::CellHit::View().

Referenced by analyze().

{
    // Loop over the hits and find the minimum seperation for each
    for(unsigned int i = 0; i < aSlice->NCell(); ++i){
      const art::Ptr<rb::CellHit>& ihit = aSlice->Cell(i);
        for (unsigned int j = (i + 1); j < (aSlice->NCell() - 1); ++j){
            const art::Ptr<rb::CellHit>& jhit = aSlice->Cell(j);
            if (ihit->Plane() == jhit->Plane()){
                int sep = fabs(ihit->Cell() - jhit->Cell());
                if (sep < min_cell_sep){
                    min_cell_sep = sep;
                    //min_cell_index = j;
                }
                if (sep < min_sep){
                    min_sep = sep;
                    //min_index = j;
                }                
            } else {
                if (ihit->View() == jhit->View()){
                    double p_s = ihit->Plane() - jhit->Plane();
                    double c_s = ihit->Cell() - jhit->Cell();
                    double sep = util::pythag(p_s, c_s);
                    if (sep < min_plane_sep){
                        min_plane_sep = sep;
                        //min_plane_index = j;
                    }
                    if (sep < min_sep){
                        min_sep = sep;
                        //min_index = j;
                    }
                }
            }
        } // End of inner loop
    } // End of outer loop
}

seed_t art::EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset, char const  key[] = "seed", seed_t const  implicit_seed = -1  )

inherited

Referenced by art::MixFilter< T >::initEngine_().

art::Handle<art::TriggerResults> art::EventObserverBase::getTriggerResults ( Event const &  e ) const

inlineinherited

Definition at line 61 of file EventObserverBase.h.

References art::detail::CachedProducts::getOneTriggerResults(), and art::EventObserverBase::selectors_.

  {
    return selectors_.getOneTriggerResults(e);
  }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it )

inherited

Definition at line 189 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ProductToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::Product,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ProductToken<T>{it};
}

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 204 of file Consumer.h.

References T.

{
  if (!moduleContext_)
    return;

  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
}

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it )

inherited

Definition at line 214 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

{
  if (!moduleContext_)
    return ViewToken<T>::invalid();

  consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                TypeID{typeid(T)},
                                it.label(),
                                it.instance(),
                                it.process());
  return ViewToken<T>{it};
}

bool art::EventObserverBase::modifiesEvent ( ) const

inlineinherited

Definition at line 25 of file EventObserverBase.h.

  {
    return false;
  }

static cet::exempt_ptr<Consumer> art::Consumer::non_module_context ( )

staticinherited

Referenced by art::RootOutput::beginSubRun(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EmptyEvent::readRun_(), and art::EmptyEvent::readSubRun_().

SupernovaAna& sn::SupernovaAna::operator= ( SupernovaAna const &  )

delete

SupernovaAna& sn::SupernovaAna::operator= ( SupernovaAna &&  )

delete

void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset )

protectedinherited

std::string const& art::EventObserverBase::processName ( ) const

inlineinherited

Definition at line 41 of file EventObserverBase.h.

References art::EventObserverBase::process_name_.

Referenced by art::FileDumperOutput::printPrincipal(), and art::RootOutput::RootOutput().

  {
    return process_name_;
  }

void art::EventObserverBase::registerProducts ( MasterProductRegistry &  , ProductDescriptions &  , ModuleDescription const &    )

inlineinherited

Definition at line 33 of file EventObserverBase.h.

References string.

  {}

fhicl::ParameterSetID art::EventObserverBase::selectorConfig ( ) const

inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

  {
    return selector_config_id_;
  }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Referenced by art::RootOutput::endJob().

void art::Consumer::validateConsumedProduct ( BranchType const  bt, ProductInfo const &  pi  )

protectedinherited

bool art::EventObserverBase::wantAllEvents ( ) const

inlineinherited

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

Referenced by art::RootOutput::RootOutput().

  {
    return wantAllEvents_;
  }

bool art::EventObserverBase::wantEvent ( Event const &  e )

inlineinherited

Definition at line 51 of file EventObserverBase.h.

References art::EventObserverBase::selectors_, and art::detail::CachedProducts::wantEvent().

  {
    return selectors_.wantEvent(e);
  }

std::string art::EDAnalyzer::workerType ( ) const

inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

    {
      return "WorkerT<EDAnalyzer>";
    }

Member Data Documentation

std::string sn::SupernovaAna::fCellHitsLabel

private

Label of module making the cell hits.

Definition at line 56 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

std::string sn::SupernovaAna::fSlicerLabel

private

Instance label making the slices.

Definition at line 57 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

bool sn::SupernovaAna::fVerbose

private

Verbose output.

Definition at line 58 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_MeanADC

private

Definition at line 82 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_Nhits

private

Definition at line 81 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_NNoiseSlices

private

Definition at line 86 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_NoiseMeanHits

private

Definition at line 88 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_NoiseNHits

private

Definition at line 87 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_NSlices

private

Definition at line 83 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_SliceMeanHits

private

Definition at line 85 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_SliceNHits

private

Definition at line 84 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::h_t_SumADC

private

Definition at line 80 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH2F* sn::SupernovaAna::hADC_vs_mb

private

Definition at line 91 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH2F* sn::SupernovaAna::hMeanADC_vs_Slice_hits

private

Definition at line 93 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hNNoiseSlices

private

Definition at line 71 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH2F* sn::SupernovaAna::hNoiseMeanADC_vs_Slice_hits

private

Definition at line 95 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH2F* sn::SupernovaAna::hNoiseSlice_hits_vs_mb

private

Definition at line 94 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hNoiseSlicesMeanADC

private

Definition at line 74 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hNoiseSlicesMinCrossPlaneSeperation

private

Definition at line 76 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hNoiseSlicesMinInPlaneSeperation

private

Definition at line 75 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hNoiseSlicesMinSeperation

private

Definition at line 77 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hNoiseSlicesNHits

private

Definition at line 72 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hNoiseSlicesSumADC

private

Definition at line 73 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hNSlices

private

Definition at line 64 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH2F* sn::SupernovaAna::hSlice_hits_vs_mb

private

Definition at line 92 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hSlicesMeanADC

private

Definition at line 67 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hSlicesMinCrossPlaneSeperation

private

Definition at line 69 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hSlicesMinInPlaneSeperation

private

Definition at line 68 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hSlicesMinSeperation

private

Definition at line 70 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hSlicesNHits

private

Definition at line 65 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

TH1F* sn::SupernovaAna::hSlicesSumADC

private

Definition at line 66 of file SupernovaAna_module.cc.

Referenced by analyze(), and SupernovaAna().

int sn::SupernovaAna::nEvent = 0

private

Definition at line 61 of file SupernovaAna_module.cc.

Referenced by analyze(), and endJob().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N20-11-28/Supernova/SupernovaAna_module.cc