# Presentation: 562138

## Results from HCAL TB2004&Inputs for P-TDR

** **
**Shuichi Kunori**
**U. of Maryland**
**20-June-2005**
**(last update: 17-June-05, 18:00)**

**For more detail, see:**
**HCAL fall meeting, Fermilab 11&13-Nov.2004  **_[http://agenda.cern.ch/fullAgenda.php?ida=a044477](http://agenda.cern.ch/fullAgenda.php?ida=a044477)_
**9****th**** RDMS Annual CMS Conference, Minsk, 01-Dec-2004  **_[http://agenda.cern.ch/fullAgenda.php?ida=a044186](http://agenda.cern.ch/fullAgenda.php?ida=a044186)_

## General

## HCAL on a Table at H2

**Pivot of table**
**       ****= IP at LHC**
**A phi slice of **
**CMS HCAL**
**HB: 2 wedges**
**        ****8 ******** segments**
**        ******** ****= 40deg.**
**HE:  4 ******** segments**
**        ******** ****= 20 deg.**
**HO: 6 ******** sector**
**        ******** ****= 30 deg.**
**ECAL: 7x7 crystals**

**ECAL**

** ****pivot**** **

**beam**

## Test Beam: 2002-2003

**Goal:**
** ****Test the integrated system**
**   ****with production modules**
** ****Verify ******** source calibration**
** ****Source/ADC vs. GeV/ADC**
** ****Operate all calibration tools**
**  ****and look for improvement.**
** ****measure basic parameters**
**  ****for MC, e.g.**
**- pulse shape**
**- signal timing**
**- attenuation**
**- noise**
**- gaps between modules**
** ****resolution and linearity**
**All achieved!  ****(repeated in 2004)**

**2 HB production wedges, 1 HE prototype wedge**
**HO layers on a movable table at CERN H2 beam line.**

**2002/03: pi- 20-300GeV, e- 20-100GeV, mu- 225GeV**

## TB2004 Goals

**Take data set for G4 validation**
**Low energy / high energy data with particle ID**
**Longitudinal/Transverse shower profile**
**+  new QIE  (time slew, noise, pulse shape)**
**Test full system ******** Slice Test with EMU **
**Clock distribution / Synchronization**
**HTR: L1 primitive generation (ADC********E****T****), Threshold **
**Operate all calibration tools and develop algorithm & code**
**Wire source/ Laser / LED    (including synchronization)**
-  detect 1% gain change
**+  demonstrate that HF calibration scheme work.  **
**Deploy prototype HCAL database   **
**Conditions DB / Construction DB / Configuration DB etc.**
**Operate remote monitoring and analysis (+debugging)**
**E-logbook, video/audio system + other tools (desk top exchange) etc. **
**DCS monitoring data.                            **
**Feed data into ORCA**
**Interface to raw data and conditions DB.      **

## Summary of TB2004 runs

**    ****May 17-Jun8 , High energy run**** **
**   ****Jun 8-14,  no beam time install HO**
**   ****Jun 14-21 25 ns run time  **
**   ****Jun 25- Jul 7: first  HF run  **
**   ****Jul 7-14  : 1****st**** week of VLE run **
**   ****Jul 14- Aug11 : second  HF run **
**   ****Aug 11- 18:  2****nd**** week of VLE run **
**   ****Aug 18- Sep 22: other experiments  in H2  **
**   ****Sep 22-Oct 4: HCAL-EMU  setup **
**   ****Oct 4-11 25ns ****HCAL-EMU**** run **
**   ****Oct 13-18: 3****rd**** week of VLE run**
**   ****Feb.’05  Wire source calibration run **

**Slice Test**

## QIE Pulse Shape: New vs Old

**88% signal collection**
**In 2 time slices.**

## QIE Time Slew

**New faster QIE for HB/HE**

## QIE Noise

## Wire Source Calibration

**Wire Source Calibration**

**HF**
**5mCi Co****60**

**HE**

**<3%**

**~5%**

## HPD Stability

**Arjan Heering, HCAL meeting @ FNAL, 12-Nov-04**

## HB

## HE (2003)

## HE Setup (TB2003)

## HE Response (2003)

**V.Ladygin, A.Volkov, RDMS Meeting @ Dubna, 05-Dec-2003**

** ****at T19-2**
**(3x3 sum)**

**8%**

**300GeV pi-**

**100GeV e-**

**104GeV**
**7.6%**

**(calibrated with 300GeV pion)**

**GeV**

**Linearilized ADC counts**

**0.163GeV/ADC**

##  Muons in HE

## TB03 HE

**V.Ladygin, A.Volkov, RDMS Meeting @ Dubna, 05-Dec-2003**

## TB03 HE

## TB03 HE

## L0 Optimization

## HE (2004)

## HE

**HE**

**Source calibration**

**(Collimated source) / (Wire source)**

**Vladimir Talvov, **
**RDMS meeting, Minsk **
**01-Dec-2004**

**21**

**28**

## HE (TB2004)

**Run  9641**
**________________________________________________________**
**|**** L ****|****phi****|****Sp****|****   ch    ****|****maximum ****|****maxreel****|****meanped ****|****meansig ****|****sig-ped****|**
**|**** 1 ****|****  3  ****|****11****|**** 1(17f) ****|****   18.380   ****|****      568  ****|****    18.151  ****|****   18.371  ****|****  0.221 **** |**
**|**** 1 ****|****  3 **** |****11****|**** 3(25f) ****|****   22.119   ****|****      682  ****|****    21.978  ****|****   22.115  ****|****  0.137  ****|**
**| ****1 ****|****  3  ****|****11****|**** 4(18f) ****|****   19.660   ****|****      579  ****|****    19.417  ****|****   19.656  ****|****  0.240  ****|**
**| ****1 ****|****  3  ****|****11****|**** 7(27f) ****| ****  19.581   ****|****      707  ****|****    19.412  ****|****   19.578  ****|****  0.166  ****|**
**|**** 1 ****|****  3  ****|****11****|****12(23f)****|****  18.736   **** | ****     654  ****|****    18.525  ****|****   18.734  ****|****  0.210  ****|**
**|**** 1 ****|****  3  ****|****11****|****14(20f)****|****  20.256    ****|****      614  ****|****    20.060  ****|****   20.255  ****|****  0.195  ****|**
**|**** 1 ****|****  3  ****|****11****|****16(19f)****|****  18.314    ****| ****     601  ****|****    18.146  ****|****   18.311  ****| **** 0.165  ****|**
**|**** 1 ****|****  3  ****|****11****|****18(21f)****|****  19.406    ****|****      629  ****|****    19.237  ****|****   19.398  ****|****  0.160  ****|**

**Vladimir Talvov, **
**RDMS meeting, Minsk **
**01-Dec-2004**

## HE (TB2004)

**Vladimir Talvov,  RDMS meeting, Minsk. 01-Dec-2004**

## Gain from Muons (TB2004)

**Selda Esen, HCAL meeting @FNAL, Nov-04**

## HE Geometry (TB2004)

**Selda Esen, HCAL meeting @FNAL, Nov-04**

## Pions in Muon beam

**Selda Esen, HCAL meeting @FNAL, Nov-04**

## Muons and Pions in electron beam

**Selda Esen, HCAL meeting @FNAL, Nov-04**

## TB2004; HE resolution

**Resolution**

**Selda Esen, HCAL meeting @FNAL, Nov-04**

## HO

## HO – Tail Catcher

**E=HB+3.6*HO**

**Seema Sharma, May’05**

**9.0%**

**11%**

## HB1+HO

**HB1+HO**
**300GeV**

## HB1+HO

**HB1+HO**
**150GeV**

## HB1+HO

**HB1+HO**
**100GeV**

## HB1+HO

**HB1+HO**
**50GeV**

## HB1+HO

**HB1+HO**
**30GeV**

## HF

## HF Response

**HF**

**ITEP Group (A.Ulyanov)  CMS Week, 17-March-05**

## HF: Long + Short

**HF: Long + Short**

******/e:   large improvement**
**  ****(0.6-0.7)****L**** ******** (0.86-1.0)****L+S**** **
**rms/E:  Minor improvement**
**  ****24% ******** 22%  for 100GeV electron**
**  ****33% ******** 30%  for 100GeV pion**

## HF: Source/Beam(100GeV e)

**1350V**

**A.Krokhotine, 17-March-05, CMS Week**

## HF Uniformity

**2x2cm****2**** scan**

## 100 GeV e- signal in HF short fiber.

**100 GeV e- signal in HF short fiber.**

## G4 Validation

## GEANT4 Validation Data

**Quantities:**
**   ********/e**
**   ****Shape of dN/dE dist.**
**       ****(resolution)**
**   ****Shower profile**
**       ****longitudinal**
**       ****transverse**
**Detectors:**
**   ****HF: **
**      ****Cerenkov based**
**   ****EC+HB**
**      ****Scintillation based**

**HB2 readout was reconfigured for longitudinal shower profile study.**

## VLE Beam Line at H2

**P-ID:    **
**     ****CK2-  electron**
**     ****CK3-  pion / kaon / proton**
**     ****V3, V6, VM – muon**

**VLE tag against punchthrough muon**

**WC**
**  ****single hit to reject**
**  ****interaction in beam line**

## Data Sets

**Electron: 5-150GeV**
**Pion:        5-300 GeV**
**Proton:    5-9 GeV**

**(2GeV and 3GeV data exist, but a lof of junks.)**

## GEANT4 Simulation

## HF

**Wedge 2.13  Long Fiber Section (EM)**

**100 GeV electrons and pions**

**Pion/elec = 0.72**

## HF (Long)

**HF (Long)**
**    ****Electrons**
**    ****Very detailed simulation with G4**
**    ****reproduces electron signal very well.**

**50GeV e-**

**100GeV e-**

**Resolution: e-**

## HF (Long)

**HF (Long)**
**       ****pi-**

**resolution**

******/e**

**50GeV**

**100GeV**

**LHEP shows good agreement.**

## EC+HB  High Energy Data (pi-)

**300 GeV**

**150 GeV**

**100 GeV**

**30 GeV**

**HB**

**ECAL**

**ECAL**

**HB**

****

**e****-**

## 9 GeV pi+ beam

**e****+**

**mip in ECAL, i.e. no-interaction in ECAL**

**Data**

**MC**
**W****/o noise**

****** ******** ******** ****decays in beam line**

**Interactions in beam line**** **

**Need a lot of clean-up !**

**Interactions in ECAL**

## Electron identification

**Electron identification**

## P-ID with Cerenkov Counter 3 (CK3)

**P-ID with Cerenkov Counter 3 (CK3)**

- The momentum thresholds
- for the range of dn are:
-   dn       P()  P()   P(p)   P(K)
- *E-6
- 2432     2.0    1.51    13.5
- 1557     2.5    1.89    16.8
- 1082     3.0    2.27    20.2
-   ***795     3.5    2.65    23.5   12.35***
-   609     4.0    3.03    26.9
-   481     4.5    3.41    30.3
-   390     5.0    3.79    33.6

**Freon Cerenkov Counter**

## Pi -

**Pi ****-**

**Pi ****+**

**(e****-**** removed)**

**(e****+**** removed)**

**4.4%**

**17.9%**

## HF Shower Library and TB2004

**Linearity**

**Resolution**

**Data vs G4:       ~10%                                                                      ~10%**

**New HF Shower Library: **
**detailed simulation from cerenkov light generation to photo electron simulation at PMT with G4.**
****** ****Test in progress.**

## 300GeV Pi-

**290**

**286**

**296**

**295**

## 150GeV Pi-

**141.4**

**141.0**

**147.7**

**149.1**

## 100GeV pi-

**94.1**

**93.7**

**99.5**

**100.2**

## 50GeV pi-

**46.2**

**45.7**

**50.0**

**50.0**

## 30GeV pi-

**26.8**

**26.3**

**29.2**

**29.2**

## 15GeV pi-

**12.5**

**12.5**

**14.4**

**14.3**

## 10GeV pi-

**7.65**

**7.82**

**9.06**

**9.10**

## 9GeV pi-

**6.72**

**6.64**

**7.97**

**7.71**

## 7GeV pi-

**5.19**

**5.17**

**6.09**

**5.92**

## 5GeV pi-

**3.63**

**3.58**

**4.29**

**4.02**

## 9GeV proton

## 7GeV proton

## 5GeV proton

## ODU (TB2004/HB2): layer-pixel maplongitudinal readout

**Pixel Numbers**

**Layer numbers**

## ODU (HB): tower-pixel map

**HB19.1**

**HB19.4**

**HB19.3**

**HB19.2**

**Tower Num.**

**Pixel Number**

## Comparison of Shapes    (Area and mean normalized.   MC: no digitization, no noise)

## Comparison of Shapes     (Area and mean normalized.   MC: no dizitization, no noise)

**Comparison of Shapes****     ****(Area and mean normalized.   MC: no dizitization, no noise)**

## EC+HC

## HB (mip in EC)

## Big PbWO4

## Big PbWO4

## Longitudinal Shower Profile

**MIP in ECAL.**
**MC- no aluminum block behind ECAL.**

## Longitudinal Shower Profile

**MIP in ECAL.**
**MC- 8cm aluminum block behind ECAL.**

## Longitudinal Shower Profile

**MIP in ECAL.**
**MC- 8cm aluminum block behind ECAL.**

## Additional Slides

## Conclusion

**We have achieved all TB2004 goal.**
**Data set (3-300GeV) for G4 validation.**
**L1 trigger primitive generation.**
**Slice test with ME.**
**Operation and verification of all calibration system.**
**Remote monitoring and data analysis.**
**Complete data flow from detector through ORCA, ****(i****ncluding access to ConditionsDB.)**
**G4 validation is in progress.**
**A lot of data to analyze.**

## HB Configuration

## GEANT4 Validation

**Test of physics models for hadron showers**
**Response**
- Resolution - shape of dN/dE (no gaussian shape)
- Linearity - pion/electron and proton/electron
**Shower profile**
- Longitudinal
- Transverse
**Calorimeter setup**
**HF**
- Cerenkov  sensitive to pi0 in shower.
**EC+HB+HO**
- Sincitillation  scintillation saturation effect (Birks’ law)?

## Data:/e (very preliminary)

**Data  HB1**

**HC only**

**EC+HC**

**EC+HC**

## First look ofLongitudinal Shower Profile

**HB2 only**

**Need**
**Clean-up**

## TB2006 Plan

**Goal:**
**Perform full system test with “production” HW and SW.**
- Generation of TPG and Raw Data.
- Operation of calibration system.
**Measure performance of EC+HC combined system.**
- Improve pi/e measurement for 2-300GeV.
- Measure material distribution between ECAL and HCAL.
- Look for any surprise, e.g. nuclear counter effects in APD, etc.
**Needs:**
**ECAL super module.**
**Improved testbeam triggering/tagging system.**
- tags – muons, beam interactions, P-ID.

## pi- :  High Energy Data

**300 GeV**

**150 GeV**

**100 GeV**

**30 GeV**

**HCAL**

**ECAL**

**ECAL**

**HCAL**

****

**e****-**

## Low Energy  Data

**Test beam Needed:**
**- Particle ID with Cerenkov   counter.**
**- Muon halo veto.**
**- Different G4 Simulations in advance - determine which parameters need to be tested.**
**Data sets in TB2004.**
**   ****Low:   2, 3, 5, 7, 9 GeV**
**   ****Intermediate:  10, 15, 20, 30**
**   ****High:  30, 50, 100, 150, 300**

**HCAL alone**

**G4 Simulation **

**It is very important to understand the calorimeter response to low energy particles for simulation and reconstruction of jets and MET.**
**Low energy data (below 30GeV) were missing in TB2002-2003.**

## G4 Validation

**Two models in G4, QGSP and LHEP show different profile.**

**5 GeV**

**10 GeV**

**20 GeV**

**100 GeV**

## G4 Prediction

**Resolution**

**Linearity for ECAL+HCAL**

## G4 vs Data: First attempt to overlay

**Data vs G4  (HB2)**

**Need more clean-up!**

## HE/ME Slice Test at H2

**With 25ns beam in Oct.2004**

## Punchthrough: EE-HE-ME

**Readout from HE and EE were successfully synchronized!**

## Trigger Primitive Generation

**HTR:  HCAL Trigger Readout**

## Example of Trigger Primitives

## Database

**HCAL Prototype for TB2004**
**      ****1 server at FNAL**
**      ****1 server at H2**
**      ****loaded- **
**          ****gain, pedestal**
**          ****HV monitoring**
**Magnet Test 2005**
**  ****- Schema for all detectors.**
**  ****- DBs filled with initial data.**
**  ****- Complete data flow chain.**
**    ****(All databases).**
**  ******** ****load initial constants for **
**  ****all HCAL channels.**
**       **
**HCAL has built a DB team at FNAL and CERN with help from FNAL/CD. Working closely with PIXEL and EMU, and Dubna DB group.**

## ORCA: HCAL Data Flow

**Raw**

**SimHit**
**   ****- G4/OSCAR**

**DAQ**

**Simulated Raw**
**   ****- pile-up with bx-id**
**   ****- pulse shape generation**
**         ****- noise, photo statistics, QIE dizitization**
**         ****- time-slew, pedestal**
**   ****- data formatting**

**RecHit**
**   ****- data validation**
**        ****- synchonization error etc.**
**        ****- error report**
**   ****- hot cell / dead cell killer (static)**
**   ****- hot cell / noisy channel killer (dynamic)**
**   ****- pulse shape analysis**
**        ****- time slew correction**
**        ****- pile-up signal subtraction**
**        ****- pedestal and gain correction**

**ConditionsDB**
**   ****- list of hot cells**
**   ****- list of dead cells**
**   ****- pedestal**
**   ****- gain**

**Zero Bias**
**    ****- Sim / Data**

**EcalPlusHcalTower**

**C.Tully, J.Mans**
**for TB2004**

**R.Wilkinson**
**M.Case**
**G.Lukhanin**
**S.Kosyakov**

**(R.Wilkinson/S.Abdullin)**

**J.Damgov,S.Piperov**

**Calibration Data**

**FroNtier**

## TB2004 data to ORCA 

**Example: reconstructing  pion events on the bare HCAL as jets using the Iterative Cone algorithm (R=0.5)**
**Full chain of creating dataframes, RecHits, and performing jet reconstruction has been demonstrated.**
**No calibration/gain correction applied (results in fC)**
**Next: **
**Access to CondtionsDB through FroNtier**

## FroNtier

**ORCA**

**CondDB**