NCExtractionCutswkNN Class Reference

A simple cuts based extraction. More...

#include <NCExtractionCutswkNN.h>

Inheritance diagram for NCExtractionCutswkNN:

List of all members.

Public Member Functions
	NCExtractionCutswkNN (NCAnalysisCuts *cuts, const Registry &r)
void	FillAnalysisInfo (NCEventInfo &evtInfo, int beamType)
	Fill the analysisInfo member of evtInfo based on the other members and the beamType (in whatever convention that is).
Protected Member Functions
virtual double	GetIdProbability (NCEventInfo &evtInfo, int)
	The probability that this event is CC.
virtual double	GetCutPosition () const
	Where to cut on the IdProbability.
Private Member Functions
bool	PassesTrackReclamation (NCEventInfo &evtInfo)
	Does the event pass the ND "good track or track reclamation" cut?
bool	PassesNCCuts (NCEventInfo &evtInfo)
	Does the event pass the NC selection algorithm?
bool	PassesCCCuts (NCEventInfo &evtInfo)
	Does the event pass the CC selection algorithm?
Private Attributes
BeamType::BeamType_t	fBeamType

---

Detailed Description

A simple cuts based extraction.

Definition at line 22 of file NCExtractionCutswkNN.h.

---

Constructor & Destructor Documentation

NCExtractionCutswkNN::NCExtractionCutswkNN	(	NCAnalysisCuts *	cuts,
		const Registry &	r
	)

Definition at line 20 of file NCExtractionCutswkNN.cxx.

References BeamType::TagToEnum().

                                                      :
  NCExtraction(cuts, r),
  fBeamType(BeamType::kUnknown)

{
  const char* tmps;
  if (r.Get("BeamType", tmps)) fBeamType=BeamType::TagToEnum(tmps);
}

---

Member Function Documentation

void NCExtractionCutswkNN::FillAnalysisInfo	(	NCEventInfo &	evtInfo,
		int	beamType
	)			[virtual]

Fill the analysisInfo member of evtInfo based on the other members and the beamType (in whatever convention that is).

The default implementation fills separationParameter with the value of GetIdProbability(), separationParameterCut with the value of GetCutPosition(), and sets isCC to 1 if separationParameter >= separationParameterCut, and isNC to 1 otherwise

Reimplemented from NCExtraction.

Definition at line 100 of file NCExtractionCutswkNN.cxx.

References NCEventInfo::analysis, NCExtraction::fCuts, GetCutPosition(), GetIdProbability(), NCAnalysisCuts::InBeamFiducialVolume(), ANtpAnalysisInfo::isCC, NCAnalysisCuts::IsGoodBeamEvent(), ANtpAnalysisInfo::isNC, ANtpAnalysisInfo::pass, PassesCCCuts(), PassesNCCuts(), NCEventInfo::reco, ANtpAnalysisInfo::Reset(), ANtpAnalysisInfo::separationParameter, ANtpAnalysisInfo::separationParameterCut, and ANtpRecoInfo::trackMomentum.

{

  ANtpAnalysisInfo* ana=evtInfo.analysis;
  double separationCut=GetCutPosition();

  ana->Reset();

  ana->separationParameterCut = separationCut;
  ana->separationParameter = GetIdProbability(evtInfo, beamType);

  // Events that pass the NC algorithm are NC...
  ana->isNC = PassesNCCuts(evtInfo);
  // ...of those that aren't, only the ones passing the CC cuts are CC
  ana->isCC = (!ana->isNC) && PassesCCCuts(evtInfo);

  if(ana->isCC) {
    ana->pass = (int)( evtInfo.reco->trackMomentum <= 0 && 
                       fCuts->IsGoodBeamEvent() && 
                       fCuts->InBeamFiducialVolume() );
  } else if(ana->isNC) {
    ana->pass = (int)fCuts->InBeamFiducialVolume();
  }
  else ana->pass=0; // neither NC nor CC
}

virtual double NCExtractionCutswkNN::GetCutPosition

(

)

const [inline, protected, virtual]

Where to cut on the IdProbability.

Implements NCExtraction.

Definition at line 32 of file NCExtractionCutswkNN.h.

Referenced by FillAnalysisInfo().

{return 0.5;}

