Online Monitoring and reconstruction for AHCAL Test Beam Campaign at CERN2017
This page is a branch of the main AHCALandHGCALTestBeamCERN2017 page.
Contents
Contents
Testbeam Monitoring
runs on flchcallab5
Tunnel to calice DB at DESY
At the moment the monitoring is set up to connect to the database at DESY (to make sure that this is the case check /home/calice/Desktop/AHCAL-Reconstruction/xml/Reco.xml, where you should see the following lines:
<parameter name="DBInit" type="string" value="127.0.0.1:calice:caliceon:Delice.1:33306"/>
(using port 33306 on 127.0.0.1)
To set up the IP tunnel to DESY, use the following command (with your user name, on the machine where the online monitoring runs):
ssh -L 33306:flccaldb02:3306 user@bastion.desy.de
The programs are in
/home/calice/Desktop/AHCAL-Reconstruction/bin
We reconstruct and apply the available calibrations to the raw events:
QtReco (reconstruction):
1. Run it by doing
./QtReco_interface
2. Edit the steering file in the folder /home/calice/Desktop/AHCAL-Reconstruction/xml (steering.xml)
Run number Input_PATH Output_PATH mapping
(Run number needs to be edited for each new Run, the rest should stay the same)
3. Then Press Configure -> Check if the lcio file is found.
4. Then Press Start, it will run the reco and the analysis.
5. Wait until it finishes. It should show Analysis Module : destroyed
Now we are ready to show the plots or event displays.
Online monitor for events after reconstruction
It is in:
/home/calice/Desktop/AHCAL-Reconstruction/bin
1. Run it by doing
./QtReco_Client_interface
QtReco-interface must be opened!
2. Click on connect then enter the Run number with a leading zero!. Check that it says connected.
3. Then choose what plots you want to look at. Careful when analysing the plots! Mis-calibrations!
Important (and implemented) distributions
Several distributions are available after analysis :
- Hit multiplicity
- Energy sum and spectra per layer
- Shower profile, beam profile
- Hit Map
- Temperature (only for raw txt!)
left plots are for 300 GeV pions in July 2017, right ones for a muon run in August 2017
Hit Profiles per layer
Number of hits per layer
Total Energy
Hit Energy (of all cells within a layer)
Shower profile, beam profile
HitMap
Calibration procedure
The QtReco program make use of the Marlin framework and the calice processors in order to reconstruct the raw data (hardware information) to physics data. The reconstruction is done in several steps and need to have access to the database containing mapping of the detector and calibration constants for each channel (gain, Pedestal, MIP).
The reconstruction is done in 4 steps :
- Event Building from raw txt or raw slcio eudaq files
- Geometry mapping of the detector
- Get Calibration constant from database/files
- Apply Calibration
This is done for each hits in each events.
Database constants (DESY)
<processor name="GeoConditions" type="ConditionsProcessor"> <parameter name="DBInit" type="string" value="127.0.0.1:calice:caliceon:Delice.1:33306"/> <parameter name="DBCondHandler" type="StringVec"> Ahc2ModuleDescription /cd_calice_Ahc2/TestbeamJuly2017/ModuleDescription HEAD Ahc2ModuleConnection /cd_calice_Ahc2/TestbeamJuly2017/ModuleConnection HEAD Ahc2ModuleLocationReference /cd_calice_Ahc2/TestbeamJuly2017/ModuleLocationReference HEAD Ahc2DetectorTransformation /cd_calice_Ahc2/TestbeamJuly2017/DetectorTransformation HEAD Ahc2HardwareConnection /cd_calice_Ahc2/TestbeamJuly2017/Ahc2HardwareConnection HEAD E4DPedestal /cd_calice_Ahc2/TestbeamJuly2017/Pedestal HEAD E4DGainConstants /cd_calice_Ahc2/TestbeamJuly2017/gain_constants HEAD E4DGainSlopes /cd_calice_Ahc2/TestbeamJuly2017/gain_slopes HEAD E4DMipConstants /cd_calice_Ahc2/TestbeamJuly2017/mip_constants HEAD E4DMipSlopes /cd_calice_Ahc2/TestbeamJuly2017/mip_slopes HEAD E4DDeadCellMap /cd_calice_Ahc2/TestbeamJuly2017/DeadCellMap HEAD E4DSaturationParameters /cd_calice_Ahc2/TestbeamJuly2017/SaturationParameters HEAD E4DIntercalibration /cd_calice_Ahc2/TestbeamJuly2017/Intercalibration HEAD E4DPhysicsCalibIntercalibration /cd_calice_Ahc2/TestbeamJuly2017/PhysicsCalibIntercalibration HEAD E4DTimeSlopes /cd_calice_Ahc2/TestbeamJuly2017/TimeSlopes HEAD E4DTimePedestal /cd_calice_Ahc2/TestbeamJuly2017/TimePedestal HEAD </parameter> </processor>
Event display
Kataj: please ignore for the moment, needs to be checked!
After applying the reconstruction and the corrections, we also can have a look on the aspect of the events in our detector. Actually, we can do it in parallel to the reconstruction, since it reads the events in the slcio file created by the QtReco program. For that we use a schematic 3D view of our detector that allows us to nicely and very intuitively discern the structure of the beam (content of pions, muons, electrons, etc), the position of the tracks, etc.
To execute the event display go in :
~/Desktop/Event_Display/
and start
glced& ./myMarlin.sh xml/CED_2017May.xml
To ship the data to the online display from a correct file, xml steering file /home/calice/Desktop/Event_Display/xml/CED_2017May.xml has to be edited. File which will be displayd is here::
<parameter name="LCIOInputFiles"> <!--/home/calice/Desktop/Reco_Data/BuiltEvents000025_tp.slcio /--> /home/calice/Desktop/Reco_Data/ROCEvents050917_tp.slcio </parameter>
and database connection has to be fixed in a same way as for reconstruction:
<processor name="GeoConditions" type="ConditionsProcessor"> <!--flccaldb02 --> <parameter name="DBInit" type="string" value="127.0.0.1:calice:caliceon:Delice.1:33306"/> <parameter name="DBCondHandler" type="StringVec"> Ahc2ModuleDescription /cd_calice_Ahc2/TestbeamMay2017/ModuleDescription HEAD Ahc2ModuleConnection /cd_calice_Ahc2/TestbeamMay2017/ModuleConnection HEAD Ahc2ModuleLocationReference /cd_calice_Ahc2/TestbeamMay2017/ModuleLocationReference HEAD Ahc2DetectorTransformation /cd_calice_Ahc2/TestbeamMay2017/DetectorTransformation HEAD Ahc2HardwareConnection /cd_calice_Ahc2/TestbeamMay2017/Ahc2HardwareConnection HEAD E4DPedestal /cd_calice_Ahc2/TestbeamMay2017/Pedestal HEAD E4DGainConstants /cd_calice_Ahc2/TestbeamMay2017/gain_constants HEAD E4DGainSlopes /cd_calice_Ahc2/TestbeamMay2017/gain_slopes HEAD E4DMipConstants /cd_calice_Ahc2/TestbeamMay2017/mip_constants HEAD E4DMipSlopes /cd_calice_Ahc2/TestbeamMay2017/mip_slopes HEAD E4DDeadCellMap /cd_calice_Ahc2/TestbeamMay2017/DeadCellMap HEAD E4DSaturationParameters /cd_calice_Ahc2/TestbeamMay2017/SaturationParameters HEAD E4DIntercalibration /cd_calice_Ahc2/TestbeamMay2017/Intercalibration HEAD E4DPhysicsCalibIntercalibration /cd_calice_Ahc2/TestbeamMay2017/PhysicsCalibIntercalibration HEAD E4DTimeSlopes /cd_calice_Ahc2/TestbeamMay2017/TimeSlopes HEAD E4DTimePedestal /cd_calice_Ahc2/TestbeamMay2017/TimePedestal HEAD </parameter> </processor>
The program has a tool to make screenshots. This is found in the "tools" tab and do screenshots in tga format.
Some standard event display are shown in the following pictures.
Electron (-) candidates (10 GeV, steel absorber)
Muon (-) candidate ( steel absorber)
Pion (-) candidate (50 GeV, steel absorber)
Pion (+) candidates (50 GeV, tungsten absorber)
Electron (+) candidates (20 GeV, tungsten absorber)
To make these nice sets of eventdisplays, we can use a few set of "Jiri's" commands:
1) Convert to png, since the saved pictures in the eventdisplay use the tga format
ls *.tga|parallel convert {} ~ {.}.png
2) Create the set, by appending individual pictures:
convert \( glced-1.png glced-2.png +append \) \( glced-3.png glced-4.png +append \) -append test.png (for 2x2) or convert \( glced-1.png glced-2.png glced-3.png +append +append \) \( glced-4.png glced-5.png glced-6.png +append +append \) \( glced-7.png glced-8.png glced-9.png \) -append -append r30079_pick_of_9events.png (for 3x3)
dCache
The procedure of copy to DCache is described here: http://flcwiki.desy.de/GridInstall
On the DCache computer, a Virtual machine is automatically started with the PC. To invoke the program, following command is needed:
ssh -X -p 2222 calice@localhost "cd dCache-gui; bin/dCache_start_gui"