Inheritance diagram for wc::CellTree:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

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

virtual	~CellTree ()

void	beginJob ()

void	endJob ()

void	beginRun (const art::Run &run)

void	analyze (const art::Event &evt)

void	reconfigure (fhicl::ParameterSet const &pset)

void	initOutput ()

void	printEvent ()

void	print_vector (ostream &out, vector< double > &v, TString desc, bool end=false)

void	processRaw (const art::Event &evt)

void	processCalib (const art::Event &evt)

void	processOpHit (const art::Event &evt)

void	processOpFlash (const art::Event &evt)

void	processSpacePoint (const art::Event &event, TString option, ostream &out=cout)

void	processSimChannel (const art::Event &evt)

void	processMC (const art::Event &evt)

void	processMCTracks ()

void	processTrigger (const art::Event &evt)

void	reset ()

void	InitProcessMap ()

bool	IsPrimary (int i)

bool	KeepMC (int i)

double	KE (float *momentum)

TString	PDGName (int pdg)

bool	DumpMCJSON (int id, ostream &out)

void	DumpMCJSON (ostream &out=cout)

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	fRawDigitLabel

std::string	fCalibLabel

std::string	fOpHitLabel

std::string	fOpFlashLabel

std::string	fTriggerLabel

std::vector< std::string >	fSpacePointLabels

std::string	fOutFileName

std::string	mcOption

int	nRawSamples

float	opMultPEThresh

bool	fSaveMCTrackPoints

bool	fSaveSimChannel

bool	fSaveRaw

bool	fSaveCalib

bool	fSaveOpHit

bool	fSaveOpFlash

bool	fSaveMC

bool	fSaveTrigger

bool	fSaveJSON

art::ServiceHandle< geo::Geometry >	fGeometry

TFile *	fOutFile

TTree *	fEventTree

std::map< std::string, int >	processMap

std::map< int, int >	savedMCTrackIdMap

int	entryNo

int	fEvent

int	fRun

int	fSubRun

double	fEventTime

unsigned int	fTriggernumber

double	fTriggertime

double	fBeamgatetime

unsigned int	fTriggerbits

int	fCalib_nChannel

std::vector< int >	fCalib_channelId

TClonesArray *	fCalib_wf

int	oh_nHits

vector< int >	oh_channel

vector< double >	oh_bgtime

vector< double >	oh_trigtime

vector< double >	oh_pe

int	of_nFlash

vector< float >	of_t

vector< float >	of_peTotal

vector< int >	of_multiplicity

TClonesArray *	fPEperOpDet

int	fRaw_nChannel

std::vector< int >	fRaw_channelId

TClonesArray *	fRaw_wf

int	fSIMIDE_size

vector< int >	fSIMIDE_channelIdY

vector< int >	fSIMIDE_trackId

vector< unsigned short >	fSIMIDE_tdc

vector< float >	fSIMIDE_x

vector< float >	fSIMIDE_y

vector< float >	fSIMIDE_z

vector< float >	fSIMIDE_numElectrons

int	mc_Ntrack

int	mc_id [MAX_TRACKS]

int	mc_pdg [MAX_TRACKS]

int	mc_process [MAX_TRACKS]

int	mc_mother [MAX_TRACKS]

float	mc_startXYZT [MAX_TRACKS][4]

float	mc_endXYZT [MAX_TRACKS][4]

float	mc_startMomentum [MAX_TRACKS][4]

float	mc_endMomentum [MAX_TRACKS][4]

std::vector< std::vector< int > >	mc_daughters

TObjArray *	fMC_trackPosition

int	mc_isnu

int	mc_nGeniePrimaries

int	mc_nu_pdg

int	mc_nu_ccnc

int	mc_nu_mode

int	mc_nu_intType

int	mc_nu_target

int	mc_hitnuc

int	mc_hitquark

double	mc_nu_Q2

double	mc_nu_W

double	mc_nu_X

double	mc_nu_Y

double	mc_nu_Pt

double	mc_nu_Theta

float	mc_nu_pos [4]

float	mc_nu_mom [4]

std::map< int, int >	trackIndex

std::vector< std::vector< int > >	trackParents

std::vector< std::vector< int > >	trackChildren

std::vector< std::vector< int > >	trackSiblings

TDatabasePDG *	dbPDG

Detailed Description

Definition at line 76 of file CellTree_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

wc::CellTree::CellTree ( fhicl::ParameterSet const & pset )

explicit

Definition at line 232 of file CellTree_module.cc.

References dbPDG, entryNo, initOutput(), InitProcessMap(), and reconfigure().

     : EDAnalyzer(parameterSet)
 {
     dbPDG = new TDatabasePDG();
 
     reconfigure(parameterSet);
     InitProcessMap();
     initOutput();
     entryNo = 0;
 }

wc::CellTree::~CellTree ( )

virtual

Definition at line 244 of file CellTree_module.cc.

245 {

246 }

Member Function Documentation

void wc::CellTree::analyze ( const art::Event & evt )

Definition at line 408 of file CellTree_module.cc.

References DumpMCJSON(), entryNo, fEvent, fEventTime, fEventTree, fRun, fSaveCalib, fSaveJSON, fSaveMC, fSaveOpFlash, fSaveOpHit, fSaveRaw, fSaveSimChannel, fSaveTrigger, fSpacePointLabels, fSubRun, processCalib(), processMC(), processMCTracks(), processOpFlash(), processOpHit(), processRaw(), processSimChannel(), processSpacePoint(), processTrigger(), reset(), art::Timestamp::timeHigh(), and art::Timestamp::timeLow().

 {
     reset();
     fEvent  = event.id().event();
     fRun    = event.run();
     fSubRun = event.subRun();
     art::Timestamp ts = event.time();
     TTimeStamp tts(ts.timeHigh(), ts.timeLow());
     fEventTime = tts.AsDouble();
 
     if (fSaveRaw) processRaw(event);
     if (fSaveCalib) processCalib(event);
     if (fSaveOpHit) processOpHit(event);
     if (fSaveOpFlash) processOpFlash(event);
     if (fSaveSimChannel) processSimChannel(event);
     if (fSaveMC) processMC(event);
     if (fSaveTrigger) processTrigger(event);
 
     if (fSaveJSON) {
         gSystem->MakeDirectory(TString::Format("data/%i", entryNo).Data());
         int nSp = fSpacePointLabels.size();
         for (int i=0; i<nSp; i++) {
             TString jsonfile;
             jsonfile.Form("data/%i/%i-%s.json", entryNo, entryNo, fSpacePointLabels[i].c_str());
             std::ofstream out(jsonfile.Data());
             processSpacePoint(event, fSpacePointLabels[i], out);
             out.close();
         }
 
         if(fSaveMC) {
             processMCTracks();
             TString jsonfile;
             jsonfile.Form("data/%i/%i-mc.json", entryNo, entryNo);
             std::ofstream out(jsonfile.Data());
             DumpMCJSON(out);
             out.close();
         }
     }
 
     // printEvent();
     fEventTree->Fill();
 
     entryNo++;
 }

void wc::CellTree::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 375 of file CellTree_module.cc.

 {
 
 }

void wc::CellTree::beginRun ( const art::Run & run )

Definition at line 402 of file CellTree_module.cc.

 {
   mf::LogInfo("CellTree") << "begin run";
 }

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

inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

Referenced by art::EDAnalyzer::doEvent(), and art::OutputModule::doWriteEvent().

   {
     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 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   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 162 of file Consumer.h.

 {
   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 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   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};
 }

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

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

 {
   return rng()->createEngine(placeholder_schedule_id(), seed);
 }

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

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make);
 }

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

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
 }

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

protectedinherited

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

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

   {
     return current_context_.get();
   }

bool wc::CellTree::DumpMCJSON	(	int	id,
		ostream &	out
	)

Definition at line 966 of file CellTree_module.cc.

References e, KE(), KeepMC(), mc_daughters, mc_endXYZT, mc_pdg, mc_startMomentum, mc_startXYZT, PDGName(), and trackIndex.

