2006 H4 ECAL testbeam

Simulation

See here

Analysis

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
* run on the RecHits
   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

Edit | Attach | Watch | Print version | History: r8 < r7 < r6 < r5 < r4 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r8 - 2006-09-13 - AndreHolzner
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback