add new cut files

Author: sdobbs

src/plugins/Analysis/p3pi_hists/DCustomAction_p3pi_Pi0Cuts_FCAL2.cc

@@ -0,0 +1,50 @@
+// $Id$
+//
+//    File: DCustomAction_p3pi_Pi0Cuts_FCAL2.cc
+// Created: Thu Jan 22 11:19:47 EST 2015
+// Creator: jrsteven (on Linux ifarm1401 2.6.32-431.el6.x86_64 x86_64)
+//
+
+#include "DCustomAction_p3pi_Pi0Cuts_FCAL2.h"
+
+void DCustomAction_p3pi_Pi0Cuts_FCAL2::Initialize(const std::shared_ptr<const JEvent>& locEvent)
+{
+	
+}
+
+bool DCustomAction_p3pi_Pi0Cuts_FCAL2::Perform_Action(const std::shared_ptr<const JEvent>& locEvent, const DParticleCombo* locParticleCombo)
+{
+	const DParticleComboStep* locParticleComboStep = locParticleCombo->Get_ParticleComboStep(2);
+        if(Get_Reaction()->Get_ReactionStep(2)->Get_InitialPID() != Pi0)
+                return false;
+	
+	// get final state particles
+   	auto locParticles = Get_UseKinFitResultsFlag() ? locParticleComboStep->Get_FinalParticles() : locParticleComboStep->Get_FinalParticles_Measured();
+
+	int nFECAL = 0;  // FCAL2 = FCAL + ECAL
+	
+        // loop final state particles
+        for(size_t loc_i = 0; loc_i < locParticles.size(); ++loc_i) {
+               if(locParticles[loc_i] == nullptr) continue;
+		
+		// get shower object
+		const DNeutralShower* locNeutralShower = static_cast<const DNeutralShower*>(locParticleComboStep->Get_FinalParticle_SourceObject(loc_i));
+		if(locNeutralShower == NULL) 
+			continue;
+		
+		// count # of FCAL photons and set separate thresholds on FCAL and BCAL energies
+		if(locNeutralShower->dDetectorSystem == SYS_FCAL) {
+			nFECAL++;
+		}
+		if(locNeutralShower->dDetectorSystem == SYS_ECAL) {
+			nFECAL++;
+		}
+	}
+
+	// require 1 or 2 photons in FCAL2 (specify parameter in DReaction setup)
+	if(nFECAL != dMinFECAL) 
+		return false;
+
+	// passed all cuts
+	return true; 
+}

src/plugins/Analysis/p3pi_hists/DCustomAction_p3pi_Pi0Cuts_FCAL2.h

@@ -0,0 +1,46 @@
+// $Id$
+//
+//    File: DCustomAction_p3pi_Pi0Cuts_FCAL2.h
+// Created: Thu Jan 22 11:19:46 EST 2015
+// Creator: jrsteven (on Linux ifarm1401 2.6.32-431.el6.x86_64 x86_64)
+//
+
+#ifndef _DCustomAction_p3pi_Pi0Cuts_FCAL2_
+#define _DCustomAction_p3pi_Pi0Cuts_FCAL2_
+
+#include <string>
+#include <iostream>
+
+#include <JANA/JEvent.h>
+#include "JANA/JApplication.h"
+
+#include "ANALYSIS/DAnalysisAction.h"
+#include "ANALYSIS/DReaction.h"
+#include "ANALYSIS/DParticleCombo.h"
+#include "ANALYSIS/DAnalysisUtilities.h"
+
+using namespace std;
+
+class DCustomAction_p3pi_Pi0Cuts_FCAL2 : public DAnalysisAction
+{
+	public:
+
+        DCustomAction_p3pi_Pi0Cuts_FCAL2(const DReaction* locReaction, bool locUseKinFitResultsFlag, double locMinFECAL, 
+        			string locActionUniqueString = "") : 
+	        DAnalysisAction(locReaction, "Custom_p3pi_Pi0Cuts", locUseKinFitResultsFlag, locActionUniqueString), 
+	        dMinFECAL(locMinFECAL) {}
+
+		void Initialize(const std::shared_ptr<const JEvent>& locEvent);
+		void Run_Update(const std::shared_ptr<const JEvent>& locEvent) {}
+
+	private:
+
+		bool Perform_Action(const std::shared_ptr<const JEvent>& locEvent, const DParticleCombo* locParticleCombo);
+
+		//Store any histograms as member variables here
+		double dMinFECAL;   // FCAL2 = FCAL+ECAL
+		
+};
+
+#endif // _DCustomAction_p3pi_Pi0Cuts_FCAL2_
+