Referenced by analyze(), and DumpMCJSON().

 {
     int i = trackIndex[id];
     if (!KeepMC(i)) return false;
 
     int e = KE(mc_startMomentum[i])*1000;
 
     int nDaughter = mc_daughters.at(i).size();
     vector<int> saved_daughters;
     for (int j=0; j<nDaughter; j++) {
         int daughter_id = mc_daughters.at(i).at(j);
         // int e_daughter = KE(mc_startMomentum[ trackIndex[daughter_id] ])*1000;
         // if (e_daughter >= thresh_KE) {
         if ( KeepMC(trackIndex[daughter_id]) ) {
             saved_daughters.push_back(daughter_id);
         }
     }
 
     out << fixed << setprecision(1);
     out << "{";
 
     out << "\"id\":" << id << ",";
     out << "\"text\":" << "\"" << PDGName(mc_pdg[i]) << "  " << e << " MeV\",";
     out << "\"data\":{";
     out << "\"start\":[" << mc_startXYZT[i][0] << ", " <<  mc_startXYZT[i][1] << ", " << mc_startXYZT[i][2] << "],";
     out << "\"end\":[" << mc_endXYZT[i][0] << ", " <<  mc_endXYZT[i][1] << ", " << mc_endXYZT[i][2] << "]";
     out << "},";
     out << "\"children\":[";
     int nSavedDaughter = saved_daughters.size();
     if (nSavedDaughter == 0) {
         out << "],";
         out << "\"icon\":" << "\"jstree-file\"";
         out << "}";
         return true;
     }
     else {
         for (int j=0; j<nSavedDaughter; j++) {
             DumpMCJSON(saved_daughters.at(j), out);
             if (j!=nSavedDaughter-1) {
                 out << ",";
             }
         }
         out << "]";
         out << "}";
         return true;
     }
 }

void wc::CellTree::DumpMCJSON ( ostream & out = cout )

Definition at line 1015 of file CellTree_module.cc.

References DumpMCJSON(), IsPrimary(), KeepMC(), mc_id, and mc_Ntrack.

 {
     out << "[";
     vector<int> primaries;
     for (int i=0; i<mc_Ntrack; i++) {
         if (IsPrimary(i)) {
             // int e = KE(mc_startMomentum[i])*1000;
             // if (e<thresh_KE) continue;
             if (KeepMC(i)) {
                 primaries.push_back(i);
             }
         }
     }
     int size = primaries.size();
     // cout << size << endl;
     for (int i=0; i<size; i++) {
         if (DumpMCJSON(mc_id[primaries[i]], out) && i!=size-1) {
             out << ", ";
         }
     }
 
     out << "]";
 }

void wc::CellTree::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 383 of file CellTree_module.cc.

References fEventTree, fOutFile, fSaveJSON, and system().

 {
     // Write fEventTree to file
     TDirectory* tmpDir = gDirectory;
     fOutFile->cd("/Event");
 
     fEventTree->Write();
 
     gDirectory = tmpDir;
 
     fOutFile->Close();
 
     if (fSaveJSON) {
         system("zip -r bee_upload data");
         gSystem->ChangeDirectory("..");
     }
 }

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

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

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

 {
   auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
   return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
 }

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_.

Referenced by art::OutputModule::doWriteEvent().

   {
     return selectors_.getOneTriggerResults(e);
   }

void wc::CellTree::initOutput ( )

Definition at line 272 of file CellTree_module.cc.

References fBeamgatetime, fCalib_channelId, fCalib_nChannel, fCalib_wf, fEvent, fEventTime, fEventTree, fMC_trackPosition, fOutFile, fOutFileName, fPEperOpDet, fRaw_channelId, fRaw_nChannel, fRaw_wf, fRun, fSaveJSON, fSIMIDE_channelIdY, fSIMIDE_numElectrons, fSIMIDE_size, fSIMIDE_tdc, fSIMIDE_trackId, fSIMIDE_x, fSIMIDE_y, fSIMIDE_z, fSubRun, fTriggerbits, fTriggernumber, fTriggertime, mc_daughters, mc_endMomentum, mc_endXYZT, mc_hitnuc, mc_hitquark, mc_id, mc_isnu, mc_mother, mc_nGeniePrimaries, mc_Ntrack, mc_nu_ccnc, mc_nu_intType, mc_nu_mode, mc_nu_mom, mc_nu_pdg, mc_nu_pos, mc_nu_Pt, mc_nu_Q2, mc_nu_target, mc_nu_Theta, mc_nu_W, mc_nu_X, mc_nu_Y, mc_pdg, mc_process, mc_startMomentum, mc_startXYZT, of_multiplicity, of_nFlash, of_peTotal, of_t, oh_bgtime, oh_channel, oh_nHits, oh_pe, oh_trigtime, and system().

Referenced by CellTree().

 {
     TDirectory* tmpDir = gDirectory;
 
     fOutFile = new TFile(fOutFileName.c_str(), "recreate");
 
     // 3.1: add mc_trackPosition
     TNamed version("version", "4.0");
     version.Write();
 
     // init Event TTree
     TDirectory* subDir = fOutFile->mkdir("Event");
     subDir->cd();
     fEventTree = new TTree("Sim", "Event Tree from Simulation");
     fEventTree->Branch("eventNo", &fEvent);
     fEventTree->Branch("runNo", &fRun);
     fEventTree->Branch("subRunNo", &fSubRun);
     fEventTree->Branch("eventTime", &fEventTime);  // timestamp
 
     fEventTree->Branch("triggerNo", &fTriggernumber);  // timestamp
     fEventTree->Branch("triggerTime", &fTriggertime);  // timestamp
     fEventTree->Branch("beamgateTime", &fBeamgatetime);  // timestamp
     fEventTree->Branch("triggerBits", &fTriggerbits);  // timestamp
 
     fEventTree->Branch("raw_nChannel", &fRaw_nChannel);  // number of hit channels above threshold
     fEventTree->Branch("raw_channelId" , &fRaw_channelId); // hit channel id; size == raw_nChannel
     fRaw_wf = new TClonesArray("TH1F");
     fEventTree->Branch("raw_wf", &fRaw_wf, 256000, 0);  // raw waveform adc of each channel
 
 
     fEventTree->Branch("calib_nChannel", &fCalib_nChannel);  // number of hit channels above threshold
     fEventTree->Branch("calib_channelId" , &fCalib_channelId); // hit channel id; size == calib_Nhit
     fCalib_wf = new TClonesArray("TH1F");
     fEventTree->Branch("calib_wf", &fCalib_wf, 256000, 0);  // calib waveform adc of each channel
     // fCalib_wf->BypassStreamer();
     // fEventTree->Branch("calib_wfTDC", &fCalib_wfTDC);  // calib waveform tdc of each channel
 
     fEventTree->Branch("oh_nHits", &oh_nHits); // number of op hits
     fEventTree->Branch("oh_channel", &oh_channel); //opchannel id; size == no ophits
     fEventTree->Branch("oh_bgtime", &oh_bgtime); // optical pulse peak time w.r.t. start of beam gate
     fEventTree->Branch("oh_trigtime", &oh_trigtime); // optical pulse peak time w.r.t. trigger
     fEventTree->Branch("oh_pe", &oh_pe); // pulse PE
 
     fEventTree->Branch("of_nFlash", &of_nFlash);
     fEventTree->Branch("of_t", &of_t); // time in us w.r.t. the trigger for each flash
     fEventTree->Branch("of_peTotal", &of_peTotal); // total PE (sum of all PMTs) for each flash
     fEventTree->Branch("of_multiplicity", &of_multiplicity); // total number of PMTs above threshold for each flash
     fPEperOpDet = new TClonesArray("TH1F");
     fEventTree->Branch("pe_opdet", &fPEperOpDet, 256000, 0);
 
     fEventTree->Branch("simide_size", &fSIMIDE_size);  // size of stored sim:IDE
     fEventTree->Branch("simide_channelIdY", &fSIMIDE_channelIdY);
     fEventTree->Branch("simide_trackId", &fSIMIDE_trackId);
     fEventTree->Branch("simide_tdc", &fSIMIDE_tdc);
     fEventTree->Branch("simide_x", &fSIMIDE_x);
     fEventTree->Branch("simide_y", &fSIMIDE_y);
     fEventTree->Branch("simide_z", &fSIMIDE_z);
     fEventTree->Branch("simide_numElectrons", &fSIMIDE_numElectrons);
 
     fEventTree->Branch("mc_Ntrack", &mc_Ntrack);  // number of tracks in MC
     fEventTree->Branch("mc_id", &mc_id, "mc_id[mc_Ntrack]/I");  // track id; size == mc_Ntrack
     fEventTree->Branch("mc_pdg", &mc_pdg, "mc_pdg[mc_Ntrack]/I");  // track particle pdg; size == mc_Ntrack
     fEventTree->Branch("mc_process", &mc_process, "mc_process[mc_Ntrack]/I");  // track generation process code; size == mc_Ntrack
     fEventTree->Branch("mc_mother", &mc_mother, "mc_mother[mc_Ntrack]/I");  // mother id of this track; size == mc_Ntrack
     fEventTree->Branch("mc_daughters", &mc_daughters);  // daughters id of this track; vector
     fEventTree->Branch("mc_startXYZT", &mc_startXYZT, "mc_startXYZT[mc_Ntrack][4]/F");  // start position of this track; size == mc_Ntrack
     fEventTree->Branch("mc_endXYZT", &mc_endXYZT, "mc_endXYZT[mc_Ntrack][4]/F");  // start position of this track; size == mc_Ntrack
     fEventTree->Branch("mc_startMomentum", &mc_startMomentum, "mc_startMomentum[mc_Ntrack][4]/F");  // start momentum of this track; size == mc_Ntrack
     fEventTree->Branch("mc_endMomentum", &mc_endMomentum, "mc_endMomentum[mc_Ntrack][4]/F");  // start momentum of this track; size == mc_Ntrack
     fMC_trackPosition = new TObjArray();
     fMC_trackPosition->SetOwner(kTRUE);
     fEventTree->Branch("mc_trackPosition", &fMC_trackPosition);
 
     fEventTree->Branch("mc_isnu", &mc_isnu);
     fEventTree->Branch("mc_nGeniePrimaries", &mc_nGeniePrimaries);
     fEventTree->Branch("mc_nu_pdg", &mc_nu_pdg);
     fEventTree->Branch("mc_nu_ccnc", &mc_nu_ccnc);
     fEventTree->Branch("mc_nu_mode", &mc_nu_mode);
     fEventTree->Branch("mc_nu_intType", &mc_nu_intType);
     fEventTree->Branch("mc_nu_target", &mc_nu_target);
     fEventTree->Branch("mc_hitnuc", &mc_hitnuc);
     fEventTree->Branch("mc_hitquark", &mc_hitquark);
     fEventTree->Branch("mc_nu_Q2", &mc_nu_Q2);
     fEventTree->Branch("mc_nu_W", &mc_nu_W);
     fEventTree->Branch("mc_nu_X", &mc_nu_X);
     fEventTree->Branch("mc_nu_Y", &mc_nu_Y);
     fEventTree->Branch("mc_nu_Pt", &mc_nu_Pt);
     fEventTree->Branch("mc_nu_Theta", &mc_nu_Theta);
     fEventTree->Branch("mc_nu_pos", &mc_nu_pos, "mc_nu_pos[4]/F");
     fEventTree->Branch("mc_nu_mom", &mc_nu_mom, "mc_nu_mom[4]/F");
 
     gDirectory = tmpDir;
 
     if (fSaveJSON) {
         system("rm -rf bee");
         gSystem->MakeDirectory("bee");
         gSystem->ChangeDirectory("bee");
         gSystem->MakeDirectory("data");
     }
 
 }

