LMU/musrfit
LMU/musrfit
5
0
Fork 0
Code Issues 6 Pull Requests Actions Packages Projects Releases 3 Activity

changed to complex Mu polarization function

Browse Source
This commit is contained in:
Thomas Prokscha 2016-02-22 22:47:58 +01:00
parent 1c5c94a86b
commit 7ae1b9ab5e
2 changed files with 30 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

37
src/tests/MuTransition/PSimulateMuTransition.cpp
View File

@ -73,6 +73,7 @@
using namespace std; using namespace std;
#include <TMath.h> #include <TMath.h>
#include <TComplex.h>
#include "PSimulateMuTransition.h" #include "PSimulateMuTransition.h"
@ -246,7 +247,7 @@ Double_t PSimulateMuTransition::PrecessionPhase(const Double_t &time, const TStr
if (chargeState == "Mu+") if (chargeState == "Mu+")
muonPhaseX = TMath::TwoPi()*fMuonPrecFreq*time; muonPhaseX = TMath::TwoPi()*fMuonPrecFreq*time;
else if (chargeState == "Mu0"){ else if (chargeState == "Mu0"){
muoniumPolX = GTFunction(time); muoniumPolX = GTFunction(time).Re();
muonPhaseX = TMath::ACos(muoniumPolX); muonPhaseX = TMath::ACos(muoniumPolX);
} }
else else
@ -263,15 +264,29 @@ Double_t PSimulateMuTransition::PrecessionPhase(const Double_t &time, const TStr
* *
* \param time (us); * \param time (us);
*/ */
Double_t PSimulateMuTransition::GTFunction(const Double_t &time) TComplex PSimulateMuTransition::GTFunction(const Double_t &time)
{ {
Double_t muoniumPolX = 0; Double_t twoPi = TMath::TwoPi();
muoniumPolX = 0.5 * TComplex complexPol = 0;
(fMuFractionState12 * (TMath::Cos(TMath::TwoPi()*fMuPrecFreq12*time) + TMath::Cos(TMath::TwoPi()*fMuPrecFreq34*time)) + complexPol =
fMuFractionState23 * (TMath::Cos(TMath::TwoPi()*fMuPrecFreq23*time) + TMath::Cos(TMath::TwoPi()*fMuPrecFreq14*time))); 0.5 * fMuFractionState12 *
(TComplex::Exp(TComplex::I()*twoPi*fMuPrecFreq12*time) +
TComplex::Exp(-TComplex::I()*twoPi*fMuPrecFreq34*time))
+
0.5 * fMuFractionState23 *
(TComplex::Exp(TComplex::I()*twoPi*fMuPrecFreq23*time) +
TComplex::Exp(TComplex::I()*twoPi*fMuPrecFreq14*time));
return complexPol;
// Double_t muoniumPolX = 0;
// muoniumPolX = 0.5 *
// (fMuFractionState12 * (TMath::Cos(twoPi*fMuPrecFreq12*time) + TMath::Cos(twoPi*fMuPrecFreq34*time)) +
// fMuFractionState23 * (TMath::Cos(twoPi*fMuPrecFreq23*time) + TMath::Cos(twoPi*fMuPrecFreq14*time)));
//
// return muoniumPolX;
return muoniumPolX;
} }
//-------------------------------------------------------------------------- //--------------------------------------------------------------------------
@ -285,6 +300,7 @@ Double_t PSimulateMuTransition::GTFunction(const Double_t &time)
*/ */
Double_t PSimulateMuTransition::GTSpinFlip(const Double_t &time) Double_t PSimulateMuTransition::GTSpinFlip(const Double_t &time)
{ {
TComplex complexPolX = 1.0;
Double_t muoniumPolX = 1.0; //initial polarization in x direction Double_t muoniumPolX = 1.0; //initial polarization in x direction
Double_t eventTime = 0; Double_t eventTime = 0;
Double_t eventDiffTime = 0; Double_t eventDiffTime = 0;
@ -292,18 +308,19 @@ Double_t PSimulateMuTransition::GTSpinFlip(const Double_t &time)
eventTime += NextEventTime(fSpinFlipRate); eventTime += NextEventTime(fSpinFlipRate);
if (eventTime >= time){ if (eventTime >= time){
muoniumPolX = GTFunction(time); muoniumPolX = GTFunction(time).Re();
} }
else{ else{
while (eventTime < time){ while (eventTime < time){
eventDiffTime = eventTime - lastEventTime; eventDiffTime = eventTime - lastEventTime;
muoniumPolX = muoniumPolX * GTFunction(eventDiffTime); complexPolX = complexPolX * GTFunction(eventDiffTime);
lastEventTime = eventTime; lastEventTime = eventTime;
eventTime += NextEventTime(fSpinFlipRate); eventTime += NextEventTime(fSpinFlipRate);
} }
// calculate for the last collision // calculate for the last collision
eventDiffTime = time - lastEventTime; eventDiffTime = time - lastEventTime;
muoniumPolX = muoniumPolX * GTFunction(eventDiffTime); complexPolX = complexPolX * GTFunction(eventDiffTime);
muoniumPolX = complexPolX.Re();
} }
return muoniumPolX; return muoniumPolX;

3
src/tests/MuTransition/PSimulateMuTransition.h
View File

@ -34,6 +34,7 @@
#include <TObject.h> #include <TObject.h>
#include <TH1F.h> #include <TH1F.h>
#include <TRandom2.h> #include <TRandom2.h>
#include <TComplex.h>
// global constants // global constants
const Double_t fMuonGyroRatio = 135.54; //!< muon gyromagnetic ratio (MHz/T) const Double_t fMuonGyroRatio = 135.54; //!< muon gyromagnetic ratio (MHz/T)
@ -95,7 +96,7 @@ class PSimulateMuTransition : public TObject
virtual Double_t NextEventTime(const Double_t &EventRate); virtual Double_t NextEventTime(const Double_t &EventRate);
// virtual Double_t PrecessionPhase(const Double_t &time, const Double_t &frequency); // virtual Double_t PrecessionPhase(const Double_t &time, const Double_t &frequency);
virtual Double_t PrecessionPhase(const Double_t &time, const TString chargeState); virtual Double_t PrecessionPhase(const Double_t &time, const TString chargeState);
virtual Double_t GTFunction(const Double_t &time); //!< transverse field polarization function of Mu0 virtual TComplex GTFunction(const Double_t &time); //!< transverse field polarization function of Mu0
virtual Double_t GTSpinFlip(const Double_t &time); //!< transverse field polarization function after spin-flip collisions virtual Double_t GTSpinFlip(const Double_t &time); //!< transverse field polarization function after spin-flip collisions
virtual void Event(const TString muonString); virtual void Event(const TString muonString);