Managed::ClusterManager Class Reference

#include <ClusterManager.h>

List of all members.

Public Member Functions
	ClusterManager ()
	~ClusterManager ()
void	MakeClusters (double min_cluster_e=0.4, double min_strip_e=0.4, double max_t_skip=0.025)
void	Reset ()
void	AddStrip (int myplane, int mystrip, double myenergy, double st, double sz, int view)
void	LoadHits ()
std::vector< int >	GetPlanes ()
std::vector< int >	GetStrips (int plane)
Managed::ManagedCluster *	GetCluster (int cid)
Managed::ManagedCluster *	GetSavedCluster (int cid)
void	SplitEnergy (Managed::ManagedCluster *cluster, double energy, int matchenergy)
int	MergeClusters (int cid1, int cid2)
void	FillClusterMap (std::map< double, std::map< double, std::pair< double, int > > > *cluster_map)
std::vector< int >	FindClustersInZ (double z, int view)
std::map< double, std::map< double, int > > *	GetClusterMap (int view=0)
std::vector< int >	GetViewIndex (int view)
void	SetHitManager (Managed::HitManager *h)
HitManager *	GetHitManager ()
void	SetClusterSaver (Managed::ClusterSaver *cs)
Managed::ClusterSaver *	GetClusterSaver ()
void	DumpClusters ()
int	SaveCluster (int cluster_id, double energy_to_save=0, int status=0)
void	ClearInUse ()
void	MarkUsed (int id)
Public Attributes
std::map< double, std::map< double, int > >	cluster_map
std::map< double, std::map< double, int > >	cluster_map_u
std::map< double, std::map< double, int > >	cluster_map_v
std::vector< Managed::ManagedCluster >	clusters
double	maxz
double	minz
double	maxu
double	minu
double	maxv
double	minv
Private Member Functions
void	EraseCluster (int cid)
void	AdjustCluster (Managed::ManagedCluster *cluster)
	ClassDef (ClusterManager, 1)
void	RebuildClusterMaps ()
Private Attributes
int	clusters_are_made
std::map< int, std::map< int, int > >	loc_map
std::vector< ManagedHit >	hits
HitManager *	hitmanager
ClusterSaver *	clustersaver
int	needMapRebuild
std::vector< int >	inuse
std::vector< int >	clusters_to_delete

---

Constructor & Destructor Documentation

ClusterManager::ClusterManager (   )

Definition at line 13 of file ClusterManager.cxx. { clusters_are_made=0; loc_map.clear(); cluster_map.clear(); cluster_map_u.clear(); cluster_map_v.clear(); hits.clear(); clusters.clear(); clusters_to_delete.clear(); minu=100000; maxu=-100000; minv=100000; maxv=-100000; minz=100000; maxz=-100000; clustersaver=0; hitmanager=0; needMapRebuild=0; }

ClusterManager::~ClusterManager (   )

Definition at line 317 of file ClusterManager.cxx. {}

---

Member Function Documentation

void ClusterManager::AddStrip (  int  myplane, int  mystrip, double  myenergy, double  st, double  sz, int  view )

Definition at line 779 of file ClusterManager.cxx. References hits, loc_map, maxu, maxv, maxz, minu, minv, and minz. { ManagedHit h(view,myplane,mystrip,sz,st,myenergy); hits.push_back(h); loc_map[myplane][mystrip]=hits.size()-1; minz=minz < sz ? minz : sz; maxz=maxz > sz ? maxz : sz; if(view==2) { minu=minu < st ? minu : st; maxu=maxu > st ? maxu : st; }else if(view==3) { minv=minv < st ? minv : st; maxv=maxv > st ? maxv : st; } }

void ClusterManager::AdjustCluster (  Managed::ManagedCluster *  cluster  )  [private]

