---

EcalDeadChannelRecoveryAlgos Class Reference

Date

2007/05/09 14:58:40

Revision

1.3

More...

#include <RecoLocalCalo/EcalDeadChannelRecoveryAlgos/src/EcalDeadChannelRecoveryAlgos.cc>

List of all members.

Public Member Functions
EcalRecHit	Correct (const EBDetId Id, const EcalRecHitCollection *hit_collection, string algo_, double Sum8Cut)
	EcalDeadChannelRecoveryAlgos ()
	EcalDeadChannelRecoveryAlgos (const CaloTopology theCaloTopology)
	~EcalDeadChannelRecoveryAlgos ()
Private Member Functions
double	MakeNxNMatrice (EBDetId itID, const EcalRecHitCollection *hit_collection, double *MNxN)
Private Attributes
CaloTopology	calotopo

---

Detailed Description

Date

2007/05/09 14:58:40

Revision

1.3

Description: <one line="" class="" summary>="">.

Implementation: <Notes on="" implementation>="">

Definition at line 37 of file EcalDeadChannelRecoveryAlgos.h.

---

Constructor & Destructor Documentation

EcalDeadChannelRecoveryAlgos::EcalDeadChannelRecoveryAlgos

(

const CaloTopology

theCaloTopology

)

Definition at line 51 of file EcalDeadChannelRecoveryAlgos.cc.

References calotopo.

{
  //now do what ever initialization is needed
  calotopo = theCaloTopology;
}

EcalDeadChannelRecoveryAlgos::EcalDeadChannelRecoveryAlgos

(

)

EcalDeadChannelRecoveryAlgos::~EcalDeadChannelRecoveryAlgos

(

)

Definition at line 59 of file EcalDeadChannelRecoveryAlgos.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

---

Member Function Documentation

EcalRecHit EcalDeadChannelRecoveryAlgos::Correct	(	const EBDetId	Id,
		const EcalRecHitCollection *	hit_collection,
		string	algo_,
		double	Sum8Cut
	)

Definition at line 73 of file EcalDeadChannelRecoveryAlgos.cc.

References CorrectDeadChannelsClassic(), CorrectDeadChannelsNN(), EBDetId::ieta(), and MakeNxNMatrice().

Referenced by EcalDeadChannelRecoveryProducers::produce().

{
  double NewEnergy=0.0;
  
  double MNxN[121];

  if(algo_=="Spline"){
     
    if(MakeNxNMatrice(Id,hit_collection,MNxN)>Sum8Cut){
      NewEnergy = CorrectDeadChannelsClassic(MNxN,Id.ieta());
    }
  }else if(algo_=="NeuralNetworks"){
   
    if(MakeNxNMatrice(Id,hit_collection,MNxN)>Sum8Cut){
      NewEnergy = CorrectDeadChannelsNN(MNxN);
    }
  }
  

  EcalRecHit NewHit(Id,NewEnergy,0);
  return NewHit;

}

double EcalDeadChannelRecoveryAlgos::MakeNxNMatrice	(	EBDetId	itID,
		const EcalRecHitCollection *	hit_collection,
		double *	MNxN
	)			[private]

Definition at line 105 of file EcalDeadChannelRecoveryAlgos.cc.

References funct::abs(), calotopo, GenMuonPlsPt100GeV_cfg::cout, DetId::Ecal, EcalBarrel, edm::SortedCollection< T, SORT >::end(), lat::endl(), edm::SortedCollection< T, SORT >::find(), CaloTopology::getSubdetectorTopology(), CaloSubdetectorTopology::getWindow(), i, EBDetId::ieta(), int, EBDetId::iphi(), DetId::null(), and size.

Referenced by Correct().

                                                                                                                        {



  //Build NxN around a given cristal
  for(int i=0; i<121;i++)MNxN[i]=0.0;
  
  //cout<<"===================================================================="<<endl;
  //cout<<" Dead Cell CENTRAL  eta = "<< itID.ieta()<<" ,  phi = "<< itID.iphi()<<endl;

  const CaloSubdetectorTopology* topology=calotopo.getSubdetectorTopology(DetId::Ecal,EcalBarrel);
  int size =5;
  std::vector<DetId> NxNaroundDC = topology->getWindow(itID,size,size);
  
  
  //cout<<"NxNaroundDC size is = "<<NxNaroundDC.size()<<endl;
  vector<DetId>::const_iterator theCells;
  
  EBDetId EBCellMax = itID;
  double EnergyMax = 0.0;

  for(theCells=NxNaroundDC.begin();theCells<NxNaroundDC.end();theCells++){
    EBDetId EBCell = EBDetId(*theCells);
      
    // We Will look for the cristal with maximum energy 
    if(!EBCell.null()){
      EcalRecHitCollection::const_iterator goS_it = hit_collection->find(EBCell);
      if( goS_it !=  hit_collection->end() && goS_it->energy()>=EnergyMax){EnergyMax=goS_it->energy(); EBCellMax = EBCell;}
    }else{
      continue; 
    }
  } 
  if(EBCellMax.null()){cout<<" Error No maximum found around dead channel, no corrections applied"<<endl;return 0;}
  //cout << " Max Cont Crystal eta phi E = " << EBCellMax.ieta() <<" "<< EBCellMax.iphi() <<" "<< EnergyMax<<endl;


  
  //NxNaroundMaxCont with N==11
  // 000 is now the maximum containement one:
  //Any modification of the following parameters will require changes in the Correction algos
  // The window is large enought to avoid modification of the number.
  int FixedSize =11;
  std::vector<DetId> NxNaroundMaxCont = topology->getWindow(EBCellMax,FixedSize,FixedSize);

  double ESUMis=0.0;
  int theIndex=0;

  vector<DetId>::const_iterator itCells;
  for(itCells=NxNaroundMaxCont.begin();itCells<NxNaroundMaxCont.end();itCells++){
    EBDetId EBitCell = EBDetId(*itCells);
    int CReta = EBitCell.ieta();
    int CRphi = EBitCell.iphi();

    double Energy = 0.0;
    if(!EBitCell.null()){
      EcalRecHitCollection::const_iterator goS_it = hit_collection->find(EBitCell);
      if( goS_it !=  hit_collection->end() ) Energy=goS_it->energy();
    }

    //cout<<"Around DC we have eta,phi,E "<<CReta<<" "<<CRphi<<" "<<Energy<<endl;

    //============
    int ietaCorr = 0;
    if((CReta * EBCellMax.ieta()) < 0 && EBCellMax.ieta()>0 )ietaCorr= -1;
    if((CReta * EBCellMax.ieta()) < 0 && EBCellMax.ieta()<0 )ietaCorr=  1;
    if((CReta * EBCellMax.ieta()) > 0 )ietaCorr= 0;
    int ieta = -( (CReta -ietaCorr) - EBCellMax.ieta() ) + int((FixedSize - 1)/2);

    int iphiCorr = 0;
    if((CRphi - EBCellMax.iphi())> 50)iphiCorr= 360;
    if((CRphi - EBCellMax.iphi())<-50)iphiCorr=-360;
    int iphi = CRphi - EBCellMax.iphi() - iphiCorr + int((FixedSize - 1)/2);

    int MIndex = ieta+iphi*FixedSize;
    if(abs(CReta)<=85)
    MNxN[MIndex]=Energy;
    
    
    //============



    
    //We add up the energy in 5x5 around the MaxCont to decide if we will correct for DCs
    if(theIndex>=36 && theIndex<=40)ESUMis +=  Energy;
    if(theIndex>=47 && theIndex<=51)ESUMis +=  Energy;
    if(theIndex>=58 && theIndex<=62)ESUMis +=  Energy;
    if(theIndex>=69 && theIndex<=73)ESUMis +=  Energy;
    if(theIndex>=80 && theIndex<=84)ESUMis +=  Energy;
    theIndex++;
  }
  //cout<<"Around MaxCont Collected Energy in 5x5 is =  "<< ESUMis  <<endl;
  //So now, we have a vector which is ordered around the Maximum containement and which contains a dead channel as:
    //Filling of the vector : NxNaroundDC with N==11 Typo are possible...
    // 000 is Maximum containement which is in +/- 5 from DC
    //
    // 120 119 118 117 116 115 114 113 112 111 110  
    // 109 108 107 106 105 104 103 102 101 100 099 
    // 098 097 096 095 094 093 092 091 090 089 088
    // 087 086 085 084 083 082 081 080 079 078 077
    // 076 075 074 073 072 071 070 069 068 067 066
    // 065 064 063 062 061 060 059 058 057 056 055
    // 054 053 052 051 050 049 048 047 046 045 044
    // 043 042 041 040 039 038 037 036 035 034 033
    // 032 031 030 029 028 027 026 025 024 023 022
    // 021 020 019 018 017 016 015 014 013 012 011
    // 010 009 008 007 006 005 004 003 002 001 000
  

  
  //================

  return ESUMis;
}

---

Member Data Documentation

CaloTopology EcalDeadChannelRecoveryAlgos::calotopo [private]

Definition at line 48 of file EcalDeadChannelRecoveryAlgos.h.

Referenced by EcalDeadChannelRecoveryAlgos(), and MakeNxNMatrice().

---

The documentation for this class was generated from the following files:

• RecoLocalCalo/EcalDeadChannelRecoveryAlgos/interface/EcalDeadChannelRecoveryAlgos.h
• RecoLocalCalo/EcalDeadChannelRecoveryAlgos/src/EcalDeadChannelRecoveryAlgos.cc

---