AnaPAW – an introduction

Leave a comment

In RIKEN Nishina Center for Accelerator-Based Science, the analysis program is called AnaPAW, which is Ana + PAW. PAW is an old analysis program developed by CERN, and no longer supported. It is going to be replaced by ROOT.

PAW is based on Fortran 77. Ana is a package for reading the speciallized data format (.ridf) in RIKEN.

some terminology for AnaPAW:

  • analyzer : a program for analyzing the raw data. some analyzers will be shared by difference detectors of same type. i.e. if there are many plastic scintillators, they may use a same program for processing the raw data.
  • analyzer ID : which is the number for each analyzer, also called analyzer flag
  • ID : this usually referring which “detectors” you want to get in an analyzer.
  • W# : what kind of data you want.
  • enc*.f : analyzer code file
  • *.fh : common value file
  • *.ana : ana code file. this file specific which analyzers will be called, and want histogram will be filled

the flow of AnaPAW is:

  1. call a data file ready to read. But NOT read yet.
  2. using an “ana code” file (.ana), calling “analyzers” ( analyzing program ), and specific which data will be plotted.
  3. the AnaPaw will read the data file according to the ana code file specification.
  4. fill histograms according to ana code file
  5. loop the data file to next data block
  6. select which histogram to plot.

one difficulty of using anapaw is that, the relationship of analyzers are hidden in the code. Some analyzers are independent, but most of them are dependent. This dependence is due to the calling of other analyzer value (common value). in order to get the common value, other analyzer must be run in advance. thus, to use an analyzer, we have to check the enc*.f file to see which analyzer is needed.

another difficulty is that, to know how to use the analyzer, we have to read the code file and most of them have comment on the top.

i always suggest they use include the analyzer dependency in the comment. However, since the dependency can be multiple levels. i.e. A depends on B, B depends on C, and so on. thus, that make anapaw is very scary for beginner.


Leave a comment

long time ago, the post on parity is quit a mess and confusing. now, i have a better understanding and try to do it once again.

here, i will clarify 2 things:

  1. what is parity.
  2. what will happen if parity is hold or violated.

the answer of 1st question is

Parity is a transform that reverse the space dimension.

it sound simple. the tricky part is, there are 2 kinds of transform: active or passive. in short, active transform is change of vector. while pass is change of coordinate. To avoid confusion and make me and reader thinking to be the same, i use active transform in here.

in 1-D space, the parity is just reverse of the direction. in 2-D space, parity is equal a rotation of 180 degree on the plane, thus, it is not really a parity but as a special case for rotation. in 3-D space, parity is just like changing our left-hand to right hand and vice versa. many people will say that, it is same as MIRROR transform, which is a plan inverse. But it is not true. the correct description is a mirror transform + 180 rotation around the normal of the mirror plan. Thus, if we only concern about the shape instead of orientation, a mirror transform is ok. Nevertheless, when we talking current loop, the direction is important that, the magnet is a mirrored inside the mirror. and we knew that all magnetic field is either by current loop or spin.  in 4-D, well, i don’t know.

To answer the 2nd question, we have to know what is parity conservation. the parity conservation is that, thing keep the same under parity transform. i.e.

P \Psi(x) = \Psi(-x) = \pm \Psi(x)

the plus sign means EVEN parity, and minus sign means ODD parity. people knowing some maths will know that, it is corresponding to even function and odd function.

Lets look at the parity of electromagnetic phenomena. lets there is a current loop, rotating at left hand. the B field will be pointing down. under MIRROR (not parity) by building an identical current loop which rotates in right hand, then the magnetic field pointing upward. Wait! is it mean magnetic field violated parity? NO, because it is the pseudovector nature of magnetic field, which, keep the direction perpendicular to the plan of mirror and change sign on parallel to the mirror plan. and as i said before, a parity is mirror + 180 rotation. thus, the parity is conserved. (polar vector has opposites properties that, under mirror, the component parallel to mirror plan is kept but the component perpendicular reversed. well this is what mirror transform was defined)

Now, we found that, all pseudovector are even parity and polar vector are odd parity. (although i don’t have a prove, i believe this is strongly related to spin1/2 and spin 1, and every pseudovector is formed by polar vector. )

Thus, a parity conservation mean that, all physics are same under parity. The other thing to say so, is, the mirrored world is same after we rotated it back.

P=MR=RM \Rightarrow RP=PR=M

in order to understood what is parity conserved or violated, we study elastics scattering. if we only have polar vector, a parity violation means the space has preference direction, i.e. the symmetry broken by a pseudovector field. but that contradict our assumption of “only polar vector”. thus, it is impossible to have parity violation for polar vectors. (i am not sure it is a prove or not)

A famous example of parity violation in nature is the experiment on Cobalt-60 nuclei, which will undergoes beta decay.


*** at the beginning, i was very confused. one major confusion is that, the active transform or passive transform and the mirror world. i found that, inside the mirror, we should NOT use “left-hand rule” instead of  “right-hand rule”.  if we do a mirror transform, and use “left-hand rule” , the transform will be cancelled. because, changing the “handiness” is same as passive transform.

another difficulty is that, a mirror transform is not parity transform, but most popular science mixed them up. a Major difference between mirror and parity is: there are 3 difference reflection plans, but parity is a “point inversion”. when you select difference plan of reflection, your current loop will keep the same or change direction, that make it sometimes keep orientation, but sometimes not.

that difficulty raised because lack of understanding the proper transform of pesudo vector under mirror transform. we the help of the picture, we can now easily imagine parity as mirror.

A few though on alpha decay

Leave a comment

I had a little chat with SASAMOTO San on her pervious research on carbon-12. In ground state, carbon-12 is a stable nucleus with 6 protons and 6 neutrons. In one of the excited state, it will use alpha decay to release energy. According to alpha decay theory, an alpha particle should be formed before the decay and oscillates or hit the column potential barrier million and million of time to get tunneling through. And theory calculation on the wavefunction show that the excited carbon-12 will form 3 alpha particles.

since the 6 protons has 15 way to combine in group of 2, thus, there are 225 ways to form 3 alpha particles. And only 1 of them will finally get out. Will these multiplicity increase the chance for tunneling and reduced the lifetime?

SASAMOTO San studied the excited state. And she knew somebody studied the alpha decay. Thus, by compare the decay rate with other nuclei, may be we can find out the role of multiplicity in decay.

The ground state of alpha particle is spin zero. And a spin 1 state of it is about 20MeV. However, adding one more proton will greatly reduced the energy. What is going on in alpha particle?

Method II (decay)

Leave a comment

there is another way to study nuclear physics, which is by observe the decay process.

there are 3 major decays, the alpha, beta and gamma. the alpha decay is an excited nucleus go to a lower energy state by emitting an Helium ( 2 protons and 2 neutrons) nucleus. This process change the nucleus constitution and make it lighter. the beta decay is an excited nucleus go to a lower energy state by emitting an electron ( or positron, or capture an electron). This process also change the nucleus constitution but the change in mass is very small. the gamma decay is an excited nucleus go to a lower energy state by emitting a photon ( or the light particle). this process does not change the nucleus constitution.

there are many other decays, like neutron, proton, and even fission can regard as decay. in general, decay is a process that a nucleus go to a lower energy state and become more stable. But other decays has too short for lifetime, so that after the earth was formed, they are almost gone. and only alpha, beta and gamma has long left time to survive.

Via a decay process, we can know the energy level, lifetime, parity of each level of a nucleus. These informations can help us to build nuclear models and theories to understand and predict nuclear properties, like the strong and weak force.

i will present a series about decay in future posts. watch out!

symmetry breaking, mass and higg field.

1 Comment

Above curie temperature, the spin of iron is isotropic. The spin can be rotated to any direction without any resistance, they like massless. Below curie temp, the iron has phase transition and all spin now point to a particular direction. And we need some force to rotate the spin direction. The spin has mass now. This is what symmetry breaking in simple manner.

When matter becomes superconduct, the magnetic field inside is decay exponentially,which is similar as Yukawa force. And we said the force carrier particle is massive. The magnetic field decay is due to the copper pair, Which respond to the magnetic field and tend to cancel it. Thus, the direction of copper pair is not isotropic and this is another symmetry breaking due to the external field and low temperature.

At very high temperature, the weak force carriers are massless. And we assign an isotropic field (scaler field) for the force carrier and call it Higg field. The Higg field quata is called Higg boson. It act like the copper pair, which respond with the force carrier. When there is a force carrier, a Higg boson will be induced. And the symmetric breaking in Higg field, the symmetry breaking makes the force carrier has mass. That we need to apply a force to change the motion.


Leave a comment

On 11st, March, a M9.0 earthquake strike east Japan around 3pm. this caused malfunction of the nuclear plant, Fukushima Daiichi.

i am not going to talk about the malfunction, but on the radioactivity.

The radioactivity is caused by nucleus decay into a more stable state. like Carbon-14 will decay to Carbon-13 by beta decay. there are 3 common types of radioactivities: alpha – a Helium nucleus, beta – a electron and gamma – high energy light.

the radioactivity was measured by Bq ( Becquerel ), which is 1 decay per second. for a typical radioactivity substance, there may be 10 billion decay per second. think about how many atom in 1 mole of matter.

each radiation carry some energy, for alpha decay, energy is about 5MeV, beta decay is 1 MeV, and for gamma decay, it is around 0.1 to 1 MeV. 1 MeV is about 1.60217653 \times 10^{-13} Joule. for boiling a water to increase 1 degree, it take 2400Joule. you can see how small it is.

1 more thing to consider is, not every radiation will be absorbed. some will just pass through and leave and no effect. the chance for being absorbed is proportion to the density of the material, or for simplicity, how much weight of the material.

for easy understanding, we multiple these 3 factors together, to give us how much energy absorbed per kg. and we call it Gy (Gray), the absolute does.  and we can have Gy per hour, Gy per second, depends on the time scale.

Thus, for 1Gy, an alpha source, typical 5MeV per decay and a total absorption, it takes 6241 billion hits. for 10 billion decay a second, it takes 10 mins for this does.

in other system, people use rad for radioactivity, for example in the game fallout. which is 1 Gy = 100 rad.

Here is a list of different does and the health effect

  • 0.25 Gy or less  –  no effect
  • 0.50 Gy  – temporary reduction of while blood cells
  • 1 Gy  –  feel nausea, vomit
  • 1.5 Gy – 50%  more chance to have cancer, cataracts or leukemia
  • 2 Gy – 5% people will die
  • 4 Gy – 50% people will die in 30 days
  • 6 Gy – 90% people will die in 14 days
  • 7 Gy – all die

However, different part of our body have different does limit. so, for 1Gy radiation, our eye will be blinded, but our skin has no effect. therefore, we introduce another does, called Effective does Sv ( Seivert ). it is produced by a weight factor.

Gy \times w = Sv

different part of out body has different weight factor, for lower weight factor, the more easy the tissue will get hurt.

here is the list of weight factor

  • whole body  – 1
  • skin — 0.01
  • lung — 0.12
  • liver — 0.05
  • breast — 0.05
  • stomach — 0.12
  • gonads – 2

we see that the weight factor for whole body is 1, therefore, in general we use Gy or Sv in equal way.

And we also observed that, when we eat some radioactive matter inside, the limit does will be decreased by 8 times in out stomach. or to say, the radioactive matter is 8 times dangerous than it is at outside. the reason is, when the matter is inside our body, we will absorb all energy.

a typical does is 3.75mSv per year. which means we received 3.75mSv in a year, 0.4uSv per hours.

for radioactive effect, there are 3 types of effect: acute ( short term ), late ( long term) and genetic effect ( permanent ) .

  • short term – Erythema, hair loss, etc
  • long term – Cateracts, Leukemia or cancer, etc

and the effects can be divided into 2 groups – deterministic and stochastic

  • deterministic mean the chance of having those effect is increase with the radiation increased.
  • stochastic is similar to deterministic, but there is a threshold. below the threshold there is no effect. – Cancer and genetic effect are belong to this group.

[Pol. p target] Found water NMR signal

Leave a comment

the DC B-field voltage is 3.100V. i found that i have have to tune the voltage slowly, to let the CRO signal stabilized.

The NMR setting is 12.8MHz, Level 200, pule length is 5uS on +Y axis. FWD is 49W, BWD is 7W. Gain is 35db, LPF is 100kHz and polarization is 90 degree.

the FID signal is about 200mV peak to peak.

Let me briefly explain the NMR circuit, for reader who want to know the theory. HERE.

Black arrows are any cable, Orange arrows are 12.8MHz cables. The direction of the arrow shows the signal flow.

the NMR pulse is a mixture of 2 signal. the RF generator provides the larmor frequency. and the pulse generator provides a square wave, the duration of the square determine the pulse time, it should be matched the requirement for π/2 pulse. and the period of the pulse should be longer than the T2, longitudinal relaxation time. since we don’t want to either distorted the RF or use 1 more RF generator, the splitter simply split the signal, 1 be a main signal and the other one be a reference signal.

The purple mixer or Double Balanced Mixer (DBM). the function of it is multiple the 2 input signal together. the circuit inside, may be i study later.

after the Mixer, the signal will be amplified and pass a λ/4 cable. i am not sure is it necessary to use a λ/4 cable in here. but anyway, we knew that the phase different from the Amplifier to the crossed-diode is λ/4.

the crossed diode has a feature that only allow signal, which larger then 0.6 V to pass. when the pulse come in, it will not see the diode. however, when the in put signal are off, the signal from the NMR coil will be generated and this signal is smaller then 0.6V, thus, the crossed diode act as a switch or duplexer.

the λ/4 cable from the crossed diode and the NMR tuner + Coil is for impedance matching. even there are little bit mismatch of the tuner, the λ/4 cable is still approximately equal to the impedance matching case. that smoothen the signal.

and i don’t know is it necessary to use  a λ/4 cable from the cross diode to the amplifier.

after the amplifier increase the signal and mixed with the reference signal at the differential mixer. The Free induction decay (FID) signal should be at the Larmor frequency. However, if the input frequency is not exactly, with just little bit of mismatch, the spin can still be rotated. by different the NMR signal with the reference signal ( which is the signal we put in the NMR coil) we can tell the different and calculate back the true Larmor frequency by Fourier transform.

and also, by different the NMR signal by the reference signal, we can more easy to tell the different then by just looking at the NMR signal. since the different usually just few kHz, compare with the Larmor frequency 12.8MHz, less then 1 % deferent.

the different mixer also decompose the signal into sine and cosine waveform. this is because, by only look at the cosine wave, we cannot know the different of the signal is positive of negative. the cos wave is natural real part of the complex wave.

the low pass filter (LPF) can be set the maximum frequency pass. this is for seeing the decay envelop of the FID. anther function is reduce the high frequency noise.

Older Entries