HGCal Test Beam Documentation

Table of Contents

• Input files:
• Reconstruction sequence
• Running the DQM plotting
• To be updated:

Download instructions: Instructions

Input files:

• channel mapping: file that contains the conversion between electronicID (SKIROC chip number) to detector ID (DetId)

  • HGCal/CondObjects/data/map_FNAL_2.txt 

• pedestals: pedestal values measured from pedestal run with dedicated analysis (PEDESTAL_ANALYSIS_ to be updated soon)

  • HGCal/CondObjects/data/Ped_HighGain_8272_Mean.txt 

  • The above file is a sample file in the format needed to be read by the RecHitProducer. This is valid for the high gain ADC and was derived from the pedestal run 8272 

• gains: ADCtoGeV conversion factors measured with dedicated analysis (ADCTOGEV_ to be updated soon)

  • N/A yet 

Pedestals and ADCtoGeV should be provided for both high gain ADC and low gain ADC in the format of the above .txt file in lines 38-41 of Reco/plugins/HGCalTBRecHitProducer.cc Now only pedestals for high gain is available. Default values if files are not provided are:

• 0 for pedestals
• 1 for ADCtoGeV

Run, event, time stamp information has been added to edm::event and can be accessed directly per event.

Reconstruction sequence

The reconstruction sequences are defined in StandardSequences/python

• RawToDigiSeq: This sequence reads the RAW output files and assigns the "decoded" adc values(both high and low gain) from individual channels of the SKIROC to the corresponding cells of the sensor.
• LocalRecoSeq: This sequence starts with the Digis above, subtracts from each cell of the sensor the pedestals as provided in the txt files mentioned above, and multiplies the resultant with the ADCtoGeV scale factor. -SimulationSeq: An example is not yet provided, but this will not have to go through the above sequences and we can directly invoke the plotting tools.

The Raw Data files: Can be found in /afs/cern.ch/user/r/rchatter/public in folders: Electron_Runs_2932016, Electron_Runs_3032016 , Electron_Runs_3132016, Proton_Runs_0142016, Electron_Runs_0142016, Proton_Runs_0242016, Proton_Runs_0442016. The run information can be found on elog http://cmsonline.cern.ch/portal/page/portal/CMS%20online%20system/Elog?_piref815_429145_815_429142_429142.strutsAction=%2FviewSubcatMessages.do%3FcatId%3D592%26subId%3D17%26page%3D1 . Please note that the folder names are a little misleading and the beam type should currently be read off from the elog. These files will be made available on CERN eos soon and a message posted on hn.

Running the DQM plotting

Currently the DQM worflow kicks in post the creation of Digis and also after the creation of RecHits.

• Post Digi we can monitor:
  • adc counts in each cell of a sensor per sample, per layer(summed over the events).
  • 1D adc distributions per cell(per layer, per sample).
• Post RecHit we can monitor:
  • RecHits(GeV) in each cell of a sensor, per layer(summed over the events).
  • The Occupancy(defined as RecHits above a certain threshold) in each cell of a sensor, per layer(summed over the events).
  • 1D RecHits distributions per cell(per layer, per sample).
  • 1D RecHits distributions summed over all cells in a layer(per layer).

On executing:

cmsRun test_cfg.py

a root file test_output.root is produced. This file has two directories namely

• hgcaltbdigisplotter: with the Digi->DQM workflow
• hgcaltbrechitsplotter: with the RecHits to DQM workflow

In the Digi->DQM workflow the following histograms may be found:

• (TH2Poly) FullLayer_Sample::M_Layer::N is the cumulative(over events) adc count in cells of the sensor in the #N layer for the #M sample.
• (TH1F) Cell_u_::P_v_::Q_Sample::M_Layer::N is the 1D adc spectrum for the sample #M in a cell in layer #N marked by u=#P and v=#Q in the local(u,v) co-ordinate system of the cells in a sensor(Not to be confused with the Sensor_U, Sensor_V used to identify sensors in a layer).

In the RecHits->DQM workflow the following histograms may be found:

• (TH2Poly) FullLayer_RecHits_Layer::N is the cumulative(over events) RecHits in cells of the sensor in the #N layer
• (TH2Poly) FullLayer_Occupancy_Layer::N is the cumulative Occupancy(as defined above) in cells of the sensor in the #N layer
• (TH1F) Cell_RecHits_u_::P_v_::Q_Layer::N is the 1D RecHits spectrum for a cell in layer #N marked by u=#P and v=#Q in the local(u,v) co-ordinate system of the cells in a sensor
• (TH1F) FullLayer_RecHits_Layer::N_Summed same as above but summed over all cells in a sensor.

Coordinate system

Git instructions: GIT SHORT instructions

To be updated:

test recHits