void wc::CellTree::InitProcessMap ( )

Definition at line 1113 of file CellTree_module.cc.

References DEFINE_ART_MODULE, and processMap.

Referenced by CellTree().

 {
     processMap["unknown"]              = 0;
     processMap["primary"]              = 1;
     processMap["compt"]                = 2;
     processMap["phot"]                 = 3;
     processMap["annihil"]              = 4;
     processMap["eIoni"]                = 5;
     processMap["eBrem"]                = 6;
     processMap["conv"]                 = 7;
     processMap["muIoni"]               = 8;
     processMap["muMinusCaptureAtRest"] = 9;
     processMap["NeutronInelastic"]     = 10;
     processMap["nCapture"]             = 11;
     processMap["hadElastic"]           = 12;
     processMap["Decay"]                = 13;
     processMap["CoulombScat"]          = 14;
     processMap["muPairProd"]           = 15;
     processMap["muBrems"]              = 16;
     processMap["muPairProd"]           = 17;
     processMap["PhotonInelastic"]      = 18;
     processMap["hIoni"]                = 19;
     processMap["ProtonInelastic"]      = 20;
     processMap["PionPlusInelastic"]    = 21;
     processMap["CHIPSNuclearCaptureAtRest"] = 22;
     processMap["PionMinusInelastic"]   = 23;
 }

bool wc::CellTree::IsPrimary ( int i )

inline

Definition at line 105 of file CellTree_module.cc.

Referenced by DumpMCJSON(), and processMCTracks().

105 { return mc_mother[i] == 0 ; }

wc::CellTree::mc_mother

int mc_mother[MAX_TRACKS]

Definition: CellTree_module.cc:195

double wc::CellTree::KE ( float * momentum )

Definition at line 1040 of file CellTree_module.cc.

Referenced by DumpMCJSON(), and KeepMC().

 {
     TLorentzVector particle(momentum);
     return particle.E()-particle.M();
 }

bool wc::CellTree::KeepMC ( int i )

Definition at line 1047 of file CellTree_module.cc.

References e, KE(), mc_pdg, and mc_startMomentum.

Referenced by DumpMCJSON().

 {
     double e = KE(mc_startMomentum[i])*1000;
     double thresh_KE_em = 5.; // MeV
     double thresh_KE_np = 10; // MeV
     if (mc_pdg[i]==22 || mc_pdg[i]==11 || mc_pdg[i]==-11) {
         if (e>=thresh_KE_em) return true;
         else return false;
     }
     else if (mc_pdg[i]==2112 || mc_pdg[i]==2212 || mc_pdg[i]>1e9) {
         if (e>=thresh_KE_np) return true;
         else return false;
     }
     return true;
 }

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 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   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 205 of file Consumer.h.

 {
   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 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   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;
   }

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

staticinherited

Definition at line 76 of file Consumer.cc.

 {
   static Consumer invalid{InvalidTag{}};
   return &invalid;
 }

TString wc::CellTree::PDGName ( int pdg )

Definition at line 1064 of file CellTree_module.cc.

References dbPDG, and z.

Referenced by DumpMCJSON().

 {
     TParticlePDG *p = dbPDG->GetParticle(pdg);
     if (p == 0) {
         if (pdg>1e9) {
             int z = (pdg - 1e9) / 10000;
             int a = (pdg - 1e9 - z*1e4) / 10;
             TString name;
             if (z == 18) name = "Ar";
 
             else if (z == 17) name = "Cl";
             else if (z == 19) name = "Ca";
             else if (z == 16) name = "S";
             else if (z == 15) name = "P";
             else if (z == 14) name = "Si";
             else if (z == 1) name = "H";
             else if (z == 2) name = "He";
 
             else return Form("%i", pdg);
             return Form("%s-%i", name.Data(), a);
         }
         return Form("%i", pdg);
     }
     else {
         return p->GetName();
     }
 }

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

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

 {
   if (!moduleContext_)
     return;
 
   pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
   for (auto& consumablesPerBranch : consumables_) {
     cet::sort_all(consumablesPerBranch);
   }
 }

void wc::CellTree::print_vector	(	ostream &	out,
		vector< double > &	v,
		TString	desc,
		bool	end = `false`
	)

Definition at line 874 of file CellTree_module.cc.

Referenced by processSpacePoint().

 {
     int N = v.size();
 
     out << '"' << desc << '"' << ":[";
     for (int i=0; i<N; i++) {
         out << v[i];
         if (i!=N-1) {
             out << ",";
         }
     }
     out << "]";
     if (!end) out << ",";
     out << endl;
 }

void wc::CellTree::printEvent ( )

Definition at line 1093 of file CellTree_module.cc.

