NuCutImps::ChairSound Class Reference

#include <NuCutImps.h>

Inheritance diagram for NuCutImps::ChairSound:

List of all members.

Public Member Functions
	ChairSound (const NuPlots *plots=0)
Static Public Member Functions
Int_t	GetRegion (const NuEvent &nu)
Protected Member Functions
void	Preselection (const NuEvent &nu)
void	Selection (const NuEvent &nu)

---

Constructor & Destructor Documentation

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

Definition at line 713 of file NuCutImps.cxx. References plots(), and NuCut::SetFidVol(). : NuCut("ChairSound", plots) // fNuCutsSelection(NuCuts::kNMB0325ChairSound, plots) { SetFidVol("cc2008"); }

---

Member Function Documentation

Int_t NuCutImps::ChairSound::GetRegion (  const NuEvent &  nu  )  [static]

Definition at line 721 of file NuCutImps.cxx. References det, NuEvent::detector, FidVol::getFarRinner(), FidVol::getFarRouter(), FidVol::getFarZMC(), FidVol::getTrkVtxZOffset(), infid(), NuEvent::simFlag, NuEvent::xTrkVtx, NuEvent::yTrkVtx, and NuEvent::zTrkVtx. Referenced by Preselection(). { // An enum for the rock detector region - this probably shouldn't go here, // but for now it will have to do (would rather use an enum than plain numbers) enum Rock_DetectorRegion { kNoTrack = -1, kFiducial = 0, kFront = 1, kBack = 2, kGap = 3, kEdge = 4, kCoilHole = 5, kCorner = 6, kUnknown = 7 }; // Placeholder to avoid nesting if's in this function Bool_t inFid = false, nearFront = false, nearEdge = false; Bool_t nearGap = false, nearCoil = false, nearBack = false; Double_t radius2 = nu.xTrkVtx*nu.xTrkVtx + nu.yTrkVtx*nu.yTrkVtx; Double_t z = nu.zTrkVtx - FidVol::getTrkVtxZOffset(); // Get the detector Detector::Detector_t det = (Detector::Detector_t)nu.detector; // If we are the near detector, we can't calculate the region - return unknown if (det == Detector::kNear) return kUnknown; // Now - examine each of the regions and determine if the vertex is in there // 0 = Fiducial inFid = infid((Detector::Detector_t)nu.detector, (SimFlag::SimFlag_t)nu.simFlag, nu.xTrkVtx, nu.yTrkVtx, z); if (inFid) return kFiducial; // 1 = Front, non-edge nearFront = z < FidVol::getFarZMC(0); nearEdge = radius2 >= (FidVol::getFarRouter()*FidVol::getFarRouter()); nearCoil = radius2 < FidVol::getFarRinner()*FidVol::getFarRinner(); if(nearFront && !nearEdge && !nearCoil) return kFront; // 2 = Back, anywhere (no corners at back) nearBack = z > FidVol::getFarZMC(3); if(nearBack) return kBack; // 3 = Supermodule gap, non-edge nearGap = z > FidVol::getFarZMC(1) && z < FidVol::getFarZMC(2); if(nearGap && !nearEdge && !nearCoil) return kGap; // 4 = edge, non-end if(nearEdge && !nearBack && !nearFront && !nearGap) return kEdge; // 5 = coil hole, non-end if(nearCoil && !nearBack && !nearFront && !nearGap) return kCoilHole; // 6 = corner cases // No other cases have been matched, it must (?) be a corner case? return kCorner; }

void NuCutImps::ChairSound::Preselection (  const NuEvent &  nu  )  [protected, virtual]

Implements NuCut. Definition at line 780 of file NuCutImps.cxx. References NuEvent::coilIsOk, NuCut::Cut_If(), NuEvent::cutOnBeamInfo, NuEvent::cutOnDataQuality, GetRegion(), NuEvent::goodBeamSntp, NuEvent::isLI, NuCut::Keep_Data_If(), NuCut::Keep_If(), NuEvent::litime, MAXMSG, and NuEvent::ntrk. { // can't use CC0325Std::Preselection because it has a fiducial cut and // an overly tight angle cut. This is copied from there except for those. Keep_If(nu.ntrk >= 1, "GoodNumberofTracks"); // Good tracks Cut_If(nu.isLI || nu.litime != -1, "LICut"); // Filter out LI events Keep_Data_If(!nu.cutOnDataQuality || nu.coilIsOk, nu, "IsCoilOkay"); Keep_Data_If(nu.goodBeamSntp && nu.coilIsOk, nu,"GoodBeamDetPOT"); Keep_Data_If(!nu.cutOnBeamInfo || CC0325Std::IsGoodBeam(nu),nu,"GoodBeam"); // Keep only front face events that are neither near the edge nor coil hole // Keep_If(nu.regionTrkVtx == 1, "front_face"); Keep_If(GetRegion(nu) == 1, "front_face"); // Per-detector cuts if (nu.detector == Detector::kFar){ Keep_If(nu.trkfitpass == 1, "trackfitpass"); Keep_Data_If(nu.isGoodDataQuality||!nu.cutOnDataQuality,nu,"DataQuality"); // [-2, 12] microseconds Keep_Data_If(NuCuts::GoodTimeToNearestSpill(nu, -2, +12) || !nu.cutOnSpillTiming, nu, "SpillTime"); // Cosmic cut? Keep_If(TMath::Abs(nu.trkvtxdcosy) < 0.4, "CosY"); } else if (nu.detector == Detector::kNear){ MAXMSG("CC0325Std",Msg::kWarning, 1) << "Front face numubar selection " " called on near detector. I doubt that this works properly.\n"; Keep_If(NuCuts::IsGoodTrackFitPassReclamation(nu), "trackfitreclamation"); Keep_Data_If(nu.coilCurrent<-1000||!nu.cutOnDataQuality,nu,"CoilCurrent"); } }

void NuCutImps::ChairSound::Selection (  const NuEvent &  nu  )  [protected, virtual]

Implements NuCut. Definition at line 815 of file NuCutImps.cxx. References NuCut::Keep_If(), NuEvent::relativeAngle, NuEvent::roID, and NuEvent::sigqp_qp. { // Other selectors have a Defer_Selection to CC0325Std here, but I // don't think it makes sense here since we don't share a preselection. // However, I honestly have no clear idea what it's even supposed to do, // so I could be wrong. Keep_If(nu.roID >= 0.25, "RoID"); // Note: >= is not the same as > Keep_If(1/nu.sigqp_qp > 3.5, "SigmaQP_QP"); Keep_If(TMath::Abs(nu.relativeAngle - TMath::Pi()) > 2.08, "RelativeAngle"); }

---

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

• NuCutImps.h
• NuCutImps.cxx

---