reco::helper::CastorJetIDHelper Class Reference

#include <CastorJetIDHelper.h>

Public Member Functions
void	calculate (const edm::Event &event, const reco::BasicJet &jet)
	CastorJetIDHelper ()
double	depth () const
double	emEnergy () const
double	fem () const
double	fhot () const
double	hadEnergy () const
void	initValues ()
int	nTowers () const
double	sigmaz () const
double	width () const
	~CastorJetIDHelper ()

Private Member Functions
double	phiangle (double testphi)

Private Attributes
double	depth_
double	emEnergy_
double	fem_
double	fhot_
double	hadEnergy_
int	nTowers_
double	sigmaz_
double	width_

Static Private Attributes
static int	sanity_checks_left_

Detailed Description

Definition at line 13 of file CastorJetIDHelper.h.

Constructor & Destructor Documentation

reco::helper::CastorJetIDHelper::CastorJetIDHelper

(

)

Definition at line 14 of file CastorJetIDHelper.cc.

References initValues().

{
 
  initValues();
}

reco::helper::CastorJetIDHelper::~CastorJetIDHelper ( )

inline

Definition at line 18 of file CastorJetIDHelper.h.

{} 

Member Function Documentation

void reco::helper::CastorJetIDHelper::calculate	(	const edm::Event &	event,
		const reco::BasicJet &	jet
	)

Definition at line 33 of file CastorJetIDHelper.cc.

References reco::CastorTower::depth(), reco::CastorTower::emEnergy(), reco::CastorTower::energy(), reco::LeafCandidate::energy(), reco::CastorTower::fhot(), reco::Jet::getJetConstituents(), reco::CastorTower::hadEnergy(), edm::Ref< C, T, F >::id(), module::module(), reco::CastorTower::phi(), reco::LeafCandidate::phi(), funct::pow(), reco::CastorTower::rechitsBegin(), reco::CastorTower::rechitsEnd(), and mathSSE::sqrt().

Referenced by CastorJetIDProducer::produce().

{
  initValues();
  
  // calculate Castor JetID properties
 
        double zmean = 0.;
        double z2mean = 0.;
    
        std::vector<CandidatePtr> ccp = jet.getJetConstituents();
        std::vector<CandidatePtr>::const_iterator itParticle;
        for (itParticle=ccp.begin();itParticle!=ccp.end();++itParticle){        
                const CastorTower* castorcand = dynamic_cast<const CastorTower*>(itParticle->get());
            emEnergy_ += castorcand->emEnergy();
            hadEnergy_ += castorcand->hadEnergy();
            depth_ += castorcand->depth()*castorcand->energy();
            width_ += pow(phiangle(castorcand->phi() - jet.phi()),2)*castorcand->energy();
                fhot_ += castorcand->fhot()*castorcand->energy();
            
            // loop over rechits
                for (edm::RefVector<edm::SortedCollection<CastorRecHit> >::iterator it = castorcand->rechitsBegin(); it != castorcand->rechitsEnd(); it++) {
                             edm::Ref<edm::SortedCollection<CastorRecHit> > rechit_p = *it;                         
                             double Erechit = rechit_p->energy();
                             HcalCastorDetId id = rechit_p->id();
                             int module = id.module();                                  
                                 double zrechit = 0;     
                                 if (module < 3) zrechit = -14390 - 24.75 - 49.5*(module-1);     
                                 if (module > 2) zrechit = -14390 - 99 - 49.5 - 99*(module-3);   
                                 zmean += Erechit*zrechit;   
                                 z2mean += Erechit*zrechit*zrechit;
                } // end loop over rechits
            
            nTowers_++;
        }
        //cout << "" << endl;
        
        depth_ = depth_/jet.energy();
        width_ = sqrt(width_/jet.energy());
        fhot_ = fhot_/jet.energy();
        fem_ = emEnergy_/jet.energy();
        
        zmean = zmean/jet.energy();
            z2mean = z2mean/jet.energy();
            double sigmaz2 = z2mean - zmean*zmean;
            if(sigmaz2 > 0) sigmaz_ = sqrt(sigmaz2);

  
}

double reco::helper::CastorJetIDHelper::depth ( void  ) const

inline

Definition at line 32 of file CastorJetIDHelper.h.

