Energy Correction with Ecal(EE) beamtest info.

2011 08 30

Basic environments

Software CMSSW4_2_5 on common tool
Framwork fwlite
Working dir cms-kr-knu-ch:~shchang/WORK/ECAL/V296/CMSSW4_2_5
Require things LazyTools, Basic Cluster, RecHit, EventSetup

Todo list (plan)

  • Figuring out and note here the structure of lazytool example codes.
     
    EcalRecHitsSorted_reducedEcalRecHitsEB__RECO
    
    EcalRecHitsSorted_reducedEcalRecHitsEE__RECO
    
    EcalRecHitsSorted_ecalRecHit_EcalRecHitsEB_RECO.
    
    EcalRecHitsSorted_ecalRecHit_EcalRecHitsEE_RECO.
    
    EcalRecHitsSorted_ecalPreshowerRecHit_EcalRecHitsES_RECO
    
    recoSuperClusters_hybridSuperClusters__RECO
    
    recoSuperClusters_hybridSuperClusters_uncleanOnlyHybridSuperClusters_RECO
    
    recoCaloClusters_multi5x5BasicClusters_multi5x5BarrelBasicClusters_RECO
    
    recoCaloClusters_multi5x5BasicClusters_multi5x5EndcapBasicClusters_RECO
    
         //showershape
         electrone1x5[nElectron]                          = jj->e1x5();
         electrone2x5Max[nElectron]                       = jj->e2x5Max();
         electrone5x5[nElectron]                          = jj->e5x5();
         electronsigmaIetaIeta[nElectron]                 = jj->sigmaIetaIeta();
         electrone3x3[nElectron] = lazyTool.e3x3(*(scRef->seed())); 
         electroneMax[nElectron] = lazyTool.eMax( *(scRef->seed())); 
         
         
          ///more variables from lazyTools
         electroneLeft[nElectron] = lazyTool.eLeft( *(scRef->seed()));
         electroneRight[nElectron] = lazyTool.eRight( *(scRef->seed()));
         electroneBottom[nElectron] = lazyTool.eBottom( *(scRef->seed()));
         electroneTop[nElectron] = lazyTool.eTop( *(scRef->seed()));
         
          electrone1x3[nElectron] = lazyTool.e1x3( *(scRef->seed()) );
          electrone3x1[nElectron] = lazyTool.e3x1( *(scRef->seed()) );
          electrone2x2[nElectron] = lazyTool.e2x2( *(scRef->seed()) );
          electrone3x2[nElectron] = lazyTool.e3x2( *(scRef->seed()) );
          electrone4x4[nElectron] = lazyTool.e4x4( *(scRef->seed()) );
          electrone2x5Right[nElectron] = lazyTool.e2x5Right( *(scRef->seed()) );
          electrone2x5Left[nElectron] = lazyTool.e2x5Left( *(scRef->seed()) );
          electrone2x5Top[nElectron] = lazyTool.e2x5Top( *(scRef->seed()) );
          electrone2x5Bottom[nElectron] = lazyTool.e2x5Bottom( *(scRef->seed()) );
          ////electrone2x5Max[nElectron] = lazyTool.e2x5Max( *(scRef->seed()) );
    
    
         ///idmax's   
         
         DetId idmax = lazyTool.getMaximum(*(scRef->seed())).first; 
         const EcalRecHitCollection & rechits  = idmax.subdetId() == EcalBarrel ? *Brechit : *Erechit; 
         EcalRecHitCollection::const_iterator it = rechits.find( idmax );
         
         ///for spike's
         electronswissCross[nElectron] = EcalTools::swissCross(idmax,rechits,1,false); //compute only if Et(seed) > 1GeV, ieta==85 as well
    
         if(idmax.subdetId() == EcalBarrel){
           EBDetId det = EBDetId(idmax);
           electronieta[nElectron] = det.ieta();
           electroniphi[nElectron] = det.iphi();
         }else{
           EEDetId det = EEDetId(idmax);
           electronieta[nElectron] = det.ix();
           electroniphi[nElectron] = det.iy();
         }
         electronE2overE9[nElectron] = mye2overe9(idmax,rechits);
         
     

  • Insert EventSetup
  • Insert RecHit and basic cluster header files.
  • Insert Lazytool header file(s).
  • Compiling without errors.

  • Implement lazytool class.(where?)
  • Compiling without errors.

  • fill tree(ntuple) with some lazytool's variables.
  • Compiling without errors.

  • Run executable and check output ntuple.

-- SunghyunChang - 30-Aug-2011

Edit | Attach | Watch | Print version | History: r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r2 - 2011-09-30 - SunghyunChang
 
    • 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