Inheritance diagram for trkf::Track3DKalman:

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 82 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 147 of file Track3DKalman_module.cc.

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

 {
 
   this->reconfigure(pset);
   
   // create a default random engine; obtain the random seed from NuRandomService,
   // unless overridden in configuration with key "Seed"
   art::ServiceHandle<rndm::NuRandomService>()->createEngine(*this, pset, "Seed");
   
   produces< std::vector<recob::Track> >();
   produces< std::vector<recob::SpacePoint>              >();
   produces< art::Assns<recob::Track, recob::Cluster>    >();
   produces< art::Assns<recob::Track, recob::SpacePoint> >();
   produces< art::Assns<recob::Track, recob::Hit>        >();
 }

trkf::Track3DKalman::~Track3DKalman ( )

Definition at line 179 of file Track3DKalman_module.cc.

 {
 
   /*
     delete stMCT ;
     delete covMCT;
     delete stREC;
     delete covREC;
     
     delete fpMCT;
     delete fpREC;
     delete fpRECL;
     delete fpRECt3D;
   */
 }

Member Function Documentation

void trkf::Track3DKalman::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 196 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.

 {
 
 
   art::ServiceHandle<art::TFileService> tfs;
   
   stMCT  = new TMatrixT<Double_t>(5,1);
   covMCT = new TMatrixT<Double_t>(5,5);
   stREC  = new TMatrixT<Double_t>(5,1);
   covREC = new TMatrixT<Double_t>(5,5);
   
   fpMCT = new Float_t[4];
   fpREC = new Float_t[4];
   fpRECL = new Float_t[4];
   fpRECt3D = new Float_t[4];
   fDimSize = 400; // if necessary will get this from pset in constructor.
 
   fshx = new Float_t[fDimSize];
   fshy = new Float_t[fDimSize];
   fshz = new Float_t[fDimSize];
 
 //   //TFile fileGENFIT("GENFITout.root","RECREATE");
 
   
   tree = tfs->make<TTree>("GENFITttree","GENFITttree");
   //tree->Branch("stMCT",&stMCT,"stMCT[5]/F"); // "TMatrixT<Double_t>"
 
   tree->Branch("stMCT","TMatrixD",&stMCT,64000,0);
   //tree->Branch("covMCT",&covMCT,"covMCT[25]/F");
   tree->Branch("covMCT","TMatrixD",&covMCT,64000,0);
   //tree->Branch("stREC",&stREC,"stREC[5]/F");
   tree->Branch("stREC","TMatrixD",&stREC,64000,0);
   //tree->Branch("covREC",&covREC,"covREC[25]/F");
   tree->Branch("covREC","TMatrixD",&covREC,64000,0);
   
   
   tree->Branch("chi2",&chi2,"chi2/F");
   tree->Branch("nfail",&nfail,"nfail/I");
   tree->Branch("ndf",&ndf,"ndf/I");
   tree->Branch("evtNo",&evtt,"evtNo/I");
   tree->Branch("chi2ndf",&chi2ndf,"chi2ndf/F");
 
   tree->Branch("trkNo",&nTrks,"trkNo/I");
   tree->Branch("ptsNo",&fptsNo,"ptsNo/I");
   tree->Branch("shx",fshx,"shx[ptsNo]/F");
   tree->Branch("shy",fshy,"shy[ptsNo]/F");
   tree->Branch("shz",fshz,"shz[ptsNo]/F");
 
   tree->Branch("pMCT",fpMCT,"pMCT[4]/F");
   tree->Branch("pRECKalF",fpREC,"pRECKalF[4]/F");
   tree->Branch("pRECKalL",fpRECL,"pRECKalL[4]/F");
   tree->Branch("pRECt3D",fpRECt3D,"pRECt3D[4]/F");
   
 
   //TGeoManager* geomGENFIT = new TGeoManager("Geometry", "Geane geometry");
   //TGeoManager::Import("config/genfitGeom.root");
   //  gROOT->Macro("config/Geane.C"); 
  
 }

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 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::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

   {
     return current_context_.get();
   }

void trkf::Track3DKalman::endJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 257 of file Track3DKalman_module.cc.

References rep, and repMC.

 {
   if (!rep) delete rep;
   if (!repMC) delete repMC;
 }

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

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

   {
     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 B, and 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 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::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

     {
       return true;
     }

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

staticinherited

Definition at line 76 of file Consumer.cc.

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

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

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

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

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

void 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 265 of file Track3DKalman_module.cc.

References genf::GFTrack::addHit(), art::PtrVector< T >::begin(), chi2, chi2ndf, recob::tracking::convertCollToPoint(), recob::tracking::convertCollToVector(), 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(), 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, and GFException::what().

 { 
 
   rep=0;
   repMC=0;
 
   // get services
   art::ServiceHandle<geo::Geometry> geom;
   // get the random number generator service and make some CLHEP generators
   art::ServiceHandle<art::RandomNumberGenerator> rng;
   CLHEP::HepRandomEngine &engine = rng->getEngine();
   CLHEP::RandGaussQ gauss(engine);
 
   // Make a std::unique_ptr<> for the thing you want to put into the event
   // because that handles the memory management for you
   std::unique_ptr<std::vector<recob::Track> > tcol(new std::vector<recob::Track>);
   std::unique_ptr< std::vector<recob::SpacePoint> >              spcol  (new std::vector<recob::SpacePoint>);
   std::unique_ptr< art::Assns<recob::Track, recob::SpacePoint> > tspassn(new art::Assns<recob::Track, recob::SpacePoint>);
   std::unique_ptr< art::Assns<recob::Track, recob::Cluster> >    tcassn (new art::Assns<recob::Track, recob::Cluster>);
   std::unique_ptr< art::Assns<recob::Track, recob::Hit> >        thassn (new art::Assns<recob::Track, recob::Hit>);
 
   // define TPC parameters
   TString tpcName = geom->GetLArTPCVolumeName();
 
 
   // get input Hit object(s).
   art::Handle< std::vector<recob::Track> > trackListHandle;
   evt.getByLabel(fSpacePtsModuleLabel,trackListHandle);
 
   art::PtrVector<simb::MCTruth> mclist;
 
   if (!evt.isRealData())
     {
 
       //      std::cout << "Track3DKalman: This is MC." << std::endl;
       // std::cout<<"Run "<<evt.run()<<" Event "<<evt.id().event()<<std::endl;
 
       art::Handle< std::vector<simb::MCTruth> > mctruthListHandle;
       evt.getByLabel(fGenieGenModuleLabel,mctruthListHandle);
 
       for (unsigned int ii = 0; ii <  mctruthListHandle->size(); ++ii)
         {
           art::Ptr<simb::MCTruth> mctparticle(mctruthListHandle,ii);
           mclist.push_back(mctparticle);
         }
 
     }
 
   //create collection of spacepoints that will be used when creating the Track object
   std::vector< art::Ptr<recob::SpacePoint> > spacepoints;
   art::PtrVector<recob::Track> trackIn;
   // std::cout<<"Run "<<evt.run()<<" Event "<<evt.id().event()<<std::endl;
   mf::LogInfo("Track3DKalman: ") << "There are " <<  trackListHandle->size() << " Track3Dreco/SpacePt tracks/groups (whichever) in this event.";
 
   art::FindManyP<recob::SpacePoint> fmsp(trackListHandle, evt, fSpacePtsModuleLabel);
   art::FindManyP<recob::Cluster>    fmc (trackListHandle, evt, fSpacePtsModuleLabel);
   art::FindManyP<recob::Hit>        fmh (trackListHandle, evt, fSpacePtsModuleLabel);
 
   for(unsigned int ii = 0; ii < trackListHandle->size(); ++ii)
     {
       art::Ptr<recob::Track> track(trackListHandle, ii);
       trackIn.push_back(track);
       
     }
       
       TVector3 MCOrigin;
       TVector3 MCMomentum;
       // TVector3 posErr(.05,.05,.05); // resolution. 0.5mm
       // TVector3 momErr(.1,.1,0.2);   // GeV
       TVector3 posErr(fPosErr[0],fPosErr[1],fPosErr[2]); // resolution. 0.5mm
       TVector3 momErr(fMomErr[0],fMomErr[1],fMomErr[2]);   // GeV
 
       // This is strictly for MC
       if (!evt.isRealData())
         {
           // Below breaks are stupid, I realize. But rather than keep all the MC
           // particles I just take the first primary, e.g., muon and only keep its
           // info in the branches of the Ttree. I could generalize later, ...
           for( unsigned int ii = 0; ii < mclist.size(); ++ii )
             {
             //art::Ptr<const simb::MCTruth> mc(mctruthListHandle,i);
             art::Ptr<simb::MCTruth> mc(mclist[ii]);
             for(int jj = 0; jj < mc->NParticles(); ++jj)
               {
                 simb::MCParticle part(mc->GetParticle(jj));
                 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 " ;
                 MCOrigin.SetXYZ(part.Vx(),part.Vy(),part.Vz()); // V for Vertex
                 MCMomentum.SetXYZ(part.Px(),part.Py(),part.Pz());
                 MCOrigin.Print();
                 MCMomentum.Print();
                 repMC = new genf::RKTrackRep(MCOrigin,
                                              MCMomentum,
                                              posErr,
                                              momErr,
                                              part.PdgCode());  
                 break;
               }
             break;
           }
           //for saving of MC truth    
           stMCT->ResizeTo(repMC->getState());
           *stMCT = repMC->getState();
           covMCT-> ResizeTo(repMC->getCov());
           *covMCT = repMC->getCov();
           mf::LogInfo("Track3DKalman: ") <<" repMC, covMC are ... " ;
           repMC->getState().Print();
           repMC->getCov().Print();
 
         } // !isRealData
       
       art::PtrVector<recob::Track>::const_iterator trackIter = trackIn.begin();
       
       nTrks=0;
       while(trackIter!=trackIn.end()) {
         spacepoints.clear();
         spacepoints = fmsp.at(nTrks);
         
         mf::LogInfo("Track3DKalman: ") << "found "<< spacepoints.size() <<" 3D spacepoint(s).";
         
         // Add the 3D track to the vector of the reconstructed tracks
         if(spacepoints.size()>0){
           // Insert the GENFIT/Kalman stuff here then make the tracks. Units are cm, GeV.
           const double resolution = 0.5; // dunno, 5 mm
           const int numIT = 3; // 3->1, EC, 6-Jan-2011. Back, 7-Jan-2011.
           
           
           //TVector3 mom(0.0,0.0,2.0);
           TVector3 mom(fMomStart[0],fMomStart[1],fMomStart[2]);
           //mom.SetMag(1.);
           TVector3 momM(mom);
           momM.SetX(gauss.fire(momM.X(),momErr.X()/* *momM.X() */));
           momM.SetY(gauss.fire(momM.Y(),momErr.Y()/* *momM.Y() */));
           momM.SetZ(gauss.fire(momM.Z(),momErr.Z()/* *momM.Z() */));
           //std::cout << "Track3DKalman: sort spacepoints by z
           
           std::sort(spacepoints.begin(), spacepoints.end(), sp_sort_3dz);
           
           //std::sort(spacepoints.begin(), spacepoints.end(), sp_sort_3dx); // Reverse sort!
           
           genf::GFFieldManager::getInstance()->init(new genf::GFConstField(0.,0.,0.0));
           genf::GFDetPlane planeG((TVector3)(spacepoints[0]->XYZ()),momM);
           
           
           //      std::cout<<"Track3DKalman about to do GAbsTrackRep."<<std::endl;
           // Initialize with 1st spacepoint location and a guess at the momentum.
           rep = new genf::RKTrackRep(//posM-.5/momM.Mag()*momM,
                                      (TVector3)(spacepoints[0]->XYZ()),
                                      momM,
                                      posErr,
                                      momErr,
                                      13);  // mu- hypothesis
           //      std::cout<<"Track3DKalman: about to do GFTrack. repDim is " << rep->getDim() <<std::endl;
           
           
           genf::GFTrack fitTrack(rep);//initialized with smeared rep
           // Gonna sort in x cuz I want to essentially transform here to volTPC coords.
           // volTPC coords, cuz that's what the Geant3/Geane stepper wants, as that's its understanding
           // from the Geant4 geometry, which it'll use. EC, 7-Jan-2011.
           int ihit = 0;
           fptsNo = 0;
           for (unsigned int point=0;point<spacepoints.size();++point)
             {
               
               TVector3 spt3(spacepoints[point]->XYZ());
                 if (point > 0 )
                   {
                     double eps (0.1);
                     TVector3 magNew(spt3[0],spt3[1],spt3[2]);
                     TVector3 magLast(spacepoints.back()->XYZ()[0],
                                      spacepoints.back()->XYZ()[1],
                                      spacepoints.back()->XYZ()[2]
                                      );
                     if (!(magNew.Mag()>=magLast.Mag()+eps || 
                           magNew.Mag()<=magLast.Mag()-eps)
                         )
                       continue;
                   }
 
               if (point%20) // Jump out of loop except on every 20th pt.
                 {
                   //continue;
                   // Icarus paper suggests we may wanna decimate our data in order to give
                   // trackfitter a better idea of multiple-scattering. EC, 7-Jan-2011.
                   //if (std::abs(spt3[0]-spacepoints.at(point-1).XYZ()[0]) < 2) continue;
                 }
               if (fptsNo<fDimSize)
                 {
                   fshx[fptsNo] = spt3[0];
                   fshy[fptsNo] = spt3[1];
                   fshz[fptsNo] = spt3[2];
                 }
                 fptsNo++;
                 
                 
                 LOG_DEBUG("Track3DKalman: ") << "ihit xyz..." << spt3[0]<<","<< spt3[1]<<","<< spt3[2];
                 fitTrack.addHit(new genf::PointHit(spt3,resolution),
                                 1,//dummy detector id
                                 ihit++
                                 );
             }
           
           //      std::cout<<"Track3DKalman about to do GFKalman."<<std::endl;
           genf::GFKalman k;
           //k.setBlowUpFactor(500); // Instead of 500 out of box. EC, 6-Jan-2011.
           //k.setInitialDirection(+1); // Instead of 1 out of box. EC, 6-Jan-2011.
           k.setNumIterations(numIT);
           //      std::cout<<"Track3DKalman back from setNumIterations."<<std::endl;
           try{
             //  std::cout<<"Track3DKalman about to processTrack."<<std::endl;
             k.processTrack(&fitTrack);
             //std::cout<<"Track3DKalman back from processTrack."<<std::endl;
           }
           catch(GFException& e){
             mf::LogError("Track3DKalman: ") << "just caught a GFException."<<std::endl;
             e.what();
             mf::LogError("Track3DKalman: ") << "Exceptions won't be further handled ->exit(1) "<<__LINE__;
             
             //  exit(1);
           }
           
           if(rep->getStatusFlag()==0) // 0 is successful completion
             {
               LOG_DEBUG("Track3DKalman: ") << __FILE__ << " " << __LINE__ ;
               LOG_DEBUG("Track3DKalman: ") << "Track3DKalman.cxx: Original plane:";
               
               if(fGenfPRINT) planeG.Print();
               LOG_DEBUG("Track3DKalman: ") << "Current (fit) reference Plane:";
               if(fGenfPRINT) rep->getReferencePlane().Print();
               
               LOG_DEBUG("Track3DKalman: ") << "Track3DKalman.cxx: Last reference Plane:";
               if(fGenfPRINT) rep->getLastPlane().Print();
               
               if(fGenfPRINT) 
                 {
                   if(planeG!=rep->getReferencePlane()) 
                     LOG_DEBUG("Track3DKalman: ") <<"Track3DKalman: Original hit plane (not surprisingly) not current reference Plane!"<<std::endl;
                 }
               
               stREC->ResizeTo(rep->getState());
               *stREC = rep->getState();
               covREC->ResizeTo(rep->getCov());
               *covREC = rep->getCov();
               if(fGenfPRINT)
                 {
                   LOG_DEBUG("Track3DKalman: ") << " Final State and Cov:";
                   stREC->Print();
                   covREC->Print();
                 }
               chi2 = rep->getChiSqu();
               ndf = rep->getNDF();
               nfail = fitTrack.getFailedHits();
               chi2ndf = chi2/ndf;
               TVector3 dircoss = (*trackIter)->VertexDirection<TVector3>();
               
               for (int ii=0;ii<3;++ii)
                 {
                   fpMCT[ii] = MCMomentum[ii]/MCMomentum.Mag();
                   fpREC[ii] = rep->getReferencePlane().getNormal()[ii];
                   fpRECL[ii] = rep->getLastPlane().getNormal()[ii];
                   fpRECt3D[ii] = dircoss[ii];
                 }
               fpMCT[3] = MCMomentum.Mag();
               fpREC[3] = -1.0/(*stREC)[0][0];
               
               evtt = (unsigned int) evt.id().event();
               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] << ".";
               
               tree->Fill();
               
               // Get the clusters associated to each track in induction and collection view
               std::vector< art::Ptr<recob::Cluster> > clusters = fmc.at(nTrks);
               std::vector< art::Ptr<recob::Hit> > hits = fmh.at(nTrks);
 
               // Use new Track constructor... EC, 21-Apr-2011.  
               std::vector<TVector3> dircos(spacepoints.size());
               dircos[0] = TVector3(fpREC);
               dircos.back() = TVector3(fpRECL);
               
               // fill a vector of TVector3 with the space points used to make this track
               std::vector<TVector3> xyz(spacepoints.size());
               size_t spStart = spcol->size();
               for(size_t tv = 0; tv < spacepoints.size(); ++tv){
                 xyz[tv] = TVector3(spacepoints[tv]->XYZ());
                 spcol->push_back(*(spacepoints[tv].get()));
               }
               size_t spEnd = spcol->size();
 
               tcol->push_back(recob::Track(recob::TrackTrajectory(recob::tracking::convertCollToPoint(xyz),
                                                                   recob::tracking::convertCollToVector(dircos),
                                                                   recob::Track::Flags_t(xyz.size()), false),
                                            0, -1., 0, recob::tracking::SMatrixSym55(), recob::tracking::SMatrixSym55(), tcol->size()-1));
 
               // associate the track with its clusters and tracks
               util::CreateAssn(*this, evt, *tcol, clusters, *tcassn);
               util::CreateAssn(*this, evt, *tcol, hits,     *thassn);
               
               // associate the track to the space points
               util::CreateAssn(*this, evt, *tcol, *spcol, *tspassn, spStart, spEnd);
 
             } // getStatusFlag 
         } // spacepoints.size()>0
         
         //
         //std::cout<<"Track3DKalman found "<< tcol->size() <<" 3D track(s)"<<std::endl;
         if(trackIter!=trackIn.end()) trackIter++;
         
         nTrks++;
         
       } // end loop over Track3Dreco/SpacePt tracks/groups (whichever) we brought into this event.
       
       evt.put(std::move(tcol));
       evt.put(std::move(spcol));
       evt.put(std::move(tcassn));
       evt.put(std::move(thassn));
       evt.put(std::move(tspassn));
 
 } // end method

void trkf::Track3DKalman::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 164 of file Track3DKalman_module.cc.

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

Referenced by Track3DKalman().

 {
 
   fSpacePtsModuleLabel   = pset.get< std::string >("SpacePtsModuleLabel");
   fGenieGenModuleLabel   = pset.get< std::string >("GenieGenModuleLabel");
   fG4ModuleLabel         = pset.get< std::string >("G4ModuleLabel");
   
   fPosErr                = pset.get< std::vector < double >  >("PosErr3");   // resolution. cm
   fMomErr                = pset.get< std::vector < double >  >("MomErr3");   // GeV
   fMomStart              = pset.get< std::vector < double >  >("MomStart3"); // Will be unit norm'd.
   fGenfPRINT             = pset.get< bool >("GenfPRINT");
   
 }

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

Member Data Documentation

Float_t trkf::Track3DKalman::chi2

private

Definition at line 109 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t trkf::Track3DKalman::chi2ndf

private

Definition at line 110 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

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

private

Definition at line 106 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

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

private

Definition at line 108 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

unsigned int trkf::Track3DKalman::evtt

private

Definition at line 118 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

unsigned int trkf::Track3DKalman::fDimSize

private

Definition at line 124 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

std::string trkf::Track3DKalman::fG4ModuleLabel

private

Definition at line 99 of file Track3DKalman_module.cc.

Referenced by reconfigure().

bool trkf::Track3DKalman::fGenfPRINT

private

Definition at line 100 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

std::string trkf::Track3DKalman::fGenieGenModuleLabel

private

Definition at line 98 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

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

private

Definition at line 127 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

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

private

Definition at line 128 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

Float_t* trkf::Track3DKalman::fpMCT

private

Definition at line 115 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

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

private

Definition at line 126 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

Float_t* trkf::Track3DKalman::fpREC

private

Definition at line 114 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fpRECL

private

Definition at line 113 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fpRECt3D

private

Definition at line 112 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

unsigned int trkf::Track3DKalman::fptsNo

private

Definition at line 120 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fshx

private

Definition at line 121 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fshy

private

Definition at line 122 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

Float_t* trkf::Track3DKalman::fshz

private

Definition at line 123 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

std::string trkf::Track3DKalman::fSpacePtsModuleLabel

private

Definition at line 97 of file Track3DKalman_module.cc.

Referenced by produce(), and reconfigure().

int trkf::Track3DKalman::ndf

private

Definition at line 117 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

int trkf::Track3DKalman::nfail

private

Definition at line 116 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

unsigned int trkf::Track3DKalman::nTrks

private

Definition at line 119 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

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

private

Definition at line 130 of file Track3DKalman_module.cc.

Referenced by endJob(), and produce().

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

private

Definition at line 129 of file Track3DKalman_module.cc.

Referenced by endJob(), and produce().

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

private

Definition at line 105 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

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

private

Definition at line 107 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

TTree* trkf::Track3DKalman::tree

private

Definition at line 103 of file Track3DKalman_module.cc.

Referenced by beginJob(), and produce().

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

/cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/TrackFinder/Track3DKalman_module.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation