2006 H4 ECAL testbeam

Simulation

Analysis

• Link to testbeam analysis mailing list archive (MMM login required)
• H4 ECAL electronic log (long page, requires CMS AFS login)
• ECALTB4 Release and hints
• Testbeam Meeting shortly after the run (2006-08-31)
• Pietro Govoni's mail concerning configuration files to reconstruct rechits
• Muon runs (from elog, click here for the individual entry, from Monday 2006-08-28, SM17):

  run	crystal	(ieta,iphi)	events
  13804	#48	(3,13)	30k
  13805	#148	(8,13)	30k
  13806	#248	(13,13)	30k
  13807	#348	(18,13)	30k
  13809	#448	(23,13)	30k
  13810	#548	(28,13)	30k
  13811	#648	(33,13)	30k
  13812	#748	(38,13)	30k
  13813	#848	(43,13)	30k
  13814	#948	(48,13)	30k
  13815	#1048	(53,13)	30k
  13816	#1148	(58,13)	30k
  13817	#1248	(63,13)	30k
  13818	#1348	(69,13)	30k
  13819	#1448	(73,13)	30k
  13820	#1548	(78,13)	30k
  13821	#1648	(83,13)	100k

The relation between crystal number and ieta,iphi seems to be the following: crystal_number = ieta * 20 + 1 - iphi

Example of how to get the reconstructed energies

Using ECALTBH4_0_2_0_pre2 (2006-08-31) which is based on CMSSW_0_8_1:

• install and build the software:

  scramv1 project -n ECALTBH4_0_2_0_pre2 CMSSW CMSSW_0_8_1
  cd ECALTBH4_0_2_0_pre2
  PackageManagement.pl --rel ECALTBH4_0_2_0_pre2
  cd src
  scramv1 b

• setup the environment variables

  eval `scramv1 runtime -sh`
  export CORAL_AUTH_PATH=/afs/cern.ch/cms/ECAL/testbeam/pedestal/2006/config_files
  export TNS_ADMIN=/afs/cern.ch/project/oracle/admin

• adapt the configuration file:

  sed \
     -e s,INPUTFOLDER,/castor/cern.ch/cms/archive/ecal/h4tb.pool, \
     -e s,INPUTFILE,h4b.00013804.A.0.0, \
     -e s,OUTPUTFILE,h4b.00013804.A.0.0, \
     -e s,OUTPUTFOLDER,.,  \
     -e s,readConfiguration2006_v1,readConfiguration2006_v1_fromDB, \
     < /afs/cern.ch/user/e/ecaltb06/public/configFiles/officialUncalib2006Production_0_3_0.cfg > test.cfg

   cmsRun -p test.cfg

• The output file has the following branches in the Events tree:

EventAux
EBDataFramesSorted_ecalTBunpack__uncalibRecHitsProd
EBDetIdedmEDCollection_ecalTBunpack_EcalIntegrityChIdErrors_uncalibRecHitsProd
EBDetIdedmEDCollection_ecalTBunpack_EcalIntegrityDCCSizeErrors_uncalibRecHitsProd
EBDetIdedmEDCollection_ecalTBunpack_EcalIntegrityGainErrors_uncalibRecHitsProd
EBDetIdedmEDCollection_ecalTBunpack_EcalIntegrityGainSwitchErrors_uncalibRecHitsProd
EBDetIdedmEDCollection_ecalTBunpack_EcalIntegrityGainSwitchStayErrors_uncalibRecHitsProd
EcalDCCHeaderBlocksSorted_ecalTBunpack__uncalibRecHitsProd
EcalElectronicsIdedmEDCollection_ecalTBunpack_EcalIntegrityMemBlockSize_uncalibRecHitsProd
EcalElectronicsIdedmEDCollection_ecalTBunpack_EcalIntegrityMemChIdErrors_uncalibRecHitsProd
EcalElectronicsIdedmEDCollection_ecalTBunpack_EcalIntegrityMemGainErrors_uncalibRecHitsProd
EcalElectronicsIdedmEDCollection_ecalTBunpack_EcalIntegrityMemTtIdErrors_uncalibRecHitsProd
EcalMatacqDigi_ecalTBunpack__uncalibRecHitsProd
EcalPnDiodeDigisSorted_ecalTBunpack__uncalibRecHitsProd
EcalTBEventHeader_ecalTBunpack__uncalibRecHitsProd
EcalTBHodoscopeRawInfo_ecalTBunpack__uncalibRecHitsProd
EcalTBHodoscopeRecInfo_ecal2006TBHodoscopeReconstructor_EcalTBHodoscopeRecInfo_uncalibRecHitsProd
EcalTBTDCRawInfo_ecalTBunpack__uncalibRecHitsProd
EcalTBTDCRecInfo_ecal2006TBTDCReconstructor_EcalTBTDCRecInfo_uncalibRecHitsProd
EcalTrigTowerDetIdedmEDCollection_ecalTBunpack_EcalIntegrityBlockSizeErrors_uncalibRecHitsProd
EcalTrigTowerDetIdedmEDCollection_ecalTBunpack_EcalIntegrityTTIdErrors_uncalibRecHitsProd
EcalUncalibratedRecHitsSorted_ecal2006TBWeightUncalibRecHit_EcalUncalibRecHitsEB_uncalibRecHitsProd

Example ROOT plots

* Useful aliases:

  Events->SetAlias("hodo","EcalTBHodoscopeRecInfo_ecal2006TBHodoscopeReconstructor_EcalTBHodoscopeRecInfo_uncalibRecHitsProd.obj")
  Events->SetAlias("header","EcalTBEventHeader_ecalTBunpack__uncalibRecHitsProd.obj")
  Events->SetAlias("urechit","EcalUncalibratedRecHitsSorted_ecal2006TBWeightUncalibRecHit_EcalUncalibRecHitsEB_uncalibRecHitsProd.obj.obj")

  Events->SetAlias("urechit_iphi","urechit.id_.rawId() & 0x1ff");
  Events->SetAlias("urechit_ieta","(urechit.id_.rawId()  >> 9) & 0x7f");

  Events->SetAlias("crystal_inbeam_iphi","header.crystalInBeam_.id_ & 0x1ff");
  Events->SetAlias("crystal_inbeam_ieta","(header.crystalInBeam_.id_ >> 9) & 0x7f");

(see http://cmslxr.fnal.gov/lxr/source/DataFormats/EcalDetId/interface/EBDetId.h for decoding the raw detector id)

* Hodoscope:

  Events->Draw("hodo.pos_y_:hodo.pos_x_","hodo.pos_y_ > -40")

• see the crystal shape (run 13804)

  Events->Draw("urechit.amplitude_:hodo.pos_x_"," urechit_iphi == 13 && urechit_ieta == 3 && urechit.amplitude_ < 100 && hodo.pos_y_ > -50");
  Events->Draw("urechit.amplitude_:hodo.pos_y_"," urechit_iphi == 13 && urechit_ieta == 3 && urechit.amplitude_ < 100 && hodo.pos_y_ > -50");

-- AndreHolzner - 01 Sep 2006