#include <SimG4HcalHitJetFinder.h>

Public Member Functions
std::vector < SimG4HcalHitCluster > *	getClusters (bool)

double	rDist (const SimG4HcalHitCluster , const CaloHit ) const

double	rDist (const double, const double, const double, const double) const

void	setCone (double)

void	setInput (std::vector< CaloHit > *)

	SimG4HcalHitJetFinder (double cone=0.5)

virtual	~SimG4HcalHitJetFinder ()

Private Attributes
std::vector< SimG4HcalHitCluster >	clusvector

std::vector< CaloHit >	input

double	jetcone

Detailed Description

Definition at line 13 of file SimG4HcalHitJetFinder.h.

Constructor & Destructor Documentation

SimG4HcalHitJetFinder::SimG4HcalHitJetFinder ( double cone = 0.5 )

Definition at line 13 of file SimG4HcalHitJetFinder.cc.

13 : jetcone(cone) {}

SimG4HcalHitJetFinder::jetcone

double jetcone

Definition: SimG4HcalHitJetFinder.h:25

SimG4HcalHitJetFinder::~SimG4HcalHitJetFinder ( )

virtual

Definition at line 15 of file SimG4HcalHitJetFinder.cc.

15 { edm::LogInfo("ValidHcal") << "SimG4HcalHitJetFinder:: Deleting"; }

edm::LogInfo

Log< level::Info, false > LogInfo

Definition: MessageLogger.h:125

Member Function Documentation

std::vector< SimG4HcalHitCluster > * SimG4HcalHitJetFinder::getClusters ( bool hcal_only )

Definition at line 21 of file SimG4HcalHitJetFinder.cc.

References haddnano::cl, clusvector, ztail::d, HcalBarrel, HcalEndcap, HcalForward, input, dqmiolumiharvest::j, jetcone, rDist(), and groupFilesInBlocks::temp.

                                                                                  {
   clusvector.erase(clusvector.begin(), clusvector.end());
   if (input.empty()) {
     return &clusvector;
   }
 
   std::vector<CaloHit>::iterator itr;
   for (itr = input.begin(); itr != input.end(); itr++) {
     edm::LogVerbatim("ValidHcal") << "HcalHitJetFinder::getClusters_1 - input :  e " << itr->e() << "  eta "
                                   << itr->eta() << "  phi " << itr->phi() << "  subdet " << itr->det();
   }
 
   sort(input.begin(), input.end());  // sort input in descending order
 
   for (itr = input.begin(); itr != input.end(); itr++) {
     edm::LogVerbatim("ValidHcal") << "HcalHitJetFinder::getClusters_2 - input :  e " << itr->e() << "  eta "
                                   << itr->eta() << "  phi " << itr->phi() << "  subdet " << itr->det();
   }
 
   std::vector<SimG4HcalHitCluster> temp;  // dummy container for clusters
 
   //  first input hit -> first cluster
 
   CaloHit hit;
   SimG4HcalHitCluster cluster;
 
   std::vector<CaloHit>::iterator itr_hits;
 
   int j, first_seed = 0;
   for (j = 0, itr_hits = input.begin(); itr_hits != input.end(); j++, itr_hits++) {
     int h_type = itr_hits->det();  // if desired HCAL hits (only) clusterfinding
     if (((h_type == static_cast<int>(HcalBarrel) || h_type == static_cast<int>(HcalEndcap) ||
           h_type == static_cast<int>(HcalForward)) &&
          hcal_only) ||
         (!hcal_only)) {
       cluster += input[j];
       edm::LogVerbatim("ValidHcal") << "HcalHitJetFinder:: First seed hit ..................\n" << (*itr_hits);
       first_seed = j;
       break;
     }
   }
 
   temp.push_back(cluster);
 
   std::vector<SimG4HcalHitCluster>::iterator itr_clus;
 
   for (j = 0, itr_hits = input.begin(); itr_hits != input.end(); j++, itr_hits++) {
     int h_type = itr_hits->det();  // if desired HCAL hits (only) clusterfinding
     if ((((h_type == static_cast<int>(HcalBarrel) || h_type == static_cast<int>(HcalEndcap) ||
            h_type == static_cast<int>(HcalForward)) &&
           hcal_only) ||
          (!hcal_only)) &&
         (j != first_seed)) {
       edm::LogVerbatim("ValidHcal") << "HcalHitJetFinder:: ........... Consider hit ..................\n"
                                     << (*itr_hits);
 
       int incl = 0;  // if the hit is included in one of clusters
 
       int iclus;
       for (itr_clus = temp.begin(), iclus = 0; itr_clus != temp.end(); itr_clus++, iclus++) {
         edm::LogVerbatim("ValidHcal") << "HcalHitJetFinder::=======> Cluster " << iclus << "\n" << (*itr_clus);
 
         double d = rDist(&(*itr_clus), &(*itr_hits));
         if (d < jetcone) {
           edm::LogVerbatim("ValidHcal") << "HcalHitJetFinder:: -> associated ... ";
           temp[iclus] += *itr_hits;
           incl = 1;
           break;
         }
       }
 
       // to here jumps "break"
       if (incl == 0) {
         SimG4HcalHitCluster cl;
         cl += *itr_hits;
         temp.push_back(cl);
         edm::LogVerbatim("ValidHcal") << "HcalHitJetFinder:: ************ NEW CLUSTER !\n" << cl;
       }
     }
   }
 
   clusvector = temp;
   return &clusvector;
 }

double SimG4HcalHitJetFinder::rDist	(	const SimG4HcalHitCluster *	cluster,
		const CaloHit *	hit
	)		const

Definition at line 106 of file SimG4HcalHitJetFinder.cc.

References SimG4HcalHitCluster::eta(), CaloHit::eta(), SimG4HcalHitCluster::phi(), and CaloHit::phi().

Referenced by getClusters().

                                                                                                 {
   double etac = cluster->eta();
   double phic = cluster->phi();
 
   double etah = hit->eta();
   double phih = hit->phi();
 
   return rDist(etac, phic, etah, phih);
 }

double SimG4HcalHitJetFinder::rDist	(	const double	etac,
		const double	phic,
		const double	etah,
		const double	phih
	)		const

Definition at line 116 of file SimG4HcalHitJetFinder.cc.

References kinem::delta_eta(), kinem::delta_phi(), M_PI, mathSSE::sqrt(), and createJobs::tmp.

                                                                                                                     {
   double delta_eta = etac - etah;
   double delta_phi = phic - phih;
 
   if (phic < phih)
     delta_phi = phih - phic;
   if (delta_phi > M_PI)
     delta_phi = 2 * M_PI - delta_phi;
 
   double tmp = sqrt(delta_eta * delta_eta + delta_phi * delta_phi);
 
   edm::LogVerbatim("ValidHcal") << "HcalHitJetFinder::rDist:\n Clus. eta, phi = " << etac << " " << phic
                                 << "\n hit   eta, phi = " << etah << " " << phih << " rDist = " << tmp;
 
   return tmp;
 }