TxCalculator Class Reference

#include <TxCalculator.h>

Public Member Functions
double	getCTx (const reco::Photon &clus, double i, double threshold, double innerDR=0, double effRatio=2)
double	getJct (const reco::Photon &cluster, double r1=0.4, double r2=0.04, double jWidth=0.015, double threshold=2)
double	getJt (const reco::Photon &cluster, double r1=0.4, double r2=0.04, double jWidth=0.015, double threshold=2)
double	getJurassicArea (double r1, double r2, double width)
double	getMPT (double ptCut=0, double etaCut=1000)
double	getTx (const reco::Photon &clus, double i, double threshold, double innerDR=0, double effRatio=2)
	TxCalculator (const edm::Event &iEvent, const edm::EventSetup &iSetup, const edm::InputTag &trackLabel)

Private Member Functions
double	calcDphi (double phi1_, double phi2_)
double	dRDistance (double eta1, double phi1, double eta2, double phi2)

Private Attributes
edm::Handle < reco::TrackCollection >	recCollection
CLHEP::RandFlat *	theDice

Detailed Description

Definition at line 27 of file TxCalculator.h.

Constructor & Destructor Documentation

TxCalculator::TxCalculator	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const edm::InputTag &	trackLabel
	)

Definition at line 24 of file TxCalculator.cc.

References edm::hlt::Exception, edm::Event::getByLabel(), edm::RandomNumberGenerator::getEngine(), and edm::Service< T >::isAvailable().

{
   iEvent.getByLabel(trackLabel, recCollection); 
   edm::Service<edm::RandomNumberGenerator> rng;
   if ( ! rng.isAvailable()) {
      throw cms::Exception("Configuration")
         << "XXXXXXX requires the RandomNumberGeneratorService\n"
         "which is not present in the configuration file.  You must add the service\n"
         "in the configuration file or remove the modules that require it.";
   }
   CLHEP::HepRandomEngine& engine = rng->getEngine();
   theDice = new CLHEP::RandFlat(engine, 0, 1);
   
} 

Member Function Documentation

double TxCalculator::calcDphi ( double  phi1_, double  phi2_  )

inlineprivate

Definition at line 54 of file TxCalculator.h.

References PI.

Referenced by dRDistance().

   {
       double dphi=phi1_-phi2_;

      if (dphi>0){
         while (dphi>2*PI) dphi-=2*PI;
         if (dphi>PI) dphi=2*PI-dphi; 
      } else {
         while (dphi<-2*PI) dphi+=2*PI;
         if (dphi<-PI) dphi=-2*PI-dphi;
      }
      return dphi;
   }

double TxCalculator::dRDistance ( double  eta1, double  phi1, double  eta2, double  phi2  )

inlineprivate

Definition at line 46 of file TxCalculator.h.

References calcDphi(), and mathSSE::sqrt().

   {
      double deta = eta1 - eta2;
      double dphi = (calcDphi(phi1, phi2));
      
      return sqrt(deta * deta + dphi * dphi);
   }

double TxCalculator::getCTx	(	const reco::Photon &	clus,
		double	i,
		double	threshold,
		double	innerDR = 0,
		double	effRatio = 2
	)

Definition at line 124 of file TxCalculator.cc.

References reco::LeafCandidate::eta(), reco::LeafCandidate::phi(), RecoTauCleanerPlugins::pt, dt_dqm_sourceclient_common_cff::reco, and x.

{
   using namespace edm;
   using namespace reco;

   double SClusterEta = cluster.eta();
   double SClusterPhi = cluster.phi();
   double TotalPt = 0;

   TotalPt = 0;

   for(reco::TrackCollection::const_iterator
      recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
      {
     double diceNum = theDice->fire();
         if ( (effRatio < 1 ) &&  ( diceNum > effRatio))
            continue;
     
     
     double pt = recTrack->pt();
     double eta2 = recTrack->eta();
     double phi2 = recTrack->phi();
     double dEta = fabs(eta2-SClusterEta);
     
      if(dEta >= 0.1 * x)
         continue;
      if(dRDistance(SClusterEta,SClusterPhi,eta2,phi2) < innerDR)
         continue;
      
      if(pt > threshold)
         TotalPt = TotalPt + pt;
   }
   
   double Tx = getTx(cluster,x,threshold,innerDR,effRatio);
   double CTx = Tx - TotalPt / 40.0 * x;

   return CTx;
}

double TxCalculator::getJct	(	const reco::Photon &	cluster,
		double	r1 = 0.4,
		double	r2 = 0.04,
		double	jWidth = 0.015,
		double	threshold = 2
	)

Definition at line 203 of file TxCalculator.cc.

References gather_cfg::cout, dPhi(), reco::LeafCandidate::eta(), eta(), phi, reco::LeafCandidate::phi(), PI, RecoTauCleanerPlugins::pt, diffTwoXMLs::r1, and dt_dqm_sourceclient_common_cff::reco.

{

   using namespace edm;
   using namespace reco;


   double SClusterEta = cluster.eta();
   double SClusterPhi = cluster.phi();
   double TotalPt = 0;

   for(reco::TrackCollection::const_iterator
          recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
      {
         double pt = recTrack->pt();
         double eta = recTrack->eta();
         double phi = recTrack->phi();
         double dEta = fabs(eta-SClusterEta);
         double dPhi = phi-SClusterPhi;
         if ( dPhi < -PI )    dPhi = dPhi + 2*PI ;
         if ( dPhi >  PI )    dPhi = dPhi - 2*PI ;
     if ( fabs(dPhi) >PI )   cout << " error!!! dphi > 2pi   : " << dPhi << endl;
     //         double dR = sqrt(dEta*dEta+dPhi*dPhi);
         

     if ( fabs(dEta) > r1 ) continue;
     if ( fabs(dPhi) <r1 ) continue;
     
     if(pt > threshold)
            TotalPt = TotalPt + pt;
      }
   
   double areaStrip = 4*PI*r1 -  4*r1*r1;
   double areaJura  = getJurassicArea(r1,r2, jWidth) ;
   double theCJ     = getJt(cluster,r1,  r2, jWidth, threshold);

   double theCCJ   = theCJ - TotalPt * areaJura / areaStrip ;
   return theCCJ;

}

double TxCalculator::getJt	(	const reco::Photon &	cluster,
		double	r1 = 0.4,
		double	r2 = 0.04,
		double	jWidth = 0.015,
		double	threshold = 2
	)

Definition at line 165 of file TxCalculator.cc.

References gather_cfg::cout, dPhi(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::LeafCandidate::eta(), eta(), phi, reco::LeafCandidate::phi(), PI, RecoTauCleanerPlugins::pt, dt_dqm_sourceclient_common_cff::reco, and mathSSE::sqrt().

{

   using namespace edm;
   using namespace reco;


   double SClusterEta = cluster.eta();
   double SClusterPhi = cluster.phi();
   double TotalPt = 0;

   for(reco::TrackCollection::const_iterator
          recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
      {
         double pt = recTrack->pt();
     double eta = recTrack->eta();
         double phi = recTrack->phi();
     double dEta = fabs(eta-SClusterEta);
     double dPhi = phi-SClusterPhi;
     if ( dPhi < -PI )    dPhi = dPhi + 2*PI ;
     if ( dPhi >  PI )    dPhi = dPhi - 2*PI ;
     if ( fabs(dPhi) >PI )   cout << " error!!! dphi > 2pi   : " << dPhi << endl;
     double dR = sqrt(dEta*dEta+dPhi*dPhi);
     // Jurassic Cone /////                                                                                                     
     if ( dR > r1 ) continue;
     if ( dR < r2 ) continue;
     if ( fabs(dEta) <  jWidth)  continue;
     // stupid bug if ( fabs(dPhi) >  jWidth)  continue;
     
     if(pt > threshold)
        TotalPt = TotalPt + pt;
      }

   return TotalPt;
}

double TxCalculator::getJurassicArea	(	double	r1,
		double	r2,
		double	width
	)

Definition at line 40 of file TxCalculator.cc.

References funct::cos(), and mathSSE::sqrt().

                                                                        {
   
   float theta1 = asin( width / r1);
   float theta2 = asin( width / r2);
   float theA   = sqrt ( r1*r1 + r2*r2 - 2 * r1 * r2 * cos ( theta1 - theta2) );
   float area1 =  0.5 * r1*r1 * ( 3.141592 - 2 * theta1 )   ;
   float area2 =  0.5 * r2*r2 * ( 3.141592 - 2 * theta2 )   ;
   float area3 =  width * theA;
   float finalArea = 2 * ( area1 - area2 - area3);
   return finalArea;
}

double TxCalculator::getMPT	(	double	ptCut = 0,
		double	etaCut = 1000
	)

Definition at line 53 of file TxCalculator.cc.

References eta(), RecoTauCleanerPlugins::pt, dt_dqm_sourceclient_common_cff::reco, and mathSSE::sqrt().

{
   using namespace edm;
   using namespace reco;
   
   double sumpx(0), sumpy(0);
   
   for(reco::TrackCollection::const_iterator
          recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
      {
     double pt = recTrack->pt();
         double eta = recTrack->eta();
     
     if(pt < ptCut ) 
        continue;
     if ( fabs( eta) > etaCut ) 
        continue;
     
     double pxTemp = recTrack->px();
     double pyTemp = recTrack->py();
     
     sumpx = sumpx + pxTemp;
     sumpy = sumpy + pyTemp;
     //  cout << " pt  = " << recTrack->pt() <<  "    and px = " << pxTemp << " and  py = " << pyTemp << endl;
      }
   //   cout << " square = " << sumpx*sumpx + sumpy*sumpy << endl;
   double theMPT = sqrt(sumpx*sumpx + sumpy*sumpy) ;
   //  cout << " mpt    = "<< theMPT << endl;
   
   return theMPT;
}

double TxCalculator::getTx	(	const reco::Photon &	clus,
		double	i,
		double	threshold,
		double	innerDR = 0,
		double	effRatio = 2
	)

Definition at line 86 of file TxCalculator.cc.

References reco::LeafCandidate::eta(), reco::LeafCandidate::phi(), RecoTauCleanerPlugins::pt, and dt_dqm_sourceclient_common_cff::reco.

{

   using namespace edm;
   using namespace reco;

   
   
   double SClusterEta = cluster.eta();
   double SClusterPhi = cluster.phi();
   double TotalPt = 0;
   
   for(reco::TrackCollection::const_iterator
      recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
      {
     double diceNum = theDice->fire();
     if ( (effRatio < 1 ) &&  ( diceNum > effRatio))
        continue;
     
     double pt = recTrack->pt();
     double eta2 = recTrack->eta();
     double phi2 = recTrack->phi();
     
     if(dRDistance(SClusterEta,SClusterPhi,eta2,phi2) >= 0.1 * x)
        continue;
     if(dRDistance(SClusterEta,SClusterPhi,eta2,phi2) < innerDR)
        continue;
     if(pt > threshold)
        TotalPt = TotalPt + pt;
      }
   
   return TotalPt;
}

Member Data Documentation

edm::Handle<reco::TrackCollection> TxCalculator::recCollection

private

Definition at line 43 of file TxCalculator.h.

CLHEP::RandFlat* TxCalculator::theDice

private

Definition at line 44 of file TxCalculator.h.