References fEvent, fRun, fSubRun, mc_daughters, mc_endMomentum, mc_endXYZT, mc_id, mc_mother, mc_Ntrack, mc_pdg, mc_startMomentum, and mc_startXYZT.

 {
     cout << " Run/SubRun/Event: " << fRun << "/" << fSubRun << "/" << fEvent << endl;
     cout << "      Ntracks:" << mc_Ntrack << endl;
 
     for (int i=0; i<mc_Ntrack; i++) {
         cout << "\n              id: " << mc_id[i];
         cout << "\n             pdg: " << mc_pdg[i];
         cout << "\n          mother: " << mc_mother[i];
         cout << "\n      Ndaughters: " << mc_daughters.at(i).size();
         cout << "\n      start XYZT: (" << mc_startXYZT[i][0] << ", " << mc_startXYZT[i][1] << ", " << mc_startXYZT[i][2] << ", " << mc_startXYZT[i][3] << ")";
         cout << "\n        end XYZT: (" << mc_endXYZT[i][0] << ", " << mc_endXYZT[i][1] << ", " << mc_endXYZT[i][2] << ", " << mc_endXYZT[i][3] << ")";
         cout << "\n  start momentum: (" << mc_startMomentum[i][0] << ", " << mc_startMomentum[i][1] << ", " << mc_startMomentum[i][2] << ", " << mc_startMomentum[i][3] << ")";
         cout << "\n    end momentum: (" << mc_endMomentum[i][0] << ", " << mc_endMomentum[i][1] << ", " << mc_endMomentum[i][2] << ", " << mc_endMomentum[i][3] << ")";
 
         cout << endl;
     }
 }

void wc::CellTree::processCalib ( const art::Event & evt )

Definition at line 559 of file CellTree_module.cc.

References fCalib_channelId, fCalib_nChannel, fCalibLabel, art::fill_ptr_vector(), art::DataViewImpl::getByLabel(), and nRawSamples.

Referenced by analyze().

 {
 
     art::Handle< std::vector<recob::Wire> > wires_handle;
     if (! event.getByLabel(fCalibLabel, wires_handle)) {
         cout << "WARNING: no label " << fCalibLabel << endl;
         return;
     }
     std::vector< art::Ptr<recob::Wire> >  wires;
     art::fill_ptr_vector(wires, wires_handle);
 
     // wires size should == Nchannels == 1992; (no hit channel has a flat 0-waveform)
     // cout << "\n wires size: " << wires.size() << endl;
     fCalib_nChannel = wires.size();
 
     int i=0;
     for (auto const& wire: wires) {
         std::vector<float> calibwf = wire->Signal();
         int chanId = wire->Channel();
         fCalib_channelId.push_back(chanId);
         TH1F *h = new((*fCalib_wf)[i]) TH1F("", "", nRawSamples, 0, nRawSamples);
         for (int j=1; j<=nRawSamples; j++) {
             h->SetBinContent(j, calibwf[j]);
         }
         // fCalib_wf.push_back(calibwf);
         // cout << chanId << ", " << nSamples << endl;
         i++;
     }
 
 }

void wc::CellTree::processMC ( const art::Event & evt )

Definition at line 681 of file CellTree_module.cc.

Referenced by analyze().

 {
     art::Handle< std::vector<simb::MCParticle> > particleHandle;
     if (! event.getByLabel("largeant", particleHandle)) return;
     std::vector< art::Ptr<simb::MCParticle> > particles;
     art::fill_ptr_vector(particles, particleHandle);
 
     art::Handle< std::vector<sim::SimChannel> > simChannelHandle;
     event.getByLabel("largeant", simChannelHandle);
 
     // art::ServiceHandle<cheat::BackTracker> bt;
     art::FindOneP<simb::MCTruth> fo(particleHandle, event, "largeant");
 
     int i=0; // track index in saved MCParticles;
     int i_all=0; // track index in all MCParticles;
     for (auto const& particle: particles ) {
         art::Ptr<simb::MCTruth> mctruth = fo.at(i_all);
         i_all++;
 
         if (mcOption == "nuOnly") {
             if ( !(mctruth->Origin() == 1 && particle->Mother() == 0) ) {
                 continue;
             }
         }
 
         // if (mctruth->Origin() == 1 || mc_mother[i] == 0) {
         //     cout << "process: " << particle->Process()
         //          << ", id: " << mc_id[i]
         //          << ", pdg: " <<  mc_pdg[i]
         //          << ", mother: " << mc_mother[i]
         //          << ", nDaughter: " << (particle->NumberDaughters())
         //          << ", truth: " << mctruth->Origin()
         //          << endl;
         // }
         // const art::Ptr<simb::MCTruth> mctruth = bt_serv->TrackIDToMCTruth(mc_id[i]);
 
         mc_process[i] = processMap[particle->Process()];
         if (mc_process[i] == 0) cout << "unknown process: " << particle->Process() << endl;
         mc_id[i] = particle->TrackId();
         mc_pdg[i] = particle->PdgCode();
         mc_mother[i] = particle->Mother();
         savedMCTrackIdMap[mc_id[i]] = mc_pdg[i];
 
         int Ndaughters = particle->NumberDaughters();
         vector<int> daughters;
         for (int i=0; i<Ndaughters; i++) {
             daughters.push_back(particle->Daughter(i));
         }
         mc_daughters.push_back(daughters);
         size_t numberTrajectoryPoints = particle->NumberTrajectoryPoints();
         int last = numberTrajectoryPoints - 1;
         const TLorentzVector& positionStart = particle->Position(0);
         const TLorentzVector& positionEnd   = particle->Position(last);
         const TLorentzVector& momentumStart = particle->Momentum(0);
         const TLorentzVector& momentumEnd   = particle->Momentum(last);
         positionStart.GetXYZT(mc_startXYZT[i]);
         positionEnd.GetXYZT(mc_endXYZT[i]);
         momentumStart.GetXYZT(mc_startMomentum[i]);
         momentumEnd.GetXYZT(mc_endMomentum[i]);
 
         if (fSaveMCTrackPoints) {
             TClonesArray *Lposition = new TClonesArray("TLorentzVector", numberTrajectoryPoints);
             // Read the position and momentum along this particle track
             for(unsigned int j=0; j<numberTrajectoryPoints; j++) {
                 new ((*Lposition)[j]) TLorentzVector(particle->Position(j));
             }
             fMC_trackPosition->Add(Lposition);
         }
 
         i++;
         if (i==MAX_TRACKS) {
             cout << "WARNING:: # tracks exceeded MAX_TRACKS " << MAX_TRACKS << endl;
             break;
         }
     } // particle loop done
     mc_Ntrack = i;
     // cout << "MC_Ntracks:" << mc_Ntrack << endl;
 
     // Generator Info
     art::Handle< std::vector<simb::MCTruth> > mctruthListHandle;
     event.getByLabel("generator",mctruthListHandle);
     std::vector<art::Ptr<simb::MCTruth> > mclist;
     art::fill_ptr_vector(mclist, mctruthListHandle);
     art::Ptr<simb::MCTruth> mctruth;
 
     if (mclist.size()>0) {
         mctruth = mclist.at(0);
         if (mctruth->NeutrinoSet()) {
             simb::MCNeutrino nu = mctruth->GetNeutrino();
             mc_isnu = 1;
             mc_nGeniePrimaries = mctruth->NParticles();
             mc_nu_pdg = nu.Nu().PdgCode();
             mc_nu_ccnc = nu.CCNC();
             mc_nu_mode = nu.Mode();
             mc_nu_intType = nu.InteractionType();
             mc_nu_target = nu.Target();
             mc_hitnuc = nu.HitNuc();
             mc_hitquark = nu.HitQuark();
             mc_nu_Q2 = nu.QSqr();
             mc_nu_W = nu.W();
             mc_nu_X = nu.X();
             mc_nu_Y = nu.Y();
             mc_nu_Pt = nu.Pt();
             mc_nu_Theta = nu.Theta();
 
             const TLorentzVector& position = nu.Nu().Position(0);
             const TLorentzVector& momentum = nu.Nu().Momentum(0);
             position.GetXYZT(mc_nu_pos);
             momentum.GetXYZT(mc_nu_mom);
             // cout << "nu: " << mc_nu_pdg << ", nPrim: " << mc_nGeniePrimaries
             //      << ", ccnc: " << mc_nu_ccnc << endl;
             // for (int i=0; i<mc_nGeniePrimaries; i++) {
             //     simb::MCParticle particle = mctruth->GetParticle(i);
             //     cout << "id: " << particle.TrackId()
             //          << ", pdg: " << particle.PdgCode()
             //          << endl;
             // }
         }
     }
 
 }

