larg4::LArG4Ana Class Reference

Inheritance diagram for larg4::LArG4Ana:

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

Public Member Functions
	LArG4Ana (fhicl::ParameterSet const &pset)
	Standard constructor and destructor for an FMWK module. More...
virtual	~LArG4Ana ()
void	analyze (const art::Event &evt)
void	beginJob ()
void	reconfigure (fhicl::ParameterSet const &pset)
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	fG4ModuleLabel
	module label for the Geant More...
std::string	fTruthModuleLabel
	module label for the Geant More...
TH1D *	fPDGCodes
TH1D *	fPi0Momentum
TH1D *	fnEnergy
TH1D *	fnDist
TH1D *	fnumChannels
	The number of channels recieving charge per event. More...
TProfile *	fnumIDEs
	Number of drift electrons per channel. More...
TH1D *	fEventCharge
	Charge collected per event. More...
TH1D *	fEventEnergy
	Energy collected per event. More...
TProfile *	fChannelCharge
	Charge per channel. More...
TProfile *	fChannelEnergy
	Energy per channel. More...
TTree *	fTree
Int_t	fTEvt
Int_t	fTSub
Int_t	fTRun
Int_t	fTPdg
Int_t	fTID
Int_t	fTNds
Int_t	fTNdsOriginal
Int_t	fTNds4
Int_t *	fTDID
Int_t *	fTDPdg
Float_t *	fTDWt
Char_t	fTProcess [35]
Char_t	fTVolume [35]
Char_t	fTTVolume [35]
Char_t	fTMaterial [35]
Char_t	fTDProcess [200][35]
Int_t	fTParentID
Int_t	fTStatus
Float_t	fTWeight
Float_t *	fT4Origin
Float_t *	fT4DOrigin
Float_t *	fT4Termination
Float_t *	fT4Momentum
Float_t *	fT4DMomentum

Detailed Description

Definition at line 56 of file LArG4Ana_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

larg4::LArG4Ana::LArG4Ana ( fhicl::ParameterSet const &  pset )

explicit

Standard constructor and destructor for an FMWK module.

Definition at line 119 of file LArG4Ana_module.cc.

References reconfigure().

    : EDAnalyzer(pset)
  {
    this->reconfigure(pset);
  }

larg4::LArG4Ana::~LArG4Ana ( )

virtual

Definition at line 127 of file LArG4Ana_module.cc.

  {
  }

Member Function Documentation

void larg4::LArG4Ana::analyze

(

const art::Event &

evt

)

Definition at line 218 of file LArG4Ana_module.cc.