Definition at line 213 of file ClusterManager.cxx. References clusters_to_delete, Managed::ManagedCluster::e, Managed::ManagedCluster::Finalize(), Managed::HitManager::FindHit(), GetCluster(), Managed::ManagedCluster::hit_id, Managed::ManagedCluster::hite, hitmanager, Managed::ManagedCluster::id, MSG, and Managed::ManagedHit::SetEnergy(). Referenced by SaveCluster(). { if(!cluster)return; ManagedCluster * current_cluster = GetCluster(cluster->id); if(!current_cluster)return; MSG("ClusterManager",Msg::kDebug)<<"cluster e "<<cluster->e<<" curcluster "<<current_cluster->e<<"\n"; //is the passed object a part of the current cluster vector, or was it copied? //if its copied, make sure it gets copied back into the version in the vector // if(fabs(cluster->e - current_cluster->e)<1e-9) // return; //now we need to update the hits! int modified_cluster=0; for(unsigned int i=0;i<cluster->hit_id.size();i++) { ManagedHit * hit = hitmanager->FindHit(cluster->hit_id[i]); double setenergy = 0; if(hit)setenergy = hit->SetEnergy(cluster->hite[i]); MSG("ClusterManager",Msg::kDebug)<<"setting energy to "<<cluster->hite[i]<<" and it set "<<setenergy<<"\n"; if(fabs(setenergy-current_cluster->hite[i])>1e-9) { MSG("ClusterManager",Msg::kDebug)<<"modifying cluster\n"; current_cluster->hite[i]=setenergy; modified_cluster=1; } } if(modified_cluster)current_cluster->Finalize(); if(current_cluster->e<0.001)clusters_to_delete.push_back(current_cluster->id); }

Managed::ClusterManager::ClassDef (  ClusterManager  , 1  )  [private]

void Managed::ClusterManager::ClearInUse (   )  [inline]

Definition at line 74 of file ClusterManager.h. Referenced by PrimaryShowerFinder::MakeChains(). {inuse.clear();};

void ClusterManager::DumpClusters (   )

Definition at line 321 of file ClusterManager.cxx. References cluster_map, cluster_map_u, cluster_map_v, Managed::ManagedCluster::dt, Managed::ManagedCluster::dz, Managed::ManagedCluster::e, GetCluster(), MSG, Managed::ManagedCluster::t, Managed::ManagedCluster::view, and Managed::ManagedCluster::z. Referenced by Finder::Process(). { ostringstream os; std::map<double, std::map<double, int> >::iterator p_iter; std::map<double, int>::iterator s_iter; os<<"dumping clusters (all)...\n"; for(p_iter=cluster_map.begin();p_iter!=cluster_map.end(); p_iter++) { std::map<double, int>::iterator s_iter; for(s_iter=p_iter->second.begin();s_iter!=p_iter->second.end(); s_iter++) { ManagedCluster *mc = GetCluster(s_iter->second); if(!mc){ os<<"missing cluster id "<<s_iter->second<<"\n"; continue; } os<<"id "<<s_iter->second<<" z "<<mc->z<<" t "<<mc->t<<" e "<<mc->e<<" dz "<<mc->dz<<" dt "<<mc->dt<<" view "<<mc->view<<"\n"; } } os<<"dumping clusters (u)...\n"; for(p_iter=cluster_map_u.begin();p_iter!=cluster_map_u.end(); p_iter++) { std::map<double, int>::iterator s_iter; for(s_iter=p_iter->second.begin();s_iter!=p_iter->second.end(); s_iter++) { ManagedCluster *mc = GetCluster(s_iter->second); if(!mc){ os<<"missing cluster id "<<s_iter->second<<"\n"; continue; } os<<"z "<<mc->z<<" t "<<mc->t<<" e "<<mc->e<<" dz "<<mc->dz<<" dt "<<mc->dt<<" view "<<mc->view<<"\n"; os<<"id "<<s_iter->second<<"\n"; } } os<<"dumping clusters (v)...\n"; for(p_iter=cluster_map_v.begin();p_iter!=cluster_map_v.end(); p_iter++) { std::map<double, int>::iterator s_iter; for(s_iter=p_iter->second.begin();s_iter!=p_iter->second.end(); s_iter++) { ManagedCluster *mc = GetCluster(s_iter->second); if(!mc){ os<<"missing cluster id "<<s_iter->second<<"\n"; continue; } os<<"z "<<mc->z<<" t "<<mc->t<<" e "<<mc->e<<" dz "<<mc->dz<<" dt "<<mc->dt<<" view "<<mc->view<<"\n"; os<<"id "<<s_iter->second<<"\n"; } } os<<"done....\n\n"; MSG("ClusterManager",Msg::kInfo)<<os; }

void ClusterManager::EraseCluster (  int  cid  )  [private]

Definition at line 907 of file ClusterManager.cxx. References cluster_map, cluster_map_u, cluster_map_v, clusters, done(), and MSG. Referenced by MergeClusters(). { MSG("ClusterManager",Msg::kDebug)<<"erasing "<<cid<<"\n"; std::map<double, std::map<double, int> >::iterator p_iterr; std::map<double, int >::iterator s_iterr; for(p_iterr=cluster_map.begin();p_iterr!=cluster_map.end(); p_iterr++) { int done =0; for(s_iterr=p_iterr->second.begin();s_iterr!=p_iterr->second.end(); s_iterr++) { if(s_iterr->second==cid) { p_iterr->second.erase(s_iterr); done=1; break; } } if(done)break; } for(p_iterr=cluster_map_u.begin();p_iterr!=cluster_map_u.end(); p_iterr++) { int done =0; for(s_iterr=p_iterr->second.begin();s_iterr!=p_iterr->second.end(); s_iterr++) { if(s_iterr->second==cid) { p_iterr->second.erase(s_iterr); done=1; break; } } if(done)break; } for(p_iterr=cluster_map_v.begin();p_iterr!=cluster_map_v.end(); p_iterr++) { int done =0; for(s_iterr=p_iterr->second.begin();s_iterr!=p_iterr->second.end(); s_iterr++) { if(s_iterr->second==cid) { p_iterr->second.erase(s_iterr); done=1; break; } } if(done)break; } std::vector<Managed::ManagedCluster>::iterator it; for(it=clusters.begin();it!=clusters.end();it++) { if(it->id==cid) { clusters.erase(it); break; } } }

void ClusterManager::FillClusterMap (  std::map< double, std::map< double, std::pair< double, int > > > *  cluster_map  )

Definition at line 975 of file ClusterManager.cxx. References clusters. Referenced by Finder::Process(). { for(unsigned int i=0;i<clusters.size();i++) { (*cluster_map)[clusters[i].z][clusters[i].t]=std::pair<double,int>(clusters[i].e, clusters[i].view); } }

std::vector< int > ClusterManager::FindClustersInZ (  double  z, int  view )

Definition at line 44 of file ClusterManager.cxx. Referenced by PrimaryShowerFinder::ExpandShowerChain(). { double planedist = 0.04; std::map<double, std::map<double, int> >*cmap = 0; if(view==2)cmap=&cluster_map_u; else if(view==3)cmap=&cluster_map_v; else cmap=&cluster_map; std::map<double, std::map<double, int> >::iterator it = cmap->lower_bound(z-planedist); std::vector<int> rec; for(;it!=cmap->end();it++) { if(it->first > z+planedist)break; std::map<double, int>::iterator s_iter; for(s_iter=it->second.begin();s_iter!=it->second.end(); s_iter++) { rec.push_back(s_iter->second); } } return rec; }

Managed::ManagedCluster * ClusterManager::GetCluster (  int  cid  )

Definition at line 833 of file ClusterManager.cxx. References clusters, clustersaver, Managed::ClusterSaver::GetCluster(), and inuse. Referenced by AdjustCluster(), LongMuonFinder::CheckChainQuality(), PrimaryShowerFinder::CleanHoughLines(), HoughView::DrawClusters(), ChainView::DrawClusters(), DumpClusters(), PrimaryShowerFinder::ExpandShowerChain(), PrimaryShowerFinder::FindHoughLineMatches(), LongMuonFinder::FindIsolationZ(), LongMuonFinder::FindMuonChain(), Finder::FindNeutrons(), PrimaryShowerFinder::FindPrimaryShower(), PrimaryShowerFinder::FindShowerChain(), PrimaryShowerFinder::LoadCloseHits(), Finder::Make3DParticle(), PrimaryShowerFinder::MakeChains(), Finder::MakeChains(), PrimaryShowerFinder::MakeHoughMap(), PrimaryShowerFinder::MakeParticle3D(), PrimaryShowerFinder::MakeShowerChain(), LongMuonFinder::MergeChainClusters(), MergeClusters(), Finder::RecordLostHits(), LongMuonFinder::RemoveNonMuonEnergy(), SaveCluster(), Finder::SetStatus(), PrimaryShowerFinder::SplitHoughLine(), and Finder::Weave(). { if(cid==0)return 0; if(cid<0) //means it is a saved cluster { if(clustersaver)return clustersaver->GetCluster(cid); else return 0; } ManagedCluster * c =0; for(unsigned int i=0;i<clusters.size();i++) { if(clusters[i].id==cid)c=&clusters[i]; } for(unsigned int i=0;i<inuse.size();i++) { if(inuse[i]==cid)return 0; } return c; }

std::map< double, std::map< double, int > > * ClusterManager::GetClusterMap (  int  view = 0  )

Definition at line 72 of file ClusterManager.cxx. References RebuildClusterMaps(). Referenced by LongMuonFinder::AbsorbMuonClusters(), HoughView::DrawClusters(), ChainView::DrawClusters(), LongMuonFinder::FindIsolationZ(), LongMuonFinder::FindMuonChain(), PrimaryShowerFinder::FindShowerChain(), PrimaryShowerFinder::LoadCloseHits(), PrimaryShowerFinder::MakeChains(), Finder::MakeChains(), PrimaryShowerFinder::MakeHoughMap(), and LongMuonFinder::RemoveNonMuonEnergy(). { //need rebuild? if(needMapRebuild)RebuildClusterMaps(); if(view==2)return &cluster_map_u; if(view==3)return &cluster_map_v; return &cluster_map; }

Managed::ClusterSaver* Managed::ClusterManager::GetClusterSaver (   )  [inline]

Definition at line 68 of file ClusterManager.h. Referenced by LongMuonFinder::MakeParticle3D(). {return clustersaver;};

HitManager* Managed::ClusterManager::GetHitManager (   )  [inline]

Definition at line 65 of file ClusterManager.h. {return hitmanager;};

std::vector<int> Managed::ClusterManager::GetPlanes (   )

Managed::ManagedCluster * ClusterManager::GetSavedCluster (  int  cid  )

Definition at line 856 of file ClusterManager.cxx. References clustersaver, and Managed::ClusterSaver::GetCluster(). Referenced by LongMuonFinder::FindLongMuon(). { return clustersaver->GetCluster(cid); }

std::vector<int> Managed::ClusterManager::GetStrips (  int  plane  )

std::vector< int > ClusterManager::GetViewIndex (  int  view  )

Definition at line 306 of file ClusterManager.cxx. References clusters. Referenced by Finder::FindNeutrons(). { std::vector<int> ret; for(unsigned int i=0;i<clusters.size();i++) if(clusters[i].view==view)ret.push_back(i); return ret; }

void ClusterManager::LoadHits (   )

Definition at line 742 of file ClusterManager.cxx. References Managed::HitManager::GetAvailableHits(), hitmanager, hits, loc_map, maxu, maxv, maxz, minu, minv, minz, and Reset(). Referenced by MakeClusters(). { if(!hitmanager)return; Reset(); std::vector<ManagedHit> avail_hits = hitmanager->GetAvailableHits(); for(unsigned int i=0;i<avail_hits.size();i++) { hits.push_back(avail_hits[i]); loc_map[avail_hits[i].GetPlane()][avail_hits[i].GetStrip()]=hits.size()-1; double sz = avail_hits[i].GetZ(); double st = avail_hits[i].GetT(); minz=minz < sz ? minz : sz; maxz=maxz > sz ? maxz : sz; int view=avail_hits[i].GetView(); if(view==2) { minu=minu < st ? minu : st; maxu=maxu > st ? maxu : st; }else if(view==3) { minv=minv < st ? minv : st; maxv=maxv > st ? maxv : st; } } }

void ClusterManager::MakeClusters (  double  min_cluster_e = 0.4, double  min_strip_e = 0.4, double  max_t_skip = 0.025 )

Definition at line 393 of file ClusterManager.cxx. References Managed::ManagedCluster::AdvanceID(), cluster_map, cluster_map_u, cluster_map_v, clusters, Managed::ManagedCluster::e, Managed::ManagedCluster::Finalize(), Managed::ManagedCluster::hit_id, Managed::ManagedCluster::hite, hitmanager, Managed::ManagedCluster::hitplane, hits, Managed::ManagedCluster::hitstrip, Managed::ManagedCluster::hitt, Managed::ManagedCluster::hitz, Managed::ManagedCluster::id, Managed::ManagedCluster::Insert(), LoadHits(), loc_map, MSG, Managed::ManagedCluster::Reset(), Managed::ManagedCluster::ResetIDCounter(), Managed::ManagedCluster::tsortmap, and Managed::ManagedCluster::view. Referenced by Finder::Process(). { //see if we have hits if(hits.size()==0 && hitmanager==0)return; LoadHits(); if(hits.size()==0)return; // double thresh=0.4; ///for a cluster // double sthresh=0.4;//for an individual strip 0.2 mips ~ 1 pe // roughly gaussian 7.361 w/ rms 1.252 pe/mip - lowest bin is 5 pe/mip //int skip=1; // double skipt=0.05; //strip width is 0.025? double thresh=min_cluster_e; double sthresh=min_strip_e; double skipt=max_t_skip; std::vector<ManagedCluster> theclusters; std::map<int, std::map<int, int> >::iterator p_iter; std::map<int, int>::iterator s_iter; double loc_z=-1; //int last_strip=-1; int numstrip=0; double sum_e_t=0; double sum_e=0; double lastt=-100000.0; int thisview=0; ManagedCluster c; c.ResetIDCounter(); c.AdvanceID(); // cout <<"idcounter 1 " <<c.idcounter<< " "<<c.id<<endl; for(p_iter=loc_map.begin();p_iter!=loc_map.end(); p_iter++) { std::map<int, int>::iterator s_iter; for(s_iter=p_iter->second.begin();s_iter!=p_iter->second.end(); s_iter++) { if(hits[s_iter->second].GetERemaining()<sthresh)continue; if(numstrip==0) lastt=hits[s_iter->second].GetT(); if((fabs(hits[s_iter->second].GetT()-lastt)>skipt|| loc_z!=hits[s_iter->second].GetZ()) && numstrip>0) { //we have finished with this cluster, so record it and reset info sum_e_t/=sum_e; //weighted average t if(thresh<sum_e) { cluster_map[loc_z][sum_e_t]=c.id; if(thisview==2)cluster_map_u[loc_z][sum_e_t]=c.id; if(thisview==3)cluster_map_v[loc_z][sum_e_t]=c.id; c.Finalize(); theclusters.push_back(c); } //printf("----cluster at %f %f %f %d id %d\n\n", loc_z, sum_e_t,sum_e, thisview,c.id); sum_e_t=0; sum_e=0; // loc_z=-1; numstrip=0; thisview=0; c.Reset(); c.AdvanceID(); // cout <<"idcounter 2 " <<c.idcounter<< " "<<c.id<<endl; if(loc_z!=hits[s_iter->second].GetZ()) { lastt=-100000.0; } //last_strip=-1; } loc_z=hits[s_iter->second].GetZ(); lastt=hits[s_iter->second].GetT(); //last_strip=s_iter->first; numstrip++; sum_e+=hits[s_iter->second].GetERemaining(); sum_e_t+=hits[s_iter->second].GetERemaining()*hits[s_iter->second].GetT(); thisview=hits[s_iter->second].GetView(); c.view=hits[s_iter->second].GetView(); c.Insert(hits[s_iter->second].GetZ(),hits[s_iter->second].GetT(),hits[s_iter->second].GetERemaining(),hits[s_iter->second].GetPlane(),hits[s_iter->second].GetStrip(),hits[s_iter->second].GetID()); //printf("%d - %d %d %f %f %f %d\n",numstrip,p_iter->first, s_iter->first, hits[s_iter->second].GetERemaining(),hits[s_iter->second].GetZ(),hits[s_iter->second].GetT(),hits[s_iter->second].GetView()); } } //get the last cluster/// if(numstrip>0) { sum_e_t/=sum_e; //weighted average t if(thresh<sum_e) { cluster_map[loc_z][sum_e_t]=c.id; if(thisview==2)cluster_map_u[loc_z][sum_e_t]=c.id; if(thisview==3)cluster_map_v[loc_z][sum_e_t]=c.id; c.Finalize(); theclusters.push_back(c); } } MSG("ClusterManager",Msg::kDebug)<<"----cluster at "<< loc_z<<" "<<sum_e_t<<" "<<sum_e<<" "<<thisview<<" id "<<c.id<<"\n"; /* //uncomment this to bypass the sharing for(unsigned int k=0;k<theclusters.size();k++) { theclusters[k].Finalize(); clusters.push_back(theclusters[k]); } return; */ c.Reset(); c.AdvanceID(); //cout <<"idcounter 3 " <<c.idcounter<< " "<<c.id<<endl; std::vector<ManagedCluster> tmponepeak; std::vector<ManagedCluster> tmpmultpeak; //go through clusters, looking for multiple peaks.... for(unsigned int k=0;k<theclusters.size();k++) { std::map<double,int>::iterator it; it=theclusters[k].tsortmap.begin(); int npeak=0; std::vector<int> maxidx; std::vector<double> maxe; std::vector<int> sortedidx; std::vector<double> sortede; MSG("ClusterManager",Msg::kDebug)<<"clu idx "<< k<<" id "<<theclusters[k].id <<" with "<<(int)theclusters[k].hite.size() <<" hits \n"; if(theclusters[k].hite.size()==1) { theclusters[k].Finalize(); if(theclusters[k].e>0.001) tmponepeak.push_back(theclusters[k]); continue; } int plastidx=-1; double plaste=0; int lastidx=it->second; double laste = theclusters[k].hite[it->second]; MSG("ClusterManager",Msg::kDebug)<<"\t z "<<theclusters[k].hitz[it->second]<<" t "<<theclusters[k].hitt[it->second]<<" e "<<theclusters[k].hite[it->second]<<"\n"; it++; for(;it!=theclusters[k].tsortmap.end();it++) { double e=theclusters[k].hite[it->second]; MSG("ClusterManager",Msg::kDebug)<<"\t z "<<theclusters[k].hitz[it->second]<<" t "<<theclusters[k].hitt[it->second]<<" e "<<theclusters[k].hite[it->second]<<"\n"; if( e < laste&&(laste>plaste || plastidx<0) ) { maxidx.push_back(lastidx); maxe.push_back(laste); MSG("ClusterManager",Msg::kDebug)<<"\t\tpeak e "<<laste<<"\n"; npeak++; }else{ //these hits needs to be shared! sortedidx.push_back(lastidx); sortede.push_back(laste); } plaste=laste; plastidx=lastidx; laste=e; lastidx=it->second; } //last hit if(plastidx>-1 && laste>plaste) { maxidx.push_back(lastidx); maxe.push_back(laste); MSG("ClusterManager",Msg::kDebug)<<"\t\tpeak e "<<laste<<"\n"; npeak++; }else{ //these hits needs to be shared! sortedidx.push_back(lastidx); sortede.push_back(laste); } MSG("ClusterManager",Msg::kDebug)<<"cluster "<<theclusters[k].id<< " npeak "<<npeak<<"\n"; if(npeak<2) //1 hit cluster will have 0 peaks! { theclusters[k].Finalize(); if(theclusters[k].e>0.001) tmponepeak.push_back(theclusters[k]); continue; } //split up the cluster! for(unsigned int i=0;i<maxidx.size();i++) { c.Reset(); c.AdvanceID(); //cout <<"idcounter 4 " <<c.idcounter<< " "<<c.id<<endl; double t=theclusters[k].hitt[maxidx[i]]; double z=theclusters[k].hitz[maxidx[i]]; double maxe=theclusters[k].hite[maxidx[i]]; int plane=theclusters[k].hitplane[maxidx[i]]; int strip=theclusters[k].hitstrip[maxidx[i]]; int hitid=theclusters[k].hit_id[maxidx[i]]; MSG("ClusterManager",Msg::kDebug)<<"filling m "<<z<<" "<<t<<" "<<maxe<<" "<<plane<<" "<<strip<<"\n"; c.Insert(z, t, maxe, plane, strip, hitid); for(unsigned int j=0;j<sortedidx.size();j++) { //calculate normalizing denominator.... double denom=0; for(unsigned int m=0;m<maxidx.size();m++) { double maxt=theclusters[k].hitt[maxidx[m]]; double sortedt=theclusters[k].hitt[sortedidx[j]]; denom+=theclusters[k].hite[maxidx[m]]/ ( (maxt-sortedt) * (maxt-sortedt) ); } double st = theclusters[k].hitt[sortedidx[j]]; double e= sortede[j];// / ( (st - t ) * (st -t) ); e*=maxe/((t-theclusters[k].hitt[sortedidx[j]]) * (t-theclusters[k].hitt[sortedidx[j]]))/denom; double t=theclusters[k].hitt[sortedidx[j]]; double z=theclusters[k].hitz[sortedidx[j]]; int plane=theclusters[k].hitplane[sortedidx[j]]; int strip=theclusters[k].hitstrip[sortedidx[j]]; int hitid=theclusters[k].hit_id[sortedidx[j]]; MSG("ClusterManager",Msg::kDebug)<<"filling s "<<z<<" "<<t<<" "<<plane<<" "<<strip<<" --- "<<st<<" --- "<<maxe/denom<<" tot "<<sortede[j]<<" used "<<e<<"\n"; if(e>0.001) { c.Insert(z, t, e, plane, strip,hitid); } } c.view = theclusters[k].view; c.Finalize(); if(c.hitt.size()>0 && c.e>0.001) { tmponepeak.push_back(c); MSG("ClusterManager",Msg::kDebug)<<"new cluster stored with size "<<c.hitt.size()<<" e "<<c.e<<" in view "<<c.view<<"\n"; } } } cluster_map.clear(); cluster_map_u.clear(); cluster_map_v.clear(); for(unsigned int k=0;k<tmponepeak.size();k++) { cluster_map[tmponepeak[k].z][tmponepeak[k].t]=tmponepeak[k].id; tmponepeak[k].Finalize(); clusters.push_back(tmponepeak[k]); if(tmponepeak[k].view==2)cluster_map_u[tmponepeak[k].z][tmponepeak[k].t]=tmponepeak[k].id; if(tmponepeak[k].view==3)cluster_map_v[tmponepeak[k].z][tmponepeak[k].t]=tmponepeak[k].id; MSG("ClusterManager",Msg::kDebug)<<"saving cluster id "<<tmponepeak[k].id<<"\n"; // minz=minz < tmponepeak[k].GetZ() ? minz : tmponepeak[k].GetZ(); // maxz=maxz > tmponepeak[k].GetZ() ? maxz : tmponepeak[k].GetZ(); // mint=mint < tmponepeak[k].GetT() ? mint : tmponepeak[k].GetT(); // maxt=maxt > tmponepeak[k].GetT() ? maxt : tmponepeak[k].GetT(); } //printf("ClusterManager dimensions z,t %f %f z,t %f %f\n",minz,mint, maxz,maxt); }

void Managed::ClusterManager::MarkUsed (  int  id  )  [inline]

Definition at line 76 of file ClusterManager.h. Referenced by PrimaryShowerFinder::MakeChains(). {inuse.push_back(id);};

int ClusterManager::MergeClusters (  int  cid1, int  cid2 )

Definition at line 863 of file ClusterManager.cxx. References cluster_map, cluster_map_u, cluster_map_v, clusters, EraseCluster(), Managed::ManagedCluster::Finalize(), GetCluster(), Managed::ManagedCluster::hit_id, Managed::ManagedCluster::hite, Managed::ManagedCluster::hitplane, Managed::ManagedCluster::hitstrip, Managed::ManagedCluster::hitt, Managed::ManagedCluster::hitz, Managed::ManagedCluster::id, Managed::ManagedCluster::Insert(), MSG, Managed::ManagedCluster::t, Managed::ManagedCluster::view, and Managed::ManagedCluster::z. Referenced by LongMuonFinder::MergeChainClusters(). { ManagedCluster *c1 = GetCluster(cid1); ManagedCluster *c2 = GetCluster(cid2); if(!c1 && c2)return cid2; if(!c2 && c1)return cid1; if(!c1 && !c2)return -1; ManagedCluster newcluster; MSG("ClusterManager",Msg::kDebug)<<"merging "<<cid1<<" "<<cid2<< " making "<<newcluster.id<<"\n"; //copy over the old hits... for(unsigned int i=0;i<c1->hit_id.size();i++) newcluster.Insert(c1->hitz[i],c1->hitt[i],c1->hite[i],c1->hitplane[i],c1->hitstrip[i],c1->hit_id[i]); for(unsigned int i=0;i<c2->hit_id.size();i++) newcluster.Insert(c2->hitz[i],c2->hitt[i],c2->hite[i],c2->hitplane[i],c2->hitstrip[i],c2->hit_id[i]); // printf("%d %d\n",c1->id,c2->id); //the pointers c1 and c2 seem to be disrupted by the assignment of the newcluster to the maps... why? //because the pointers are into a vector... which has its elements shifted! EraseCluster(cid1); EraseCluster(cid2); cluster_map[newcluster.z][newcluster.t]=newcluster.id; newcluster.Finalize(); clusters.push_back(newcluster); if(newcluster.view==2)cluster_map_u[newcluster.z][newcluster.t]=newcluster.id; if(newcluster.view==3)cluster_map_v[newcluster.z][newcluster.t]=newcluster.id; return newcluster.id; }

void ClusterManager::RebuildClusterMaps (   )  [private]

Definition at line 88 of file ClusterManager.cxx. References cluster_map, cluster_map_u, cluster_map_v, clusters, clusters_to_delete, and needMapRebuild. Referenced by GetClusterMap(). { if(clusters_to_delete.size()>0) { std::vector<Managed::ManagedCluster> cluster_temp; for(unsigned int i=0;i<clusters.size();i++) { int keep=1; for(unsigned int j=0;j<clusters_to_delete.size();j++) { if(clusters_to_delete[j]==clusters[i].id) { keep=0; break; } } if(!keep)continue; cluster_temp.push_back(clusters[i]); } clusters=cluster_temp; cluster_temp.clear(); clusters_to_delete.clear(); } cluster_map.clear(); cluster_map_u.clear(); cluster_map_v.clear(); for(unsigned int i=0;i<clusters.size();i++) { cluster_map[clusters[i].z][clusters[i].t]=clusters[i].id; if(clusters[i].view==2)cluster_map_u[clusters[i].z][clusters[i].t]=clusters[i].id; if(clusters[i].view==3)cluster_map_v[clusters[i].z][clusters[i].t]=clusters[i].id; } needMapRebuild=0; }

void ClusterManager::Reset (   )

Definition at line 807 of file ClusterManager.cxx. References cluster_map, cluster_map_u, cluster_map_v, clusters, clusters_are_made, clusters_to_delete, hits, inuse, loc_map, maxu, maxv, maxz, minu, minv, minz, and needMapRebuild. Referenced by LoadHits(). { clusters_are_made=0; loc_map.clear(); cluster_map.clear(); cluster_map_u.clear(); cluster_map_v.clear(); hits.clear(); clusters.clear(); minu=100000; maxu=-100000; minv=100000; maxv=-100000; minz=100000; maxz=-100000; needMapRebuild=0; clusters_to_delete.clear(); inuse.clear(); }

int ClusterManager::SaveCluster (  int  cluster_id, double  energy_to_save = 0, int  status = 0 )

Definition at line 132 of file ClusterManager.cxx. References AdjustCluster(), clustersaver, Managed::ManagedCluster::e, GetCluster(), needMapRebuild, Managed::ClusterSaver::SaveCluster(), Managed::ManagedCluster::SetStatus(), and SplitEnergy(). Referenced by LongMuonFinder::FindLongMuon(), PrimaryShowerFinder::FindPrimaryShower(), and Finder::Weave(). { if(cluster_id<0) { // printf("!!!! attempt to save an already saved cluster!\n"); return cluster_id;//its already a saved cluster! } //tell them we need a new map because we may be removing clusters.... needMapRebuild=1; ManagedCluster *cluster = GetCluster(cluster_id); if(!cluster)return 0; //allow the ability to fake a save cluster for tempory work if(!clustersaver) { // printf("SAVING WITHOUT SAVER!\n"); //can only allocation the full energy in the cluster... if(energy_to_save-cluster->e >= 0) { // printf("setting e to 0, was %f\n",cluster->e); SplitEnergy(cluster,cluster->e,0); if(hitmanager)AdjustCluster(cluster); return 0; } SplitEnergy(cluster,energy_to_save,0); if(hitmanager)AdjustCluster(cluster); return 0; } //can only allocation the full energy in the cluster... if(energy_to_save-cluster->e >= 0 || energy_to_save==0) { energy_to_save=cluster->e; int newid= clustersaver->SaveCluster(cluster); ManagedCluster *newcluster=GetCluster(newid); if(newcluster) { newcluster->SetStatus(status); SplitEnergy(cluster,newcluster->e,0); }else{ SplitEnergy(cluster,0,0);//throw away a cluster since the new cluster can't be made... probably because the old cluster had 0 e anyways. } AdjustCluster(cluster); // if(newcluster)printf("saving full cluster with e %f old has %f\n",newcluster->e,cluster->e); return newid; } //if we are here we want to save only part of the cluster! ManagedCluster newcluster = *cluster; SplitEnergy(&newcluster,energy_to_save,1); SplitEnergy(cluster,energy_to_save,0); newcluster.SetStatus(status); int newid = clustersaver->SaveCluster(&newcluster); AdjustCluster(cluster); // printf("saving split save has %f old has %f\n",GetCluster(newid)->e,cluster->e); return newid; }

void Managed::ClusterManager::SetClusterSaver (  Managed::ClusterSaver *  cs  )  [inline]

Definition at line 67 of file ClusterManager.h. Referenced by PrimaryShowerFinder::MakeChains(), and ParticleFinder::Reco(). {clustersaver=cs;};

void Managed::ClusterManager::SetHitManager (  Managed::HitManager *  h  )  [inline]

Definition at line 64 of file ClusterManager.h. Referenced by PrimaryShowerFinder::MakeChains(), and ParticleFinder::Reco(). {hitmanager=h;};

void ClusterManager::SplitEnergy (  Managed::ManagedCluster *  cluster, double  energy, int  matchenergy )

Definition at line 252 of file ClusterManager.cxx. References Managed::ManagedCluster::e, Managed::ManagedCluster::Finalize(), Managed::ManagedCluster::hite, and MSG. Referenced by SaveCluster(). { MSG("ClusterManager",Msg::kDebug)<<"splitting cluster with "<<cluster->e<<" to have "<<energy <<" match "<<matchenergy<<"\n"; double totale=cluster->e; if(totale<1e-9)return; for(unsigned int i=0;i<cluster->hite.size();i++) { double newe =0; if(matchenergy)newe=cluster->hite[i]/totale*energy; else newe=cluster->hite[i]/totale*(totale-energy); /* printf("hit with %f should have %f\n",cluster->hite[i],newe); if(matchenergy) { //needs to also adjust the hit stored in the hit manager here..... Managed::ManagedHit *hit = hitmanager->FindHit(cluster->hit_id[i]); double sete = hit->SetEnergy(newe); cluster->hite[i]=sete; printf("setting to %f\n",sete); }else{ //add save the rest as a new hit if(newe>0) { printf("newhit energy %f\n",newe); Managed::ManagedHit *hit = hitmanager->FindHit(cluster->hit_id[i]); int newid = hitmanager->InsertHit(hit->GetView(), hit->GetPlane(), hit->GetStrip(), hit->GetZ(), hit->GetT(), newe); cluster->hite[i]=newe; cluster->hit_id[i]=newid; printf("setting to %f\n",newe); }else{ cluster->hite[i]=0; cluster->hit_id[i]=-1; printf("setting to 0"); } } */ cluster->hite[i]=newe; // if(matchenergy)hite[i]=hite[i]/totale*energy; // else hite[i]=hite[i]/totale*(totale-energy); } cluster->Finalize(); }

---

Member Data Documentation

std::map<double, std::map<double, int> > Managed::ClusterManager::cluster_map

Definition at line 46 of file ClusterManager.h. Referenced by DumpClusters(), EraseCluster(), MakeClusters(), MergeClusters(), RebuildClusterMaps(), and Reset().

std::map<double, std::map<double, int> > Managed::ClusterManager::cluster_map_u

Definition at line 47 of file ClusterManager.h. Referenced by DumpClusters(), EraseCluster(), MakeClusters(), MergeClusters(), RebuildClusterMaps(), and Reset().

std::map<double, std::map<double, int> > Managed::ClusterManager::cluster_map_v

Definition at line 48 of file ClusterManager.h. Referenced by DumpClusters(), EraseCluster(), MakeClusters(), MergeClusters(), RebuildClusterMaps(), and Reset().

std::vector<Managed::ManagedCluster> Managed::ClusterManager::clusters

Definition at line 52 of file ClusterManager.h. Referenced by EraseCluster(), FillClusterMap(), Finder::FindNeutrons(), GetCluster(), GetViewIndex(), MakeClusters(), MergeClusters(), Finder::Process(), RebuildClusterMaps(), Finder::RecordLostHits(), and Reset().

int Managed::ClusterManager::clusters_are_made [private]

Definition at line 79 of file ClusterManager.h. Referenced by Reset().

std::vector<int> Managed::ClusterManager::clusters_to_delete [private]

Definition at line 96 of file ClusterManager.h. Referenced by AdjustCluster(), RebuildClusterMaps(), and Reset().

ClusterSaver* Managed::ClusterManager::clustersaver [private]

Definition at line 83 of file ClusterManager.h. Referenced by GetCluster(), GetSavedCluster(), and SaveCluster().

HitManager* Managed::ClusterManager::hitmanager [private]

Definition at line 82 of file ClusterManager.h. Referenced by AdjustCluster(), LoadHits(), and MakeClusters().

std::vector<ManagedHit> Managed::ClusterManager::hits [private]

Definition at line 81 of file ClusterManager.h. Referenced by AddStrip(), LoadHits(), MakeClusters(), and Reset().

std::vector<int> Managed::ClusterManager::inuse [private]

Definition at line 94 of file ClusterManager.h. Referenced by Finder::FindNeutrons(), GetCluster(), Finder::RecordLostHits(), and Reset().

std::map<int, std::map<int,int> > Managed::ClusterManager::loc_map [private]

Definition at line 80 of file ClusterManager.h. Referenced by AddStrip(), LoadHits(), MakeClusters(), and Reset().

double Managed::ClusterManager::maxu

Definition at line 59 of file ClusterManager.h. Referenced by AddStrip(), PrimaryShowerFinder::FindFirstIntersection(), LoadHits(), and Reset().

double Managed::ClusterManager::maxv

Definition at line 61 of file ClusterManager.h. Referenced by AddStrip(), PrimaryShowerFinder::FindFirstIntersection(), LoadHits(), and Reset().

double Managed::ClusterManager::maxz

Definition at line 57 of file ClusterManager.h. Referenced by AddStrip(), PrimaryShowerFinder::FindFirstIntersection(), LoadHits(), and Reset().

double Managed::ClusterManager::minu

Definition at line 60 of file ClusterManager.h. Referenced by AddStrip(), PrimaryShowerFinder::FindFirstIntersection(), LoadHits(), and Reset().

double Managed::ClusterManager::minv

Definition at line 62 of file ClusterManager.h. Referenced by AddStrip(), PrimaryShowerFinder::FindFirstIntersection(), LoadHits(), and Reset().

double Managed::ClusterManager::minz

Definition at line 58 of file ClusterManager.h. Referenced by AddStrip(), PrimaryShowerFinder::FindFirstIntersection(), LoadHits(), and Reset().

int Managed::ClusterManager::needMapRebuild [private]

Definition at line 91 of file ClusterManager.h. Referenced by RebuildClusterMaps(), Reset(), and SaveCluster().

---

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

• ClusterManager.h
• ClusterManager.cxx

---