void wc::CellTree::processMCTracks ( )

Definition at line 891 of file CellTree_module.cc.

References IsPrimary(), mc_daughters, mc_id, mc_mother, mc_Ntrack, trackChildren, trackIndex, trackParents, and trackSiblings.

Referenced by analyze().

 {
     // map track id to track index in the array
     for (int i=0; i<mc_Ntrack; i++) {
         trackIndex[mc_id[i]] = i;
     }
 
     // in trackParents, trackChildren, trackSiblings vectors, store track index (not track id)
     for (int i=0; i<mc_Ntrack; i++) {
         // currently, parent size == 1;
         // for primary particle, parent id = 0;
         vector<int> parents;
         if ( !IsPrimary(i) ) {
             parents.push_back(trackIndex[mc_mother[i]]);
         }
         trackParents.push_back(parents); // primary track will have 0 parents
 
         vector<int> children;
         int nChildren = mc_daughters.at(i).size();
         for (int j=0; j<nChildren; j++) {
             children.push_back(trackIndex[mc_daughters.at(i).at(j)]);
         }
         trackChildren.push_back(children);
 
     }
 
     // siblings
     for (int i=0; i<mc_Ntrack; i++) {
         vector<int> siblings;
         if ( IsPrimary(i) ) {
             for (int j=0; j<mc_Ntrack; j++) {
                 if( IsPrimary(j) ) {
                     siblings.push_back(j);
                 }
             }
         }
         else {
             // siblings are simply children of the mother
             int mother = trackIndex[mc_mother[i]];
             int nSiblings = trackChildren.at(mother).size();
             for (int j=0; j<nSiblings; j++) {
                 siblings.push_back(trackChildren.at(mother).at(j));
             }
         }
         trackSiblings.push_back(siblings);
     }
 
 }

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 wc::CellTree::processOpFlash ( const art::Event & evt )

Definition at line 612 of file CellTree_module.cc.

References fGeometry, art::fill_ptr_vector(), fOpFlashLabel, art::DataViewImpl::getByLabel(), geo::GeometryCore::NOpDets(), of_multiplicity, of_nFlash, of_peTotal, of_t, and opMultPEThresh.

Referenced by analyze().

   {
     art::Handle<std::vector<recob::OpFlash> > flash_handle;
     if(! event.getByLabel(fOpFlashLabel, flash_handle)){
       cout << "WARNING: no label " << fOpFlashLabel << endl;
       return;
     }
     std::vector<art::Ptr<recob::OpFlash> > flashes;
     art::fill_ptr_vector(flashes, flash_handle);
     of_nFlash = (int)flashes.size();
 
     int a=0;
     int nOpDet = fGeometry->NOpDets();
 
     for(auto const& flash: flashes){
       of_t.push_back(flash->Time());
       of_peTotal.push_back(flash->TotalPE());
       TH1F *h = new ((*fPEperOpDet)[a]) TH1F("","",nOpDet,0,nOpDet);
 
       int mult = 0;
       for(int i=0; i<nOpDet; ++i){
         if(flash->PE(i) >= opMultPEThresh){
           mult++;
         }
         h->SetBinContent(i, flash->PE(i));
       }
       of_multiplicity.push_back(mult);
       a++;
     }
   }

void wc::CellTree::processOpHit ( const art::Event & evt )

Definition at line 591 of file CellTree_module.cc.

References art::fill_ptr_vector(), fOpHitLabel, art::DataViewImpl::getByLabel(), oh_bgtime, oh_channel, oh_nHits, oh_pe, and oh_trigtime.

Referenced by analyze().

 {
   art::Handle<std::vector<recob::OpHit> > ophit_handle;
   if(! event.getByLabel(fOpHitLabel, ophit_handle)){
     cout << "WARNING: no label " << fOpHitLabel << endl;
     return;
   }
   std::vector<art::Ptr<recob::OpHit> > ophits;
   art::fill_ptr_vector(ophits, ophit_handle);
   oh_nHits = (int)ophits.size();
 
   for(auto const& oh : ophits){
     oh_channel.push_back(oh->OpChannel());
     oh_bgtime.push_back(oh->PeakTime());
     oh_trigtime.push_back(oh->PeakTimeAbs());
     oh_pe.push_back(oh->PE());
   }
 
 }

void wc::CellTree::processRaw ( const art::Event & evt )

Definition at line 525 of file CellTree_module.cc.

References art::fill_ptr_vector(), fRaw_channelId, fRaw_nChannel, fRawDigitLabel, art::DataViewImpl::getByLabel(), nRawSamples, and raw::Uncompress().

Referenced by analyze().

 {
     art::Handle< std::vector<raw::RawDigit> > rawdigit;
     if (! event.getByLabel(fRawDigitLabel, rawdigit)) {
         cout << "WARNING: no label " << fRawDigitLabel << endl;
         return;
     }
     std::vector< art::Ptr<raw::RawDigit> >  wires;
     art::fill_ptr_vector(wires, rawdigit);
 
     fRaw_nChannel = wires.size();
 
     int i=0;
     for (auto const& wire: wires) {
         int chanId = wire->Channel();
         fRaw_channelId.push_back(chanId);
 
         int nSamples = wire->Samples();
         std::vector<short> uncompressed(nSamples);
         raw::Uncompress(wire->ADCs(), uncompressed, wire->Compression());
 
         TH1F *h = new((*fRaw_wf)[i]) TH1F("", "", nRawSamples, 0, nRawSamples);
         for (int j=1; j<=nSamples; j++) {
             h->SetBinContent(j, uncompressed[j-1]);
         }
         i++;
         if (i==1) {
           cout << nSamples << " samples expanding to " << nRawSamples << endl;
         }
     }
 
 }

void wc::CellTree::processSimChannel ( const art::Event & evt )

Definition at line 644 of file CellTree_module.cc.

References fSIMIDE_channelIdY, fSIMIDE_numElectrons, fSIMIDE_size, fSIMIDE_tdc, fSIMIDE_trackId, fSIMIDE_x, fSIMIDE_y, and fSIMIDE_z.

Referenced by analyze().

 {
     art::Handle< std::vector<sim::SimChannel> > simChannelHandle;
     event.getByLabel("largeant", simChannelHandle);
     // cout << "total simChannel: " << (*simChannelHandle).size() << endl;
     fSIMIDE_size = 0;
     for ( auto const& channel : (*simChannelHandle) ) {
         auto channelNumber = channel.Channel();
         // cout << channelNumber << endl;
         // if (! (fGeometry->SignalType( channelNumber ) == geo::kCollection) ) {
         //     continue;
         // }
         auto const& timeSlices = channel.TDCIDEMap();
         for ( auto const& timeSlice : timeSlices ) {
             auto const& energyDeposits = timeSlice.second;
             for ( auto const& energyDeposit : energyDeposits ) {
                 fSIMIDE_size++;
                 fSIMIDE_channelIdY.push_back(channelNumber);
                 fSIMIDE_tdc.push_back(timeSlice.first);
                 fSIMIDE_trackId.push_back(energyDeposit.trackID);
                 fSIMIDE_x.push_back(energyDeposit.x);
                 fSIMIDE_y.push_back(energyDeposit.y);
                 fSIMIDE_z.push_back(energyDeposit.z);
                 fSIMIDE_numElectrons.push_back(energyDeposit.numElectrons);
                 // cout << channelNumber << ": " << energyDeposit.trackID << ": " << timeSlice.first << ": "
                 //      << energyDeposit.x << ", " << energyDeposit.y << ", " << energyDeposit.z << ", "
                 //      << energyDeposit.numElectrons << endl;
             }
         }
 
     }
     cout << "total IDEs: " << fSIMIDE_size << endl;
 
 
 }

void wc::CellTree::processSpacePoint	(	const art::Event &	event,
		TString	option,
		ostream &	out = `cout`
	)

Definition at line 805 of file CellTree_module.cc.