References d, DEFINE_ART_MODULE, E, sim::ParticleList::end(), art::EventID::event(), fChannelCharge, fChannelEnergy, fEventCharge, fEventEnergy, fG4ModuleLabel, sim::ParticleList::find(), fnDist, fnEnergy, fnumChannels, fnumIDEs, fPDGCodes, fPi0Momentum, fT4DMomentum, fT4DOrigin, fT4Momentum, fT4Origin, fT4Termination, fTDID, fTDPdg, fTDProcess, fTDWt, fTEvt, fTID, fTMaterial, fTNds, fTNdsOriginal, fTParentID, fTPdg, fTProcess, fTree, fTRun, fTStatus, fTSub, fTTVolume, fTVolume, fTWeight, art::DataViewImpl::getView(), art::Event::id(), geo::GeometryCore::MaterialName(), sim::NoParticleId, cheat::ParticleInventoryService::ParticleList(), art::Event::run(), s, sim::ParticleList::size(), art::Event::subRun(), and geo::GeometryCore::VolumeName().

  {

    //get the list of particles from this event
    art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
    const sim::ParticleList& plist = pi_serv->ParticleList();
    art::ServiceHandle<geo::Geometry> geom;

    // loop over all sim::SimChannels in the event and make sure there are no
    // sim::IDEs with trackID values that are not in the sim::ParticleList
    std::vector<const sim::SimChannel*> sccol;
    evt.getView(fG4ModuleLabel, sccol);

    double totalCharge=0.0;
    double totalEnergy=0.0;
    fnumChannels->Fill(sccol.size());
    for(size_t sc = 0; sc < sccol.size(); ++sc){
      double numIDEs=0.0;
      double scCharge=0.0;
      double scEnergy=0.0;
      const auto & tdcidemap = sccol[sc]->TDCIDEMap();
      for(auto mapitr = tdcidemap.begin(); mapitr != tdcidemap.end(); mapitr++){
        const std::vector<sim::IDE> idevec = (*mapitr).second;
        numIDEs += idevec.size();
        for(size_t iv = 0; iv < idevec.size(); ++iv){
          if(plist.find( idevec[iv].trackID ) == plist.end()
             && idevec[iv].trackID != sim::NoParticleId) 
          mf::LogWarning("LArG4Ana") << idevec[iv].trackID << " is not in particle list"; 
          totalCharge +=idevec[iv].numElectrons;
          scCharge += idevec[iv].numElectrons;
          totalEnergy +=idevec[iv].energy;
          scEnergy += idevec[iv].energy;
        }
      }
      fnumIDEs->Fill(sc,numIDEs);
      fChannelCharge->Fill(sc,scCharge);
      fChannelEnergy->Fill(sc,scEnergy);
    }
    fEventCharge->Fill(totalCharge);
    fEventEnergy->Fill(totalEnergy);

    // get the particles from the back tracker
    const sim::ParticleList& Particles = pi_serv->ParticleList();
    std::vector<const simb::MCParticle*> pvec;
    pvec.reserve(Particles.size());
    for (const auto& PartPair: Particles) {
      pvec.push_back(PartPair.second);
      fPDGCodes->Fill(PartPair.second->PdgCode());
    }

    // now look for pi0's that decay to 2 gammas
    int pi0loc = -1;
    int numpi0gamma = 0;
    for(unsigned int i = 0; i < pvec.size(); ++i){
      if(pvec[i]->PdgCode() == 111) pi0loc = i;
      if(pvec[i]->Mother() == pi0loc+1 && 
         pi0loc > 0 && 
         pvec[i]->PdgCode() == 22){
        mf::LogInfo("LArG4Ana") << pvec[i]->E() << " gamma energy ";
        ++numpi0gamma;
      }
      
      // n,Lambda,K0s,K0L,K0
      if (pvec[i]->PdgCode() == 2112 || 
          pvec[i]->PdgCode() == 3122 || 
          pvec[i]->PdgCode() == 130  || 
          pvec[i]->PdgCode() == 310  || 
          pvec[i]->PdgCode() == 311 ) {
        fnEnergy->Fill(pvec[i]->E(),pvec[i]->Weight());
        fnDist->Fill(pvec[i]->Vx(),pvec[i]->Weight());
      }

      fTPdg = pvec[i]->PdgCode();
      fTID = pvec[i]->TrackId();
      // 0 out strings, else there may be cruft in here from prev evt.
      for (unsigned int s = 0; s < 35; ++s){
          *(fTProcess+s)  = 0; 
          *(fTProcess+s)  = 0;
          *(fTMaterial+s) = 0; 
          *(fTMaterial+s) = 0;
          *(fTVolume+s)   = 0;
          *(fTVolume+s)   = 0;
          *(fTTVolume+s)  = 0; 
          *(fTTVolume+s)  = 0;
      }
      
      for(unsigned int s = 0; s < pvec[i]->Process().length(); ++s) *(fTProcess+s) = pvec[i]->Process()[s];

      TVector3 dum = pvec[i]->Position().Vect();

      for (unsigned int s = 0; s < geom->MaterialName(pvec[i]->Position().Vect()).length(); ++s) 
        *(fTMaterial+s) = geom->MaterialName(pvec[i]->Position().Vect())[s];

      for (unsigned int s = 0; s < geom->VolumeName(pvec[i]->Position().Vect()).length(); ++s) 
        *(fTVolume+s) = geom->VolumeName(pvec[i]->Position().Vect())[s];

      for (unsigned int s = 0; s < geom->VolumeName(pvec[i]->EndPosition().Vect()).length(); ++s) 
        *(fTTVolume+s) = geom->VolumeName(pvec[i]->EndPosition().Vect())[s];

      fTEvt = evt.id().event(); 
      fTSub = evt.subRun();
      fTRun = evt.run();
      fTParentID = pvec[i]->Mother();
      fTStatus = pvec[i]->StatusCode();
      int daughter = 9999;
      fTNds = TMath::Min(pvec[i]->NumberDaughters(),fTNdsOriginal);
      for( int d = 0; d < fTNds; d++ ){
        daughter = pvec[i]->Daughter(d);
        fTDID[d] = daughter; 
        // zero it out.
        for (unsigned int s = 0; s < 35; ++s) *(fTDProcess[d]+s) = 0; 
        
        for(unsigned int jj = i; jj < pvec.size(); ++jj){ // Don't look below i.
       
          if (fTDID[d] == pvec[jj]->TrackId()){
            fTDPdg[d] = pvec[jj]->PdgCode(); // get the pointer,  
            fTDWt[d]  = pvec[jj]->Weight();
            
            for (unsigned int s = 0; s < pvec[jj]->Process().length(); ++s) 
              *(fTDProcess[d]+s) = pvec[jj]->Process()[s];
                
            for (unsigned int kk = 0; kk < 4; ++kk){
              fT4DOrigin[d*4+kk] = pvec[jj]->Position()[kk];
              fT4DMomentum[d*4+kk] = pvec[jj]->Momentum()[kk];
            }
            break;
          }           
        }
      }//end loop over d
      
      for (unsigned int ii = 0; ii < 4; ++ii){
        fT4Termination[ii] = 1e9;
        fT4Origin[ii] = pvec[i]->Position()[ii];
        if (ii!=3) fT4Termination[ii] = pvec[i]->EndPosition()[ii];
        if (ii==4) fT4Termination[ii] = pvec[i]->Momentum()[ii]; // yes, odd
        fT4Momentum[ii] = pvec[i]->Momentum()[ii];
      }

      fTWeight = pvec[i]->Weight();
      fTree->Fill();
      
    } // end loop on particles in list 
    if(numpi0gamma == 2 && pi0loc > 0){
      mf::LogInfo("LArG4Ana") << pvec[pi0loc]->E();
      fPi0Momentum->Fill(pvec[pi0loc]->E());
    }

    return;
  }