References depth_.

Referenced by CastorJetIDProducer::produce().

{ return    depth_;}  

double reco::helper::CastorJetIDHelper::emEnergy ( ) const

inline

Definition at line 28 of file CastorJetIDHelper.h.

References emEnergy_.

Referenced by CastorJetIDProducer::produce().

{ return    emEnergy_;}           

double reco::helper::CastorJetIDHelper::fem ( ) const

inline

Definition at line 30 of file CastorJetIDHelper.h.

References fem_.

Referenced by CastorJetIDProducer::produce().

{ return    fem_;}

double reco::helper::CastorJetIDHelper::fhot ( ) const

inline

Definition at line 33 of file CastorJetIDHelper.h.

References fhot_.

Referenced by CastorJetIDProducer::produce().

{ return    fhot_;}  

double reco::helper::CastorJetIDHelper::hadEnergy ( ) const

inline

Definition at line 29 of file CastorJetIDHelper.h.

References hadEnergy_.

Referenced by CastorJetIDProducer::produce().

{ return    hadEnergy_;}           

void reco::helper::CastorJetIDHelper::initValues

(

)

Definition at line 20 of file CastorJetIDHelper.cc.

Referenced by CastorJetIDHelper().

{
      emEnergy_ = 0.0;
      hadEnergy_ = 0.0; 
      fem_ = 0.0;
      width_ = 0.0;
      depth_ = 0.0;
      fhot_ = 0.0;
      sigmaz_ = 0.0;
      nTowers_ = 0;
}

int reco::helper::CastorJetIDHelper::nTowers ( ) const

inline

Definition at line 35 of file CastorJetIDHelper.h.

References nTowers_.

Referenced by CastorJetIDProducer::produce().

{ return    nTowers_;}

double reco::helper::CastorJetIDHelper::phiangle ( double  testphi )

private

Definition at line 83 of file CastorJetIDHelper.cc.

References M_PI, and phi.

                                                            {
  double phi = testphi;
  while (phi>M_PI) phi -= (2*M_PI);
  while (phi<-M_PI) phi += (2*M_PI);
  return phi;
}

double reco::helper::CastorJetIDHelper::sigmaz ( ) const

inline

Definition at line 34 of file CastorJetIDHelper.h.

References sigmaz_.

Referenced by CastorJetIDProducer::produce().

{ return    sigmaz_;}  

double reco::helper::CastorJetIDHelper::width ( ) const

inline

Definition at line 31 of file CastorJetIDHelper.h.

References width_.

Referenced by Vispa.Main.MainWindow.MainWindow::_saveIni(), CastorJetIDProducer::produce(), Vispa.Views.PropertyView.PropertyView::resizeEvent(), Vispa.Views.PropertyView.PropertyView::sectionResized(), and Vispa.Main.MainWindow.MainWindow::updateStartupScreenGeometry().

{ return    width_;}  

Member Data Documentation

double reco::helper::CastorJetIDHelper::depth_

private

Definition at line 48 of file CastorJetIDHelper.h.

Referenced by depth().

double reco::helper::CastorJetIDHelper::emEnergy_

private

Definition at line 44 of file CastorJetIDHelper.h.

Referenced by emEnergy().

double reco::helper::CastorJetIDHelper::fem_

private

Definition at line 46 of file CastorJetIDHelper.h.

Referenced by fem().

double reco::helper::CastorJetIDHelper::fhot_

private

Definition at line 49 of file CastorJetIDHelper.h.

Referenced by fhot().

double reco::helper::CastorJetIDHelper::hadEnergy_

private

Definition at line 45 of file CastorJetIDHelper.h.

Referenced by hadEnergy().

int reco::helper::CastorJetIDHelper::nTowers_

private

Definition at line 51 of file CastorJetIDHelper.h.

Referenced by nTowers().

int reco::helper::CastorJetIDHelper::sanity_checks_left_

staticprivate

Definition at line 53 of file CastorJetIDHelper.h.

double reco::helper::CastorJetIDHelper::sigmaz_

private

Definition at line 50 of file CastorJetIDHelper.h.

Referenced by sigmaz().

double reco::helper::CastorJetIDHelper::width_

private

Definition at line 47 of file CastorJetIDHelper.h.

Referenced by width().