diff --git a/run/macros/NewSpec.C b/run/macros/NewSpec.C index 510c6eb..377c497 100644 --- a/run/macros/NewSpec.C +++ b/run/macros/NewSpec.C @@ -464,7 +464,7 @@ void NewSpec::CoinIOUD( Bool_t FigFlag, Double_t eCut ) } -void NewSpec::TBLRCoinIO( Bool_t FigFlag, Double_t eCut ) +void NewSpec::TBLRCoinIO(Bool_t FigFlag, Double_t eCut, Bool_t MDecayNotInSample) { // Double_t eCut = 0.7; //minimum energy deposition in scintillators [MeV] Bool_t tofFlag = 0; @@ -501,9 +501,11 @@ void NewSpec::TBLRCoinIO( Bool_t FigFlag, Double_t eCut ) // if (ientry < 0) break; nb = fChain->GetEntry(jentry); nbytes += nb; // fChain->GetEntry(jentry); - -// if (muDecayPosZ>16.9&&muDecayPosZ<17.1) continue; //muons decaying outside of sample plane (in the radiation shield etc) - + + if (MDecayNotInSample){ + if (muDecayPosZ>16.9&&muDecayPosZ<17.1) continue; //muons decaying outside of sample plane (in the radiation shield etc) + } + // FILL IN HISTOGRAMS hEdepositCF->Fill(save_ke[0]/1000.); hBeamSpot->Fill(save_x[0], save_y[0]); diff --git a/run/macros/NewSpec.h b/run/macros/NewSpec.h index 0f80a75..199d240 100644 --- a/run/macros/NewSpec.h +++ b/run/macros/NewSpec.h @@ -130,7 +130,7 @@ public : virtual void CreateIO(Bool_t FigFlag, Double_t eCut); virtual void CoinIO(Bool_t FigFlag, Double_t eCut); virtual void CoinIOUD(Bool_t FigFlag, Double_t eCut); - virtual void TBLRCoinIO(Bool_t FigFlag, Double_t eCut); + virtual void TBLRCoinIO(Bool_t FigFlag, Double_t eCut, Bool_t MDecayNotInSample); virtual void TBLRCoinDown(Bool_t FigFlag, Double_t eCut); virtual void TBLRCoinUp(Bool_t FigFlag, Double_t eCut); virtual void SCoinIO(Bool_t FigFlag, Double_t eCut);