References geo::GeometryCore::DetectorName(), fEvent, fGeometry, art::fill_ptr_vector(), fRun, fSubRun, print_vector(), x, y, and z.

Referenced by analyze().

 {
 
     art::Handle< std::vector<recob::SpacePoint> > sp_handle;
     art::Handle< std::vector<recob::PointCharge> > pc_handle;
     bool sp_exists = event.getByLabel(option.Data(), sp_handle);
     bool pc_exists = event.getByLabel(option.Data(), pc_handle);
     if (! sp_exists) {
         cout << "WARNING: no label " << option << endl;
         return;
     }
     std::vector< art::Ptr<recob::SpacePoint> >  sps;
     std::vector< art::Ptr<recob::PointCharge> >  pcs;
     art::fill_ptr_vector(sps, sp_handle);
     if (pc_exists) {
         art::fill_ptr_vector(pcs, pc_handle);
         if (sps.size() != pcs.size()) {
             cout << "WARNING: SpacePoint and PointCharge length mismatch" << endl;
             return;
         }
     }
     double x=0, y=0, z=0, q=0, nq=1;
     vector<double> vx, vy, vz, vq, vnq;
 
     for (uint i=0; i < sps.size(); i++ ) {
         // cout << sp->XYZ()[0] << ", " << sp->XYZ()[1] << ", " << sp->XYZ()[2] << endl;
         x = sps[i]->XYZ()[0];
         y = sps[i]->XYZ()[1];
         z = sps[i]->XYZ()[2];
         if (pc_exists && pcs[i]->hasCharge()) {
           q = pcs[i]->charge();
         } else {
           q = 0;
         }
         vx.push_back(x);
         vy.push_back(y);
         vz.push_back(z);
         vq.push_back(q);
         vnq.push_back(nq);
     }
 
     out << fixed << setprecision(1);
     out << "{" << endl;
 
     out << '"' << "runNo" << '"' << ":" << '"' << fRun << '"' << "," << endl;
     out << '"' << "subRunNo" << '"' << ":" << '"' << fSubRun << '"' << "," << endl;
     out << '"' << "eventNo" << '"' << ":" << '"' << fEvent << '"' << "," << endl;
 
     TString geomName(fGeometry->DetectorName().c_str());
     if (geomName.Contains("35t")) { geomName = "dune35t"; }
     else if (geomName.Contains("protodune")) { geomName = "protodune"; }
     else if (geomName.Contains("workspace")) { geomName = "dune10kt_workspace"; }
     else { geomName = "uboone"; } // use uboone as default
     out << '"' << "geom" << '"' << ":" << '"' << geomName << '"' << "," << endl;
 
 
     print_vector(out, vx, "x");
     print_vector(out, vy, "y");
     print_vector(out, vz, "z");
 
     out << fixed << setprecision(0);
     print_vector(out, vq, "q");
     print_vector(out, vnq, "nq");
 
     out << '"' << "type" << '"' << ":" << '"' << option << '"' << endl;
     out << "}" << endl;
 }

void wc::CellTree::processTrigger ( const art::Event & evt )

Definition at line 941 of file CellTree_module.cc.

References fBeamgatetime, art::fill_ptr_vector(), fTriggerbits, fTriggerLabel, fTriggernumber, fTriggertime, and art::DataViewImpl::getByLabel().

Referenced by analyze().

 {
   art::Handle< std::vector<raw::Trigger>> triggerListHandle;
   std::vector<art::Ptr<raw::Trigger>> triggerlist;
   if (event.getByLabel(fTriggerLabel, triggerListHandle)) {
     art::fill_ptr_vector(triggerlist, triggerListHandle);
   }
   else {
     cout << "WARNING: no label " << fTriggerLabel << endl;
   }
   if (triggerlist.size()){
     fTriggernumber = triggerlist[0]->TriggerNumber();
     fTriggertime   = triggerlist[0]->TriggerTime();
     fBeamgatetime  = triggerlist[0]->BeamGateTime();
     fTriggerbits   = triggerlist[0]->TriggerBits();
   }
   else {
     fTriggernumber = 0;
     fTriggertime   = 0;
     fBeamgatetime  = 0;
     fTriggerbits   = 0;
   }
 }

void wc::CellTree::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 249 of file CellTree_module.cc.

References fCalibLabel, fOpFlashLabel, fOpHitLabel, fOutFileName, fRawDigitLabel, fSaveCalib, fSaveJSON, fSaveMC, fSaveMCTrackPoints, fSaveOpFlash, fSaveOpHit, fSaveRaw, fSaveSimChannel, fSaveTrigger, fSpacePointLabels, fTriggerLabel, fhicl::ParameterSet::get(), mcOption, nRawSamples, and opMultPEThresh.

Referenced by CellTree().

                                                     {
     fRawDigitLabel   = p.get<std::string>("RawDigitLabel");
     fCalibLabel      = p.get<std::string>("CalibLabel");
     fOpHitLabel      = p.get<std::string>("OpHitLabel");
     fOpFlashLabel    = p.get<std::string>("OpFlashLabel");
     fTriggerLabel    = p.get<std::string>("TriggerLabel");
     fSpacePointLabels= p.get<std::vector<std::string> >("SpacePointLabels");
     fOutFileName     = p.get<std::string>("outFile");
     mcOption        = p.get<std::string>("mcOption");
     fSaveMCTrackPoints = p.get<bool>("saveMCTrackPoints");
     fSaveRaw         = p.get<bool>("saveRaw");
     fSaveCalib       = p.get<bool>("saveCalib");
     fSaveOpHit       = p.get<bool>("saveOpHit");
     fSaveOpFlash     = p.get<bool>("saveOpFlash");
     fSaveMC          = p.get<bool>("saveMC");
     fSaveSimChannel  = p.get<bool>("saveSimChannel");
     fSaveTrigger     = p.get<bool>("saveTrigger");
     fSaveJSON        = p.get<bool>("saveJSON");
     opMultPEThresh   = p.get<float>("opMultPEThresh");
     nRawSamples      = p.get<int>("nRawSamples");
 }

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

inlineinherited

Definition at line 33 of file EventObserverBase.h.

36 {}

void wc::CellTree::reset ( void )

Definition at line 454 of file CellTree_module.cc.

References fCalib_channelId, fCalib_wf, fMC_trackPosition, fPEperOpDet, fRaw_channelId, fRaw_wf, fSIMIDE_channelIdY, fSIMIDE_numElectrons, fSIMIDE_tdc, fSIMIDE_trackId, fSIMIDE_x, fSIMIDE_y, fSIMIDE_z, MAX_TRACKS, mc_daughters, mc_endMomentum, mc_endXYZT, mc_hitnuc, mc_hitquark, mc_id, mc_isnu, mc_mother, mc_nGeniePrimaries, mc_Ntrack, mc_nu_ccnc, mc_nu_intType, mc_nu_mode, mc_nu_mom, mc_nu_pdg, mc_nu_pos, mc_nu_Pt, mc_nu_Q2, mc_nu_target, mc_nu_Theta, mc_nu_W, mc_nu_X, mc_nu_Y, mc_pdg, mc_startMomentum, mc_startXYZT, of_multiplicity, of_peTotal, of_t, oh_bgtime, oh_channel, oh_pe, oh_trigtime, savedMCTrackIdMap, trackChildren, trackIndex, trackParents, and trackSiblings.

