add separate sum for FCAL2 with gains for trigger equation

Author: jrstevenjlab

libraries/DSelector/DSelector.cc

@@ -97,6 +97,10 @@ void DSelector::Setup_Branches(void)
 	dL1BCALEnergy = (locL1BCALEnergyBranch != NULL) ? (Double_t*)locL1BCALEnergyBranch->GetAddress() : NULL;
 	TBranch* locL1FCALEnergyBranch = dTreeInterface->Get_Branch("L1FCALEnergy");
 	dL1FCALEnergy = (locL1FCALEnergyBranch != NULL) ? (Double_t*)locL1FCALEnergyBranch->GetAddress() : NULL;
+	TBranch* locL1ECALEnergyBranch = dTreeInterface->Get_Branch("L1ECALEnergy");
+	dL1ECALEnergy = (locL1ECALEnergyBranch != NULL) ? (Double_t*)locL1ECALEnergyBranch->GetAddress() : NULL;
+	TBranch* locL1FCAL2EnergyBranch = dTreeInterface->Get_Branch("L1FCAL2Energy");
+	dL1FCAL2Energy = (locL1FCAL2EnergyBranch != NULL) ? (Double_t*)locL1FCAL2EnergyBranch->GetAddress() : NULL;
 
 	// MC
 	if(dIsMCFlag)

libraries/DSelector/DSelector.h

@@ -92,6 +92,8 @@ class DSelector : public TSelector
 		UInt_t Get_L1TriggerBits(void) const;
 		Double_t Get_L1BCALEnergy(void) const;
 		Double_t Get_L1FCALEnergy(void) const;
+	        Double_t Get_L1ECALEnergy(void) const;
+		Double_t Get_L1FCAL2Energy(void) const;
 		Bool_t Get_IsThrownTopology(void) const;
 		Float_t Get_MCWeight(void) const;
 		Float_t Get_GeneratedEnergy(void) const;
@@ -155,6 +157,8 @@ class DSelector : public TSelector
 		UInt_t* dL1TriggerBits;
 		Double_t* dL1BCALEnergy;
 		Double_t* dL1FCALEnergy;
+	        Double_t* dL1ECALEnergy;
+		Double_t* dL1FCAL2Energy;
 		Float_t* dMCWeight; //only present if simulated data
 		Float_t* dGeneratedEnergy; //only present if simulated data
 		Bool_t* dIsThrownTopology; //only present if simulated data
@@ -207,6 +211,7 @@ inline DSelector::DSelector(TTree* locTree) :
 		dBeamWrapper(NULL), dComboWrapper(NULL), dAnalysisActions(vector<DAnalysisAction*>()),
 		dFile(NULL), dOutputFlatTreeFile(NULL), dProofFile(NULL), dOutputFlatTreeProofFile(NULL),
 		dTreeNumber(0), dRunNumber(NULL), dEventNumber(NULL), dL1TriggerBits(NULL), dMCWeight(NULL), dGeneratedEnergy(NULL), dIsThrownTopology(NULL), dX4_Production(NULL),
+		dL1BCALEnergy(NULL), dL1FCALEnergy(NULL), dL1ECALEnergy(NULL), dL1FCAL2Energy(NULL), 
 		dNumBeam(NULL), dNumChargedHypos(NULL), dNumNeutralHypos(NULL), dNumCombos(NULL), dNumThrown(NULL),
 		dNumPIDThrown_FinalState(NULL), dPIDThrown_Decaying(NULL) {}
 
@@ -238,6 +243,16 @@ inline Double_t DSelector::Get_L1FCALEnergy(void) const
 	return *dL1FCALEnergy;
 }
 
+inline Double_t DSelector::Get_L1ECALEnergy(void) const
+{
+	return *dL1ECALEnergy;
+}
+
+inline Double_t DSelector::Get_L1FCAL2Energy(void) const
+{
+	return *dL1FCAL2Energy;
+}
+
 inline Bool_t DSelector::Get_IsThrownTopology(void) const
 {
 	return ((dIsThrownTopology != NULL) ? *dIsThrownTopology : false);