CMS 3D CMS Logo
muonisolation::CaloExtractor

muonisolation::CaloExtractor Class Reference

#include <RecoMuon/MuonIsolation/plugins/CaloExtractor.h>

Inheritance diagram for muonisolation::CaloExtractor:

reco::isodeposit::IsoDepositExtractor

List of all members.

Public Member Functions

 CaloExtractor (const edm::ParameterSet &par)
 CaloExtractor ()
virtual reco::IsoDeposit deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const
 make single IsoDeposit based on track as input purely virtual: have to implement in concrete implementations
virtual void fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &tracks)
 fill vetoes: to exclude deposits at IsoDeposit creation stage check concrete extractors if it's no-op !
virtual ~CaloExtractor ()

Static Public Member Functions

static GlobalPoint MuonAtCaloPosition (const reco::Track &muon, const double bz, const GlobalPoint &endpos, bool fixVxy=false, bool fixVz=false)
 Extrapolate muons to calorimeter-object positions.

Private Member Functions

double noiseEcal (const CaloTower &tower) const
double noiseHcal (const CaloTower &tower) const

Private Attributes

edm::InputTag theCaloTowerCollectionLabel
std::string theDepositLabel
double theDR_Max
double theDR_Veto_E
double theDR_Veto_H
double theThreshold_E
double theThreshold_H
std::vector< DetId > theVetoCollection
double theWeight_E
double theWeight_H
bool vertexConstraintFlag_XY
bool vertexConstraintFlag_Z

Detailed Description

Definition at line 18 of file CaloExtractor.h.

Constructor & Destructor Documentation

muonisolation::CaloExtractor::CaloExtractor (  )  [inline]

Definition at line 22 of file CaloExtractor.h.

00022 {};

CaloExtractor::CaloExtractor ( const edm::ParameterSet &  par  ) 

Definition at line 25 of file CaloExtractor.cc.

00025                                                     :
00026   theCaloTowerCollectionLabel(par.getParameter<edm::InputTag>("CaloTowerCollectionLabel")),
00027   theDepositLabel(par.getUntrackedParameter<string>("DepositLabel")),
00028   theWeight_E(par.getParameter<double>("Weight_E")),
00029   theWeight_H(par.getParameter<double>("Weight_H")),
00030   theThreshold_E(par.getParameter<double>("Threshold_E")),
00031   theThreshold_H(par.getParameter<double>("Threshold_H")),
00032   theDR_Veto_E(par.getParameter<double>("DR_Veto_E")),
00033   theDR_Veto_H(par.getParameter<double>("DR_Veto_H")),
00034   theDR_Max(par.getParameter<double>("DR_Max")),
00035   vertexConstraintFlag_XY(par.getParameter<bool>("Vertex_Constraint_XY")),
00036   vertexConstraintFlag_Z(par.getParameter<bool>("Vertex_Constraint_Z"))
00037 {
00038 }

virtual muonisolation::CaloExtractor::~CaloExtractor (  )  [inline, virtual]

Definition at line 25 of file CaloExtractor.h.

00025 {}

Member Function Documentation

IsoDeposit CaloExtractor::deposit ( const edm::Event &  ev,
const edm::EventSetup &  evSetup,
const reco::Track &  track 
) const [virtual]

make single IsoDeposit based on track as input purely virtual: have to implement in concrete implementations

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 92 of file CaloExtractor.cc.

References reco::deltaR(), PV3DBase< T, PVType, FrameType >::eta(), reco::TrackBase::eta(), find(), edm::EventSetup::get(), LogDebug, MuonAtCaloPosition(), noiseEcal(), noiseHcal(), PV3DBase< T, PVType, FrameType >::phi(), reco::TrackBase::phi(), reco::TrackBase::pt(), theCaloTowerCollectionLabel, theDR_Max, theDR_Veto_E, theDR_Veto_H, theThreshold_E, theThreshold_H, theVetoCollection, theWeight_E, theWeight_H, vertexConstraintFlag_XY, and vertexConstraintFlag_Z.

00093 {
00094   IsoDeposit dep(muon.eta(), muon.phi() );
00095   LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00096           << " >>> Muon: pt " << muon.pt()
00097           << " eta " << muon.eta()
00098           << " phi " << muon.phi();
00099 
00100   Handle<CaloTowerCollection> towers;
00101   event.getByLabel(theCaloTowerCollectionLabel,towers);
00102 
00103   edm::ESHandle<CaloGeometry> caloGeom;
00104   eventSetup.get<CaloGeometryRecord>().get(caloGeom);
00105 
00106   edm::ESHandle<MagneticField> bField;
00107   eventSetup.get<IdealMagneticFieldRecord>().get(bField);
00108   double bz = bField->inInverseGeV(GlobalPoint(0.,0.,0.)).z();
00109 
00110   CaloTowerCollection::const_iterator cal;
00111   CaloTowerCollection::const_iterator calEnd(towers->end());
00112   for ( cal = towers->begin(); cal != calEnd; ++cal ) {
00114       double dEta = fabs(muon.eta()-cal->eta());
00115       if (fabs(dEta) > theDR_Max) continue;
00116     
00117       double deltar0 = reco::deltaR(muon,*cal);
00118       if (deltar0>theDR_Max) continue;
00119 
00120       double etecal = cal->emEt();
00121       double eecal = cal->emEnergy();
00122       bool doEcal = theWeight_E>0 && etecal>theThreshold_E && eecal>3*noiseEcal(*cal);
00123       double ethcal = cal->hadEt();
00124       double ehcal = cal->hadEnergy();
00125       bool doHcal = theWeight_H>0 && ethcal>theThreshold_H && ehcal>3*noiseHcal(*cal);
00126       if ((!doEcal) && (!doHcal)) continue;
00127 
00128       DetId calId = cal->id();
00129       GlobalPoint endpos = caloGeom->getPosition(calId);
00130       GlobalPoint muatcal = MuonAtCaloPosition(muon,bz,endpos,vertexConstraintFlag_XY, vertexConstraintFlag_Z);
00131       double deltar = reco::deltaR(muatcal,endpos);
00132 
00133       if (deltar<theDR_Veto_H) {
00134               dep.setVeto(IsoDeposit::Veto(reco::isodeposit::Direction(muatcal.eta(), muatcal.phi()), theDR_Veto_H));
00135       }
00136 
00137       if (doEcal) {
00138             if (deltar<theDR_Veto_E) { 
00139                   double calodep = theWeight_E*etecal;
00140                   if (doHcal) calodep += theWeight_H*ethcal;
00141                   dep.addCandEnergy(calodep);
00142                     LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00143                     << " >>> Calo deposit inside veto (with ECAL): deltar " << deltar
00144                     << " calodep " << calodep
00145                     << " ecaldep " << etecal
00146                     << " hcaldep " << ethcal
00147                     << " eta " << cal->eta()
00148                     << " phi " << cal->phi();
00149                   continue;
00150             }
00151       } else {
00152             if (deltar<theDR_Veto_H) { 
00153                   dep.addCandEnergy(theWeight_H*ethcal);
00154                     LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00155                     << " >>> Calo deposit inside veto (no ECAL): deltar " << deltar
00156                     << " calodep " << theWeight_H*ethcal
00157                     << " eta " << cal->eta()
00158                     << " phi " << cal->phi();
00159                   continue;
00160             }
00161       }
00162 
00163       if (std::find(theVetoCollection.begin(), theVetoCollection.end()
00164                   , calId)!=theVetoCollection.end()) {
00165             LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00166             << " >>> Deposits belongs to other track: deltar, etecal, ethcal= " 
00167             << deltar << ", " << etecal << ", " << ethcal;
00168             continue;
00169       }
00170 
00171       if (doEcal) {
00172             if (deltar>theDR_Veto_E) { 
00173                   double calodep = theWeight_E*etecal;
00174                   if (doHcal) calodep += theWeight_H*ethcal;
00175                   dep.addDeposit(reco::isodeposit::Direction(endpos.eta(), endpos.phi()),calodep);
00176                     LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00177                     << " >>> Calo deposit (with ECAL): deltar " << deltar
00178                     << " calodep " << calodep
00179                     << " ecaldep " << etecal
00180                     << " hcaldep " << ethcal
00181                     << " eta " << cal->eta()
00182                     << " phi " << cal->phi();
00183             }
00184       } else {
00185             if (deltar>theDR_Veto_H) { 
00186                     dep.addDeposit(reco::isodeposit::Direction(endpos.eta(), endpos.phi()),theWeight_H*ethcal);
00187                     LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00188                     << " >>> Calo deposit (no ECAL): deltar " << deltar
00189                     << " calodep " << theWeight_H*ethcal
00190                     << " eta " << cal->eta()
00191                     << " phi " << cal->phi();
00192             }
00193       }
00194   }
00195 
00196   return dep;
00197 
00198 }

void CaloExtractor::fillVetos ( const edm::Event &  ev,
const edm::EventSetup &  evSetup,
const reco::TrackCollection &  tracks 
) [virtual]

fill vetoes: to exclude deposits at IsoDeposit creation stage check concrete extractors if it's no-op !

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 40 of file CaloExtractor.cc.

References reco::deltaR(), edm::EventSetup::get(), MuonAtCaloPosition(), noiseEcal(), noiseHcal(), theCaloTowerCollectionLabel, theDR_Max, theDR_Veto_E, theDR_Veto_H, theThreshold_E, theThreshold_H, theVetoCollection, theWeight_E, theWeight_H, vertexConstraintFlag_XY, and vertexConstraintFlag_Z.

00041 {
00042   theVetoCollection.clear();
00043 
00044   Handle<CaloTowerCollection> towers;
00045   event.getByLabel(theCaloTowerCollectionLabel,towers);
00046 
00047   edm::ESHandle<CaloGeometry> caloGeom;
00048   eventSetup.get<CaloGeometryRecord>().get(caloGeom);
00049 
00050   edm::ESHandle<MagneticField> bField;
00051   eventSetup.get<IdealMagneticFieldRecord>().get(bField);
00052   double bz = bField->inInverseGeV(GlobalPoint(0.,0.,0.)).z();
00053 
00054   TrackCollection::const_iterator mu;
00055   TrackCollection::const_iterator muEnd(muons.end());
00056   
00057   CaloTowerCollection::const_iterator cal;
00058   CaloTowerCollection::const_iterator calEnd(towers->end());
00059   
00060   for ( mu = muons.begin(); mu != muEnd; ++mu ) {
00061     for ( cal = towers->begin(); cal != calEnd; ++cal ) {
00063       double dEta = fabs(mu->eta()-cal->eta());
00064       if (fabs(dEta) > theDR_Max) continue;
00065 
00066             double deltar0 = reco::deltaR(*mu,*cal);
00067             if (deltar0>theDR_Max) continue;
00068 
00069             double etecal = cal->emEt();
00070             double eecal = cal->emEnergy();
00071             bool doEcal = theWeight_E>0 && etecal>theThreshold_E && eecal>3*noiseEcal(*cal);
00072             double ethcal = cal->hadEt();
00073             double ehcal = cal->hadEnergy();
00074             bool doHcal = theWeight_H>0 && ethcal>theThreshold_H && ehcal>3*noiseHcal(*cal);
00075             if ((!doEcal) && (!doHcal)) continue;
00076 
00077             DetId calId = cal->id();
00078             GlobalPoint endpos = caloGeom->getPosition(calId);
00079             GlobalPoint muatcal = MuonAtCaloPosition(*mu,bz,endpos, vertexConstraintFlag_XY, vertexConstraintFlag_Z);
00080             double deltar = reco::deltaR(muatcal,endpos);
00081 
00082             if (doEcal) {
00083                   if (deltar<theDR_Veto_E) theVetoCollection.push_back(calId);
00084             } else {
00085                   if (deltar<theDR_Veto_H) theVetoCollection.push_back(calId);
00086             }
00087       }
00088   }
00089      
00090 }

GlobalPoint CaloExtractor::MuonAtCaloPosition ( const reco::Track &  muon,
const double  bz,
const GlobalPoint &  endpos,
bool  fixVxy = false,
bool  fixVz = false 
) [static]

Extrapolate muons to calorimeter-object positions.

Definition at line 200 of file CaloExtractor.cc.

References reco::TrackBase::charge(), funct::cos(), reco::TrackBase::covariance(), reco::TrackBase::dsz(), reco::TrackBase::dxy(), reco::TrackBase::dz(), reco::TrackBase::i_dsz, reco::TrackBase::i_dxy, reco::TrackBase::i_lambda, reco::TrackBase::i_phi, reco::TrackBase::i_qoverp, normalizedPhi(), reco::TrackBase::p(), reco::TrackBase::phi(), funct::pow(), reco::TrackBase::pt(), reco::TrackBase::qoverp(), funct::sin(), funct::tan(), theta, reco::TrackBase::theta(), reco::TrackBase::vx(), reco::TrackBase::vy(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), y, PV3DBase< T, PVType, FrameType >::z(), and z.

Referenced by deposit(), and fillVetos().

00200                                                                                                                                     {
00201       double qoverp= muon.qoverp();
00202       double cur = bz*muon.charge()/muon.pt();
00203       double phi0 = muon.phi();
00204       double dca = muon.dxy();
00205       double theta = muon.theta();
00206       double dz = muon.dz();
00207 
00208       //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00209           //<< " Pt(GeV): " <<  muon.pt()
00210           //<< ", phi0 " <<  muon.phi0()
00211              //<< ", eta " <<  muon.eta();
00212       //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00213           //<< " d0 " <<  muon.d0()
00214           //<< ", dz " <<  muon.dz();
00215       //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
00216           //<< " rhocal " <<  endpos.perp()
00217           //<< ", zcal " <<  endpos.z();
00218 
00219       if (fixVxy && fixVz) {
00220             // Note that here we assume no correlation between XY and Z projections
00221             // This should be a reasonable approximation for our purposes
00222             double errd02 = muon.covariance(muon.i_dxy,muon.i_dxy);
00223             if (pow(muon.dxy(),2)<4*errd02) {
00224                   phi0 -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_phi)
00225                                      /errd02;
00226                   cur -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_qoverp)
00227                                      /errd02 * (cur/qoverp);
00228                   dca = 0;
00229             } 
00230             double errdsz2 = muon.covariance(muon.i_dsz,muon.i_dsz);
00231             if (pow(muon.dsz(),2)<4*errdsz2) {
00232                   theta += muon.dsz()*muon.covariance(muon.i_dsz,muon.i_lambda)
00233                                      /errdsz2;
00234                   dz = 0;
00235             } 
00236       } else if (fixVxy) {
00237             double errd02 = muon.covariance(muon.i_dxy,muon.i_dxy);
00238             if (pow(muon.dxy(),2)<4*errd02) {
00239                   phi0  -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_phi)
00240                                      /errd02;
00241                   cur -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_qoverp)
00242                                      /errd02 * (cur/qoverp);
00243                   theta += muon.dxy()*muon.covariance(muon.i_dxy,muon.i_lambda)
00244                                      /errd02;
00245                   dz    -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_dsz)
00246                                      /errd02 * muon.p()/muon.pt();
00247                   dca = 0;
00248             } 
00249       } else if (fixVz) {
00250             double errdsz2 = muon.covariance(muon.i_dsz,muon.i_dsz);
00251             if (pow(muon.dsz(),2)<4*errdsz2) {
00252                   theta += muon.dsz()*muon.covariance(muon.i_dsz,muon.i_lambda)
00253                                      /errdsz2;
00254                   phi0  -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_phi)
00255                                      /errdsz2;
00256                   cur -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_qoverp)
00257                                      /errdsz2 * (cur/qoverp);
00258                   dca   -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_dxy)
00259                                      /errdsz2;
00260                   dz = 0;
00261             } 
00262       }
00263 
00264       double sphi0 = sin(phi0);
00265       double cphi0 = cos(phi0);
00266 
00267       double xsin =  endpos.x()*sphi0 - endpos.y()*cphi0;
00268       double xcos =  endpos.x()*cphi0 + endpos.y()*sphi0;
00269       double fcdca = fabs(1-cur*dca);
00270       double phif = atan2( fcdca*sphi0-cur*endpos.x()
00271                          , fcdca*cphi0+cur*endpos.y());
00272       double tphif2 = tan(0.5*(phif-phi0));
00273       double dcaf = dca + xsin + xcos*tphif2;
00274 
00275       double x = endpos.x() - dcaf*sin(phif);
00276       double y = endpos.y() + dcaf*cos(phif);
00277 
00278       double deltas =  (x-muon.vx())*cphi0 + (y-muon.vy())*sphi0;
00279       double deltaphi = normalizedPhi(phif-phi0);
00280       if (deltaphi!=0) deltas = deltas*deltaphi/sin(deltaphi);
00281 
00282       double z =dz;
00283       double tantheta = tan(theta);
00284       if (tantheta!=0) {
00285             z += deltas/tan(theta);
00286       } else {
00287             z = endpos.z();
00288       }
00289 
00290       return GlobalPoint(x,y,z);
00291 }

double CaloExtractor::noiseEcal ( const CaloTower &  tower  )  const [private]

Definition at line 293 of file CaloExtractor.cc.

References reco::Particle::eta(), and eta.

Referenced by deposit(), and fillVetos().

00293                                                             {
00294       double noise = 0.04;
00295       double eta = tower.eta();
00296       if (fabs(eta)>1.479) noise = 0.15;
00297       return noise;
00298 }

double CaloExtractor::noiseHcal ( const CaloTower &  tower  )  const [private]

Definition at line 300 of file CaloExtractor.cc.

Referenced by deposit(), and fillVetos().

00300                                                             {
00301       double noise = 0.2;
00302       return noise;
00303 }

Member Data Documentation

edm::InputTag muonisolation::CaloExtractor::theCaloTowerCollectionLabel [private]

Definition at line 35 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

std::string muonisolation::CaloExtractor::theDepositLabel [private]

Definition at line 38 of file CaloExtractor.h.

double muonisolation::CaloExtractor::theDR_Max [private]

Definition at line 47 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theDR_Veto_E [private]

Definition at line 45 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theDR_Veto_H [private]

Definition at line 46 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theThreshold_E [private]

Definition at line 43 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theThreshold_H [private]

Definition at line 44 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

std::vector<DetId> muonisolation::CaloExtractor::theVetoCollection [private]

Definition at line 52 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theWeight_E [private]

Definition at line 41 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theWeight_H [private]

Definition at line 42 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

bool muonisolation::CaloExtractor::vertexConstraintFlag_XY [private]

Definition at line 48 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

bool muonisolation::CaloExtractor::vertexConstraintFlag_Z [private]

Definition at line 49 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:49:13 2009 for CMSSW by  doxygen 1.5.4