void larg4::LArG4Ana::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 132 of file LArG4Ana_module.cc.

References fChannelCharge, fChannelEnergy, fEventCharge, fEventEnergy, fnDist, fnEnergy, fnumChannels, fnumIDEs, fPDGCodes, fPi0Momentum, fT4DMomentum, fT4DOrigin, fT4Momentum, fT4Origin, fT4Termination, fTDID, fTDPdg, fTDProcess, fTDWt, fTEvt, fTID, fTMaterial, fTNds, fTNds4, fTParentID, fTPdg, fTProcess, fTree, fTRun, fTStatus, fTSub, fTTVolume, fTVolume, fTWeight, art::TFileDirectory::make(), and geo::GeometryCore::Nchannels().

  {
    art::ServiceHandle<art::TFileService> tfs;
    art::ServiceHandle<geo::Geometry> geo;
    
    fPDGCodes    = tfs->make<TH1D>("pdgcodes", ";PDG Code;",               5000, -2500, 2500);
    fPi0Momentum = tfs->make<TH1D>("pi0mom",   ";#pi^{0} Momentum (GeV);", 1000, 0.,    1000.);

    fTree = tfs->make<TTree>("MCTTree","MCTTree");
    fnEnergy = tfs->make<TH1D>("nEnergy", ";n,#Lambda^{0},K^{0} Momentum (GeV);", 100, 0., 10.);
    fnDist = tfs->make<TH1D>("nDistance", ";n,#Lambda^{0},K^{0} Distance (m);", 200, -30000.0, +30000.);


    // Some histograms relating to drift electrons, active detector
    // channels and charge/energy on channels
    fnumChannels = tfs->make<TH1D>("fnumChannels", 
                                   "Active channels;Active channels;# events",
                                   256, 0, geo->Nchannels()); 
    fnumIDEs = tfs->make<TProfile>("fnumIDEs", 
                                   "Drift Electrons per channel;Channel;Drift electrons",
                                   geo->Nchannels()+1, 0, geo->Nchannels(),
                                   0, 1e4); 
    fEventCharge = tfs->make<TH1D>("fEventCharge", 
                                   "Charge in event;Total charge per event;# events",
                                   100, 0, 2.5e8); 
    fEventEnergy = tfs->make<TH1D>("fEventEnergy", 
                                   "Energy in event;Total energy per event;# events",
                                   100, 0, 1e4); 
    fChannelCharge = tfs->make<TProfile>("fChannelCharge", 
                                         "Charge on channel;Channel;Total charge per channel",
                                         geo->Nchannels()+1,0,geo->Nchannels(),
                                         0,1e5); 
    fChannelEnergy = tfs->make<TProfile>("fChannelEnergy", 
                                         "Energy on channel;Channel;Total energy per channel",
                                         geo->Nchannels()+1,0,geo->Nchannels(),
                                         0, 1e3); 


    fT4Origin = new Float_t[4];
    fT4DOrigin = new Float_t[fTNds*4];
    fT4Termination = new Float_t[4];
    fT4Momentum = new Float_t[4];
    fT4DMomentum = new Float_t[fTNds*4];
    fTDID = new Int_t[fTNds];
    fTDPdg = new Int_t[fTNds];
    fTDWt  = new Float_t[fTNds];
    fTNds4 = fTNds*4; //  TTree/Branch requirement to store this.

    fTree->Branch("MCEvt", &fTEvt, "MCEvt/I");
    fTree->Branch("MCSub", &fTSub, "MCSub/I");
    fTree->Branch("MCRun", &fTRun, "MCRun/I");
    fTree->Branch("MCWt", &fTWeight, "MCWt/F");
    fTree->Branch("MCPdg", &fTPdg, "MCPdg/I");
    fTree->Branch("MCID", &fTID, "MCID/I");
    fTree->Branch("MCParentID", &fTParentID, "MCParentID/I");
    fTree->Branch("MCNumDs", &fTNds, "MCNumDs/I");
    fTree->Branch("MCNumDs4", &fTNds4, "MCNumDs4/I");
    fTree->Branch("MCDID", fTDID, "MCDID[MCNumDs]/I");
    fTree->Branch("MCDPdg", fTDPdg, "MCDPdg[MCNumDs]/I");
    fTree->Branch("MCDWt", fTDWt, "MCDWt[MCNumDs]/I");
    fTree->Branch("MCProcess", fTProcess, "MCProcess/C");
    fTree->Branch("MCVolume", fTVolume, "MCVolume/C");
    fTree->Branch("MCTVolume", fTTVolume, "MCTVolume/C");
    fTree->Branch("MCMaterial", fTMaterial, "MCMaterial/C");
    fTree->Branch("MCDProcess", fTDProcess, "MCDProcess[MCNumDs]/C");
    fTree->Branch("MCStatus", &fTStatus, "MCStatus/I");
    fTree->Branch("MCOrigin", fT4Origin, "MCOrigin[4]/F");
    fTree->Branch("MCDOrigin", fT4DOrigin, "MCDOrigin[MCNumDs4]/F");
    fTree->Branch("MCTermination", fT4Termination, "MCTermination[4]/F");
    fTree->Branch("MCMomentum", fT4Momentum, "MCMomentum[4]/F");
    fTree->Branch("MCDMomentum", fT4DMomentum, "MCDMomentum[MCNumDs4]/F");
  
  }

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();
  }

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

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.

Referenced by art::RootOutput::beginSubRun(), art::OutputModule::doBeginRun(), art::OutputModule::doBeginSubRun(), art::OutputModule::doEndRun(), art::OutputModule::doEndSubRun(), art::ProducingService::doPostReadEvent(), art::ProducingService::doPostReadRun(), art::ProducingService::doPostReadSubRun(), art::OutputModule::doWriteEvent(), 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::EventProcessor::readEvent(), art::EventProcessor::readRun(), art::EmptyEvent::readRun_(), art::EventProcessor::readSubRun(), and art::EmptyEvent::readSubRun_().

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

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

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 larg4::LArG4Ana::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 207 of file LArG4Ana_module.cc.

References fG4ModuleLabel, fTNds, fTNdsOriginal, fTruthModuleLabel, and fhicl::ParameterSet::get().

Referenced by LArG4Ana().

  {
    fG4ModuleLabel    = p.get< std::string >("GeantModuleLabel");
    fTNdsOriginal     = p.get< int         >("Ndaughters"      );
    fTruthModuleLabel = p.get< std::string >("TruthModuleLabel"); 
    fTNds = fTNdsOriginal;
    
    return;
  }

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

inlineinherited

Definition at line 33 of file EventObserverBase.h.

  {}

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

TProfile* larg4::LArG4Ana::fChannelCharge

