Energy Correction with Ecal(EE) beamtest info.

2011 08 30

Basic environments

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