CmsMTDLevelBuilder::ExtractPhiModule Struct Reference

#include <CmsMTDLevelBuilder.h>

Inheritance diagram for CmsMTDLevelBuilder::ExtractPhiModule:

Public Member Functions
double	operator() (const GeometricTimingDet *a) const

Detailed Description

Definition at line 70 of file CmsMTDLevelBuilder.h.

Member Function Documentation

double CmsMTDLevelBuilder::ExtractPhiModule::operator() ( const GeometricTimingDet *  a ) const

inline

Definition at line 71 of file CmsMTDLevelBuilder.h.

References AlCaHLTBitMon_QueryRunRegistry::comp, GeometricTimingDet::components(), objects.autophobj::float, mps_fire::i, phi, pi, and groupFilesInBlocks::temp.

                                                       {
      const double pi = 3.141592653592;
      std::vector<const GeometricTimingDet*> const & comp = a->components().back()->components();
      float phi = 0.;
      bool sum = true;
      
      for(auto i : comp){
    if(fabs(i->phi())>pi/2.) { 
      sum = false;
      break;
    }
      }
      
      if(sum){
    for(auto i : comp){
      phi+= i->phi();
    }
    
    double temp = phi/float(comp.size()) < 0. ? 
      2*pi + phi/float(comp.size()):
      phi/float(comp.size());
    return temp;
    
      }else{
    for(auto i : comp){
      double phi1 = i->phi() >= 0 ? i->phi(): 
        i->phi()+2*pi; 
      phi+= phi1;
    }
    
    double com = comp.front()->phi() >= 0 ? comp.front()->phi():
      2*pi + comp.front()->phi();
    double temp = fabs(phi/float(comp.size()) - com) > 2. ? 
      pi - phi/float(comp.size()):
      phi/float(comp.size());
    temp = temp >= 0? temp:2*pi+temp;
    return temp;
      }
    }