LArSoft  v06_85_00
Liquid Argon Software toolkit - http://larsoft.org/
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Attributes | List of all members
trkf::Track3DKalman Class Reference
Inheritance diagram for trkf::Track3DKalman:
art::EDProducer art::ProducerBase art::Consumer art::EngineCreator art::ProductRegistryHelper

Public Types

using ModuleType = EDProducer
 
using WorkerType = WorkerT< EDProducer >
 
template<typename UserConfig , typename KeysToIgnore = void>
using Table = ProducerBase::Table< UserConfig, KeysToIgnore >
 

Public Member Functions

 Track3DKalman (fhicl::ParameterSet const &pset)
 
 ~Track3DKalman ()
 
void produce (art::Event &evt)
 
void beginJob ()
 
void endJob ()
 
void reconfigure (fhicl::ParameterSet const &p)
 
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 Attributes

std::string fSpacePtsModuleLabel
 
std::string fGenieGenModuleLabel
 
std::string fG4ModuleLabel
 
bool fGenfPRINT
 
TTree * tree
 
TMatrixT< Double_t > * stMCT
 
TMatrixT< Double_t > * covMCT
 
TMatrixT< Double_t > * stREC
 
TMatrixT< Double_t > * covREC
 
Float_t chi2
 
Float_t chi2ndf
 
Float_t * fpRECt3D
 
Float_t * fpRECL
 
Float_t * fpREC
 
Float_t * fpMCT
 
int nfail
 
int ndf
 
unsigned int evtt
 
unsigned int nTrks
 
unsigned int fptsNo
 
Float_t * fshx
 
Float_t * fshy
 
Float_t * fshz
 
unsigned int fDimSize
 
std::vector< double > fPosErr
 
std::vector< double > fMomErr
 
std::vector< double > fMomStart
 
genf::GFAbsTrackRep * repMC
 
genf::GFAbsTrackRep * rep
 

Detailed Description

Definition at line 83 of file Track3DKalman_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

trkf::Track3DKalman::Track3DKalman ( fhicl::ParameterSet const &  pset)
explicit

Definition at line 148 of file Track3DKalman_module.cc.

References art::EngineCreator::createEngine(), and reconfigure().

149 {
150 
151  this->reconfigure(pset);
152 
153  // create a default random engine; obtain the random seed from NuRandomService,
154  // unless overridden in configuration with key "Seed"
155  art::ServiceHandle<rndm::NuRandomService>()->createEngine(*this, pset, "Seed");
156 
157  produces< std::vector<recob::Track> >();
158  produces< std::vector<recob::SpacePoint> >();
159  produces< art::Assns<recob::Track, recob::Cluster> >();
160  produces< art::Assns<recob::Track, recob::SpacePoint> >();
161  produces< art::Assns<recob::Track, recob::Hit> >();
162 }
art::EngineCreator::createEngine
base_engine_t & createEngine(seed_t seed)
Definition: EngineCreator.cc:26
trkf::Track3DKalman::reconfigure
void reconfigure(fhicl::ParameterSet const &p)
Definition: Track3DKalman_module.cc:165
trkf::Track3DKalman::~Track3DKalman ( )

Definition at line 180 of file Track3DKalman_module.cc.

181 {
182 
183  /*
184  delete stMCT ;
185  delete covMCT;
186  delete stREC;
187  delete covREC;
188 
189  delete fpMCT;
190  delete fpREC;
191  delete fpRECL;
192  delete fpRECt3D;
193  */
194 }

Member Function Documentation

void trkf::Track3DKalman::beginJob ( )
virtual

Reimplemented from art::EDProducer.

Definition at line 197 of file Track3DKalman_module.cc.

References chi2, chi2ndf, covMCT, covREC, evtt, fDimSize, fpMCT, fpREC, fpRECL, fpRECt3D, fptsNo, fshx, fshy, fshz, art::TFileDirectory::make(), ndf, nfail, nTrks, stMCT, stREC, and tree.

