NuCutImps::CCA_NC Class Reference

#include <NuCutImps.h>

Inheritance diagram for NuCutImps::CCA_NC:

List of all members.

Public Member Functions
	CCA_NC (const NuPlots *plots=0)
Bool_t	InFidVol (const NuEvent &nu) const
void	Preselection (const NuEvent &nu)
void	Selection (const NuEvent &nu)
Protected Member Functions
void	FDCleaning (const NuEvent &nu)
void	FDCosmics (const NuEvent &nu)
void	NDCleaning (const NuEvent &nu)
Private Attributes
CCAPresel	fCCAPresel

---

Constructor & Destructor Documentation

NuCutImps::CCA_NC::CCA_NC (  const NuPlots *  plots = 0  )

Definition at line 582 of file NuCutImps.cxx. References plots(), and NuCut::SetFidVol(). : NuCut("CCA_NC", plots) { // Match CC fiducial volume for simplicity // TODO CJB - will this try to use the track vertex? Is that bad? SetFidVol("cc2008"); }

---

Member Function Documentation

void NuCutImps::CCA_NC::FDCleaning (  const NuEvent &  nu  )  [protected]

Definition at line 632 of file NuCutImps.cxx. References NuCut::Cut_If(), NuEvent::evtphsigcor, FDCosmics(), max, NuEvent::nearestSpillNanosec, NuEvent::nevt, NuEvent::nstripEvt, NuEvent::snarlPulseHeight, and NuEvent::trigtime. Referenced by Preselection(). { // See NCUtils/Cuts/NCAnalysisCutsNC for original version Cut_If(nu.nevt == 0, "NumEvents"); const double phFrac = nu.evtphsigcor/max(nu.snarlPulseHeight, 1e-10); Cut_If(nu.nevt > 2 || (nu.nevt == 2 && phFrac < .75), "PHfrac"); FDCosmics(nu); // Cut out cosmics Cut_If(nu.evtphsigcor < 5000 && nu.nstripEvt <= 4, "FibreNoise"); // Timing cut. Data only const double dt_spill = nu.trigtime*Munits::s - nu.nearestSpillNanosec*Munits::ns; Cut_Data_If(dt_spill < -2*Munits::microsecond || dt_spill > 12*Munits::microsecond, nu, "Timing"); }

void NuCutImps::CCA_NC::FDCosmics (  const NuEvent &  nu  )  [protected]

Definition at line 656 of file NuCutImps.cxx. Referenced by FDCleaning(). { //Ultra-steep showers Cut_If(nu.planeEvtN == 0 || nu.nstripEvt >= SQR(nu.planeEvtN), "UltraSteep"); //Steep showers const double tRMS = TMath::Sqrt(SQR(nu.transverseRMSU)+ SQR(nu.transverseRMSV)); Cut_If(nu.nshw > 0 && tRMS > .3+.1901*TMath::Log10(nu.planeEvtN), "Steep"); //Track angle Cut_If(nu.ntrk > 0 && nu.trkvtxdcosz < 0.4, "TrkDCosZ"); const bool upExitTrack = ((nu.zTrkEnd > 28.78 || nu.endMetersToCloseEdge < .5) && nu.yTrkEnd > -1.657); Cut_If(nu.ntrk > 0 && upExitTrack && nu.dtdz < -.5 && nu.trkvtxdcosy < -.4, "UpExitTrack"); }

Bool_t NuCutImps::CCA_NC::InFidVol (  const NuEvent &  nu  )  const [virtual]

Fiducial volume calculation function. This is the basic 'infid' derived version, and should be overridden/passed to for any more complicated evaluations. Reimplemented from NuCut. Definition at line 591 of file NuCutImps.cxx. References NuCut::InFidVol(), NuCut::InFidVolEvt(), NuEvent::nplaneShw, NuEvent::ntrk, NuEvent::trknplane, and NuEvent::zTrkVtx. Referenced by Preselection(). { // If there's a track and it's longer than the shower then use its // vertex, otherwise use the event vertex. This is what is done in // the NC analysis. // TODO CJB - is this maximally consistent with the CC analysis? if(nu.ntrk > 0 && nu.trknplane > nu.nplaneShw){ // Move trk vtx upstream by 3.92cm from scintillator to steel NuEvent nuc = nu; nuc.zTrkVtx = nu.zTrkVtx - (0.0392*Munits::m); return NuCut::InFidVol(nuc); } return NuCut::InFidVolEvt(nu); }

void NuCutImps::CCA_NC::NDCleaning (  const NuEvent &  nu  )  [protected]

Definition at line 679 of file NuCutImps.cxx. References NuCut::Cut_If(). Referenced by Preselection(). { // See NCUtils/Cuts/NCAnalysisCuts::IsMultiCutsClean for original version // 35ns timing cut Cut_If(TMath::Abs(nu.minTimeSeparation) < 35e-9, "TimeSep"); // tiny events are mostly junk Cut_If(nu.nstripEvt < 5, "TotalStrips"); // this cuts very steep showers, leaking in Cut_If(nu.nstripEvt > 1.15*SQR(nu.planeEvtN), "Steep"); // this cuts leakage which leaves activity in partially instrumented // region Cut_If(nu.edgeActivityStrips > 2 && nu.edgeActivityPH > 1000 && nu.energy < 5 && nu.nplaneShw > nu.trknplane, "Leakage"); Cut_If(nu.oppEdgeStrips > 2 && nu.energy < 5.0 && nu.nplaneShw > nu.trknplane, "OppositeEdge"); // make additional deltaZ cuts if (|minDeltaT|<175ns) Cut_If(TMath::Abs(nu.closeTimeDeltaZ) < 1.0 && TMath::Abs(nu.minTimeSeparation) < 175e-9, "dz&dt"); }

void NuCutImps::CCA_NC::Preselection (  const NuEvent &  nu  )  [virtual]

Implements NuCut. Definition at line 607 of file NuCutImps.cxx. References NuCut::Defer_Preselection(), NuEvent::detector, fCCAPresel, FDCleaning(), InFidVol(), NuCut::Keep_If(), and NDCleaning(). { Defer_Preselection(fCCAPresel, nu); // Now do the rest of the preselection // No good track cut // Is it in the fiducial volume? Keep_If(InFidVol(nu), "FidVol"); // No track fit pass or direction cosines or track reclamation // TODO CJB - check CC cuts don't overcut these if(nu.detector == Detector::kFar ) FDCleaning(nu); if(nu.detector == Detector::kNear) NDCleaning(nu); }

void NuCutImps::CCA_NC::Selection (  const NuEvent &  nu  )  [virtual]

Implements NuCut. Definition at line 626 of file NuCutImps.cxx. References NuCut::Keep_If(), and NuEvent::roID. { Keep_If(nu.roID < 0.3, "RoID"); }

---

Member Data Documentation

CCAPresel NuCutImps::CCA_NC::fCCAPresel [private]

Definition at line 138 of file NuCutImps.h. Referenced by Preselection().

---

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

• NuCutImps.h
• NuCutImps.cxx

---