JiveXML::SiSpacePointRetrieverHelpers Namespace Reference

Functions
unsigned int	getTruthBarcodes (const Identifier idFirst, const Identifier idSecond, const PRD_MultiTruthCollection *truthColl, DataVect &barcodes)
	Get the barcodes of associated truth particles for a SpacePoint ( we require both clusters, in case there are)

Function Documentation

getTruthBarcodes()

unsigned int JiveXML::SiSpacePointRetrieverHelpers::getTruthBarcodes	(	const Identifier	idFirst,
		const Identifier	idSecond,
		const PRD_MultiTruthCollection *	truthColl,
		DataVect &	barcodes )

Get the barcodes of associated truth particles for a SpacePoint ( we require both clusters, in case there are)

Parameters

idFirst	the SpacePoint identifier of first cluster
idSecond	the SpacePoint identifier of second cluster (only for SCT)
truthColl	the truth collection
barcodes	the date vector to which to append the barcodes

Returns

number of associated truth particles

NOTE: This function could be simplified if we could directly insert into the DataVect, instead of barcodesCommon. At the moment, this doesnt work as we are lacking operator< and operator== for DataType

Definition at line 55 of file SiSpacePointRetriever.cxx.

                                                                                                                                                       {
      
      //Make life easier
      using PRDTruthIter = PRD_MultiTruthCollection::const_iterator;
 

 
      //Sets of barcodes associated with first, second and both cluster
      std::set<int> barcodesFirst;
      std::set<int> barcodesSecond;
      std::set<int> barcodesCommon;
 
      //Get the set of particle barcodes associated with the first cluster identifier
      std::pair<PRDTruthIter,PRDTruthIter> equalRangeFirst = truthColl->equal_range(idFirst);
      for(PRDTruthIter TruthCollItr=equalRangeFirst.first; TruthCollItr!= equalRangeFirst.second; ++TruthCollItr)
        barcodesFirst.insert(TruthCollItr->second.barcode());
 
      //Check if we have only have one valid cluster identifier
      if (! idSecond.is_valid()){
 
        //Make this our list of barcodes (swap is O(1))
        barcodesCommon.swap(barcodesFirst);
 
      //otherwise only store barcodes associated with both identifiers
      } else {
 
        //Get the set of particle barcodes associated with the second cluster identifier
        std::pair<PRDTruthIter,PRDTruthIter> equalRangeSecond = truthColl->equal_range(idSecond);
        for(PRDTruthIter TruthCollItr=equalRangeSecond.first; TruthCollItr!= equalRangeSecond.second; ++TruthCollItr)
          barcodesSecond.insert(TruthCollItr->second.barcode());
 
        //Copy the list of particle barcodes that are associated with both clusters
        std::set_intersection(barcodesFirst.begin(), barcodesFirst.end(), barcodesSecond.begin(), barcodesSecond.end(), 
                              std::insert_iterator< std::set<int> >(barcodesCommon,barcodesCommon.begin()) );
      }
 
      //Finally add the list of barcodes to our DataVect
      for (const auto barcodeCommon : barcodesCommon) barcodes.push_back(DataType(barcodeCommon));
 
      //return the number of added barcodes
      return barcodesCommon.size();
    }