evd::GraphCluster Class Reference

Inheritance diagram for evd::GraphCluster:

Public Member Functions
	GraphCluster (fhicl::ParameterSet const &)
void	produce (art::Event &evt)

Private Member Functions
void	GetStartEndHits (unsigned int plane, recob::Hit *starthit, recob::Hit *endhit)
void	GetStartEndHits (unsigned int plane)
void	GetHitList (unsigned int plane, art::PtrVector< recob::Hit > &ptrhitlist)
std::vector< util::PxLine >	GetSeedLines ()

Private Attributes
GraphClusterAlg	fGClAlg
unsigned int	fNPlanes
int	TestFlag
int	fRun
int	fSubRun
int	fEvent
std::vector< recob::Hit * >	starthit
std::vector< recob::Hit * >	endhit
std::vector< util::PxLine >	startendpoints

Detailed Description

Definition at line 83 of file GraphCluster_module.cc.

Constructor & Destructor Documentation

evd::GraphCluster::GraphCluster ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 138 of file GraphCluster_module.cc.

                                                          :
    EDProducer{pset},
  fGClAlg(pset.get< fhicl::ParameterSet >("GraphClusterAlg"))
  {
    art::ServiceHandle<geo::Geometry const>  geo;


    produces< std::vector<recob::Cluster>                >();
    produces< art::Assns<recob::Cluster, recob::Hit>     >();
    produces< std::vector < art::PtrVector <recob::Cluster> >  >();


    fNPlanes = geo->Nplanes();
    starthit.resize(fNPlanes);
    endhit.resize(fNPlanes);


    startendpoints.resize(fNPlanes);
//     swire.resize(fNPlanes);
//     ewire.resize(fNPlanes);
//     stime.resize(fNPlanes);
//     etime.resize(fNPlanes);
  }

Member Function Documentation

void evd::GraphCluster::GetHitList ( unsigned int  plane, art::PtrVector< recob::Hit > &  ptrhitlist  )

private

std::vector< util::PxLine > evd::GraphCluster::GetSeedLines ( )

private

void evd::GraphCluster::GetStartEndHits ( unsigned int  plane, recob::Hit *  starthit, recob::Hit *  endhit  )

private

void evd::GraphCluster::GetStartEndHits ( unsigned int  plane )

private

void evd::GraphCluster::produce

(

art::Event &

evt

)

Todo:

This method appears to produce a recob::Cluster really as it is

a collection of 2D clusters from single planes

Definition at line 167 of file GraphCluster_module.cc.

  {

    std::unique_ptr<std::vector<recob::Cluster> > Graphcol(new std::vector<recob::Cluster>);
    std::unique_ptr< art::Assns<recob::Cluster, recob::Hit>     > hassn(new art::Assns<recob::Cluster, recob::Hit>);
    //  std::unique_ptr< art::Assns<recob::Cluster, recob::Cluster>     > classn(new art::Assns<recob::Cluster, recob::Cluster>);
    std::unique_ptr< std::vector < art::PtrVector < recob::Cluster > > > classn(new std::vector < art::PtrVector < recob::Cluster > >);




    art::ServiceHandle<geo::Geometry const>  geo;

    // check if evt and run numbers check out, etc...
    if(fGClAlg.CheckValidity(  evt ) == -1)
      {
        return;
      }


    for(unsigned int ip=0;ip<fNPlanes;ip++) {
      startendpoints[ip]=util::PxLine();  //assign empty PxLine

    }

    std::vector < art::PtrVector < recob::Hit > > hitlist;
    hitlist.resize(fNPlanes);

    for(unsigned int ip=0;ip<fNPlanes;ip++) {

        fGClAlg.GetHitListAndEndPoints(ip,hitlist[ip],startendpoints[ip]);
      // Read in the Hit List object(s).
     //fGClAlg.GetHitList(ip,hitlist[ip]);

      if(hitlist[ip].size()==0)
        continue;
      //Read in the starthit:
      // GetStartEndHits(ip, starthit[ip],endhit[ip]);




      //fGClAlg.GetStartEndHits(&startendpoints[ip]);

      if(hitlist[ip].size()>0 && !(TestFlag==-1 ) )   //old event or transfer not ready
        {
          double swterror=0.,ewterror=0.;

          if(startendpoints[ip].w0==0 )
            swterror=999;

          if(startendpoints[ip].t1==0 )
            ewterror=999;

          std::cout << " clustering @ " <<startendpoints[ip].w0 << " +/- "<< swterror
                    <<" " <<  startendpoints[ip].t0<< " +/- "<< swterror
                    <<" " <<  startendpoints[ip].w1<< " +/- "<< ewterror
                    <<" " << startendpoints[ip].t1<< " +/- "<< ewterror << std::endl;

          // this is all we can do easily without getting the full-blown
          // ClusterParamsAlg (that means lareventdisplay has to depend on larreco!)
          lar::util::StatCollector<float> integral, summedADC;
          for (art::Ptr<recob::Hit> const& hit: hitlist[ip]) {
            integral.add(hit->Integral());
            summedADC.add(hit->SummedADC());
          } // for

          // get the plane ID from the first hit
          geo::PlaneID planeID = hitlist[ip].front()->WireID().planeID();
          Graphcol->emplace_back(
            startendpoints[ip].w0,
            swterror,
            startendpoints[ip].t0,
            swterror,
            0.,                         // start_charge
            0.,                         // start_angle
            0.,                         // start_opening
            startendpoints[ip].w1,
            ewterror,
            startendpoints[ip].t1,
            ewterror,
            0.,                         // end_charge
            0.,                         // end_angle
            0.,                         // end_opening
            integral.Sum(),             // integral
            integral.RMS(),             // integral_stddev
            summedADC.Sum(),            // summedADC
            summedADC.RMS(),            // summedADC_stddev
            hitlist[ip].size(),         // n_hits
            0.,                         // multiple_hit_density
            0.,                         // width
            ip,
            geo->Plane(ip,planeID.TPC,planeID.Cryostat).View(),
            planeID,
            recob::Cluster::Sentry
            );

          // associate the hits to this cluster
          util::CreateAssn(*this, evt, *Graphcol, hitlist[ip], *hassn);
        }

    }// end of loop on planes

    art::PtrVector < recob::Cluster > cvec;
    cvec.reserve(fNPlanes);

    for(unsigned int ip=0;ip<fNPlanes;ip++)  {
      art::ProductID aid = evt.getProductID< std::vector < recob::Cluster > >();
      art::Ptr< recob::Cluster > aptr(aid, ip, evt.productGetter(aid));
      cvec.push_back(aptr);
    }

    classn->push_back(cvec);

    //  for(unsigned int ip=0;ip<fNPlanes;ip++) {
    //    for(unsigned int jp=ip+1;jp<fNPlanes;jp++) {
    //      util::CreateSameAssn(*this, evt, *Graphcol, *Graphcol, *classn,ip,ip+1,jp );
    //    //  std::cout << "associating cluster" << ip <<" with cluster " << jp << std::endl;
    //    }
    //  }
    //

    evt.put(std::move(Graphcol));
    evt.put(std::move(hassn));
    evt.put(std::move(classn));

    return;
  } // end of produce

Member Data Documentation

std::vector< recob::Hit * > evd::GraphCluster::endhit

private

Definition at line 120 of file GraphCluster_module.cc.

int evd::GraphCluster::fEvent

private

Definition at line 114 of file GraphCluster_module.cc.

GraphClusterAlg evd::GraphCluster::fGClAlg

private

Definition at line 93 of file GraphCluster_module.cc.

unsigned int evd::GraphCluster::fNPlanes

private

Definition at line 109 of file GraphCluster_module.cc.

int evd::GraphCluster::fRun

private

Definition at line 112 of file GraphCluster_module.cc.

int evd::GraphCluster::fSubRun

private

Definition at line 113 of file GraphCluster_module.cc.

std::vector< util::PxLine > evd::GraphCluster::startendpoints

private

Definition at line 124 of file GraphCluster_module.cc.

std::vector< recob::Hit * > evd::GraphCluster::starthit

private

Definition at line 119 of file GraphCluster_module.cc.

int evd::GraphCluster::TestFlag

private

Definition at line 111 of file GraphCluster_module.cc.

---

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

• GraphCluster_module.cc