remid Namespace Reference

A PID for muons. More...

Classes
struct	Counts
class	RecoMuon
class	ReMId
class	ReMIdDedx
class	ReMIdDedxFD
class	ReMIdDedxRock
class	ReMIdDedxStudies
class	ReMIdTrain
class	ReMIdValidate
	A module to analyze remid objects produced by the RecoMuon module. More...

Functions
unsigned int	HighestPIDTrack (const std::vector< art::Ptr< rb::Track > > &sliceTracks, const std::string &remidModuleLabel, const art::Event &e)
unsigned int	HighestPIDTrack (std::vector< art::Ptr< rb::Track > > const &sliceTracks, art::InputTag const &remidTag, art::Event const &e)

Detailed Description

A PID for muons.

A module for analyzing the ReMId objects in neutrino interactions.

Author

Nicholas Raddatz - radda.nosp@m.tz@p.nosp@m.hysic.nosp@m.s.um.nosp@m.n.edu

Date

September 2012

Function Documentation

unsigned int remid::HighestPIDTrack	(	const std::vector< art::Ptr< rb::Track > > &	sliceTracks,
		const std::string &	remidModuleLabel,
		const art::Event &	e
	)

Function which finds index of the best ReMId track among a vector of tracks. Requires the event to fetch associated ReMId objects and the module label

Definition at line 249 of file ReMId.cxx.

References getGoodRuns4SAM::tag.

Referenced by numue::NumuEAna::analyze(), numue::NumuEAlg::FDEnergy(), skim::SkimmerAna::FillNumuHistograms(), skim::CosmicBeamComparison::FillNumuHistograms(), numue::NumuEAlg::GetUCMuonENonSingle(), numue::NumuEAlg::GetUCMuonESingle(), skim::EvaluatorNumuCCpi::Initialize(), skim::EvaluatorNumu::Initialize(), numue::NumuEAlg::MCTruthEnergyVariables(), numue::NumuEAlg::NDEnergy(), qeef::QeFinder::produce(), wsnumu::WSNumu::produce(), cosrej::CosRej::produce(), xsrec::MakeXSecCCPi0Inc::produce(), caf::CAFMaker::produce(), and numue::NumuEAlg::QEFormulaEnergy().

  {
    
    art::InputTag tag(remidModuleLabel);
    return remid::HighestPIDTrack(sliceTracks, tag, e);
  }

unsigned int remid::HighestPIDTrack	(	std::vector< art::Ptr< rb::Track > > const &	sliceTracks,
		art::InputTag const &	remidTag,
		art::Event const &	e
	)

Definition at line 259 of file ReMId.cxx.

References plot_validation_datamc::c, allTimeWatchdog::endl, geo::kXorY, cet::sqlite::max(), rb::PID::Value(), and POTSpillRate::view.

  {
    int    bestTrack        = -1;
    double highestPID       = -2.0;
    double bestTrackLength  = 0.0;
    double bestStartZ       = std::numeric_limits<double>::max();
    double bestADC          = std::numeric_limits<double>::max();
    
    art::FindOneP<remid::ReMId> tracksToReMId(sliceTracks, e, remidTag);
    
    for(size_t c = 0; c < sliceTracks.size(); ++c){
      
      geo::View_t view = sliceTracks[c]->View();
      
        //Only use 3D tracks
      if (view != geo::kXorY){
        continue;
      } //End of loop over 2d tracks
      
      if(tracksToReMId.isValid()){
        art::Ptr<remid::ReMId> trackReMId = tracksToReMId.at(c);
        
          // If the PID is lower, we can bail, nothing to update
        if(trackReMId->Value() < highestPID) continue;
        
        if(trackReMId->Value() == highestPID)
        { // If the new guy is shorter, we can bail
          if (sliceTracks[c]->TotalLength() < bestTrackLength) continue;
          
            // If equal, need to fall to the next level, low start Z
          if (sliceTracks[c]->TotalLength() == bestTrackLength){
            
              // If start Z is higher, we can bail
            if (sliceTracks[c]->Start().Z() > bestStartZ) continue;
            
              // If equal, need to fall to the next level, total ADC
              // You might think this is stupid, but it's equivalent to flipping
              // a coin which always lands on the same side.  This is good, this
              // will happen so rarely that you should just pee up a rope.
              // Use lower ADC
            if (sliceTracks[c]->Start().Z() == bestStartZ){
              
              if(sliceTracks[c]->TotalADC() > bestADC) continue;
              
              if(sliceTracks[c]->TotalADC() == bestADC)
              throw cet::exception("ReMId Matched Track Exception")<<std::endl
              << "Can't break a tie between two ReMId/tracks.  "
              << "This should have never happend, but it did.  "
              << "The numu group is going to have to reflect on this "
              << "and find a way to fix it." << std::endl;
              
            }
            
            
          }
        }
        
        
        
        highestPID = trackReMId->Value();
        bestTrackLength = sliceTracks[c]->TotalLength();
        bestStartZ = sliceTracks[c]->Start().Z();
        bestADC    = sliceTracks[c]->TotalADC();
        bestTrack = c;
        
        
      }//End of valid ReMId object
      
    }//End of for loop over tracks
    
    
    if (bestTrack < 0)
    return 999;
    
    return (unsigned int)bestTrack;
    
    
  }