double NCExtractionCutswkNN::GetIdProbability	(	NCEventInfo &	evtInfo,
		int
	)			[protected, virtual]

The probability that this event is CC.

Implements NCExtraction.

Definition at line 31 of file NCExtractionCutswkNN.cxx.

References Msg::kDebug, MAXMSG, PassesCCCuts(), and PassesNCCuts().

Referenced by FillAnalysisInfo().

{
  MAXMSG("NCExtractionCuts", Msg::kDebug,20) << "GetIdProbabilty" << endl;

  bool passesNCCuts=PassesNCCuts(evtInfo);
  bool passesCCCuts=PassesCCCuts(evtInfo);

  // Mini truth table for return value:
  //     |  CC  | !CC
  //  ----------------
  //  NC |  -1  |   0 
  // !NC |   1  |  -2

  if     ( passesCCCuts && !passesNCCuts) return 1;
  else if(!passesCCCuts &&  passesNCCuts) return 0;
  else if( passesCCCuts &&  passesNCCuts) return -1;
  else if(!passesCCCuts && !passesNCCuts) return -2;

  assert(0 && "This should never happen");
}

bool NCExtractionCutswkNN::PassesCCCuts

(

NCEventInfo &

evtInfo

)

[private]

Does the event pass the CC selection algorithm?

Definition at line 85 of file NCExtractionCutswkNN.cxx.

References ANtpTrackInfoNC::kNN, ANtpTrackInfoNC::kNNShort, and NCEventInfo::track.

Referenced by FillAnalysisInfo(), and GetIdProbability().

{
  return (evtInfo.track->kNN>0.25625 || evtInfo.track->kNNShort > 0.50625 );
}

bool NCExtractionCutswkNN::PassesNCCuts

(

NCEventInfo &

evtInfo

)

[private]

Does the event pass the NC selection algorithm?

Definition at line 53 of file NCExtractionCutswkNN.cxx.

References det, ANtpHeaderInfo::detector, NCEventInfo::event, NCEventInfo::header, Detector::kNear, ANtpEventInfoNC::lengthInPlanes, ANtpTrackInfo::planes, ANtpShowerInfo::planes, ANtpEventInfo::planes, NCEventInfo::shower, NCEventInfo::track, and ANtpEventInfo::tracks.

Referenced by FillAnalysisInfo(), and GetIdProbability().

{
  const Detector::Detector_t det = (Detector::Detector_t)evtInfo.header->detector;

  int tracks = evtInfo.event->tracks;

  int planes = det==Detector::kNear ? evtInfo.event->lengthInPlanes :
                                      evtInfo.event->planes;

  int trackExtension = evtInfo.shower->planes - evtInfo.track->planes;

  //NC Algorithm:

  // 1. Long events are definitely not NC
  if(planes >= 47) return false;
  else{
    // 2. Events without tracks definitely *are* NC
    if(tracks == 0) {
      return true;
    } else {
      // 3. Short events with tracks need to pass the track extension cut
      if(trackExtension > -6) {
        return true;
      } else {
        return false;
      }
    }
  }

}

bool NCExtractionCutswkNN::PassesTrackReclamation

(

NCEventInfo &

evtInfo

)

[private]

Does the event pass the ND "good track or track reclamation" cut?

Definition at line 91 of file NCExtractionCutswkNN.cxx.

References ANtpTrackInfo::begPlaneU, ANtpTrackInfo::begPlaneV, ANtpTrackInfo::endPlane, ANtpTrackInfo::endPlaneU, ANtpTrackInfo::endPlaneV, ANtpTrackInfo::passedFit, and NCEventInfo::track.

{
  return evtInfo.track->passedFit>0 ||
    (TMath::Abs(evtInfo.track->begPlaneU - evtInfo.track->begPlaneV) < 6
     && TMath::Abs(evtInfo.track->endPlaneU - evtInfo.track->endPlaneV) < 41
     && evtInfo.track->endPlane < 270);
}

---

Member Data Documentation

BeamType::BeamType_t NCExtractionCutswkNN::fBeamType [private]

Definition at line 45 of file NCExtractionCutswkNN.h.

---

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

• NCExtractionCutswkNN.h
• NCExtractionCutswkNN.cxx

---