private

Charge per channel.

Definition at line 81 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TProfile* larg4::LArG4Ana::fChannelEnergy

private

Energy per channel.

Definition at line 82 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TH1D* larg4::LArG4Ana::fEventCharge

private

Charge collected per event.

Definition at line 79 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TH1D* larg4::LArG4Ana::fEventEnergy

private

Energy collected per event.

Definition at line 80 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

std::string larg4::LArG4Ana::fG4ModuleLabel

private

module label for the Geant

Definition at line 69 of file LArG4Ana_module.cc.

Referenced by analyze(), and reconfigure().

TH1D* larg4::LArG4Ana::fnDist

private

Definition at line 76 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TH1D* larg4::LArG4Ana::fnEnergy

private

Definition at line 75 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TH1D* larg4::LArG4Ana::fnumChannels

private

The number of channels recieving charge per event.

Definition at line 77 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TProfile* larg4::LArG4Ana::fnumIDEs

private

Number of drift electrons per channel.

Definition at line 78 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TH1D* larg4::LArG4Ana::fPDGCodes

private

Definition at line 73 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TH1D* larg4::LArG4Ana::fPi0Momentum

private

Definition at line 74 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Float_t* larg4::LArG4Ana::fT4DMomentum

private

Definition at line 110 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Float_t* larg4::LArG4Ana::fT4DOrigin

private

Definition at line 107 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Float_t* larg4::LArG4Ana::fT4Momentum

private

Definition at line 109 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Float_t* larg4::LArG4Ana::fT4Origin

private

Definition at line 106 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Float_t* larg4::LArG4Ana::fT4Termination

private

Definition at line 108 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Int_t* larg4::LArG4Ana::fTDID

private

Definition at line 95 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Int_t* larg4::LArG4Ana::fTDPdg

private

Definition at line 96 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Char_t larg4::LArG4Ana::fTDProcess[200][35]

private

Definition at line 102 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Float_t* larg4::LArG4Ana::fTDWt

private

Definition at line 97 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Int_t larg4::LArG4Ana::fTEvt

private

Definition at line 87 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Int_t larg4::LArG4Ana::fTID

private

Definition at line 91 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Char_t larg4::LArG4Ana::fTMaterial[35]

private

Definition at line 101 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Int_t larg4::LArG4Ana::fTNds

private

Definition at line 92 of file LArG4Ana_module.cc.

Referenced by analyze(), beginJob(), and reconfigure().

Int_t larg4::LArG4Ana::fTNds4

private

Definition at line 94 of file LArG4Ana_module.cc.

Referenced by beginJob().

Int_t larg4::LArG4Ana::fTNdsOriginal

private

Definition at line 93 of file LArG4Ana_module.cc.

Referenced by analyze(), and reconfigure().

Int_t larg4::LArG4Ana::fTParentID

private

Definition at line 103 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Int_t larg4::LArG4Ana::fTPdg

private

Definition at line 90 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Char_t larg4::LArG4Ana::fTProcess[35]

private

Definition at line 98 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

TTree* larg4::LArG4Ana::fTree

private

Definition at line 86 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Int_t larg4::LArG4Ana::fTRun

private

Definition at line 89 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

std::string larg4::LArG4Ana::fTruthModuleLabel

private

module label for the Geant

Definition at line 70 of file LArG4Ana_module.cc.

Referenced by reconfigure().

Int_t larg4::LArG4Ana::fTStatus

private

Definition at line 104 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Int_t larg4::LArG4Ana::fTSub

private

Definition at line 88 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Char_t larg4::LArG4Ana::fTTVolume[35]

private

Definition at line 100 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Char_t larg4::LArG4Ana::fTVolume[35]

private

Definition at line 99 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

Float_t larg4::LArG4Ana::fTWeight

private

Definition at line 105 of file LArG4Ana_module.cc.

Referenced by analyze(), and beginJob().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larsim/v07_09_00/source/larsim/LArG4/LArG4Ana_module.cc