LMU/musrsim
LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Option to generate decay e+ spectra for muons which do not decay in the sample plane

Browse Source
This commit is contained in:
prokscha 2017-10-11 16:35:08 +02:00
parent b37eb3b12f
commit d7b33f07ae
2 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
run/macros/NewSpec.C
View File

@ -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] // Double_t eCut = 0.7; //minimum energy deposition in scintillators [MeV]
Bool_t tofFlag = 0; Bool_t tofFlag = 0;
@ -502,7 +502,9 @@ void NewSpec::TBLRCoinIO( Bool_t FigFlag, Double_t eCut )
nb = fChain->GetEntry(jentry); nbytes += nb; nb = fChain->GetEntry(jentry); nbytes += nb;
// fChain->GetEntry(jentry); // 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 // FILL IN HISTOGRAMS
hEdepositCF->Fill(save_ke[0]/1000.); hEdepositCF->Fill(save_ke[0]/1000.);

2
run/macros/NewSpec.h
View File

@ -130,7 +130,7 @@ public :
virtual void CreateIO(Bool_t FigFlag, Double_t eCut); virtual void CreateIO(Bool_t FigFlag, Double_t eCut);
virtual void CoinIO(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 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 TBLRCoinDown(Bool_t FigFlag, Double_t eCut);
virtual void TBLRCoinUp(Bool_t FigFlag, Double_t eCut); virtual void TBLRCoinUp(Bool_t FigFlag, Double_t eCut);
virtual void SCoinIO(Bool_t FigFlag, Double_t eCut); virtual void SCoinIO(Bool_t FigFlag, Double_t eCut);