198 {
199 
200 
201  art::ServiceHandle<art::TFileService> tfs;
202 
203  stMCT = new TMatrixT<Double_t>(5,1);
204  covMCT = new TMatrixT<Double_t>(5,5);
205  stREC = new TMatrixT<Double_t>(5,1);
206  covREC = new TMatrixT<Double_t>(5,5);
207 
208  fpMCT = new Float_t[4];
209  fpREC = new Float_t[4];
210  fpRECL = new Float_t[4];
211  fpRECt3D = new Float_t[4];
212  fDimSize = 400; // if necessary will get this from pset in constructor.
213 
214  fshx = new Float_t[fDimSize];
215  fshy = new Float_t[fDimSize];
216  fshz = new Float_t[fDimSize];
217 
218 // //TFile fileGENFIT("GENFITout.root","RECREATE");
219 
220 
221  tree = tfs->make<TTree>("GENFITttree","GENFITttree");
222  //tree->Branch("stMCT",&stMCT,"stMCT[5]/F"); // "TMatrixT<Double_t>"
223 
224  tree->Branch("stMCT","TMatrixD",&stMCT,64000,0);
225  //tree->Branch("covMCT",&covMCT,"covMCT[25]/F");
226  tree->Branch("covMCT","TMatrixD",&covMCT,64000,0);
227  //tree->Branch("stREC",&stREC,"stREC[5]/F");
228  tree->Branch("stREC","TMatrixD",&stREC,64000,0);
229  //tree->Branch("covREC",&covREC,"covREC[25]/F");
230  tree->Branch("covREC","TMatrixD",&covREC,64000,0);
231 
232 
233  tree->Branch("chi2",&chi2,"chi2/F");
234  tree->Branch("nfail",&nfail,"nfail/I");
235  tree->Branch("ndf",&ndf,"ndf/I");
236  tree->Branch("evtNo",&evtt,"evtNo/I");
237  tree->Branch("chi2ndf",&chi2ndf,"chi2ndf/F");
238 
239  tree->Branch("trkNo",&nTrks,"trkNo/I");
240  tree->Branch("ptsNo",&fptsNo,"ptsNo/I");
241  tree->Branch("shx",fshx,"shx[ptsNo]/F");
242  tree->Branch("shy",fshy,"shy[ptsNo]/F");
243  tree->Branch("shz",fshz,"shz[ptsNo]/F");
244 
245  tree->Branch("pMCT",fpMCT,"pMCT[4]/F");
246  tree->Branch("pRECKalF",fpREC,"pRECKalF[4]/F");
247  tree->Branch("pRECKalL",fpRECL,"pRECKalL[4]/F");
248  tree->Branch("pRECt3D",fpRECt3D,"pRECt3D[4]/F");
249 
250 
251  //TGeoManager* geomGENFIT = new TGeoManager("Geometry", "Geane geometry");
252  //TGeoManager::Import("config/genfitGeom.root");
253  // gROOT->Macro("config/Geane.C");
254 
255 }
art::TFileDirectory::make
T * make(ARGS...args) const
Definition: TFileDirectory.h:82
trkf::Track3DKalman::stMCT
TMatrixT< Double_t > * stMCT
Definition: Track3DKalman_module.cc:106
trkf::Track3DKalman::stREC
TMatrixT< Double_t > * stREC
Definition: Track3DKalman_module.cc:108
trkf::Track3DKalman::covMCT
TMatrixT< Double_t > * covMCT
Definition: Track3DKalman_module.cc:107
trkf::Track3DKalman::covREC
TMatrixT< Double_t > * covREC
Definition: Track3DKalman_module.cc:109
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().

148 {
149  if (!moduleContext_)
150  return ProductToken<T>::invalid();
151 
152  consumables_[BT].emplace_back(ConsumableType::Product,
153  TypeID{typeid(T)},
154  it.label(),
155  it.instance(),
156  it.process());
157  return ProductToken<T>{it};
158 }
art::ProductToken::invalid
static ProductToken< T > invalid()
Definition: ProductToken.h:47
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
template<typename T , art::BranchType BT>
void art::Consumer::consumesMany ( )
inherited

Definition at line 162 of file Consumer.h.

163 {
164  if (!moduleContext_)
165  return;
166 
167  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
168 }
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
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().

173 {
174  if (!moduleContext_)
175  return ViewToken<T>::invalid();
176 
177  consumables_[BT].emplace_back(ConsumableType::ViewElement,
178  TypeID{typeid(T)},
179  it.label(),
180  it.instance(),
181  it.process());
182  return ViewToken<T>{it};
183 }
art::ViewToken::invalid
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
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 Track3DKalman().

27 {
28  return rng()->createEngine(placeholder_schedule_id(), seed);
29 }
seed
long seed
Definition: chem4.cc:68
art::EngineCreator::rng
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
Definition: EngineCreator.cc:58
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().

34 {
35  return rng()->createEngine(
36  placeholder_schedule_id(), seed, kind_of_engine_to_make);
37 }
seed
long seed
Definition: chem4.cc:68
art::EngineCreator::rng
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
Definition: EngineCreator.cc:58
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().

43 {
44  return rng()->createEngine(
45  placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
46 }
seed
long seed
Definition: chem4.cc:68
art::EngineCreator::rng
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
Definition: EngineCreator.cc:58
CurrentProcessingContext const * art::EDProducer::currentContext ( ) const
protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

121  {
122  return current_context_.get();
123  }
art::EDProducer::current_context_
CPC_exempt_ptr current_context_
Definition: EDProducer.h:116
void trkf::Track3DKalman::endJob ( )
virtual

Reimplemented from art::EDProducer.

Definition at line 258 of file Track3DKalman_module.cc.

References rep, and repMC.

259 {
260  if (!rep) delete rep;
261  if (!repMC) delete repMC;
262 }
trkf::Track3DKalman::rep
genf::GFAbsTrackRep * rep
Definition: Track3DKalman_module.cc:131
trkf::Track3DKalman::repMC
genf::GFAbsTrackRep * repMC
Definition: Track3DKalman_module.cc:130
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_().

52 {
53  auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
54  return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
55 }
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_.

124  {
125  return ProducerBase::getProductID<PROD, B>(moduleDescription_,
126  instanceName);
127  }
art::EDProducer::moduleDescription_
ModuleDescription moduleDescription_
Definition: EDProducer.h:115
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 B, and art::ModuleDescription::moduleLabel().

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

58  {
59  auto const& pd =
60  get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
61  return pd.productID();
62  }
B
Int_t B
Definition: plot.C:25
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().

191 {
192  if (!moduleContext_)
193  return ProductToken<T>::invalid();
194 
195  consumables_[BT].emplace_back(ConsumableType::Product,
196  TypeID{typeid(T)},
197  it.label(),
198  it.instance(),
199  it.process());
200  return ProductToken<T>{it};
201 }
art::ProductToken::invalid
static ProductToken< T > invalid()
Definition: ProductToken.h:47
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
template<typename T , art::BranchType BT>
void art::Consumer::mayConsumeMany ( )
inherited

Definition at line 205 of file Consumer.h.

206 {
207  if (!moduleContext_)
208  return;
209 
210  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
211 }
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
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().

216 {
217  if (!moduleContext_)
218  return ViewToken<T>::invalid();
219 
220  consumables_[BT].emplace_back(ConsumableType::ViewElement,
221  TypeID{typeid(T)},
222  it.label(),
223  it.instance(),
224  it.process());
225  return ViewToken<T>{it};
226 }
art::ViewToken::invalid
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
bool art::ProducerBase::modifiesEvent ( ) const
inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

41  {
42  return true;
43  }
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_().

77 {
78  static Consumer invalid{InvalidTag{}};
79  return &invalid;
80 }
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().

90 {
91  if (!moduleContext_)
92  return;
93 
94  pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
95  for (auto& consumablesPerBranch : consumables_) {
96  cet::sort_all(consumablesPerBranch);
97  }
98 }
art::Consumer::requireConsumes_
bool requireConsumes_
Definition: Consumer.h:137
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
void trkf::Track3DKalman::produce ( art::Event &  evt)
virtual
Todo:
: remove this statement, there is no place for checks on isRealData in reconstruction code

Implements art::EDProducer.

Definition at line 266 of file Track3DKalman_module.cc.

References genf::GFTrack::addHit(), art::PtrVector< T >::begin(), chi2, chi2ndf, covMCT, covREC, util::CreateAssn(), DEFINE_ART_MODULE, e, art::PtrVector< T >::end(), art::EventID::event(), evtt, fDimSize, fGenfPRINT, fGenieGenModuleLabel, fMomErr, fMomStart, fpMCT, fPosErr, fpREC, fpRECL, fpRECt3D, fptsNo, fshx, fshy, fshz, fSpacePtsModuleLabel, art::DataViewImpl::getByLabel(), genf::GFAbsTrackRep::getChiSqu(), genf::GFAbsTrackRep::getCov(), art::RandomNumberGenerator::getEngine(), genf::GFTrack::getFailedHits(), genf::GFFieldManager::getInstance(), geo::GeometryCore::GetLArTPCVolumeName(), genf::GFAbsTrackRep::getLastPlane(), genf::GFAbsTrackRep::getNDF(), genf::GFDetPlane::getNormal(), simb::MCTruth::GetParticle(), genf::GFAbsTrackRep::getReferencePlane(), genf::GFAbsTrackRep::getState(), genf::GFAbsTrackRep::getStatusFlag(), hits(), art::Event::id(), genf::GFFieldManager::init(), art::Event::isRealData(), util::kBogusD, LOG_DEBUG, ndf, nfail, simb::MCTruth::NParticles(), nTrks, part, genf::GFDetPlane::Print(), genf::GFKalman::processTrack(), art::PtrVector< T >::push_back(), art::Event::put(), rep, repMC, art::EngineCreator::rng(), genf::GFKalman::setNumIterations(), art::PtrVector< T >::size(), sp_sort_3dz(), stMCT, stREC, track, tree, lar::dump::vector(), and GFException::what().

267 {
268 
269  rep=0;
270  repMC=0;
271 
272  // get services
273  art::ServiceHandle<geo::Geometry> geom;
274  // get the random number generator service and make some CLHEP generators
275  art::ServiceHandle<art::RandomNumberGenerator> rng;
276  CLHEP::HepRandomEngine &engine = rng->getEngine();
277  CLHEP::RandGaussQ gauss(engine);
278 
280  // Make a std::unique_ptr<> for the thing you want to put into the event
281  // because that handles the memory management for you
283  std::unique_ptr<std::vector<recob::Track> > tcol(new std::vector<recob::Track>);
284  std::unique_ptr< std::vector<recob::SpacePoint> > spcol (new std::vector<recob::SpacePoint>);
285  std::unique_ptr< art::Assns<recob::Track, recob::SpacePoint> > tspassn(new art::Assns<recob::Track, recob::SpacePoint>);
286  std::unique_ptr< art::Assns<recob::Track, recob::Cluster> > tcassn (new art::Assns<recob::Track, recob::Cluster>);
287  std::unique_ptr< art::Assns<recob::Track, recob::Hit> > thassn (new art::Assns<recob::Track, recob::Hit>);
288 
289  // define TPC parameters
290  TString tpcName = geom->GetLArTPCVolumeName();
291 
292 
293  // get input Hit object(s).
294  art::Handle< std::vector<recob::Track> > trackListHandle;
295  evt.getByLabel(fSpacePtsModuleLabel,trackListHandle);
296 
297  art::PtrVector<simb::MCTruth> mclist;
298 
299  if (!evt.isRealData())
300  {
301 
302  // std::cout << "Track3DKalman: This is MC." << std::endl;
303  // std::cout<<"Run "<<evt.run()<<" Event "<<evt.id().event()<<std::endl;
304 
305  art::Handle< std::vector<simb::MCTruth> > mctruthListHandle;
306  evt.getByLabel(fGenieGenModuleLabel,mctruthListHandle);
307 
308  for (unsigned int ii = 0; ii < mctruthListHandle->size(); ++ii)
309  {
310  art::Ptr<simb::MCTruth> mctparticle(mctruthListHandle,ii);
311  mclist.push_back(mctparticle);
312  }
313 
314  }
315 
316  //create collection of spacepoints that will be used when creating the Track object
317  std::vector< art::Ptr<recob::SpacePoint> > spacepoints;
318  art::PtrVector<recob::Track> trackIn;
319  // std::cout<<"Run "<<evt.run()<<" Event "<<evt.id().event()<<std::endl;
320  mf::LogInfo("Track3DKalman: ") << "There are " << trackListHandle->size() << " Track3Dreco/SpacePt tracks/groups (whichever) in this event.";
321 
322  art::FindManyP<recob::SpacePoint> fmsp(trackListHandle, evt, fSpacePtsModuleLabel);
323  art::FindManyP<recob::Cluster> fmc (trackListHandle, evt, fSpacePtsModuleLabel);
324  art::FindManyP<recob::Hit> fmh (trackListHandle, evt, fSpacePtsModuleLabel);
325 
326  for(unsigned int ii = 0; ii < trackListHandle->size(); ++ii)
327  {
328  art::Ptr<recob::Track> track(trackListHandle, ii);
329  trackIn.push_back(track);
330 
331  }
332 
333  TVector3 MCOrigin;
334  TVector3 MCMomentum;
335  // TVector3 posErr(.05,.05,.05); // resolution. 0.5mm
336  // TVector3 momErr(.1,.1,0.2); // GeV
337  TVector3 posErr(fPosErr[0],fPosErr[1],fPosErr[2]); // resolution. 0.5mm
338  TVector3 momErr(fMomErr[0],fMomErr[1],fMomErr[2]); // GeV
339 
340  // This is strictly for MC
342  if (!evt.isRealData())
343  {
344  // Below breaks are stupid, I realize. But rather than keep all the MC
345  // particles I just take the first primary, e.g., muon and only keep its
346  // info in the branches of the Ttree. I could generalize later, ...
347  for( unsigned int ii = 0; ii < mclist.size(); ++ii )
348  {
349  //art::Ptr<const simb::MCTruth> mc(mctruthListHandle,i);
350  art::Ptr<simb::MCTruth> mc(mclist[ii]);
351  for(int jj = 0; jj < mc->NParticles(); ++jj)
352  {
353  simb::MCParticle part(mc->GetParticle(jj));
354  mf::LogInfo("Track3DKalman: ") << "FROM MC TRUTH, the particle's pdg code is: "<<part.PdgCode()<< " with energy = "<<part.E() <<", with energy = "<<part.E()<< " and vtx and momentum in Global (not volTPC) coords are " ;
355  MCOrigin.SetXYZ(part.Vx(),part.Vy(),part.Vz()); // V for Vertex
356  MCMomentum.SetXYZ(part.Px(),part.Py(),part.Pz());
357  MCOrigin.Print();
358  MCMomentum.Print();
359  repMC = new genf::RKTrackRep(MCOrigin,
360  MCMomentum,
361  posErr,
362  momErr,
363  part.PdgCode());
364  break;
365  }
366  break;
367  }
368  //for saving of MC truth
369  stMCT->ResizeTo(repMC->getState());
370  *stMCT = repMC->getState();
371  covMCT-> ResizeTo(repMC->getCov());
372  *covMCT = repMC->getCov();
373  mf::LogInfo("Track3DKalman: ") <<" repMC, covMC are ... " ;
374  repMC->getState().Print();
375  repMC->getCov().Print();
376 
377  } // !isRealData
378 
379  art::PtrVector<recob::Track>::const_iterator trackIter = trackIn.begin();
380 
381  nTrks=0;
382  while(trackIter!=trackIn.end()) {
383  spacepoints.clear();
384  spacepoints = fmsp.at(nTrks);
385 
386  mf::LogInfo("Track3DKalman: ") << "found "<< spacepoints.size() <<" 3D spacepoint(s).";
387 
388  // Add the 3D track to the vector of the reconstructed tracks
389  if(spacepoints.size()>0){
390  // Insert the GENFIT/Kalman stuff here then make the tracks. Units are cm, GeV.
391  const double resolution = 0.5; // dunno, 5 mm
392  const int numIT = 3; // 3->1, EC, 6-Jan-2011. Back, 7-Jan-2011.
393 
394 
395  //TVector3 mom(0.0,0.0,2.0);
396  TVector3 mom(fMomStart[0],fMomStart[1],fMomStart[2]);
397  //mom.SetMag(1.);
398  TVector3 momM(mom);
399  momM.SetX(gauss.fire(momM.X(),momErr.X()/* *momM.X() */));
400  momM.SetY(gauss.fire(momM.Y(),momErr.Y()/* *momM.Y() */));
401  momM.SetZ(gauss.fire(momM.Z(),momErr.Z()/* *momM.Z() */));
402  //std::cout << "Track3DKalman: sort spacepoints by z
403 
404  std::sort(spacepoints.begin(), spacepoints.end(), sp_sort_3dz);
405 
406  //std::sort(spacepoints.begin(), spacepoints.end(), sp_sort_3dx); // Reverse sort!
407 
408  genf::GFFieldManager::getInstance()->init(new genf::GFConstField(0.,0.,0.0));
409  genf::GFDetPlane planeG((TVector3)(spacepoints[0]->XYZ()),momM);
410 
411 
412  // std::cout<<"Track3DKalman about to do GAbsTrackRep."<<std::endl;
413  // Initialize with 1st spacepoint location and a guess at the momentum.
414  rep = new genf::RKTrackRep(//posM-.5/momM.Mag()*momM,
415  (TVector3)(spacepoints[0]->XYZ()),
416  momM,
417  posErr,
418  momErr,
419  13); // mu- hypothesis
420  // std::cout<<"Track3DKalman: about to do GFTrack. repDim is " << rep->getDim() <<std::endl;
421 
422 
423  genf::GFTrack fitTrack(rep);//initialized with smeared rep
424  // Gonna sort in x cuz I want to essentially transform here to volTPC coords.
425  // volTPC coords, cuz that's what the Geant3/Geane stepper wants, as that's its understanding
426  // from the Geant4 geometry, which it'll use. EC, 7-Jan-2011.
427  int ihit = 0;
428  fptsNo = 0;
429  for (unsigned int point=0;point<spacepoints.size();++point)
430  {
431 
432  TVector3 spt3(spacepoints[point]->XYZ());
433  if (point > 0 )
434  {
435  double eps (0.1);
436  TVector3 magNew(spt3[0],spt3[1],spt3[2]);
437  TVector3 magLast(spacepoints.back()->XYZ()[0],
438  spacepoints.back()->XYZ()[1],
439  spacepoints.back()->XYZ()[2]
440  );
441  if (!(magNew.Mag()>=magLast.Mag()+eps ||
442  magNew.Mag()<=magLast.Mag()-eps)
443  )
444  continue;
445  }
446 
447  if (point%20) // Jump out of loop except on every 20th pt.
448  {
449  //continue;
450  // Icarus paper suggests we may wanna decimate our data in order to give
451  // trackfitter a better idea of multiple-scattering. EC, 7-Jan-2011.
452  //if (std::abs(spt3[0]-spacepoints.at(point-1).XYZ()[0]) < 2) continue;
453  }
454  if (fptsNo<fDimSize)
455  {
456  fshx[fptsNo] = spt3[0];
457  fshy[fptsNo] = spt3[1];
458  fshz[fptsNo] = spt3[2];
459  }
460  fptsNo++;
461 
462 
463  LOG_DEBUG("Track3DKalman: ") << "ihit xyz..." << spt3[0]<<","<< spt3[1]<<","<< spt3[2];
464  fitTrack.addHit(new genf::PointHit(spt3,resolution),
465  1,//dummy detector id
466  ihit++
467  );
468  }
469 
470  // std::cout<<"Track3DKalman about to do GFKalman."<<std::endl;
471  genf::GFKalman k;
472  //k.setBlowUpFactor(500); // Instead of 500 out of box. EC, 6-Jan-2011.
473  //k.setInitialDirection(+1); // Instead of 1 out of box. EC, 6-Jan-2011.
474  k.setNumIterations(numIT);
475  // std::cout<<"Track3DKalman back from setNumIterations."<<std::endl;
476  try{
477  // std::cout<<"Track3DKalman about to processTrack."<<std::endl;
478  k.processTrack(&fitTrack);
479  //std::cout<<"Track3DKalman back from processTrack."<<std::endl;
480  }
481  catch(GFException& e){
482  mf::LogError("Track3DKalman: ") << "just caught a GFException."<<std::endl;
483  e.what();
484  mf::LogError("Track3DKalman: ") << "Exceptions won't be further handled ->exit(1) "<<__LINE__;
485 
486  // exit(1);
487  }
488 
489  if(rep->getStatusFlag()==0) // 0 is successful completion
490  {
491  LOG_DEBUG("Track3DKalman: ") << __FILE__ << " " << __LINE__ ;
492  LOG_DEBUG("Track3DKalman: ") << "Track3DKalman.cxx: Original plane:";
493 
494  if(fGenfPRINT) planeG.Print();
495  LOG_DEBUG("Track3DKalman: ") << "Current (fit) reference Plane:";
496  if(fGenfPRINT) rep->getReferencePlane().Print();
497 
498  LOG_DEBUG("Track3DKalman: ") << "Track3DKalman.cxx: Last reference Plane:";
499  if(fGenfPRINT) rep->getLastPlane().Print();
500 
501  if(fGenfPRINT)
502  {
503  if(planeG!=rep->getReferencePlane())
504  LOG_DEBUG("Track3DKalman: ") <<"Track3DKalman: Original hit plane (not surprisingly) not current reference Plane!"<<std::endl;
505  }
506 
507  stREC->ResizeTo(rep->getState());
508  *stREC = rep->getState();
509  covREC->ResizeTo(rep->getCov());
510  *covREC = rep->getCov();
511  if(fGenfPRINT)
512  {
513  LOG_DEBUG("Track3DKalman: ") << " Final State and Cov:";
514  stREC->Print();
515  covREC->Print();
516  }
517  chi2 = rep->getChiSqu();
518  ndf = rep->getNDF();
519  nfail = fitTrack.getFailedHits();
520  chi2ndf = chi2/ndf;
521  double dircoss[3],dircose[3];
522  (*trackIter)->Direction(dircoss,dircose);
523 
524  for (int ii=0;ii<3;++ii)
525  {
526  fpMCT[ii] = MCMomentum[ii]/MCMomentum.Mag();
527  fpREC[ii] = rep->getReferencePlane().getNormal()[ii];
528  fpRECL[ii] = rep->getLastPlane().getNormal()[ii];
529  fpRECt3D[ii] = dircoss[ii];
530  }
531  fpMCT[3] = MCMomentum.Mag();
532  fpREC[3] = -1.0/(*stREC)[0][0];
533 
534  evtt = (unsigned int) evt.id().event();
535  mf::LogInfo("Track3DKalman: ") << "Track3DKalman about to do tree->Fill(). Chi2/ndf is " << chi2/ndf << ". All in volTPC coords .... pMCT[0-3] is " << fpMCT[0] << ", " << fpMCT[1] << ", " << fpMCT[2] << ", " << fpMCT[3] << ". pREC[0-3] is " << fpREC[0] << ", "<< fpREC[1] << ", " << fpREC[2] << ", " << fpREC[3] << ".";
536 
537  tree->Fill();
538 
539  // Get the clusters associated to each track in induction and collection view
540  std::vector< art::Ptr<recob::Cluster> > clusters = fmc.at(nTrks);
541  std::vector< art::Ptr<recob::Hit> > hits = fmh.at(nTrks);
542 
543  // Use new Track constructor... EC, 21-Apr-2011.
544  std::vector<TVector3> dircos(spacepoints.size());
545  dircos[0] = TVector3(fpREC);
546  dircos.back() = TVector3(fpRECL);
547 
548  // fill a vector of TVector3 with the space points used to make this track
549  std::vector<TVector3> xyz(spacepoints.size());
550  size_t spStart = spcol->size();
551  for(size_t tv = 0; tv < spacepoints.size(); ++tv){
552  xyz[tv] = TVector3(spacepoints[tv]->XYZ());
553  spcol->push_back(*(spacepoints[tv].get()));
554  }
555  size_t spEnd = spcol->size();
556 
557  tcol->push_back(recob::Track(xyz, dircos, std::vector< std::vector<double> >(0),
558  std::vector<double>(2, util::kBogusD), tcol->size()-1));
559 
560  // associate the track with its clusters and tracks
561  util::CreateAssn(*this, evt, *tcol, clusters, *tcassn);
562  util::CreateAssn(*this, evt, *tcol, hits, *thassn);
563 
564  // associate the track to the space points
565  util::CreateAssn(*this, evt, *tcol, *spcol, *tspassn, spStart, spEnd);
566 
567  } // getStatusFlag
568  } // spacepoints.size()>0
569 
570  //
571  //std::cout<<"Track3DKalman found "<< tcol->size() <<" 3D track(s)"<<std::endl;
572  if(trackIter!=trackIn.end()) trackIter++;
573 
574  nTrks++;
575 
576  } // end loop over Track3Dreco/SpacePt tracks/groups (whichever) we brought into this event.
577 
578  evt.put(std::move(tcol));
579  evt.put(std::move(spcol));
580  evt.put(std::move(tcassn));
581  evt.put(std::move(thassn));
582  evt.put(std::move(tspassn));
583 
584 } // end method
genf::GFAbsTrackRep::getNDF
unsigned int getNDF() const
Definition: GFAbsTrackRep.h:253
sp_sort_3dz
static bool sp_sort_3dz(const art::Ptr< recob::SpacePoint > &h1, const art::Ptr< recob::SpacePoint > &h2)
Definition: Track3DKalman_module.cc:138
genf::GFAbsTrackRep::getReferencePlane
const GFDetPlane & getReferencePlane() const
Definition: GFAbsTrackRep.h:285
genf::GFAbsTrackRep::getLastPlane
GFDetPlane getLastPlane() const
Definition: GFAbsTrackRep.h:242
mf::LogInfo
MaybeLogger_< ELseverityLevel::ELsev_info, false > LogInfo
Definition: MessageLogger.h:212
art::PtrVector::begin
iterator begin()
Definition: PtrVector.h:223
trkf::Track3DKalman::fMomStart
std::vector< double > fMomStart
Definition: Track3DKalman_module.cc:129
GFException::what
virtual const char * what() const
standard error message handling for exceptions. use like "std::cerr << e.what();" ...
Definition: GFException.cxx:46
genf::GFDetPlane::Print
void Print(std::ostream &out=std::cout) const
Definition: GFDetPlane.cxx:242
genf::GFAbsTrackRep::getState
const TMatrixT< Double_t > & getState() const
Definition: GFAbsTrackRep.h:198
mf::LogError
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Definition: MessageLogger.h:214
art::Event::isRealData
bool isRealData() const
Definition: Event.h:83
trkf::Track3DKalman::fMomErr
std::vector< double > fMomErr
Definition: Track3DKalman_module.cc:128
trkf::Track3DKalman::rep
genf::GFAbsTrackRep * rep
Definition: Track3DKalman_module.cc:131
art::Event::put
ProductID put(std::unique_ptr< PROD > &&product)
Definition: Event.h:102
genf::GFDetPlane::getNormal
TVector3 getNormal() const
Definition: GFDetPlane.cxx:140
lar::dump::vector
auto vector(Vector const &v)
Returns a manipulator which will print the specified array.
Definition: DumpUtils.h:265
art::RandomNumberGenerator::getEngine
base_engine_t & getEngine() const
hits
void hits()
Definition: readHits.C:15
art::PtrVector::push_back
void push_back(Ptr< U > const &p)
Definition: PtrVector.h:441
part
TString part[npart]
Definition: Style.C:32
art::PtrVector::end
iterator end()
Definition: PtrVector.h:237
genf::GFKalman::setNumIterations
void setNumIterations(Int_t i)
Set number of iterations for Kalman Filter.
Definition: GFKalman.h:90
genf::GFKalman::processTrack
void processTrack(GFTrack *)
Performs fit on a GFTrack.
Definition: GFKalman.cxx:48
genf::GFFieldManager::getInstance
static GFFieldManager * getInstance()
Definition: GFFieldManager.h:70
trkf::Track3DKalman::repMC
genf::GFAbsTrackRep * repMC
Definition: Track3DKalman_module.cc:130
util::CreateAssn
bool CreateAssn(PRODUCER const &prod, art::Event &evt, std::vector< T > const &a, art::Ptr< U > const &b, art::Assns< U, T > &assn, std::string a_instance, size_t indx=UINT_MAX)
Creates a single one-to-one association.
Definition: AssociationUtil.h:591
genf::GFAbsTrackRep::getCov
const TMatrixT< Double_t > & getCov() const
Definition: GFAbsTrackRep.h:199
art::PtrVector::size
size_type size() const
Definition: PtrVector.h:308
art::PtrVector::const_iterator
data_t::const_iterator const_iterator
Definition: PtrVector.h:61
art::EngineCreator::rng
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
Definition: EngineCreator.cc:58
GFException
Exception class for error handling in GENFIT (provides storage for diagnostic information) ...
Definition: GFException.h:50
trkf::Track3DKalman::stMCT
TMatrixT< Double_t > * stMCT
Definition: Track3DKalman_module.cc:106
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
Definition: DataViewImpl.h:344
trkf::Track3DKalman::stREC
TMatrixT< Double_t > * stREC
Definition: Track3DKalman_module.cc:108
genf::GFFieldManager::init
void init(GFAbsBField *b)
set the magntic field here. Magnetic field classes must be derived from GFAbsBField ...
Definition: GFFieldManager.h:66
LOG_DEBUG
#define LOG_DEBUG(id)
Definition: MessageLogger.h:263
art::EventID::event
EventNumber_t event() const
Definition: EventID.h:117
util::kBogusD
constexpr double kBogusD
obviously bogus double value
Definition: PhysicalConstants.h:71
genf::GFAbsTrackRep::getChiSqu
double getChiSqu() const
Definition: GFAbsTrackRep.h:245
e
Float_t e
Definition: plot.C:34
track
Float_t track
Definition: plot.C:34
trkf::Track3DKalman::fPosErr
std::vector< double > fPosErr
Definition: Track3DKalman_module.cc:127
art::Event::id
EventID id() const
Definition: Event.h:56
recob::Track
Track from a non-cascading particle.A recob::Track consists of a recob::TrackTrajectory, plus additional members relevant for a "fitted" track:
Definition: Track.h:51
trkf::Track3DKalman::covMCT
TMatrixT< Double_t > * covMCT
Definition: Track3DKalman_module.cc:107
geo::GeometryCore::GetLArTPCVolumeName
std::string GetLArTPCVolumeName(geo::TPCID const &tpcid) const
Return the name of specified LAr TPC volume.
Definition: GeometryCore.cxx:667
trkf::Track3DKalman::covREC
TMatrixT< Double_t > * covREC
Definition: Track3DKalman_module.cc:109
void trkf::Track3DKalman::reconfigure ( fhicl::ParameterSet const &  p)

Definition at line 165 of file Track3DKalman_module.cc.

References fG4ModuleLabel, fGenfPRINT, fGenieGenModuleLabel, fMomErr, fMomStart, fPosErr, fSpacePtsModuleLabel, and fhicl::ParameterSet::get().

Referenced by Track3DKalman().

166 {
167 
168  fSpacePtsModuleLabel = pset.get< std::string >("SpacePtsModuleLabel");
169  fGenieGenModuleLabel = pset.get< std::string >("GenieGenModuleLabel");
170  fG4ModuleLabel = pset.get< std::string >("G4ModuleLabel");
171 
172  fPosErr = pset.get< std::vector < double > >("PosErr3"); // resolution. cm
173  fMomErr = pset.get< std::vector < double > >("MomErr3"); // GeV
174  fMomStart = pset.get< std::vector < double > >("MomStart3"); // Will be unit norm'd.
175  fGenfPRINT = pset.get< bool >("GenfPRINT");
176 
177 }
trkf::Track3DKalman::fMomStart
std::vector< double > fMomStart
Definition: Track3DKalman_module.cc:129
trkf::Track3DKalman::fMomErr
std::vector< double > fMomErr
Definition: Track3DKalman_module.cc:128
trkf::Track3DKalman::fPosErr
std::vector< double > fPosErr
Definition: Track3DKalman_module.cc:127
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().

126 {
127  if (!moduleContext_)
128  return;
129 
130  // If none of the branches have missing consumes statements, exit early.
131  if (std::all_of(cbegin(missingConsumes_),
132  cend(missingConsumes_),
133  [](auto const& perBranch) { return perBranch.empty(); }))
134  return;
135 
136  constexpr cet::HorizontalRule rule{60};
137  mf::LogPrint log{"MTdiagnostics"};
138  log << '\n'
139  << rule('=') << '\n'
140  << "The following consumes (or mayConsume) statements are missing from\n"
141  << module_context(moduleDescription_) << '\n'
142  << rule('-') << '\n';
143 
144  cet::for_all_with_index(
145  missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
146  for (auto const& pi : perBranch) {
147  log << " "
148  << assemble_consumes_statement(static_cast<BranchType>(i), pi)
149  << '\n';
150  }
151  });
152  log << rule('=');
153 }
mf::MaybeLogger_
Definition: MessageLogger.h:114
art::Consumer::moduleDescription_
cet::exempt_ptr< ModuleDescription const > moduleDescription_
Definition: Consumer.h:140
util::pi
constexpr T pi()
Returns the constant pi (up to 35 decimal digits of precision)
Definition: PhysicalConstants.h:80
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
art::Consumer::missingConsumes_
ConsumableProductSets missingConsumes_
Definition: Consumer.h:139
void art::Consumer::validateConsumedProduct ( BranchType const  bt,
ProductInfo const &  pi 
)
protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

103 {
104  // Early exits if consumes tracking has been disabled or if the
105  // consumed product is an allowed consumable.
106  if (!moduleContext_)
107  return;
108 
109  if (cet::binary_search_all(consumables_[bt], pi))
110  return;
111 
112  if (requireConsumes_) {
113  throw Exception(errors::ProductRegistrationFailure,
114  "Consumer: an error occurred during validation of a "
115  "retrieved product\n\n")
116  << "The following consumes (or mayConsume) statement is missing from\n"
117  << module_context(moduleDescription_) << ":\n\n"
118  << " " << assemble_consumes_statement(bt, pi) << "\n\n";
119  }
120 
121  missingConsumes_[bt].insert(pi);
122 }
art::Consumer::moduleDescription_
cet::exempt_ptr< ModuleDescription const > moduleDescription_
Definition: Consumer.h:140
art::Consumer::requireConsumes_
bool requireConsumes_
Definition: Consumer.h:137
art::Exception
cet::coded_exception< errors::ErrorCodes, ExceptionDetail::translate > Exception
Definition: Exception.h:66
util::pi
constexpr T pi()
Returns the constant pi (up to 35 decimal digits of precision)
Definition: PhysicalConstants.h:80
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
art::Consumer::missingConsumes_
ConsumableProductSets missingConsumes_
Definition: Consumer.h:139

Member Data Documentation

Float_t trkf::Track3DKalman::chi2
private

Definition at line 110 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t trkf::Track3DKalman::chi2ndf
private

Definition at line 111 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

TMatrixT<Double_t>* trkf::Track3DKalman::covMCT
private

Definition at line 107 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

TMatrixT<Double_t>* trkf::Track3DKalman::covREC
private

Definition at line 109 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

unsigned int trkf::Track3DKalman::evtt
private

Definition at line 119 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

unsigned int trkf::Track3DKalman::fDimSize
private

Definition at line 125 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

std::string trkf::Track3DKalman::fG4ModuleLabel
private

Definition at line 100 of file Track3DKalman_module.cc.

Referenced by reconfigure().

bool trkf::Track3DKalman::fGenfPRINT
private

Definition at line 101 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

std::string trkf::Track3DKalman::fGenieGenModuleLabel
private

Definition at line 99 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

std::vector<double> trkf::Track3DKalman::fMomErr
private

Definition at line 128 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

std::vector<double> trkf::Track3DKalman::fMomStart
private

Definition at line 129 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

Float_t* trkf::Track3DKalman::fpMCT
private

Definition at line 116 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

std::vector<double> trkf::Track3DKalman::fPosErr
private

Definition at line 127 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

Float_t* trkf::Track3DKalman::fpREC
private

Definition at line 115 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fpRECL
private

Definition at line 114 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fpRECt3D
private

Definition at line 113 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

unsigned int trkf::Track3DKalman::fptsNo
private

Definition at line 121 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fshx
private

Definition at line 122 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fshy
private

Definition at line 123 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fshz
private

Definition at line 124 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

std::string trkf::Track3DKalman::fSpacePtsModuleLabel
private

Definition at line 98 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

int trkf::Track3DKalman::ndf
private

Definition at line 118 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

int trkf::Track3DKalman::nfail
private

Definition at line 117 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

unsigned int trkf::Track3DKalman::nTrks
private

Definition at line 120 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

genf::GFAbsTrackRep* trkf::Track3DKalman::rep
private

Definition at line 131 of file Track3DKalman_module.cc.

Referenced by endJob(), and produce().

genf::GFAbsTrackRep* trkf::Track3DKalman::repMC
private

Definition at line 130 of file Track3DKalman_module.cc.

Referenced by endJob(), and produce().

TMatrixT<Double_t>* trkf::Track3DKalman::stMCT
private

Definition at line 106 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

TMatrixT<Double_t>* trkf::Track3DKalman::stREC
private

Definition at line 108 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

TTree* trkf::Track3DKalman::tree
private

Definition at line 104 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

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