Referenced by analyze().

 {
 
     fRaw_channelId.clear();
     // fRaw_wf->Clear();
     fRaw_wf->Delete();
 
     fCalib_channelId.clear();
     fCalib_wf->Clear();
 
     oh_channel.clear();
     oh_bgtime.clear();
     oh_trigtime.clear();
     oh_pe.clear();
 
     of_t.clear();
     of_peTotal.clear();
     of_multiplicity.clear();
     fPEperOpDet->Delete();
 
     fSIMIDE_channelIdY.clear();
     fSIMIDE_trackId.clear();
     fSIMIDE_tdc.clear();
     fSIMIDE_x.clear();
     fSIMIDE_y.clear();
     fSIMIDE_z.clear();
     fSIMIDE_numElectrons.clear();
 
     mc_Ntrack = 0;
     for (int i=0; i<MAX_TRACKS; i++) {
         mc_id[i] = 0;
         mc_pdg[i] = 0;
         mc_mother[i] = 0;
         for (int j=0; j<4; j++) {
             mc_startXYZT[i][j]      = 0;
             mc_endXYZT[i][j]        = 0;
             mc_startMomentum[i][j] = 0;
             mc_endMomentum[i][j]   = 0;
         }
     }
     mc_daughters.clear();
     savedMCTrackIdMap.clear();
     fMC_trackPosition->Clear();
 
     mc_isnu = 0;
     mc_nGeniePrimaries = -1;
     mc_nu_pdg = -1;
     mc_nu_ccnc = -1;
     mc_nu_mode = -1;
     mc_nu_intType = -1;
     mc_nu_target = -1;
     mc_hitnuc = -1;
     mc_hitquark = -1;
     mc_nu_Q2 = -1;
     mc_nu_W = -1;
     mc_nu_X = -1;
     mc_nu_Y = -1;
     mc_nu_Pt = -1;
     mc_nu_Theta = -1;
     for (int i=0; i<4; i++) {
         mc_nu_pos[i] = 0;
         mc_nu_mom[i] = 0;
     }
 
     trackIndex.clear();
     trackParents.clear();
     trackChildren.clear();
     trackSiblings.clear();
 }

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

inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

Referenced by art::RootOutputFile::writeOne().

   {
     return selector_config_id_;
   }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

 {
   if (!moduleContext_)
     return;
 
   // If none of the branches have missing consumes statements, exit early.
   if (std::all_of(cbegin(missingConsumes_),
                   cend(missingConsumes_),
                   [](auto const& perBranch) { return perBranch.empty(); }))
     return;
 
   constexpr cet::HorizontalRule rule{60};
   mf::LogPrint log{"MTdiagnostics"};
   log << '\n'
       << rule('=') << '\n'
       << "The following consumes (or mayConsume) statements are missing from\n"
       << module_context(moduleDescription_) << '\n'
       << rule('-') << '\n';
 
   cet::for_all_with_index(
     missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
       for (auto const& pi : perBranch) {
         log << "  "
             << assemble_consumes_statement(static_cast<BranchType>(i), pi)
             << '\n';
       }
     });
   log << rule('=');
 }

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

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

 {
   // Early exits if consumes tracking has been disabled or if the
   // consumed product is an allowed consumable.
   if (!moduleContext_)
     return;
 
   if (cet::binary_search_all(consumables_[bt], pi))
     return;
 
   if (requireConsumes_) {
     throw Exception(errors::ProductRegistrationFailure,
                     "Consumer: an error occurred during validation of a "
                     "retrieved product\n\n")
       << "The following consumes (or mayConsume) statement is missing from\n"
       << module_context(moduleDescription_) << ":\n\n"
       << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
   }
 
   missingConsumes_[bt].insert(pi);
 }

bool art::EventObserverBase::wantAllEvents ( ) const

inlineinherited

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), art::RootOutput::RootOutput(), and art::OutputWorker::wantAllEvents().

   {
     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().

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), and art::OutputModule::doWriteEvent().

   {
     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

TDatabasePDG* wc::CellTree::dbPDG

private

Definition at line 226 of file CellTree_module.cc.

Referenced by CellTree(), and PDGName().

int wc::CellTree::entryNo

private

Definition at line 145 of file CellTree_module.cc.

Referenced by analyze(), and CellTree().

double wc::CellTree::fBeamgatetime

private

Definition at line 155 of file CellTree_module.cc.

Referenced by initOutput(), and processTrigger().

std::vector<int> wc::CellTree::fCalib_channelId

private

Definition at line 161 of file CellTree_module.cc.

Referenced by initOutput(), processCalib(), and reset().

int wc::CellTree::fCalib_nChannel

private

Definition at line 158 of file CellTree_module.cc.

Referenced by initOutput(), and processCalib().

TClonesArray* wc::CellTree::fCalib_wf

private

Definition at line 163 of file CellTree_module.cc.

Referenced by initOutput(), and reset().

std::string wc::CellTree::fCalibLabel

private

Definition at line 117 of file CellTree_module.cc.

Referenced by processCalib(), and reconfigure().

int wc::CellTree::fEvent

private

Definition at line 148 of file CellTree_module.cc.

Referenced by analyze(), initOutput(), printEvent(), and processSpacePoint().

double wc::CellTree::fEventTime

private

Definition at line 151 of file CellTree_module.cc.

Referenced by analyze(), and initOutput().

TTree* wc::CellTree::fEventTree

private

Definition at line 141 of file CellTree_module.cc.

Referenced by analyze(), endJob(), and initOutput().

art::ServiceHandle<geo::Geometry> wc::CellTree::fGeometry

private

Definition at line 135 of file CellTree_module.cc.

Referenced by processOpFlash(), and processSpacePoint().

TObjArray* wc::CellTree::fMC_trackPosition

private

Definition at line 201 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

std::string wc::CellTree::fOpFlashLabel

private

Definition at line 119 of file CellTree_module.cc.

Referenced by processOpFlash(), and reconfigure().

std::string wc::CellTree::fOpHitLabel

private

Definition at line 118 of file CellTree_module.cc.

Referenced by processOpHit(), and reconfigure().

TFile* wc::CellTree::fOutFile

private

Definition at line 140 of file CellTree_module.cc.

Referenced by endJob(), and initOutput().

std::string wc::CellTree::fOutFileName

private

Definition at line 122 of file CellTree_module.cc.

Referenced by initOutput(), and reconfigure().

TClonesArray* wc::CellTree::fPEperOpDet

private

Definition at line 176 of file CellTree_module.cc.

Referenced by initOutput(), and reset().

std::vector<int> wc::CellTree::fRaw_channelId

private

Definition at line 179 of file CellTree_module.cc.

Referenced by initOutput(), processRaw(), and reset().

int wc::CellTree::fRaw_nChannel

private

Definition at line 178 of file CellTree_module.cc.

Referenced by initOutput(), and processRaw().

TClonesArray* wc::CellTree::fRaw_wf

private

Definition at line 180 of file CellTree_module.cc.

Referenced by initOutput(), and reset().

std::string wc::CellTree::fRawDigitLabel

private

Definition at line 116 of file CellTree_module.cc.

Referenced by processRaw(), and reconfigure().

int wc::CellTree::fRun

private

Definition at line 149 of file CellTree_module.cc.

Referenced by analyze(), initOutput(), printEvent(), and processSpacePoint().

bool wc::CellTree::fSaveCalib

private

Definition at line 129 of file CellTree_module.cc.

Referenced by analyze(), and reconfigure().

bool wc::CellTree::fSaveJSON

private

Definition at line 134 of file CellTree_module.cc.

Referenced by analyze(), endJob(), initOutput(), and reconfigure().

bool wc::CellTree::fSaveMC

private

Definition at line 132 of file CellTree_module.cc.

Referenced by analyze(), and reconfigure().

bool wc::CellTree::fSaveMCTrackPoints

private

Definition at line 126 of file CellTree_module.cc.

Referenced by processMC(), and reconfigure().

bool wc::CellTree::fSaveOpFlash

private

Definition at line 131 of file CellTree_module.cc.

Referenced by analyze(), and reconfigure().

bool wc::CellTree::fSaveOpHit

private

Definition at line 130 of file CellTree_module.cc.

Referenced by analyze(), and reconfigure().

bool wc::CellTree::fSaveRaw

private

Definition at line 128 of file CellTree_module.cc.

Referenced by analyze(), and reconfigure().

bool wc::CellTree::fSaveSimChannel

private

Definition at line 127 of file CellTree_module.cc.

Referenced by analyze(), and reconfigure().

bool wc::CellTree::fSaveTrigger

private

Definition at line 133 of file CellTree_module.cc.

Referenced by analyze(), and reconfigure().

vector<int> wc::CellTree::fSIMIDE_channelIdY

private

Definition at line 183 of file CellTree_module.cc.

Referenced by initOutput(), processSimChannel(), and reset().

vector<float> wc::CellTree::fSIMIDE_numElectrons

private

Definition at line 189 of file CellTree_module.cc.

Referenced by initOutput(), processSimChannel(), and reset().

int wc::CellTree::fSIMIDE_size

private

Definition at line 182 of file CellTree_module.cc.

Referenced by initOutput(), and processSimChannel().

vector<unsigned short> wc::CellTree::fSIMIDE_tdc

private

Definition at line 185 of file CellTree_module.cc.

Referenced by initOutput(), processSimChannel(), and reset().

vector<int> wc::CellTree::fSIMIDE_trackId

private

Definition at line 184 of file CellTree_module.cc.

Referenced by initOutput(), processSimChannel(), and reset().

vector<float> wc::CellTree::fSIMIDE_x

private

Definition at line 186 of file CellTree_module.cc.

Referenced by initOutput(), processSimChannel(), and reset().

vector<float> wc::CellTree::fSIMIDE_y

private

Definition at line 187 of file CellTree_module.cc.

Referenced by initOutput(), processSimChannel(), and reset().

vector<float> wc::CellTree::fSIMIDE_z

private

Definition at line 188 of file CellTree_module.cc.

Referenced by initOutput(), processSimChannel(), and reset().

std::vector<std::string> wc::CellTree::fSpacePointLabels

private

Definition at line 121 of file CellTree_module.cc.

Referenced by analyze(), and reconfigure().

int wc::CellTree::fSubRun

private

Definition at line 150 of file CellTree_module.cc.

Referenced by analyze(), initOutput(), printEvent(), and processSpacePoint().

unsigned int wc::CellTree::fTriggerbits

private

Definition at line 156 of file CellTree_module.cc.

Referenced by initOutput(), and processTrigger().

std::string wc::CellTree::fTriggerLabel

private

Definition at line 120 of file CellTree_module.cc.

Referenced by processTrigger(), and reconfigure().

unsigned int wc::CellTree::fTriggernumber

private

Definition at line 153 of file CellTree_module.cc.

Referenced by initOutput(), and processTrigger().

double wc::CellTree::fTriggertime

private

Definition at line 154 of file CellTree_module.cc.

Referenced by initOutput(), and processTrigger().

std::vector<std::vector<int> > wc::CellTree::mc_daughters

private

Definition at line 200 of file CellTree_module.cc.

Referenced by DumpMCJSON(), initOutput(), printEvent(), processMC(), processMCTracks(), and reset().

float wc::CellTree::mc_endMomentum[MAX_TRACKS][4]

private

Definition at line 199 of file CellTree_module.cc.

Referenced by initOutput(), printEvent(), processMC(), and reset().

float wc::CellTree::mc_endXYZT[MAX_TRACKS][4]

private

Definition at line 197 of file CellTree_module.cc.

Referenced by DumpMCJSON(), initOutput(), printEvent(), processMC(), and reset().

int wc::CellTree::mc_hitnuc

private

Definition at line 210 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_hitquark

private

Definition at line 211 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_id[MAX_TRACKS]

private

Definition at line 192 of file CellTree_module.cc.

Referenced by DumpMCJSON(), initOutput(), printEvent(), processMC(), processMCTracks(), and reset().

int wc::CellTree::mc_isnu

private

Definition at line 203 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_mother[MAX_TRACKS]

private

Definition at line 195 of file CellTree_module.cc.

Referenced by initOutput(), printEvent(), processMC(), processMCTracks(), and reset().

int wc::CellTree::mc_nGeniePrimaries

private

Definition at line 204 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_Ntrack

private

Definition at line 191 of file CellTree_module.cc.

Referenced by DumpMCJSON(), initOutput(), printEvent(), processMC(), processMCTracks(), and reset().

int wc::CellTree::mc_nu_ccnc

private

Definition at line 206 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_nu_intType

private

Definition at line 208 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_nu_mode

private

Definition at line 207 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

float wc::CellTree::mc_nu_mom[4]

private

Definition at line 219 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_nu_pdg

private

Definition at line 205 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

float wc::CellTree::mc_nu_pos[4]

private

Definition at line 218 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

double wc::CellTree::mc_nu_Pt

private

Definition at line 216 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

double wc::CellTree::mc_nu_Q2

private

Definition at line 212 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_nu_target

private

Definition at line 209 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

double wc::CellTree::mc_nu_Theta

private

Definition at line 217 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

double wc::CellTree::mc_nu_W

private

Definition at line 213 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

double wc::CellTree::mc_nu_X

private

Definition at line 214 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

double wc::CellTree::mc_nu_Y

private

Definition at line 215 of file CellTree_module.cc.

Referenced by initOutput(), processMC(), and reset().

int wc::CellTree::mc_pdg[MAX_TRACKS]

private

Definition at line 193 of file CellTree_module.cc.

Referenced by DumpMCJSON(), initOutput(), KeepMC(), printEvent(), processMC(), and reset().

int wc::CellTree::mc_process[MAX_TRACKS]

private

Definition at line 194 of file CellTree_module.cc.

Referenced by initOutput(), and processMC().

float wc::CellTree::mc_startMomentum[MAX_TRACKS][4]

private

Definition at line 198 of file CellTree_module.cc.

Referenced by DumpMCJSON(), initOutput(), KeepMC(), printEvent(), processMC(), and reset().

float wc::CellTree::mc_startXYZT[MAX_TRACKS][4]

private

Definition at line 196 of file CellTree_module.cc.

Referenced by DumpMCJSON(), initOutput(), printEvent(), processMC(), and reset().

std::string wc::CellTree::mcOption

private

Definition at line 123 of file CellTree_module.cc.

Referenced by processMC(), and reconfigure().

int wc::CellTree::nRawSamples

private

Definition at line 124 of file CellTree_module.cc.

Referenced by processCalib(), processRaw(), and reconfigure().

vector<int> wc::CellTree::of_multiplicity

private

Definition at line 175 of file CellTree_module.cc.

Referenced by initOutput(), processOpFlash(), and reset().

int wc::CellTree::of_nFlash

private

Definition at line 172 of file CellTree_module.cc.

Referenced by initOutput(), and processOpFlash().

vector<float> wc::CellTree::of_peTotal

private

Definition at line 174 of file CellTree_module.cc.

Referenced by initOutput(), processOpFlash(), and reset().

vector<float> wc::CellTree::of_t

private

Definition at line 173 of file CellTree_module.cc.

Referenced by initOutput(), processOpFlash(), and reset().

vector<double> wc::CellTree::oh_bgtime

private

Definition at line 168 of file CellTree_module.cc.

Referenced by initOutput(), processOpHit(), and reset().

vector<int> wc::CellTree::oh_channel

private

Definition at line 167 of file CellTree_module.cc.

Referenced by initOutput(), processOpHit(), and reset().

int wc::CellTree::oh_nHits

private

Definition at line 166 of file CellTree_module.cc.

Referenced by initOutput(), and processOpHit().

vector<double> wc::CellTree::oh_pe

private

Definition at line 170 of file CellTree_module.cc.

Referenced by initOutput(), processOpHit(), and reset().

vector<double> wc::CellTree::oh_trigtime

private

Definition at line 169 of file CellTree_module.cc.

Referenced by initOutput(), processOpHit(), and reset().

float wc::CellTree::opMultPEThresh

private

Definition at line 125 of file CellTree_module.cc.

Referenced by processOpFlash(), and reconfigure().

std::map<std::string, int> wc::CellTree::processMap

private

Definition at line 142 of file CellTree_module.cc.

Referenced by InitProcessMap(), and processMC().

std::map<int, int> wc::CellTree::savedMCTrackIdMap

private

Definition at line 143 of file CellTree_module.cc.

Referenced by processMC(), and reset().

std::vector<std::vector<int> > wc::CellTree::trackChildren

private

Definition at line 224 of file CellTree_module.cc.

Referenced by processMCTracks(), and reset().

std::map<int, int> wc::CellTree::trackIndex

private

Definition at line 222 of file CellTree_module.cc.

Referenced by DumpMCJSON(), processMCTracks(), and reset().

std::vector<std::vector<int> > wc::CellTree::trackParents

private

Definition at line 223 of file CellTree_module.cc.

Referenced by processMCTracks(), and reset().

std::vector<std::vector<int> > wc::CellTree::trackSiblings

private

Definition at line 225 of file CellTree_module.cc.

Referenced by processMCTracks(), and reset().

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

/cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/WireCell/CellTree_module.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation