Inheritance diagram for ncid::MakeNCCosRej:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

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

	~MakeNCCosRej ()

void	beginRun (art::Run &evt) override

void	reconfigure (const fhicl::ParameterSet &p)

void	produce (art::Event &e) override

template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const

template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const

bool	modifiesEvent () 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

void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)

void	prepareForJob (fhicl::ParameterSet const &pset)

void	showMissingConsumes () const

Private Member Functions
double	TransMomFraction (const std::vector< art::Ptr< rb::Prong > > &prongs) const

Private Attributes
std::string	fSliceLabel

std::string	fProngLabel

std::string	fProngInstance

std::string	fCVNLabel

std::string	fShowerLIDLabel

std::string	fShowerLabel

std::string	fCosmicTrackLabel

std::string	fCellHitLabel

std::string	fCosRejLabel

std::string	fPIDLibPath

std::string	fNCCosRejPIDFile

std::string	fNCCosRejPIDGFile

TMVA::Reader *	fReaderDT

TMVA::Reader *	fReaderGBDT

float	TMVAvars [13]

fhicl::ParameterSet	fParticleIDAlgPSet
	Particle ID alorithm (loglikelihoods and dE/dx) More...

Detailed Description

Definition at line 53 of file MakeNCCosRej_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

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

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

ncid::MakeNCCosRej::MakeNCCosRej ( fhicl::ParameterSet const & p )

explicit

Definition at line 104 of file MakeNCCosRej_module.cc.

References reconfigure().

                                                         :
     fReaderDT(nullptr),
     fReaderGBDT(nullptr),
     //fParticleIDAlg(0),
     fParticleIDAlgPSet(p.get< fhicl::ParameterSet >("ParticleIDAlgPSet"))
   {
     produces< std::vector<ncid::NCCosRej> > ();
     produces< art::Assns<ncid::NCCosRej, rb::Cluster> > ();
     this->reconfigure(p);
   }

ncid::MakeNCCosRej::~MakeNCCosRej ( )

Definition at line 115 of file MakeNCCosRej_module.cc.

References fReaderDT, and fReaderGBDT.

   {
     if(fReaderDT)   delete fReaderDT;
     if(fReaderGBDT) delete fReaderGBDT;
   }

Member Function Documentation

void ncid::MakeNCCosRej::beginRun ( art::Run & evt )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 121 of file MakeNCCosRej_module.cc.

References allTimeWatchdog::endl, util::EnvExpansion(), cet::file_exists(), fNCCosRejPIDFile, fNCCosRejPIDGFile, fPIDLibPath, fReaderDT, fReaderGBDT, geom(), string, and TMVAvars.

   {
     std::string pidlibpath = util::EnvExpansion(fPIDLibPath);
     
     // Make sure we can find the Real Adaptive Boosted Decision Trees weight file before we set up MVA
     if(!cet::file_exists(pidlibpath+fNCCosRejPIDFile)){
       throw cet::exception("MakeNCCosRej")
         << "Couldn't find file " << pidlibpath+fNCCosRejPIDFile <<std::endl;
     }
     // Make sure we can find the Gradient Boosted Decision Trees weight file before we set up MVA
     if(!cet::file_exists(pidlibpath+fNCCosRejPIDGFile)){
       throw cet::exception("MakeNCCosRej")
         << "Couldn't find file " << pidlibpath+fNCCosRejPIDGFile <<std::endl;
     }
     
     art::ServiceHandle<geo::Geometry>  geom;
     
     fReaderDT = new TMVA::Reader;
     fReaderDT->AddVariable("cosmicid",   &TMVAvars[0]);  //cosmicid
     fReaderDT->AddVariable("partptp",    &TMVAvars[1]);  //partptp
     fReaderDT->AddVariable("shwnhit",    &TMVAvars[2]);  //shwnhit
     fReaderDT->AddVariable("shwxminusy", &TMVAvars[3]);  //shwminusy
     fReaderDT->AddVariable("shwxplusy",  &TMVAvars[4]);  //shwxplusy
     fReaderDT->AddVariable("shwxovery",  &TMVAvars[5]);  //shwxpvery
     fReaderDT->AddVariable("shwcalE",    &TMVAvars[6]);  //shwcalE
     fReaderDT->AddVariable("shwdirY",    &TMVAvars[7]);  //shwdirY
     fReaderDT->AddVariable("shwlen",     &TMVAvars[8]);  //shwlen
     fReaderDT->AddVariable("shwwwidth",  &TMVAvars[9]);  //shwwwidth
     fReaderDT->AddVariable("shwGap",     &TMVAvars[10]); //shwGap
     fReaderDT->AddVariable("nshwlid",    &TMVAvars[11]); //nshwlid
     fReaderDT->AddVariable("nmiphit",    &TMVAvars[12]); //nmiphit    
     fReaderDT->BookMVA("BDTA",pidlibpath+fNCCosRejPIDFile);
 
     fReaderGBDT = new TMVA::Reader;
     fReaderGBDT->AddVariable("cosmicid",   &TMVAvars[0]);  //cosmicid
     fReaderGBDT->AddVariable("partptp",    &TMVAvars[1]);  //partptp
     fReaderGBDT->AddVariable("shwnhit",    &TMVAvars[2]);  //shwnhit
     fReaderGBDT->AddVariable("shwxminusy", &TMVAvars[3]);  //shwminusy
     fReaderGBDT->AddVariable("shwxplusy",  &TMVAvars[4]);  //shwxplusy
     fReaderGBDT->AddVariable("shwxovery",  &TMVAvars[5]);  //shwxpvery
     fReaderGBDT->AddVariable("shwcalE",    &TMVAvars[6]);  //shwcalE
     fReaderGBDT->AddVariable("shwdirY",    &TMVAvars[7]);  //shwdirY
     fReaderGBDT->AddVariable("shwlen",     &TMVAvars[8]);  //shwlen
     fReaderGBDT->AddVariable("shwwwidth",  &TMVAvars[9]);  //shwwwidth
     fReaderGBDT->AddVariable("shwGap",     &TMVAvars[10]); //shwGap
     fReaderGBDT->AddVariable("nshwlid",    &TMVAvars[11]); //nshwlid
     fReaderGBDT->AddVariable("nmiphit",    &TMVAvars[12]); //nmiphit    
     fReaderGBDT->BookMVA("BDTG", pidlibpath+fNCCosRejPIDGFile);
 }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

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

inherited

Definition at line 146 of file Consumer.h.

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

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

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 161 of file Consumer.h.

References T.

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

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

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

inherited

Definition at line 171 of file Consumer.h.

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

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

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

inherited

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

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

inherited

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

inherited

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

protectedinherited

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

inherited

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

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const & instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

Referenced by skim::NueSkimmer::CopyMichelSlice(), and skim::NueSkimmer::CopyMichelTrack().

   {
     return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                                instanceName);
   }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID	(	ModuleDescription const &	moduleDescription,
		std::string const &	instanceName
	)		const

inherited

Definition at line 56 of file ProducerBase.h.

References art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

   {
     auto const& pd =
       get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
     return pd.productID();
   }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

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

inherited

Definition at line 189 of file Consumer.h.

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

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

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 204 of file Consumer.h.

References T.

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

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

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

inherited

Definition at line 214 of file Consumer.h.

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

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

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID(), and string.

     {
       return true;
     }

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

staticinherited

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

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

protectedinherited

void ncid::MakeNCCosRej::produce ( art::Event & e )

overridevirtual

Implements art::EDProducer.

Definition at line 189 of file MakeNCCosRej_module.cc.

References rb::CellHit::Cell(), slid::CompareByEnergy(), ana::cosmicid, util::CreateAssn(), fCellHitLabel, fCosRejLabel, fCVNLabel, art::fill_ptr_vector(), check_time_usage::float, fProngInstance, fProngLabel, fReaderDT, fReaderGBDT, fShowerLIDLabel, fSliceLabel, geom(), art::DataViewImpl::getByLabel(), rb::RecoHit::IsCalibrated(), rb::Cluster::IsNoise(), m, ana::nmiphit, ana::nshwlid, partptp, rb::RecoHit::PECorr(), art::PtrVector< T >::push_back(), art::Event::put(), rb::RecoHit::RecoHit(), ncid::NCCosRej::SetCosPIDDT(), ncid::NCCosRej::SetCosPIDDTG(), ncid::NCCosRej::SetProngTransMom(), ana::shwcalE, ana::shwdirY, ana::shwGap, ana::shwlen, ana::shwnhit, ana::shwnhitx, ana::shwnhity, ana::shwxminusy, ana::shwxovery, ana::shwxplusy, TMVAvars, and TransMomFraction().

   {
     art::ServiceHandle<geo::Geometry>  geom;
     art::ServiceHandle<cheat::BackTracker> bt;
     
     std::unique_ptr< std::vector<ncid::NCCosRej> >
       ncCosRejjies(new std::vector< ncid::NCCosRej >);
     std::unique_ptr< art::Assns<ncid::NCCosRej, rb::Cluster> >
       sliceAssn(new art::Assns<ncid::NCCosRej, rb::Cluster>);
     // Get the slices in the event
     art::Handle<std::vector<rb::Cluster> > sliceHandle;
     evt.getByLabel( fSliceLabel, sliceHandle);
     
     std::vector< art::Ptr<rb::Cluster> > slices;
     art::fill_ptr_vector(slices, sliceHandle);
     
     art::Handle< std::vector<rb::CellHit> > htHandle;
     evt.getByLabel(fCellHitLabel, htHandle);
     
     art::PtrVector<rb::CellHit> allHits;
     
     for(unsigned int m = 0; m < htHandle->size(); ++m){
       allHits.push_back(art::Ptr<rb::CellHit>(htHandle,m));
     }
     // Get things associated with slices
     art::FindManyP<rb::Prong> prongAssn(sliceHandle, evt, art::InputTag(fProngLabel, fProngInstance));
     art::FindMany<slid::ShowerLID> showerLidAssn(sliceHandle, evt, fShowerLIDLabel);    
     art::FindManyP<cosrej::NueCosRej> fmcr(slices, evt, fCosRejLabel);
 
     const int nslices = slices.size();     
     for(int isli = 0; isli < nslices; ++isli){
       // Loop over slices
       art::Ptr<rb::Cluster> thisSlice = slices.at(isli);
       
       if(thisSlice->IsNoise())  continue;
       if(!prongAssn.isValid())  continue;
       
       ncid::NCCosRej nccosrej;
       std::vector< art::Ptr<rb::Prong> > prongs = prongAssn.at(isli);
       
       if(prongs.empty()) continue;
       
       nccosrej.SetProngTransMom(TransMomFraction(prongs));
       std::vector< art::Ptr<cosrej::NueCosRej> > cosrejs = fmcr.at(isli);
       
       if(cosrejs.empty()) continue;
       
       // Need CVN, input to Real Adaptive BDT
       art::FindManyP<cvn::Result> fmCVN(slices, evt, fCVNLabel);
 
       if(!fmCVN.isValid()) continue;
       
       std::vector< art::Ptr<cvn::Result> > cvns = fmCVN.at(isli);   
       // CVN values
       float cosmicid  = 0.;
       if(!cvns.empty()) cosmicid = cvns[0]->fOutput[14];
       
       // LIDShower values 
       float shwGap    = -5.;
       float shwnhit   = -5.; 
       float shwnhitx  = -5.;
       float shwnhity  = -5.;
       float shwwwidth = -5.;
       float shwcalE   = -5.;    
       float shwdirY   = -5.;
       float shwlen    = -5.;
       float partptp   = -5.;
       /////////////////////// Fill the PIDs ////////////////////////////  
       // Only execute when there are showers associated with the shwlids
       if(!showerLidAssn.isValid()) continue;
       
       std::vector<const slid::ShowerLID*> shwlids = showerLidAssn.at(isli);
       if(shwlids.empty()) continue;
       
       sort(shwlids.begin(), shwlids.end(), slid::CompareByEnergy);
       art::FindOneP<rb::Shower> foSh(shwlids, evt, fShowerLIDLabel);
       float nshwlid = shwlids.size();
       
       if(!foSh.isValid()) continue;
       
       cet::maybe_ref< art::Ptr<rb::Shower> const > roSh(foSh.at(0));
       art::Ptr<rb::Shower> shower = roSh.ref();
       std::vector< art::Ptr<rb::Shower> > showers;
       showers.push_back(shower);
       art::FindOneP<slid::ShowerLID> foShLID(showers, evt,fShowerLIDLabel);
 
       if(!foShLID.isValid()) continue;
       
       cet::maybe_ref<art::Ptr<slid::ShowerLID> const> roShLID(foShLID.at(0));
       art::Ptr<slid::ShowerLID> showerLID = roShLID.ref();
 
       int nMipHit = 0;
       for(size_t hitIdx = 0; hitIdx < slices[isli]->NCell(); ++hitIdx){
         const art::Ptr<rb::CellHit>&  chptr = slices[isli]->Cell(hitIdx);
         const rb::RecoHit rhit = slices[isli]->RecoHit(chptr);
         if(!rhit.IsCalibrated()) continue;
         if((rhit.PECorr() > 100.0) && (rhit.PECorr() < 245.0)) nMipHit++;
       } // end of hits loop 
    
       shwGap    = showerLID->Gap();
       shwnhit   = shower->NCell();
       shwwwidth = showerLID->Radius();
       shwcalE   = showerLID->ShowerEnergy();
       shwnhitx  = shower->NXCell();
       shwdirY   = (float)shower->Dir().Y();
       shwnhity  = shower->NYCell();
       shwlen    = shower->TotalLength();
       partptp   = cosrejs[0]->ParticleShowerTransMom();
       float shwxminusy = shwnhitx - shwnhity; 
       float shwxplusy  = shwnhitx + shwnhity; 
       float shwxovery  = shwxminusy/shwxplusy;        
       float nmiphit    = nMipHit;
 
 
       TMVAvars[0]  = cosmicid;                     //ncid
       TMVAvars[1]  = partptp;                      //partptp
       //TMVAvars[1]  = nccosrej.ProngTransMom();     //partptp
       TMVAvars[2]  = shwnhit;                      //shwnhit;
       TMVAvars[3]  = shwxminusy;
       TMVAvars[4]  = shwxplusy;
       TMVAvars[5]  = shwxovery;
       TMVAvars[6]  = shwcalE;                     //shwcalE
       TMVAvars[7]  = shwdirY;                     //shwdirY
       TMVAvars[8]  = shwlen;                      //shwlen
       TMVAvars[9]  = shwwwidth;                   //shwwidth
       TMVAvars[10] = shwGap;                      //shwgap
       TMVAvars[11] = nshwlid;
       TMVAvars[12] = nmiphit;                     //nmiphit;
       
       nccosrej.SetCosPIDDT(fReaderDT->EvaluateMVA("BDTA"));
       nccosrej.SetCosPIDDTG(fReaderGBDT->EvaluateMVA("BDTG"));    
       ncCosRejjies->push_back(nccosrej);
       util::CreateAssn(*this, evt, *ncCosRejjies, thisSlice ,*sliceAssn);
     } // end loop over slices
     
     evt.put(std::move(sliceAssn));
     evt.put(std::move(ncCosRejjies));
   } // end of producer

void ncid::MakeNCCosRej::reconfigure ( const fhicl::ParameterSet & p )

Definition at line 171 of file MakeNCCosRej_module.cc.

References fCellHitLabel, fCosRejLabel, fCVNLabel, fNCCosRejPIDFile, fNCCosRejPIDGFile, fParticleIDAlgPSet, fPIDLibPath, fProngInstance, fProngLabel, fShowerLabel, fShowerLIDLabel, fSliceLabel, fhicl::ParameterSet::get(), and string.

Referenced by MakeNCCosRej().

   {
     //fMCTruthModuleLabel   = p.get< std::string >("MCTruthModuleLabel");
     fSliceLabel           = p.get< std::string >("SliceLabel");
     fCellHitLabel         = p.get< std::string >("CellHitLabel");
     fProngLabel           = p.get< std::string >("ProngLabel");
     fProngInstance        = p.get< std::string >("ProngInstance");
     fCVNLabel             = p.get< std::string >("CVNLabel");
     fShowerLabel          = p.get< std::string >("ShowerLabel");
     fShowerLIDLabel       = p.get< std::string >("ShowerLIDLabel");
     fPIDLibPath           = p.get< std::string >("PIDLibPath");
     fCosRejLabel          = p.get< std::string >("CosRejLabel");
     fNCCosRejPIDFile      = p.get< std::string >("NCCosRejPIDFile");
     fNCCosRejPIDGFile     = p.get< std::string >("NCCosRejPIDGFile");
     fParticleIDAlgPSet    = p.get< fhicl::ParameterSet >("ParticleIDAlgPSet");  
   }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

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

double ncid::MakeNCCosRej::TransMomFraction ( const std::vector< art::Ptr< rb::Prong > > & prongs ) const

private

Return transverse momentum fraction of the event. Calculation based on reconstructed prongs.

Definition at line 329 of file MakeNCCosRej_module.cc.

References geom(), geo::GeometryBase::NuMIBeamDirection(), runNovaSAM::ret, and w.

Referenced by produce().

   {
     if(prongs.empty()) return 0;
     
     TVector3 ret;
     for(const art::Ptr<rb::Prong>& prong: prongs){
       if(!prong->Is3D()) continue;
       const double w = prong->CalorimetricEnergy();
       ret += w*prong->Dir();
     }
     
     ret = ret.Unit();
     
     art::ServiceHandle<geo::Geometry> geom;
     const TVector3 beamDir = geom->NuMIBeamDirection();
     const double  beamProj = ret.Dot(beamDir);
     
     return (ret-beamProj*beamDir).Mag();
   }