Differences between revisions 4 and 6 (spanning 2 versions)
 ⇤ ← Revision 4 as of 2017-09-15 14:57:31 → Size: 12274 Editor: KatjaKrueger Comment: ← Revision 6 as of 2017-09-15 15:07:53 → ⇥ Size: 11870 Editor: KatjaKrueger Comment: Deletions are marked like this. Additions are marked like this. Line 7: Line 7: '''runs on flchcallab5''' Line 92: Line 94: ==== Hit Profiles per layer ==== ==== Number of hits per layer ==== Line 96: Line 98: ==== Energy ====[[attachment:Energy_sum.png|{{attachment:Energy_sum.png|Energy Sum|height="250"}}]] [[attachment:Energy_per_layer.png|{{attachment:Energy_per_layer.png|Energy Per Layer|height="250"}}]] ==== Total Energy ====[[attachment:Energy_sum_pion.png|{{attachment:Energy_sum_pion.png|Energy Sum Pions|height="400"}}]][[attachment:Energy_sum_muon.png|{{attachment:Energy_sum_muon.png|Energy Sum Muons|height="400"}}]]==== Hit Energy (of all cells within a layer) ====[[attachment:HitEnergy_per_layer_pion.png|{{attachment:HitEnergy_per_layer_pion.png|HitEnergy Per Layer Pions|height="400"}}]][[attachment:HitEnergy_per_layer_muon.png|{{attachment:HitEnergy_per_layer_muon.png|HitEnergy Per Layer Muons|height="400"}}]] Line 100: Line 108: Muons +120 GeV, tungsten absorber[[attachment:shower_muon120_tungsten.png|{{attachment:shower_muon120_tungsten.png|shower_muon120_tungsten|height="250"}}]]Pions +50 GeV, tungsten absorber[[attachment:shower_pion50_tungsten.png|{{attachment:shower_pion50_tungsten.png|shower_pion50_tungsten|height="250"}}]]Electrons +20 GeV, tungsten absorber[[attachment:shower_electron20_tungsten.png|{{attachment:shower_electron20_tungsten.png|shower_electro20_tungsten|height="250"}}]] [[attachment:Shower_pion.png|{{attachment:Shower_pion.png|Shower variables Pions|height="400"}}]][[attachment:Shower_muon.png|{{attachment:Shower_muon.png|Shower variables Muons|height="400"}}]] Line 113: Line 112: You can see some reference hitmaps from August 2015 data taking period.For Muons, +120 GeV, tungsten absorber.[[attachment:muon120GeV_tungsten_hitmap.png|{{attachment:muon120GeV_tungsten_hitmap.png|Hit Map Muon 120 GeV, tungsten absorber|height="250"}}]]For Pions +50 GeV, tungsten absorber.[[attachment:pion50GeV_tungsten_hitmap.png|{{attachment:pion50GeV_tungsten_hitmap.png|Hit Map Pion 50 GeV, tungsten absorber|height="250"}}]]For Electrons +20 GeV, tungsten absorber.[[attachment:electron20GeV_tungsten_hitmap.png|{{attachment:electron20GeV_tungsten_hitmap.png|Hit Map electron 20 GeV, tungsten absorber|height="250"}}]] [[attachment:Hitmap_pion.png|{{attachment:Hitmap_pion.png|Hit Map Pions|height="400"}}]][[attachment:Hitmap_muon.png|{{attachment:Hitmap_muon.png|Hit Map Muons|height="400"}}]]

# Online Monitoring and reconstruction for AHCAL Test Beam Campaign at CERN2017

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, right ones for a muon run

### 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">
</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">
</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"`

AHCALandHGCALTestBeamCERN2017/Monitoring (last edited 2017-09-15 15:10:18 by KatjaKrueger)