EIS Team Science Meeting Agenda#
18th March
This page refers to one day of a closed meeting of the EIS consortium
1 Technical#
1.1 Report from Technical Splinter Meeting (previous day)#
Khalid reported on the instrument status, w.r.t. PZT non-linearity.
Also discussed the statistics of observations carried out to date.
Remote operations initial testing and model also reported on.
Another instrumental issue is the optimisation of the focus within the instrument.
Harry Warren was called on to talk about the slit and grating focus.
Google was apparently able to show that the grating focus happened in Aug 2008.
- Post-SLA manoeuvre, the slit images from each channel overlap to 0.25 arcsec. An explosive event observed would help to prove the alignment.
KFJ: Just to clarify there was an offset between the two channels, and now they appear simultaneously on the two detectors.
HPW: True. We haven't quite seen dynamic events
GAD: What about images on the same detector. Wasn't there a slight offset between Fe VIII and Fe XII.
HPW: in my analysis, the offset is pretty well taken care of with PRY's formula. Extremely well correlated once you account for that.
- Grating focus. The grating was driven back and forth to find the best focus position, and the images moved a lot, and consequently the dispersion relation constant changed a lot(!). But then we moved it back. HPW took histograms of the line width. We thought it might improve. At the new position, though, it's unchanged. At the extreme position, there was no usable data because the grating had moved the image so far away :-)
Linear trend in the increase of the number of hot pixels on the detector with time. These are in the CCD detectors, approximately 1% of pixels.
- Warm pixels have mean + (5 sigma) to mean + (50 sigma) pedestal values. The rate of increase has had multiple level-offs. After the levelling off in 2008, there was an increase again in July, then levelled off in December. Seems to be a temperature dependence throughout Earth's orbit about the Sun (January is hottest). By 2010, 16% of pixels may be affected; 2012, 26% might be affected. Call on PRY to present happier news about this, though...
Sounds like a huge fraction of the CCD.
What's the best way to treat warm pixels for scientific analysis. My thought from the start of the mission was to ignore them, and neglect them in the analysis. But HPW and PRY talked and found that interpolating actually gave good results.
Wrote a document to be posted on this Wiki after the meeting presenting some of the results.
Used the standard EIS_PREP processing.
Then artificially inserted 30% bad pixels. So how badly were the fits degraded by ignoring these fake "bad pixels".
HEM: Interpolatin after fitting?
PRY: This is done by EIS_PREP. And is done in the solar_Y direction
KPD: How do you know they're better fits?
PRY: I have the original data, and identify the places where there never were bad pixels. Then I process these data by inserting fake bad pixels, interpolate, and then compare the line fit parameters.
KPD: I think you're inventing data...
PRY: That was my worry, too, but it seems to work!
HEM: What about clusters of missing pixels?
PRY: There are different methods of tackling these.
You could reproduce the original data 97%, but the bad interpolations tend to happen in the high-intensity areas. My recommendation is to use interpolated data rather than just ignoring the missing pixels. I was amazed how well you could do even with 30% warm pixels, but it seems we can continue longer than expected by a couple, even if we don't remove them.
It may be related to EIS's oversampling w.r.t. its resolution, so there is correlation of data which allows some justification of the interpolation: there's a contribution from neighbouring pixels anyway.
- Sensitivity of the instrument over time (KFJ)
We use flat-field LEDs inside the instrument, and we also measure the solar radiation variation over time. We don't seem to be seeing substantial degradation in our sensitivity.
Contamination is also measured with Quartz Crystal Mircobalance sensors, which show no appreciable increase in contamination so far.
PRY: how about the event flags? Are they implemented?
KFJ: Implemented in GSW, but not tested yet. They need testing on the Sun. BP trigger could be tested soon.
PRY: Is there a problem with remote planning? We'll be more cautious about doing that?
KFJ: If we sit together remotely, it may be easier.
JTM: are there problems with Data Volume
DRW: It may be easier in the current system.
JLC: How about use of the XRT flare flag? Do we have plans to use it?
HH: it is posible for us to use it
HPW: He II may be better than XRT
HEM: What does the flag give us? Is it coordinates?
HH: Yes, coordinates.
Discussion followed about how to test this flag.
1.2 Analysis Software#
1.2.1 Status of the reduction package#
- 10:00 EIS_PREP (Alessandro Gardini)
AG went through the ordering of calibration steps.
- Dark current subtraction and options for this
- Marking hot pixels affected by CRs
- Hot pixel flagging
- Warm pixel flagging
- Dusty pixel flagging
- Absolute calibration (specific intensity units; photons)
The 1st step returns missing pixels, but they aren't filled in by other values. To do this, you have to allow the CR option (step 2?).
PRY: If I don't use /DEFAULT what happens?
AG: If you don't select default, it should ask you to choose the DC file.
Computation of moments is available in XFILES -> XCONTROL -> Line Fit
Can't construct moments in anything other than absolute units. Also doesn't work with negative values retained.
Can then, in XFILES, use moments OR the Gaussian fitting routine supplied by PRY.
HEM: This is assuming no blends and a single line per window?
AG: When you select a line, you have to then specify the line and continuum interactively. If you have two lines very close to each other, you can choose only a few pixels plus continuum to choose just that line. But if there is a blend then it won't work.
After the line and contin. are selected, the 0th, 1st, and 2nd order moments are computed. The last term 2*ln(2) is necessary to turn sigma into the FWHM. That's it.
XFILES also gives you the errors on these moments for intensity and velocity. The error on the line width hasn't been implemented yet.
HEM: is there documentation for this yet? What you said is really useful, but is there a written version of it.
Brief discussion followed about the desire for tutorials and guides.
1.2.2 Useful additional software#
(With descriptions or demonstrations)
- 10:15 Gaussian fitting routines for EIS (Peter Young)
EIS_AUTO_FIT and EIS_AUTO_FIT_GEN (newer) both available.
When you want to fit multiple gaussians the latter is better (?)
You might want to add many spatial pixels together to make a total or mean spectrum.
You can select sub-regions of the array.
eis_auto_fit, windata, fitdata, refwvl=195.12
- First choose the FITS file I want to look at.
- I then choose a pre-calibrated file
- extract the window data with EIS_GET_WINDATA
- choose one at 264 Å
- eis_mk_fit_template averages the spectrum over space
- go to the window, select the lines roughly and the continuum
- wanted to speed up the fitting for this talk, so...
- have a routine that takes the windata and rebins it by some factor which gives you that factor^2 less pixels to fit
- then read fit template eis_read_fit_template
- then start the fitting process, which takes a while, meantime...
The fit template stores the initial guesses as a text file which you can modify if you want to.
GAD: How accurate does it have to be?
PRY: Not too accurate. In this case the lines are quite widely spaced,
GAD: But how about if they're close together
PB: in the case where there is a blend, then in some cases I found I had to restrict the parameters, e.g. fixing the widths to be the same
PRY: You could try the lines having a narrow range of width?
PB: what I did was to force the widths to be equal, which meant I had to change the code a bit
PRY: sometimes you have to write a more specialised routine, because it's difficult to write something completely general?
YKK: Is the continuum fixed as linear?
PRY: yes, linear or constant
okay, the fit has finished:
- when yo do the orbital correction, have to specify a line -- in this case line 1.
- you can then view the quality of the fit with FIT_VIEWER
- see the I, V and dl maps. Click on a position and you see the spectrum AND fit.
It's available in SSW, and there's a tutorial on the web to show you how to use it.
That was the Gaussian fitting.
I also mentioned that sometimes you see a BP, for example, and you want to add all the pixels' spectra together. Arbitrary spatial pixels. It's complicated by the grating tilt and inter-channel offset, but that can be taken care of with this software.
I'll use the same dataset as before. Have a simple routine that makes an image in the specified line. eis_make_image
eis_pixel_mask to make the mask...
how is the pixel mask stored?
This performs the averaging and gives you a Channel A (long) and Channel B (short) spectrum.
In spec_gauss you can see the data quality with the yellow line (something to do with number of missing pixels?)
Just wanted to show the Gaussian fitting routines I wrote.
- Auto fits single
- Auto fits multiple
- Spatially averaged spectra
uses EIS_CCD_OFFSET for each wavelength in the spectrum to compensate for intra-channel vertical offsets.
- 10:35 A quick display method (David Williams)
Coffee#
10:40 to 11:00
1.3 Choosing data#
- 11:00 Discussion on Data access
There are huge volumes of EIS data, but it would be good to have a discussion on how to go about accessing the kind of data you're looking for.
It seems like it's possible to cross-match the EIS and RHESSI catalogues already. RM has done this. And HPW will talk to SDO AIA people to discuss their requirements for integrating our information with the Heliospheric knowledge base.
Ultimately, though, it might be good to set up ways to make event catalogues; it could be possible for people to construct and contribute custom cataloguing/data-mining routines.
Seems like we have issues more with the Wiki than with Data Access per se, though.
- HEM suggested Nuggets à la Hudson for RHESSI.
- KPD felt that we weren't advertising or communicating enough.
- DRW felt that it was easy and good to put information on the EISWiki.
2 Science#
2.1 EIS Science Achievements#
Each national PI will give their perspective on the scientific achievements of EIS in its first 2.5 years. These are personal perspectives on where we stand on two things:
- phenomena and scientific issues that were known before launch
- the discoveries Hinode has made with EIS
( they won't necessarily agree :-) )
Correlation of Doppler speed and non-thermal width implies unresolved flows in acrive regions.
Magnetic reconnection, if it at all takes place for coronal heating, is far below the height of the reconnection point that causes major flares.
Recent results from DHB: even if we look at 1 arcsec pixel of Quiet Sun, we usually have the same shape of EM distributions. Unresolved very common nature of a plasma temperature distribution that causes this universal shape of DEM even in a size less than 1".
Triggered by Sakao et al.'s observations, we observe high-velocity outflows, although they're slow in the solar wind sense. We've actually found plasma upflows taking place close to coronal holes and quiet Sun beside active regions. The same sort of mechanism took place the data analysed by Imada after the largest 20061213 flare. In quiet regions, the kinetic energy is roughly enough to heat the QS corona.
20061217 LDE was analysed by Hara: Reconnection inflows seen in LOS velocity by HH;
20060519 flare: EIS counterpart of a Masuda source? The EIS slit was very lucily at the flaring loop top when its initial phase started. TRACE shows a clear two-ribbon pattern; RHESSI image contours for 12 - 25 keV which show a high-temperature RHESSI structure. Fe XXIII emission corresponding to it shows a slight red-shift.
- 11:50 Louise Harra
- 12:10 George Doschek
Lunch#
12:30 to 13:30
2.2 EIS Science Results#
Speakers are encouraged to put special emphasis on unsolved components in the context of what they talk about, and how they can be addressed using all three instruments on Hinode.
2.2.1 Active Regions#
- 13:30 Cambridge active region studies (Helen Mason)
- 13:45 The thermal structure of active regions (Harry Warren)
- 14:00 Multi-component active region flows (Paul Bryans)
- 14:15 Flows & motions in moss in the core of a flaring active region (David Brooks)
2.2.2 Flares#
- 14:30 Evaporated Plasma (Ryan Milligan)
Science with other, non-Hinode instruments#
I’d also like to solicit talks on any of the wealth of cross-observatory data taken over Hinode’s life so far (Hinode-SUMER campaigns, regular CDS co-observing, TRACE co-observing, WHI 2008, etc.) What projects do we have (under-way or recently completed) with these multi-instrument observations that couldn’t otherwise be tackled?
2.3 Where to next for EIS science?#
14:45 to 15:15
Not only having listened to the contributions given this afternoon, but also having considered our progress in achieving EIS's main science goals — both pre- and post-launch — it would be good to have a discussion before and after coffee on what burning questions we think are obviously being left unanswered.
I'd like to budget 30 minutes for this.
Afternoon Tea#
15:15
2.4 Analysis Techniques#
- 15:30 Line identifications with EIS and CHIANTI (Peter Young)
- 15:50 Diagnostics with the slot (Ignacio Ugarte)
2.4 Observing Techniques & Strategies#
What we’ve learned to date on observing with the S-band antenna
- 16:10 Some techniques for observing (David Williams)
- 16:20 Raster stitching (David Brooks)
- Problems
- Some examples of studies that have proven useful in overcoming these problems
- Prospects & challenges for the future:
- How do we address observing challenging targets like flares?
- 16:25 Flares (Ryan Milligan)
- 16:45 Very active regions (Harry Warren / John Mariska)
- 17:05 Observation Planning - Core Team proposals and HOPS (Len Culhane)
4 EIS Website & Wiki#
17:15
These are our public face for information. It would be good to talk about:
- Changes made since last year
- How to increase participation
- Good practice
- Keeping the contents up-to-date
- New types of content
5 Action Items and Wrap-up#
17:45
6 終わり#
18:00 END
Confirmed Attendees#
- Khalid al-Janabi (MSSL)
- Danielle Bewsher (RAL)
- David Brooks (NRL)
- Charlie Brown (NRL)
- Paul Bryans (NRL)
- Len Culhane (MSSL)
- Ken Dere (GMU)
- George Doschek (US PI; NRL)
- Alessandro Gardini (UiO)
- Hirohisa Hara (NAOJ)
- Louise Harra (PI; MSSL)
- Shinsuke Imada (NAOJ)
- Yuan-Kuen Ko (NRL)
- John Mariska (NRL)
- Helen Mason (DAMTP)
- Keiichi Matsuzaki (ISAS)
- Ryan Milligan (GSFC)
- Karin Muglach (NRL)
- Steve Myers (NRL)
- John Rainnie (RAL)
- Toshifumi Shimizu (ISAS)
- Ignacio Ugarte Urra (NRL)
- Harry Warren (NRL)
- Tetsuya Watanabe (Japan PI; NAOJ)
- David Williams (MSSL)
- Peter Young (NRL)
Back to the
EIS Team Meeting top page.