Inheritance diagram for pdsp::WireChgAna:

Public Member Functions
	WireChgAna (fhicl::ParameterSet const &p)

	WireChgAna (WireChgAna const &)=delete

	WireChgAna (WireChgAna &&)=delete

WireChgAna &	operator= (WireChgAna const &)=delete

WireChgAna &	operator= (WireChgAna &&)=delete

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

void	beginJob () override

void	endJob () override

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
TTree *	ftree

TH1D *	rawwf [3][800]

TH1D *	deconwf [3][800]

double	charge [3][800]

double	deconchg [3][800]

double	hitchg [3][800]

double	hitsumadc [3][800]

double	nelec [3][800]

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 35 of file WireChgAna_module.cc.

Constructor & Destructor Documentation

pdsp::WireChgAna::WireChgAna ( fhicl::ParameterSet const & p )

explicit

Definition at line 69 of file WireChgAna_module.cc.

   : EDAnalyzer{p}  // ,
   // More initializers here.
 {
   // Call appropriate consumes<>() for any products to be retrieved by this module.
 }

pdsp::WireChgAna::WireChgAna ( WireChgAna const & )

delete

pdsp::WireChgAna::WireChgAna ( WireChgAna && )

delete

Member Function Documentation

void pdsp::WireChgAna::analyze ( art::Event const & e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 76 of file WireChgAna_module.cc.

 {
 
   for (int i = 0; i<3; ++i){
     for (int j = 0; j<800; ++j){
       charge[i][j] = 0;
       nelec[i][j] = 0;
       deconchg[i][j] = 0;
       hitchg[i][j] = 0;
       hitsumadc[i][j] = 0;
     }
   }
 
   // Implementation of required member function here.
   std::vector<art::Ptr<raw::RawDigit> > rawlist;
   art::InputTag itag("tpcrawdecoder","daq");
   auto rawListHandle = e.getHandle< std::vector<raw::RawDigit> >(itag);
   if (rawListHandle) art::fill_ptr_vector(rawlist, rawListHandle);
 
   std::vector < art::Ptr < recob::Wire > > wires;
   art::InputTag itag2("wclsdatanfsp","gauss");
   auto wireListHandle = e.getHandle< std::vector < recob::Wire > >(itag2);
   if (wireListHandle) {
     art::fill_ptr_vector(wires, wireListHandle);
   }
 
   std::vector <art::Ptr<recob::Hit> > hits;
   auto hitListHandle = e.getHandle<std::vector<recob::Hit> >("gaushit");
   if (hitListHandle) {
     art::fill_ptr_vector(hits, hitListHandle);
   }
 
   auto const* geo = lar::providerFrom<geo::Geometry>();
 
   auto simChannelHandle = e.getValidHandle<std::vector<sim::SimChannel>>("tpcrawdecoder:simpleSC");
 
   for (unsigned int ich = 0; ich < rawlist.size(); ++ich){
     const auto & digitVec = rawlist[ich];
     const auto & wireid = geo->ChannelToWire(digitVec->Channel());
     if (wireid[0].TPC != 1) continue;
     //if (wireid[0].Plane != 2) continue;
     int wire = wireid[0].Wire;
     int plane = wireid[0].Plane;
     double this_charge = 0;
     std::vector<short> rawadc(6000);
     raw::Uncompress(digitVec->ADCs(), rawadc, digitVec->GetPedestal(), digitVec->Compression());
     for (size_t itck = 0; itck < rawadc.size(); ++itck){
       this_charge += rawadc[itck] - digitVec->GetPedestal();
       rawwf[plane][wire]->SetBinContent(itck, rawadc[itck] - digitVec->GetPedestal());
     }
     charge[plane][wire] = this_charge;
   }
 
   for ( auto const& channel : (*simChannelHandle) ){
     const auto & wireid = geo->ChannelToWire(channel.Channel());
     if (wireid[0].TPC != 1) continue;
     //if (wireid[0].Plane != 2) continue;
     int wire = wireid[0].Wire;
     int plane = wireid[0].Plane;
     double this_nelec = 0;
     auto const& timeSlices = channel.TDCIDEMap();
     for ( auto const& timeSlice : timeSlices ){
       auto const& energyDeposits = timeSlice.second;
       for ( auto const& energyDeposit : energyDeposits ){
         this_nelec += energyDeposit.numElectrons;
       }
     }
     nelec[plane][wire] = this_nelec;
   }
 
   for (unsigned int ich = 0; ich < wires.size(); ++ich){
     const auto & wireid = geo->ChannelToWire(wires[ich]->Channel());
     if (wireid[0].TPC != 1) continue;
     //if (wireid[0].Plane != 2) continue;
     int wire = wireid[0].Wire;
     int plane = wireid[0].Plane;
     double this_charge = 0;
     const auto & signal = wires[ich]->Signal();
 
     for (size_t itck = 0; itck < signal.size(); ++itck){
       this_charge += signal[itck];
       deconwf[plane][wire]->SetBinContent(itck, signal[itck]);
     }
     deconchg[plane][wire] = this_charge;
   }
 
   for (auto const & hit : hits){
     if (hit->WireID().TPC != 1) continue;
     int wire = hit->WireID().Wire;
     int plane = hit->WireID().Plane;
     hitchg[plane][wire] += hit->Integral();
     hitsumadc[plane][wire] += hit->SummedADC();
   }
 
   ftree->Fill();
 }

void pdsp::WireChgAna::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 173 of file WireChgAna_module.cc.

 {
   art::ServiceHandle<art::TFileService> fileServiceHandle;
   ftree = fileServiceHandle->make<TTree>("ftree", "raw digit info");
   ftree->Branch("charge", charge, "charge[3][800]/D");
   ftree->Branch("nelec", nelec, "nelec[3][800]/D");
   ftree->Branch("deconchg", deconchg, "deconchg[3][800]/D");
   ftree->Branch("hitchg", hitchg, "hitchg[3][800]/D");
   ftree->Branch("hitsumadc", hitsumadc, "hitsumadc[3][800]/D");
 
   for (int i = 0; i<3; ++i){
     for (int j = 0; j<800; ++j){
       rawwf[i][j] = fileServiceHandle->make<TH1D>(Form("rawwf_%d_%d",i,j),Form("Plane %d Wire %d",i,j),6000,0,6000);
       deconwf[i][j] = fileServiceHandle->make<TH1D>(Form("deconwf_%d_%d",i,j),Form("Plane %d Wire %d",i,j),6000,0,6000);
     }
   }
 
 }

void pdsp::WireChgAna::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 192 of file WireChgAna_module.cc.

 {
   // Implementation of optional member function here.
 }

WireChgAna& pdsp::WireChgAna::operator= ( WireChgAna const & )

delete

WireChgAna& pdsp::WireChgAna::operator= ( WireChgAna && )

delete

Member Data Documentation

double pdsp::WireChgAna::charge[3][800]

private

Definition at line 60 of file WireChgAna_module.cc.

double pdsp::WireChgAna::deconchg[3][800]

private

Definition at line 61 of file WireChgAna_module.cc.

TH1D* pdsp::WireChgAna::deconwf[3][800]

private

Definition at line 59 of file WireChgAna_module.cc.

TTree* pdsp::WireChgAna::ftree

private

Definition at line 57 of file WireChgAna_module.cc.

double pdsp::WireChgAna::hitchg[3][800]

private

Definition at line 62 of file WireChgAna_module.cc.

double pdsp::WireChgAna::hitsumadc[3][800]

private

Definition at line 63 of file WireChgAna_module.cc.

double pdsp::WireChgAna::nelec[3][800]

private

Definition at line 64 of file WireChgAna_module.cc.

TH1D* pdsp::WireChgAna::rawwf[3][800]

private

Definition at line 58 of file WireChgAna_module.cc.

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

protoduneana/protoduneana/singlephase/truepion/WireChgAna_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation