Anp::RunkNN Class Reference

#include <RunkNN.h>

Inheritance diagram for Anp::RunkNN:

List of all members.

Public Types
typedef std::map< std::string, AlgThread * >	ThreadMap
Public Member Functions
	RunkNN ()
virtual	~RunkNN ()
bool	Run (Record &record)
void	Config (const Registry &reg)
bool	Init (const Header &header)
void	End (const DataBlock &)
Private Member Functions
bool	SaveParts (const Header &header)
bool	MakeNames (const Header &header)
void	MakeAlg (const Header &header)
void	MakePid (const Header &header, const std::string &file, int abase)
void	MakeKin (const Header &header, const std::string &file, int abase)
bool	Alg (const std::string &name, const Registry &reg, int base, int k1, int k2, int k3, int k4)
bool	Add (const std::string &name, int abase)
void	AddData (Record &record)
void	CutData (Record &record)
Private Attributes
bool	fCreateNew
bool	fEraseData
bool	fMakeAllBeam
bool	fMultiThread
bool	fPrint
int	fMuonKey
int	fKinBase
int	fPidBase
double	fComCut
double	fPurCut
double	fQoPCut
Registry	fConfig
std::string	fPath
std::map< std::string, int >	fFile
std::string	PREFIX
std::string	DETECTOR
std::string	MCRLTYPE
std::string	RECOTYPE
std::string	BEAMTYPE
std::string	fSelect
ThreadMap	fMap

---

Member Typedef Documentation

typedef std::map<std::string, AlgThread *> Anp::RunkNN::ThreadMap

Definition at line 47 of file RunkNN.h.

---

Constructor & Destructor Documentation

Anp::RunkNN::RunkNN (   )

Definition at line 26 of file RunkNN.cxx. :fCreateNew(false), fEraseData(true), fMakeAllBeam(false), fMultiThread(false), fPrint(false), fMuonKey(7000), fKinBase(1100), fPidBase(4000), fComCut(0.90), fPurCut(0.90), fQoPCut(-0.005), fConfig(false), fPath(""), fFile(), PREFIX("prototype"), DETECTOR(""), MCRLTYPE(""), RECOTYPE(""), BEAMTYPE(""), fSelect("physics") { }

Anp::RunkNN::~RunkNN (   )  [virtual]

Definition at line 51 of file RunkNN.cxx. { }

---

Member Function Documentation

bool Anp::RunkNN::Add (  const std::string &  name, int  abase )  [private]

Definition at line 623 of file RunkNN.cxx. References fFile. Referenced by MakeNames(). { // // If file name and key are unique then add values // if(fFile.find(name) != fFile.end()) { cerr << "RunkNN::Add - file name already exists: " << name << endl; return false; } for(map<string, int>::const_iterator it = fFile.begin(); it != fFile.end(); ++it) { if(it -> second == kbase) { cerr << "RunkNN::Add - key already exists: " << kbase << endl; return false; } } fFile[name] = kbase; return true; }

void Anp::RunkNN::AddData (  Record &  record  )  [private]

Definition at line 650 of file RunkNN.cxx. References Anp::Event::Add(), Anp::RecoNu::CosN(), Anp::Record::EventBegIterator(), Anp::Record::EventEndIterator(), Anp::EventIterator, Lit::Find(), Anp::Record::FindTruth(), fKinBase, fMuonKey, fPurCut, Anp::Record::GetHeader(), Anp::Event::GetNu(), Anp::Header::IsData(), Anp::LongestTrack(), Anp::RecoNu::PLep(), Anp::RecoNu::Q2(), Anp::Record::TrackEnd(), Anp::TrackIter, Anp::Record::TruthBeg(), Anp::Record::TruthEnd(), Anp::TruthIter, Anp::TruthRecoIter, Anp::RecoNu::W2(), Anp::RecoNu::X(), and Anp::RecoNu::Y(). Referenced by Run(). { // // Add temporary data required for kinematic kNN selection // for(EventIterator ievent = record.EventBegIterator(); ievent != record.EventEndIterator(); ++ievent) { Event &event = *ievent; // // Find longest reconstructed track // TrackIter itrack = Anp::LongestTrack(event, record); // // Event must have a track // if(itrack == record.TrackEnd()) continue; // // Track must pass fit // if(!(itrack -> PassFit())) continue; // // Add data keys // const double cosn = event.GetNu().CosN(); if(cosn < -1.0 || cosn > +1.0) { continue; } else { event.Add(fKinBase + 10, std::acos(cosn)); } event.Add(fKinBase + 4, event.GetNu().PLep()); event.Add(fKinBase + 5, event.GetNu().Q2()); event.Add(fKinBase + 6, event.GetNu().Y()); event.Add(fKinBase + 7, event.GetNu().W2()); event.Add(fKinBase + 8, event.GetNu().X()); event.Add(fKinBase + 9, cosn); // // Pass this point only if running on Monte-Carlo with truth records filled // if(record.GetHeader().IsData() || record.TruthBeg() == record.TruthEnd()) continue; const TruthIter itruth = record.FindTruth(event); if(itruth == record.TruthEnd()) { cerr << "RunkNN::AddData - failed to find Truth" << endl; continue; } TruthRecoIter ireco = itruth -> Find(event); if(ireco == itruth -> RecoEnd()) { cerr << "RunkNN::AddData - failed to find reconstructed truth info" << endl; continue; } else { // // Add purity and completeness variables for debugging // event.Add(fKinBase + 2, ireco -> Purity()); event.Add(fKinBase + 3, ireco -> ComplAll()); } // // Use tracks with reconstructed negative curvature // if(itrack -> QP() > fQoPCut) continue; // // Use true muon neutrino interactions // if(itruth -> Particle() != Particle::kNuM) continue; // // Use true QES, RES or DIS charged current events // if(itruth -> Interaction() != 1 || (itruth -> Resonance() != 1001 && itruth -> Resonance() != 1002 && itruth -> Resonance() != 1003)) continue; // // Use well reconstructed events // if(ireco -> Purity() < fPurCut || ireco -> ComplAll() < fComCut) continue; // // Find muon truth key which is filled by FillMuonId // if(!(itrack -> KeyExists(fMuonKey))) { cerr << "RunkNN::AddData - failed to find muon truth key: " << fMuonKey << endl; continue; } // // Reconstucted track must match true muon !!!hardcoded numbers!!! // if(itrack -> DataAt(fMuonKey) < 0.999 || itrack -> DataAt(fMuonKey) > 1.001) { continue; } // // Add resonance code for QES, RES or DIS and use this event for training // event.Add(fKinBase + 1, float(itruth -> Resonance() % 1000)); } }

bool Anp::RunkNN::Alg (  const std::string &  name, const Registry &  reg, int  base, int  k1, int  k2, int  k3, int  k4 )  [private]

Definition at line 580 of file RunkNN.cxx. References base, Config(), fMap, reg, Registry::Set(), and Anp::Handle< T >::valid(). Referenced by MakeKin(), and MakePid(). { // // Create one FillkNN algorithm // Registry kreg(reg); stringstream keys; if(k1 > 0) keys << std::setw(5) << std::setfill(' ') << std::right << k1 << " "; if(k2 > 0) keys << std::setw(5) << std::setfill(' ') << std::right << k2 << " "; if(k3 > 0) keys << std::setw(5) << std::setfill(' ') << std::right << k3 << " "; if(k4 > 0) keys << std::setw(5) << std::setfill(' ') << std::right << k4 << " "; kreg.Set("FillkNNKeyBase", base); kreg.Set("FillkNNKeyList", keys.str().c_str()); const string name = file + ":" + keys.str(); if(fMap.find(name) != fMap.end()) { cerr << "RunkNN::Alg - name already exists: " << name << endl; return false; } Handle<AlgSnarl> alg = Factory<AlgSnarl>::Instance().Create("FillkNN"); if(!alg.valid()) { cerr << "RunkNN::Alg - failed to create FillkNN algorithm!" << endl; return false; } // // Configure algorithm - this is the only direct access // to FillkNN algorithm from this class. // alg -> Config(kreg); return fMap.insert(ThreadMap::value_type(name, new AlgThread(alg, name))).second; }

void Anp::RunkNN::Config (  const Registry &  reg  )  [virtual]

Reimplemented from Anp::AlgSnarl. Definition at line 171 of file RunkNN.cxx. References fComCut, fConfig, fCreateNew, fEraseData, fKinBase, fMakeAllBeam, fMultiThread, fMuonKey, fPath, fPidBase, fPrint, fPurCut, fQoPCut, fSelect, Registry::Get(), Registry::KeyExists(), Registry::Merge(), PREFIX, Anp::Read(), and reg. Referenced by Alg(). { // // Configure self and save Registry for children // Anp::Read(reg, "RunkNNCreateNew", fCreateNew); Anp::Read(reg, "RunkNNEraseData", fEraseData); Anp::Read(reg, "RunkNNMultiThread", fMultiThread); Anp::Read(reg, "RunkNNMakeAllBeam", fMakeAllBeam); Anp::Read(reg, "RunkNNPrint", fPrint); // // Do not multi-thread when making new knn prototype(s) // if(fCreateNew) fMultiThread = false; const char *value_path = 0; if(reg.Get("RunkNNPath", value_path) && value_path) { fPath = value_path; } if(!fPath.empty() && fPath[fPath.size()-1] != '/') fPath += "/"; const char *value_pref = 0; if(reg.Get("RunkNNPrefix", value_pref) && value_pref) { PREFIX = value_pref; } const char *value_selc = 0; if(reg.Get("RunkNNSelect", value_selc) && value_selc) { fSelect = value_selc; } reg.Get("RunkNNComCut", fComCut); reg.Get("RunkNNPurCut", fPurCut); reg.Get("RunkNNQoPCut", fQoPCut); if(fSelect=="QPall") fQoPCut = 10000; if(reg.KeyExists("PrintConfig")) { cout << "RunkNN::Config" << endl << " CreateNew = " << fCreateNew << endl << " EraseData = " << fEraseData << endl << " MultiThread = " << fMultiThread << endl << " MuonKey = " << fMuonKey << endl << " KinBase = " << fKinBase << endl << " PidBase = " << fPidBase << endl << " ComCut = " << fComCut << endl << " PurCut = " << fPurCut << endl << " QoPCut = " << fQoPCut << endl << " Print = " << fPrint << endl << " Path = " << fPath << endl << " PREFIX = " << PREFIX << endl << " Select = " << fSelect << endl; } // // Configure generic Registry for kNN // fConfig.Merge(reg); }

void Anp::RunkNN::CutData (  Record &  record  )  [private]

Definition at line 769 of file RunkNN.cxx. References Anp::Record::EventBegIterator(), Anp::Record::EventEndIterator(), Anp::EventIterator, and fKinBase. Referenced by Run(). { // // Erase temporary data required for kinematic kNN selection // vector<short> kvec; kvec.push_back(fKinBase + 1); kvec.push_back(fKinBase + 2); kvec.push_back(fKinBase + 3); kvec.push_back(fKinBase + 4); kvec.push_back(fKinBase + 5); kvec.push_back(fKinBase + 6); kvec.push_back(fKinBase + 7); kvec.push_back(fKinBase + 8); kvec.push_back(fKinBase + 9); kvec.push_back(fKinBase + 10); for(EventIterator ievent = record.EventBegIterator(); ievent != record.EventEndIterator(); ++ievent) { ievent -> Erase(std::remove_if(ievent -> DataBegIterator(), ievent -> DataEndIterator(), Anp::PrevDataKey(kvec)), ievent -> DataEndIterator()); } }

void Anp::RunkNN::End (  const DataBlock &   )  [virtual]

Reimplemented from Anp::AlgSnarl. Definition at line 258 of file RunkNN.cxx. References fMap. { // // Destroy AlgThread pointers // for(ThreadMap::iterator it = fMap.begin(); it != fMap.end(); ++it) { AlgThread *alg = it -> second; assert(alg && "logic error"); cout << " " << it -> first << endl; alg -> End(block); delete alg; } fMap.clear(); }

bool Anp::RunkNN::Init (  const Header &  header  )  [virtual]

Reimplemented from Anp::AlgSnarl. Definition at line 237 of file RunkNN.cxx. References fMap, MakeAlg(), MakeNames(), and SaveParts(). { // // Initiliaze self and AlgThread algorithm(s) // assert(fMap.empty() && "logic error - Init() function is probably called twice"); SaveParts(header); MakeNames(header); MakeAlg(header); for(ThreadMap::iterator ialg = fMap.begin(); ialg != fMap.end(); ++ialg) { ialg -> second -> Init(header); } return true; }

void Anp::RunkNN::MakeAlg (  const Header &  header  )  [private]

Definition at line 397 of file RunkNN.cxx. References fFile, fPath, infile, MakeKin(), and MakePid(). Referenced by Init(). { // // Make kNN algorithms... lots of hardcoded numbers and decisions... // if(fFile.empty()) { cerr << "RunkNN::MakeAlg - no kNN file names available" << endl; return; } for(map<string, int>::const_iterator ifile = fFile.begin(); ifile != fFile.end(); ++ifile) { // // If not creating new kNN algorithm then check if prototype file exists // const string file = fPath + ifile -> first; ifstream infile(file.c_str(), std::ios_base::in); if(fCreateNew) { if(infile.is_open()) { cerr << "RunkNN::MakeAlg - file already exists:\n " << file << endl; continue; } } else { if(!infile.is_open()) { cerr << "RunkNN::MakeAlg - file does not exist:\n " << file << endl; continue; } } // // Make FillkNN algorithms for current file and key // MakePid(header, ifile -> first, ifile -> second); MakeKin(header, ifile -> first, ifile -> second); } }

void Anp::RunkNN::MakeKin (  const Header &  header, const std::string &  file, int  abase )  [private]

Definition at line 535 of file RunkNN.cxx. References Alg(), fConfig, fKinBase, fPath, fSelect, reg, and Registry::Set(). Referenced by MakeAlg(). { // // Lots of hardcoded numbers, will overwrite some FillkNN keys // Registry reg(fConfig); reg.Set("FillkNNFilePath", (fPath + file).c_str()); reg.Set("FillkNNKeyTruth", int(fKinBase + 1)); reg.Set("FillkNNKNeighbor", int(122)); reg.Set("FillkNNKNeighborMod", int(40)); reg.Set("FillkNNAddAll", "yes"); reg.Set("FillkNNAddEvent", "yes"); reg.Set("FillkNNUseTrack", "no"); reg.Set("FillkNNTrim", "yes"); reg.Set("FillkNNTrimDelta", int(1000)); if(fCreateNew) { reg.Set("FillkNN", "yes"); } else { reg.Set("FillkNN", "no"); } // // Input variables for 1st selection: // if(fCreateNew) { Alg(file, reg, 0, fKinBase+5, fKinBase+6, fKinBase+ 7, -1); // Q^2, Y, W^2 Alg(file, reg, 0, fKinBase+4, fKinBase+6, fKinBase+10, -1); // P_{\mu}, Y, \theta_{\mu} } else { if(fSelect == "kinem") { Alg(file, reg, fKinBase + kbase, fKinBase+5, fKinBase+6, fKinBase+ 7, -1); Alg(file, reg, fKinBase + kbase + 300, fKinBase+4, fKinBase+6, fKinBase+10, -1); } } }

bool Anp::RunkNN::MakeNames (  const Header &  header  )  [private]

Definition at line 303 of file RunkNN.cxx. References Add(), BEAMTYPE, DETECTOR, fMakeAllBeam, fSelect, Anp::Header::IsData(), MCRLTYPE, PREFIX, and RECOTYPE. Referenced by Init(). { // // Make file names for kNN algorithms. This function has // some hardcoded names and decisions. // // File name is ideally created as follows: // // PREFIX.DETECTOR.MCRLTYPE.RECOTYPE.BEAMTYPE.root // // // Case when creating new prototype // if(fCreateNew) { // // Can not create prototype from data events // if(header.IsData()) return true; // // Create single file, base key is 0 because no keys should be computed // return Add(PREFIX+"."+DETECTOR+"."+MCRLTYPE+"."+RECOTYPE+"."+BEAMTYPE+".root", 0); } // // !!!KLUDGE!!! Only know how to deal with cedar-like releases // if(RECOTYPE.find("cedar") == string::npos && RECOTYPE.find("dogwood1")==string::npos) { cerr << "RunkNN::MakeNames - unknown reconstruction release: " << RECOTYPE << endl; return false; } // // !!!KLUDGE!!! Use 3 different Monte-Carlo daikon/cedar releases // if(fSelect == "physics" && !(RECOTYPE=="dogwood1")) { Add(PREFIX + "." + DETECTOR + ".daikon_04.cedar_phy_bhcurv.L010z185i.root", 0); Add(PREFIX + "." + DETECTOR + ".daikon_00.cedar_phy.L010z185i.root", 1000); } else if(fSelect=="physics") { Add(PREFIX + "." + DETECTOR + ".daikon_04.dogwood1.L010z185i.root", 0); } else if(RECOTYPE == "cedar_phy_bhcurv") { Add(PREFIX + "." + DETECTOR + ".daikon_04.cedar_phy_bhcurv.L010z185i.root", 0); } else if(RECOTYPE == "cedar_phy") { Add(PREFIX + "." + DETECTOR + ".daikon_00.cedar_phy.L010z185i.root", 1000); } else if(RECOTYPE == "dogwood1") { Add(PREFIX + "." + DETECTOR + ".daikon_04.dogwood1.L010z185i.root", 0); } else if(RECOTYPE == "dogwood1" && fSelect=="QPall") { Add(PREFIX + "." + DETECTOR + ".daikon_04.dogwood1.L010z185i.QPall.root", 1000); } return true; // // Native beam // if(BEAMTYPE != "L010z185i" && (RECOTYPE == "cedar" || RECOTYPE == "cedar_phy" || RECOTYPE == "cedar_phy_bhcurv" || RECOTYPE == "dogwood1" )) { Add(PREFIX+"."+DETECTOR+"."+MCRLTYPE+"."+RECOTYPE+"."+BEAMTYPE+".root", 10000); } // // Many beams combined // if(fMakeAllBeam && (RECOTYPE == "cedar" || RECOTYPE == "cedar_phy" || RECOTYPE == "cedar_phy_bhcurv" || RECOTYPE == "dogwood1" )) { Add(PREFIX+"."+DETECTOR+"."+MCRLTYPE+"."+RECOTYPE+"."+"ALLBEAM.root", 11000); } return true; }

void Anp::RunkNN::MakePid (  const Header &  header, const std::string &  file, int  abase )  [private]

Definition at line 443 of file RunkNN.cxx. References Alg(), fConfig, fPath, fPidBase, fSelect, reg, and Registry::Set(). Referenced by MakeAlg(). { // // Lots of hardcoded numbers, will overwrite some FillkNN keys // Registry reg(fConfig); reg.Set("FillkNNFilePath", (fPath + file).c_str()); reg.Set("FillkNNKeySignal", int(1)); // for historical reasons muon have type=1 reg.Set("FillkNNKeyTruth", int(7000)); // value for this key is set by FillMuonId reg.Set("FillkNNKNeighbor", int(82)); reg.Set("FillkNNKNeighborMod", int(20)); reg.Set("FillkNNAddEvent", "no"); reg.Set("FillkNNUseTrack", "yes"); reg.Set("FillkNNTrim", "yes"); reg.Set("FillkNNTrimDelta", int(1000)); if(fCreateNew) { reg.Set("FillkNN", "yes"); } else { reg.Set("FillkNN", "no"); } // // !!!KLUDGE!!! Multiple selections... // if(fCreateNew) { if(fSelect == "physics" ) { Alg(file, reg, 0, 7001, 7010, 7020, 7040); Alg(file, reg, 0, 7001, 7210, 7220, 7240); Alg(file, reg, 0, 7001, 7310, 7320, 7340); Alg(file, reg, 0, 7001, 7410, 7420, 7440); Alg(file, reg, 0, 11001, 11210, 11220, 11240); Alg(file, reg, 0, 11001, 11410, 11420, 11440); } else { Alg(file, reg, 0, 7001, 7010, 7020, 7040); Alg(file, reg, 0, 7001, 7210, 7220, 7240); Alg(file, reg, 0, 7001, 7310, 7320, 7340); Alg(file, reg, 0, 7001, 7410, 7420, 7440); Alg(file, reg, 0, 9001, 9210, 9220, 9240); Alg(file, reg, 0, 9001, 9410, 9420, 9440); Alg(file, reg, 0, 10001, 10210, 10220, 10240); Alg(file, reg, 0, 10001, 10410, 10420, 10440); Alg(file, reg, 0, 11001, 11210, 11220, 11240); Alg(file, reg, 0, 11001, 11410, 11420, 11440); Alg(file, reg, 0, 12001, 12210, 12220, 12240); Alg(file, reg, 0, 12001, 12410, 12420, 12440); } } else { if(fSelect == "physics" || fSelect=="QPall") { Alg(file, reg, fPidBase + kbase, 7001, 7010, 7020, 7040); Alg(file, reg, fPidBase + kbase + 200, 11001, 11410, 11420, 11440); } else if(fSelect == "study-sigall") { Alg(file, reg, fPidBase + kbase, 7001, 7010, 7020, 7040); Alg(file, reg, fPidBase + kbase + 100, 7001, 7210, 7220, 7240); Alg(file, reg, fPidBase + kbase + 200, 7001, 7310, 7320, 7340); Alg(file, reg, fPidBase + kbase + 300, 7001, 7410, 7420, 7440); } else if(fSelect == "study-sigcor") { Alg(file, reg, fPidBase + kbase, 7001, 7010, 7020, 7040); Alg(file, reg, fPidBase + kbase + 100, 7001, 7210, 7220, 7240); Alg(file, reg, fPidBase + kbase + 200, 9001, 9210, 9220, 9240); Alg(file, reg, fPidBase + kbase + 300, 10001, 10210, 10220, 10240); Alg(file, reg, fPidBase + kbase + 400, 11001, 11210, 11220, 11240); Alg(file, reg, fPidBase + kbase + 500, 12001, 12210, 12220, 12240); } else if(fSelect == "study-sigmip") { Alg(file, reg, fPidBase + kbase, 7001, 7310, 7320, 7340); Alg(file, reg, fPidBase + kbase + 100, 7001, 7410, 7420, 7440); Alg(file, reg, fPidBase + kbase + 200, 9001, 9410, 9420, 9440); Alg(file, reg, fPidBase + kbase + 300, 10001, 10410, 10420, 10440); Alg(file, reg, fPidBase + kbase + 400, 11001, 11410, 11420, 11440); Alg(file, reg, fPidBase + kbase + 500, 12001, 12410, 12420, 12440); } } }

bool Anp::RunkNN::Run (  Record &  record  )  [virtual]

Implements Anp::AlgSnarl. Definition at line 56 of file RunkNN.cxx. References AddData(), copy(), CutData(), Anp::Record::EventBeg(), Anp::Record::EventBegIterator(), Anp::Record::EventEnd(), Anp::Record::EventEndIterator(), Anp::EventIter, Anp::EventIterator, fMap, fMultiThread, Anp::Record::GetStdHep(), Anp::Record::GetStrip(), Anp::Record::TrackBeg(), Anp::Record::TrackBegIterator(), Anp::Record::TrackEnd(), Anp::Record::TrackEndIterator(), Anp::TrackIter, and Anp::TrackIterator. { // // Add data to record, run multiple kNN algorithm, the collect and save knn ouput // if(fPrint) { for(EventIter it = record.EventBeg(); it != record.EventEnd(); ++it) { cout << "beg: event index " << it -> EventIndex() << endl; std::copy(it -> DataBeg(), it -> DataEnd(), ostream_iterator<Data>(cout, "\n")); } for(TrackIter it = record.TrackBeg(); it != record.TrackEnd(); ++it) { cout << "beg: track index " << it -> TrackIndex() << endl; std::copy(it -> DataBeg(), it -> DataEnd(), ostream_iterator<Data>(cout, "\n")); } } // // Add data needed by kinematic kNN selection // AddData(record); // // Make lightweight copy of Record for AlgThread(s) // Record copy = record; copy.GetStrip().clear(); copy.GetStdHep().clear(); // // First pass, run all kNN algorithms // for(ThreadMap::iterator ialg = fMap.begin(); ialg != fMap.end(); ++ialg) { ialg -> second -> Run(copy, fMultiThread); } // // Second pass, collect data // for(ThreadMap::iterator ialg = fMap.begin(); ialg != fMap.end(); ++ialg) { AlgThread *alg = ialg -> second; if(fMultiThread && alg -> Join() != 0) { cerr << "RunkNN::Run - thread join function failed" << endl; continue; } // // Collect data from current AlgThread algorithm // map<int, PrevDataKey> event_map, track_map; for(EventIter it = alg -> GetRecord().EventBeg(); it != alg -> GetRecord().EventEnd(); ++it) { event_map[it -> EventIndex()] = PrevDataKey(it -> DataBeg(), it -> DataEnd()); } for(TrackIter it = alg -> GetRecord().TrackBeg(); it != alg -> GetRecord().TrackEnd(); ++it) { track_map[it -> TrackIndex()] = PrevDataKey(it -> DataBeg(), it -> DataEnd()); } // // Add data from current AlgThread algorithm to events and tracks // for(EventIterator it = record.EventBegIterator(); it != record.EventEndIterator(); ++it) { event_map[it -> EventIndex()](*it); } for(TrackIterator it = record.TrackBegIterator(); it != record.TrackEndIterator(); ++it) { track_map[it -> TrackIndex()](*it); } } // // Sort event and track data by keys // for(EventIterator it = record.EventBegIterator(); it != record.EventEndIterator(); ++it) { std::sort(it -> DataBegIterator(), it -> DataEndIterator()); } for(TrackIterator it = record.TrackBegIterator(); it != record.TrackEndIterator(); ++it) { std::sort(it -> DataBegIterator(), it -> DataEndIterator()); } // // Erase temporary data added by this algorithm // if(fEraseData) CutData(record); if(fPrint) { for(EventIter it = record.EventBeg(); it != record.EventEnd(); ++it) { cout << "end: event index " << it -> EventIndex() << endl; std::copy(it -> DataBeg(), it -> DataEnd(), ostream_iterator<Data>(cout, "\n")); } for(TrackIter it = record.TrackBeg(); it != record.TrackEnd(); ++it) { cout << "end: track index " << it -> TrackIndex() << endl; std::copy(it -> DataBeg(), it -> DataEnd(), ostream_iterator<Data>(cout, "\n")); } } return true; }

bool Anp::RunkNN::SaveParts (  const Header &  header  )  [private]

Definition at line 277 of file RunkNN.cxx. References BEAMTYPE, DETECTOR, Anp::Header::GetBEAMTYPE(), Anp::Header::GetDETECTOR(), Anp::Header::GetMCRLTYPE(), Anp::Header::GetRECOTYPE(), MCRLTYPE, and RECOTYPE. Referenced by Init(). { // // Collect parts of information for knn prototype file name. // DETECTOR = header.GetDETECTOR(); MCRLTYPE = header.GetMCRLTYPE(); RECOTYPE = header.GetRECOTYPE(); // // Do not specify BEAMTYPE if fMakeAllBeam is set to true. // if(fMakeAllBeam) { BEAMTYPE = "ALLBEAM"; } else { BEAMTYPE = header.GetBEAMTYPE(); } return false; }

---

Member Data Documentation

std::string Anp::RunkNN::BEAMTYPE [private]

Definition at line 90 of file RunkNN.h. Referenced by MakeNames(), and SaveParts().

std::string Anp::RunkNN::DETECTOR [private]

Definition at line 87 of file RunkNN.h. Referenced by MakeNames(), and SaveParts().

double Anp::RunkNN::fComCut [private]

Definition at line 77 of file RunkNN.h. Referenced by Config().

Registry Anp::RunkNN::fConfig [private]

Definition at line 81 of file RunkNN.h. Referenced by Config(), MakeKin(), and MakePid().

bool Anp::RunkNN::fCreateNew [private]

Definition at line 67 of file RunkNN.h. Referenced by Config().

bool Anp::RunkNN::fEraseData [private]

Definition at line 68 of file RunkNN.h. Referenced by Config().

std::map<std::string, int> Anp::RunkNN::fFile [private]

Definition at line 84 of file RunkNN.h. Referenced by Add(), and MakeAlg().

int Anp::RunkNN::fKinBase [private]

Definition at line 74 of file RunkNN.h. Referenced by AddData(), Config(), CutData(), and MakeKin().

bool Anp::RunkNN::fMakeAllBeam [private]

Definition at line 69 of file RunkNN.h. Referenced by Config(), and MakeNames().

ThreadMap Anp::RunkNN::fMap [private]

Definition at line 93 of file RunkNN.h. Referenced by Alg(), End(), Init(), and Run().

bool Anp::RunkNN::fMultiThread [private]

Definition at line 70 of file RunkNN.h. Referenced by Config(), and Run().

int Anp::RunkNN::fMuonKey [private]

Definition at line 73 of file RunkNN.h. Referenced by AddData(), and Config().

std::string Anp::RunkNN::fPath [private]

Definition at line 83 of file RunkNN.h. Referenced by Config(), MakeAlg(), MakeKin(), and MakePid().

int Anp::RunkNN::fPidBase [private]

Definition at line 75 of file RunkNN.h. Referenced by Config(), and MakePid().

bool Anp::RunkNN::fPrint [private]

Definition at line 71 of file RunkNN.h. Referenced by Config().

double Anp::RunkNN::fPurCut [private]

Definition at line 78 of file RunkNN.h. Referenced by AddData(), and Config().

double Anp::RunkNN::fQoPCut [private]

Definition at line 79 of file RunkNN.h. Referenced by Config().

std::string Anp::RunkNN::fSelect [private]

Definition at line 91 of file RunkNN.h. Referenced by Config(), MakeKin(), MakeNames(), and MakePid().

std::string Anp::RunkNN::MCRLTYPE [private]

Definition at line 88 of file RunkNN.h. Referenced by MakeNames(), and SaveParts().

std::string Anp::RunkNN::PREFIX [private]

Definition at line 86 of file RunkNN.h. Referenced by Config(), and MakeNames().

std::string Anp::RunkNN::RECOTYPE [private]

Definition at line 89 of file RunkNN.h. Referenced by MakeNames(), and SaveParts().

---

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

• RunkNN.h
• RunkNN.cxx

---