---

CSCTFPtMethods Class Reference

#include <L1Trigger/CSCTrackFinder/interface/CSCTFPtMethods.h>

List of all members.

Public Types
enum	{ kMaxParameters = 4 }
enum	{   kME1andME2 = 1, kME1andME3, kME2andME3, kME2andME4,   kME3andME4, kME1andME2ovr, kME2andMB1, kME2andMB2 }
	Allowed station combinations for extrapolation units. More...
Public Member Functions
bool	chargeValid (unsigned Pt, unsigned Quality, unsigned Eta, unsigned method) const
	The hybrid method may be changing soon to: 1st Calculate PT with Darin's method 2nd if BELOW a certain cut call Cathy's method 3rd if Cathy's < Darin's use Cathy's otherwise return Darin's A study needs to be performed to determine any gains from this procedure.
	CSCTFPtMethods (const L1MuTriggerPtScale *ptScale=0)
float	Pt2Stn (int type, float eta, float dphi, int fr=-1) const
	2-station Pt measurement for types (see SP class for 2-stn types)
float	Pt2StnChiSq (int type, float eta, int dphi, int fr) const
float	Pt2StnHybrid (int type, float eta, int dphi, int fr) const
	Third is the hybrid method.
float	Pt3Stn (int type, float eta, float dphi1, float dphi2, int fr=-1) const
	3-station Pt measurement for types (see SP class for 3-stn types)
float	Pt3StnChiSq (int type, float eta, int dphi1, int dphi2, int fr) const
float	Pt3StnHybrid (int type, float eta, int dphi1, int dphi2, int fr) const
float	PtEff90 (float pt, float eta, int mode) const
	Legacy Pt90 calculation function.
Static Public Attributes
static const float	AkHighEta_Fit1 [kME2andMB2][kMaxParameters]
static const float	AkHighEta_Fit2 [kME2andMB2][kMaxParameters]
static const float	AkLowEta_Fit1 [kME2andMB2][kMaxParameters]
	The three station pt measument only needs a one constant fit, but the dependence on eta is still there.
static const float	AkLowEta_Fit2 [kME2andMB2][kMaxParameters]
	First is the parameterizations of Acosta/McDonald.
static const float	BkHighEta_Fit2 [kME2andMB2][kMaxParameters]
static const float	BkLowEta_Fit2 [kME2andMB2][kMaxParameters]
static const float	dphifr0 [4][15][28]
static const float	dphifr1 [4][15][28]
static const float	etabins [16]
static const float	FRCorrHighEta [kME2andMB2][2]
static const float	FRCorrLowEta [kME2andMB2][2]
	Corrections for ME1 F/R bit.
static const float	kGlobalScaleFactor = 1.36
static const float	ptbins [29]
	Second are the parameterizations of Acosta/Yeh.
static const float	sigmafr0 [4][15][28]
static const float	sigmafr1 [4][15][28]
Private Attributes
const L1MuTriggerPtScale *	trigger_scale

---

Detailed Description

Definition at line 8 of file CSCTFPtMethods.h.

---

Member Enumeration Documentation

anonymous enum

Enumerator:

kMaxParameters

Definition at line 11 of file CSCTFPtMethods.h.

{kMaxParameters = 4};  

anonymous enum

Allowed station combinations for extrapolation units.

Enumerator:

kME1andME2
kME1andME3
kME2andME3
kME2andME4
kME3andME4
kME1andME2ovr
kME2andMB1
kME2andMB2

Definition at line 13 of file CSCTFPtMethods.h.

       {kME1andME2=1, kME1andME3, kME2andME3, kME2andME4, 
        kME3andME4, kME1andME2ovr, kME2andMB1, kME2andMB2};

---

Constructor & Destructor Documentation

CSCTFPtMethods::CSCTFPtMethods

(

const L1MuTriggerPtScale *

ptScale = 0

)

Definition at line 130 of file CSCTFPtMethods.cc.

  : trigger_scale( ptScale )
{
}

---

Member Function Documentation

bool CSCTFPtMethods::chargeValid	(	unsigned	Pt,
		unsigned	Quality,
		unsigned	Eta,
		unsigned	method
	)			const

The hybrid method may be changing soon to: 1st Calculate PT with Darin's method 2nd if BELOW a certain cut call Cathy's method 3rd if Cathy's < Darin's use Cathy's otherwise return Darin's A study needs to be performed to determine any gains from this procedure.

A method to calculate the charge valid bit Regions where this bit is were determined via simulation

Definition at line 1023 of file CSCTFPtMethods.cc.

References HLT_VtxMuL3::result.

{
  bool result = false;

  switch(method)
    {
    case 1:
      break;
    case 2:
      break;
    case 3:
      if(quality != 1)
        {
          if(pT <= 19)
            {
              if(eta >=4 && eta <=9)
                result = true;         
            }
          if(pT <= 18 && pT >= 10)
            {
              if(eta > 0 && eta < 4)
                result = true;
            }
        }
    };  
  
  return result;
}

float CSCTFPtMethods::Pt2Stn	(	int	type,
		float	eta,
		float	dphi,
		int	fr = -1
	)			const

2-station Pt measurement for types (see SP class for 2-stn types)

DEA: kluge to get verilog and original model to agree

Definition at line 139 of file CSCTFPtMethods.cc.

References funct::A, AkHighEta_Fit2, AkLowEta_Fit2, BkHighEta_Fit2, BkLowEta_Fit2, FRCorrHighEta, FRCorrLowEta, L1MuScale::getLowEdge(), L1MuTriggerPtScale::getPtScale(), kGlobalScaleFactor, kME1andME2, kME2andMB2, funct::sqrt(), and trigger_scale.

Referenced by CSCTFPtLUT::calcPt(), Pt2StnChiSq(), Pt2StnHybrid(), and Pt3Stn().

{
  float A = 0; 
  float B = 0;
  if (dphi == 0.0) dphi = 1.e-6;
  if (eta  < 0.0) eta  = static_cast<float>(fabs(static_cast<double>(eta)));
  if (dphi < 0.0) dphi = static_cast<float>(fabs(static_cast<double>(dphi)));
#ifdef L1CSC_STANDALONE
  if (type == kME1andME2 && eta <= 1.25) fr = 1; 
#else
  if (type == kME1andME2 && eta < 1.2) fr = 1; 
#endif
  if (type >= kME1andME2 && 
      type <= kME2andMB2 && eta < 2.5) 
    {
      if (eta >= 0.0 && eta < 1.6)
        {
          A = AkLowEta_Fit2[type-1][0] + AkLowEta_Fit2[type-1][1]*eta 
            + AkLowEta_Fit2[type-1][2]*eta*eta + AkLowEta_Fit2[type-1][3]*eta*eta*eta;
          B = BkLowEta_Fit2[type-1][0] + BkLowEta_Fit2[type-1][1]*eta 
            + BkLowEta_Fit2[type-1][2]*eta*eta + BkLowEta_Fit2[type-1][3]*eta*eta*eta;
          if (fr >= 0 && fr <= 1) 
            {     
              A *= FRCorrLowEta[type-1][fr];
              B *= FRCorrLowEta[type-1][fr];
            } 
        }       
      
      if (eta >= 1.6 && eta < 2.5)
        {
          A = AkHighEta_Fit2[type-1][0] + AkHighEta_Fit2[type-1][1]*eta 
            + AkHighEta_Fit2[type-1][2]*eta*eta + AkHighEta_Fit2[type-1][3]*eta*eta*eta;
          B = BkHighEta_Fit2[type-1][0] + BkHighEta_Fit2[type-1][1]*eta 
            + BkHighEta_Fit2[type-1][2]*eta*eta + BkHighEta_Fit2[type-1][3]*eta*eta*eta;
        if (fr >= 0 && fr <= 1) 
          {
            A *= FRCorrHighEta[type-1][fr];
            B *= FRCorrHighEta[type-1][fr];
          } 
        }
      
      A *= kGlobalScaleFactor;
      B *= kGlobalScaleFactor;
      float Pt = (A + sqrt(A*A + 4.*dphi*B))/(2.*dphi);

      //      return (Pt>0.0) ? Pt : 0.0;
      return (Pt>trigger_scale->getPtScale()->getLowEdge(1)) ? Pt 
        : trigger_scale->getPtScale()->getLowEdge(1);
    }
    return 0.0;
}

float CSCTFPtMethods::Pt2StnChiSq	(	int	type,
		float	eta,
		int	dphi,
		int	fr
	)			const

Definition at line 618 of file CSCTFPtMethods.cc.

References diff, dphifr0, dphifr1, etabins, cmsRelvalreport::exit, i, j, k, min, Pt2Stn(), ptbins, ptmin, CSCTFConstants::SECTOR_RAD, sigmafr0, and sigmafr1.

Referenced by CSCTFPtLUT::calcPt(), and Pt2StnHybrid().

{

  float diff, min, ptmin, ptmax;
  float mypt = 0.0;

  int dphicopy = dphi;

  if (type<0 || type>3) 
    {
//      std::cout << "CSCTFPtMethods: illegal track type for Chi-square method" << std::endl;
      edm::LogError("CSCTFPtMethods::Pt2StnChiSq()")<<"Illegal track type for Chi-square method";
      return 0.;
    }
  
  // flip sign
  //dphi = -dphi;
  
  // force positive
  eta = fabs(eta);
  
  //determine which eta bin muon falls into
  int i = 0;
  for(i = 0; i < 15; i++) 
    {
      if (eta >= etabins[i] && eta < etabins[i+1]) break;
    }
  if ( i == 15 ) 
    {
//      std::cout<<" CSCTFPtMethods: muon not within any eta range"<< std::endl;
      edm::LogWarning("CSCTFPtMethods::Pt2StnChiSq()")<<"Muon not within any eta range";
      if (eta<etabins[0]) 
        { 
          eta = etabins[0];
          i = 0;
        }
      else if (eta>=etabins[15]) 
        {
          eta = etabins[15];
          i = 15;
        }
      else
        exit(0);
    }

  bool bCallOldMethod = false;
  
  if (fr == 0) 
    {
      //advance past bins in array w/ default values
      //default is: dphifr0[x][y][z] = -1, sigmafr0[x][y][z] = 1
      int j = 0;  //start from 1st column, j=0
      while (dphifr0[type][i][j] == -1 && sigmafr0[type][i][j] == 1 && j != 28) j++;
      if ( j == 28 ) 
        {
//        std::cout <<" CSCTFPtMethods: every entry in row is default"<< std::endl;
          edm::LogInfo("CSCTFPtMethods::Pt2StnChiSq()")<<"Every entry in row is default";
          //      exit(0); //normal termination
          bCallOldMethod = true;
        }

      if (!bCallOldMethod) 
        {
          //set min, ptmin, ptmax to first bin
          //in row that is not default
          min = fabs(dphi - dphifr0[type][i][j]) / sigmafr0[type][i][j];
          ptmin = ptbins[j];
          ptmax = ptbins[j+1];
          
          //loop through all pt bins.  assign pt of bin w/ least diff
          //do not include default bins in calculation
          for ( int k = j; k < 28; k++ ) 
            {
              if(dphifr0[type][i][k] != -1 || sigmafr0[type][i][k] != 1) 
                {
                  diff = fabs(dphi - dphifr0[type][i][k]) / sigmafr0[type][i][k];
                  if (diff < min) 
                    {
                      min = diff;
                      ptmin = ptbins[k];
                      ptmax = ptbins[k+1];
                    }
                }
            }
          mypt = (ptmin + ptmax)/2;
        }
    }

  if (fr == 1) 
    {
      //advance past bins in array w/ default values
      //default is: dphifr1[x][y][z] = -1, sigmafr1[x][y][z] = 1
      int j = 0;  //start from 1st column, j=0
      
      while (dphifr1[type][i][j] == -1 && sigmafr1[type][i][j] == 1 && j != 28) j++;
      
      if ( j == 28 ) 
        {
//        std::cout <<" CSCTFPtMethods: every entry in row is default"<< std::endl;
          edm::LogInfo("CSCTFPtMethods::Pt2StnChiSq()")<<"Every entry in row is default";
          //      exit(0); //normal termination
          bCallOldMethod = true;
        }

    if (!bCallOldMethod) 
      {
        //set min, ptmin, ptmax to first bin
        //in row that is not default
        min = fabs(dphi - dphifr1[type][i][j]) / sigmafr1[type][i][j];
        ptmin = ptbins[j];
        ptmax = ptbins[j+1];
        
        //loop through all pt bins.  assign pt of bin w/ least diff
        //do not include default bins in calculation
        for ( int k = j; k < 28; k++ ) 
          {
            if(dphifr1[type][i][k] != -1 || sigmafr1[type][i][k] != 1) 
              {
                diff = fabs(dphi - dphifr1[type][i][k]) / sigmafr1[type][i][k];
                if (diff < min) 
                  {
                    min = diff;
                    ptmin = ptbins[k];
                    ptmax = ptbins[k+1];
                  }
              }
          }
        mypt = (ptmin + ptmax)/2;
      }
    }

  // hybrid approach:
  if (bCallOldMethod)
    {
      float dphiR = static_cast<float>(dphicopy) / static_cast<float>(1<<12) * CSCTFConstants::SECTOR_RAD;
      // must change type definition, just add one
      mypt =  Pt2Stn(type+1, eta, dphiR, fr);
    }
  
  return mypt;
}

float CSCTFPtMethods::Pt2StnHybrid	(	int	type,
		float	eta,
		int	dphi,
		int	fr
	)			const

Third is the hybrid method.

Definition at line 993 of file CSCTFPtMethods.cc.

References Pt2Stn(), Pt2StnChiSq(), and CSCTFConstants::SECTOR_RAD.

Referenced by CSCTFPtLUT::calcPt().

{
  float mypt = 0.0;

  mypt = Pt2StnChiSq(type, eta, dphi, fr);
  if(mypt >= 8.0)
    {
      float dphiR = static_cast<float>(dphi) / static_cast<float>(1<<12) * CSCTFConstants::SECTOR_RAD;
      mypt = Pt2Stn(type+1, eta, dphiR, fr);
    }

  return mypt;
}

float CSCTFPtMethods::Pt3Stn	(	int	type,
		float	eta,
		float	dphi1,
		float	dphi2,
		int	fr = -1
	)			const

3-station Pt measurement for types (see SP class for 3-stn types)

Definition at line 192 of file CSCTFPtMethods.cc.

References AkHighEta_Fit1, AkLowEta_Fit1, c1, c2, FRCorrHighEta, FRCorrLowEta, L1MuScale::getLowEdge(), L1MuTriggerPtScale::getPtScale(), kGlobalScaleFactor, kME1andME2, kME1andME3, kME2andME3, kME2andME4, kME3andME4, Pt2Stn(), r, funct::sqrt(), and trigger_scale.

Referenced by CSCTFPtLUT::calcPt(), Pt3StnChiSq(), and Pt3StnHybrid().

{
    int ty1 = 0, ty2 =0;

    switch (type)
      {
      case 1 :
        ty1 = kME1andME2 -1; // subtype sets the right position for array A
        ty2 = kME2andME3 -1;
        break;
      case 2 :
        ty1 = kME1andME2 -1;
        ty2 = kME2andME4 -1;
        break;
      case 3 :
        ty1 = kME1andME3 -1;
        ty2 = kME3andME4 -1;
        break;
      case 4 :
        ty1 = kME2andME3 -1;
        ty2 = kME3andME4 -1;
        break;
      default:
        return 0.0;
      }

    // Switch to 2-Station measurement if dphi is too small
    // box cut around Pt of 10 GeV
    float Pt;
//    if ( (fabs(static_cast<double>(dphi2))<0.004) &&
//         (fabs(static_cast<double>(dphi1))<0.030)) {
    if ( fabs(static_cast<double>(dphi2))<0.004 ) 
      {
        Pt = Pt2Stn((ty1+1), eta, dphi1, fr);
      }
    else 
      {
        float c1=(-.2999 * (eta*eta*eta) + 2.030 * (eta*eta) - 4.235 * eta + 3.024) + 0.02;
        float c2=(-2.484 * (eta*eta*eta) + 14.43 * (eta*eta) - 27.66 * eta + 18.47)*.72;
        float r = 0.6; //correlation
        float A1 = 0; 
        float A2 = 0;
        
        if (dphi1 == 0.0) dphi1 = 1.e-6;
        if (dphi2 == 0.0) dphi2 = 1.e-6;    
        if (eta  < 0.0) eta  = static_cast<float>(fabs(static_cast<double>(eta)));
        if (eta >= 0.0 && eta < 1.6) 
          {
            A1 = AkLowEta_Fit1[ty1][0]         + AkLowEta_Fit1[ty1][1]*eta 
              + AkLowEta_Fit1[ty1][2]*eta*eta + AkLowEta_Fit1[ty1][3]*eta*eta*eta;
            A2 = AkLowEta_Fit1[ty2][0]         + AkLowEta_Fit1[ty2][1]*eta 
              + AkLowEta_Fit1[ty2][2]*eta*eta + AkLowEta_Fit1[ty2][3]*eta*eta*eta;
            if (fr >= 0 && fr <= 1) 
              {
                A1 *= FRCorrLowEta[ty1][fr];
              }
          }
        
      if (eta >= 1.6 && eta < 2.5)
        {       
          A1 = AkHighEta_Fit1[ty1][0]         + AkHighEta_Fit1[ty1][1]*eta
            + AkHighEta_Fit1[ty1][2]*eta*eta + AkHighEta_Fit1[ty1][3]*eta*eta*eta;
          A2 = AkHighEta_Fit1[ty2][0]         + AkHighEta_Fit1[ty2][1]*eta
            + AkHighEta_Fit1[ty2][2]*eta*eta + AkHighEta_Fit1[ty2][3]*eta*eta*eta;
          if (fr >= 0 && fr <= 1) 
            {
              A1 *= FRCorrHighEta[ty1][fr];
            }
        }
      A1 *= kGlobalScaleFactor;
      A2 *= kGlobalScaleFactor;
      Pt = 1/((( -dphi1/c1/c1/A1+r*dphi2/c1/c2/A2+dphi1*r/c1/c2/A1-dphi2/c2/c2/A2)
               + sqrt( (dphi1/c1/c1/A1-r*dphi2/c1/c2/A2-dphi1*r/c1/A1/c2+dphi2/c2/c2/A2)
                       *(dphi1/c1/c1/A1-r*dphi2/c1/c2/A2-dphi1*r/c1/A1/c2+dphi2/c2/c2/A2)
                       + 8*(1-r*r)*(dphi1*dphi1/c1/c1/A1/A1-2*dphi1*dphi2*r/c1/A1/c2/A2
                                    +dphi2*dphi2/c2/c2/A2/A2)))/(4*(1-r*r)));
      
      } // end 2 or 3 station method
    //    if (Pt<0.0) Pt = 0.0;
    if (Pt<trigger_scale->getPtScale()->getLowEdge(1)) 
      Pt = trigger_scale->getPtScale()->getLowEdge(1);

    // switch to 2-stn Pt above 10 GeV:
    /*
      if (Pt>10.0) {
      switch (type)
      {
      case 1 :
      case 2 :
      Pt = Pt2Stn(1, eta, dphi1);
      break;
      case 3:
      Pt = Pt2Stn(2, eta, dphi1);
      break;
      case 4:
      Pt = Pt2Stn(3, eta, dphi1);
      }
      }
    */
    
    return Pt;
}

float CSCTFPtMethods::Pt3StnChiSq	(	int	type,
		float	eta,
		int	dphi1,
		int	dphi2,
		int	fr
	)			const

Definition at line 761 of file CSCTFPtMethods.cc.

References funct::abs(), diff, dphifr0, dphifr1, etabins, cmsRelvalreport::exit, i, j, k, min, funct::pow(), Pt3Stn(), ptbins, ptmin, CSCTFConstants::SECTOR_RAD, sigmafr0, and sigmafr1.

Referenced by CSCTFPtLUT::calcPt(), and Pt3StnHybrid().

{

  float diff, min, ptmin, ptmax;
  float mypt = 0.0;

  int dphi1copy = dphi1, dphi2copy = dphi2;

  if (type<4 || type>5) 
    {
//      std::cout << "PtParams: illegal track type for Chi-square method" << std::endl;
      edm::LogError("CSCTFPtMethods::Pt3StnChiSq()")<<"Illegal track type for Chi-square method";
      return 0.;
    }

  // flip sign
  //dphi1 = -dphi1; 
  //dphi2 = -dphi2; 


  //determine which eta bin muon falls into
  int i = 0;
  for( i = 0; i < 15; i++ ) 
    {
      if ( eta >= etabins[i] && eta < etabins[i + 1] ) break;
    }
  if ( i == 15 ) 
    {
      //std::cout<<"muon not within any eta range";
      edm::LogWarning("CSCTFPtMethods::Pt3StnChiSq()")<<"Muon not within any eta range";
      if (eta<etabins[0]) 
        { 
          eta = etabins[0];
          i = 0;
        }
      else if (eta>=etabins[15]) 
        {
          eta = etabins[15];
          i = 15;
        }
      else
        exit(0);
    }

  bool bCallOldMethod = false;

  if(abs(dphi2) < 32) dphi2 = 0; //dphi2=dphi23 or dphi24

  //sta1-2-3, dphi1=dphi12, dphi2=dphi23
  if (type == 4) 
    {
      if (fr == 0) 
        {
          //advance past bins in array w/ default values
          //default is: dphifr0[x][y][z] = -1, sigmafr0[x][y][z] = 1
          int j = 0;  //start from 1st column, j = 0
          while ( ((dphifr0[0][i][j] == -1 && sigmafr0[0][i][j] == 1) || (dphifr0[2][i][j] == -1 && sigmafr0[2][i][j] == 1)) && j != 28 ) j++;
          if ( j == 28 ) 
            {
//            std::cout<<" L1MuCSCPtParams: every entry in row is default"<<std::endl;  
              edm::LogInfo("CSCTFPtMethods::Pt3StnChiSq()")<<"Every entry in row is default";
              //        exit(0); //normal termination
              bCallOldMethod = true;
            }
          
          if (!bCallOldMethod) 
            {
              //set min, ptmin, ptmax to first bin
              //in row that is not default
              min = 1.5625 * ( (pow((dphi1 - dphifr0[0][i][j]),2) / pow(sigmafr0[0][i][j],2)) + (pow((dphi2 - dphifr0[2][i][j]),2) / pow(sigmafr0[2][i][j],2)) - (1.2 * (dphi1 - dphifr0[0][i][j]) * (dphi2 - dphifr0[2][i][j]) / (sigmafr0[0][i][j] * sigmafr0[2][i][j])) ); //calculate chi square
              ptmin = ptbins[j];
              ptmax = ptbins[j + 1];
              
              //loop through all pt bins.  assign pt of bin w/ least diff
              //do not include default bins in calculation
              for ( int k = j; k < 28; k++ ) 
                {
                  if( (dphifr0[0][i][k] != -1 || sigmafr0[0][i][k] != 1) && (dphifr0[2][i][k] != -1 || sigmafr0[2][i][k] != 1) ) 
                    {
                      diff = 1.5625 * ( (pow((dphi1 - dphifr0[0][i][k]),2) / pow(sigmafr0[0][i][k],2)) + (pow((dphi2 - dphifr0[2][i][k]),2) / pow(sigmafr0[2][i][k],2)) - (1.2 * (dphi1 - dphifr0[0][i][k]) * (dphi2 - dphifr0[2][i][k]) / (sigmafr0[0][i][k] * sigmafr0[2][i][k])) );
                      if (diff < min) 
                        {
                          min = diff;
                          ptmin = ptbins[k];
                          ptmax = ptbins[k + 1];
                        }
                    }
                }
              mypt = (ptmin + ptmax) / 2;
            }
        }
      
      if (fr == 1) 
        {
          //advance past bins in array w/ default values
          //default is: dphifr1[x][y][z] = -1, sigmafr1[x][y][z] = 1
          int j = 0;  //start from 1st column, j = 0
          while ( ((dphifr1[0][i][j] == -1 && sigmafr1[0][i][j] == 1) || (dphifr1[2][i][j] == -1 && sigmafr1[2][i][j] == 1)) && j != 28 ) j++;
          if ( j == 28 ) 
            {
//            std::cout<<" L1MuCSCPtParams: every entry in row is default"<<std::endl;  
              edm::LogInfo("CSCTFPtMethods::Pt3StnChiSq()")<<"Every entry in row is default";
              //        exit(0); //normal termination
              bCallOldMethod = true;
            }

          if (!bCallOldMethod) 
            {
              //set min, ptmin, ptmax to first bin
              //in row that is not default
              min = 1.5625 * ( (pow((dphi1 - dphifr1[0][i][j]),2) / pow(sigmafr1[0][i][j],2)) + (pow((dphi2 - dphifr1[2][i][j]),2) / pow(sigmafr1[2][i][j],2)) - (1.2 * (dphi1 - dphifr1[0][i][j]) * (dphi2 - dphifr1[2][i][j]) / (sigmafr1[0][i][j] * sigmafr1[2][i][j])) ); //calculate chi square
              ptmin = ptbins[j];
              ptmax = ptbins[j + 1];
              
              //loop through all pt bins.  assign pt of bin w/ least diff
              //do not include default bins in calculation
              for ( int k = j; k  < 28; k++ ) 
                {
                  if( (dphifr1[0][i][k] != -1 || sigmafr1[0][i][k] != 1) && (dphifr1[2][i][k] != -1 || sigmafr1[2][i][k] != 1) ) 
                    {
                      diff = 1.5625 * ( (pow((dphi1 - dphifr1[0][i][k]),2) / pow(sigmafr1[0][i][k],2)) + (pow((dphi2 - dphifr1[2][i][k]),2) / pow(sigmafr1[2][i][k],2)) - (1.2 * (dphi1 - dphifr1[0][i][k]) * (dphi2 - dphifr1[2][i][k]) / (sigmafr1[0][i][k] * sigmafr1[2][i][k])) );
                      if (diff < min) 
                        {
                          min = diff;
                          ptmin = ptbins[k];
                          ptmax = ptbins[k + 1];
                        }
                    }
                  mypt = (ptmin + ptmax) / 2;
                }
            }
        }
    }
  
  //sta1-2-4, dphi1=dphi12, dphi2=dphi24
  if (type == 5) 
    {
      
      if (fr == 0) 
        {
          //advance past bins in array w/ default values
          //default is: dphifr0[x][y][z] = -1, sigmafr0[x][y][z] = 1
          int j = 0;  //start from 1st column, j = 0
          while ( ((dphifr0[0][i][j] == -1 && sigmafr0[0][i][j] == 1) || (dphifr0[3][i][j] == -1 && sigmafr0[3][i][j] == 1)) && j != 28 ) j++;
          if ( j == 28 ) 
            {
//            std::cout<<" L1MuCSCPtParams: every entry in row is default"<<std::endl;  
              edm::LogInfo("CSCTFPtMethods::Pt3StnChiSq()")<<"Every entry in row is default";
              //        exit(0); //normal termination
              bCallOldMethod = true;
            }

          if (!bCallOldMethod) 
            {
              //set min, ptmin, ptmax to first bin
              //in row that is not default
              min = 1.5625 * ( (pow((dphi1 - dphifr0[0][i][j]),2) / pow(sigmafr0[0][i][j],2)) + (pow((dphi2 - dphifr0[3][i][j]),2) / pow(sigmafr0[3][i][j],2)) - (1.2 * (dphi1 - dphifr0[0][i][j]) * (dphi2 - dphifr0[3][i][j]) / (sigmafr0[0][i][j] * sigmafr0[3][i][j])) ); //calculate chi square
              ptmin = ptbins[j];
              ptmax = ptbins[j + 1];
              
              //loop through all pt bins.  assign pt of bin w/ least diff
      //do not include default bins in calculation
              for ( int k = j; k < 28; k++ ) 
                {
                  if( (dphifr0[0][i][k] != -1 || sigmafr0[0][i][k] != 1) && (dphifr0[3][i][k] != -1 || sigmafr0[3][i][k] != 1) ) 
                    {
                      diff = 1.5625 * ( (pow((dphi1 - dphifr0[0][i][k]),2) / pow(sigmafr0[0][i][k],2)) + (pow((dphi2 - dphifr0[3][i][k]),2) / pow(sigmafr0[3][i][k],2)) - (1.2 * (dphi1 - dphifr0[0][i][k]) * (dphi2 - dphifr0[3][i][k]) / (sigmafr0[0][i][k] * sigmafr0[3][i][k])) );
                      if (diff < min) 
                        {
                          min = diff;
                          ptmin = ptbins[k];
                          ptmax = ptbins[k + 1];
                        }
                    }
                }
              mypt = (ptmin + ptmax) / 2;
            }
        }
      
      if (fr == 1) 
        {
          //advance past bins in array w/ default values
          //default is: dphifr1[x][y][z] = -1, sigmafr1[x][y][z] = 1
          int j = 0;  //start from 1st column, j = 0
          while ( ((dphifr1[0][i][j] == -1 && sigmafr1[0][i][j] == 1) || (dphifr1[3][i][j] == -1 && sigmafr1[3][i][j] == 1)) && j != 28 ) j++;
          if ( j == 28 ) 
            {
//            std::cout<<" L1MuCSCPtParams: every entry in row is default"<<std::endl;  
              edm::LogInfo("CSCTFPtMethods::Pt3StnChiSq()")<<"Every entry in row is default";
              //        exit(0); //normal termination
              bCallOldMethod = true;
            }
          
          if (!bCallOldMethod) 
            {
              //set min, ptmin, ptmax to first bin
              //in row that is not default
              min = 1.5625 * ( (pow((dphi1 - dphifr1[0][i][j]),2) / pow(sigmafr1[0][i][j],2)) + (pow((dphi2 - dphifr1[3][i][j]),2) / pow(sigmafr1[3][i][j],2)) - (1.2 * (dphi1 - dphifr1[0][i][j]) * (dphi2 - dphifr1[3][i][j]) / (sigmafr1[0][i][j] * sigmafr1[3][i][j])) ); //calculate chi square
              ptmin = ptbins[j];
              ptmax = ptbins[j + 1];
              
              //loop through all pt bins.  assign pt of bin w/ least diff
              //do not include default bins in calculation
              for ( int k = j; k < 28; k++ ) 
                {
                  if( (dphifr1[0][i][k] != -1 || sigmafr1[0][i][k] != 1) && (dphifr1[3][i][k] != -1 || sigmafr1[3][i][k] != 1) ) 
                    {
                      diff = 1.5625 * ( (pow((dphi1 - dphifr1[0][i][k]),2) / pow(sigmafr1[0][i][k],2)) + (pow((dphi2 - dphifr1[3][i][k]),2) / pow(sigmafr1[3][i][k],2)) - (1.2 * (dphi1 - dphifr1[0][i][k]) * (dphi2 - dphifr1[3][i][k]) / (sigmafr1[0][i][k] * sigmafr1[3][i][k])) );
                      if (diff < min) 
                        {
                          min = diff;
                          ptmin = ptbins[k];
                          ptmax = ptbins[k + 1];
                        }
                    }
                }
              mypt = (ptmin + ptmax) / 2;
            }
        }
    }
    
  if (bCallOldMethod) 
    {
      float dphi12R = (static_cast<float>(dphi1copy)) / static_cast<float>(1<<12) * CSCTFConstants::SECTOR_RAD;
      float dphi23R = (static_cast<float>(dphi2copy)) / static_cast<float>(1<<12) * CSCTFConstants::SECTOR_RAD;
      // change defintion of track type to old method
      mypt =  Pt3Stn(type-3, eta, dphi12R, dphi23R, fr);
    }
  return mypt;

}

float CSCTFPtMethods::Pt3StnHybrid	(	int	type,
		float	eta,
		int	dphi1,
		int	dphi2,
		int	fr
	)			const

Definition at line 1007 of file CSCTFPtMethods.cc.

References Pt3Stn(), Pt3StnChiSq(), and CSCTFConstants::SECTOR_RAD.

Referenced by CSCTFPtLUT::calcPt().

{
  float mypt = 0.0;

  mypt = Pt3StnChiSq(type, eta, dphi1, dphi2, fr);
  if(mypt >= 8.0)
    {
      float dphi12R = (static_cast<float>(dphi1)) / static_cast<float>(1<<12) * CSCTFConstants::SECTOR_RAD;
      float dphi23R = (static_cast<float>(dphi2)) / static_cast<float>(1<<12) * CSCTFConstants::SECTOR_RAD;
      // change defintion of track type to old method
      mypt =  Pt3Stn(type-3, eta, dphi12R, dphi23R, fr);      
    }

  return mypt;
}

float CSCTFPtMethods::PtEff90	(	float	pt,
		float	eta,
		int	mode
	)			const

Legacy Pt90 calculation function.

Definition at line 1052 of file CSCTFPtMethods.cc.

References c.

{
  // set pT resolution 
  float c;
  switch (mode) 
    {
      // 3-stn with ME1
    case 2:
    case 3:
    case 4:
      // leave as 2 Stn for now:
      //        if (eta < 2.0)
      //            c = 0.2;
      //        else
      //            c = 0.2 + (eta-2.0)*0.5;
      //        break;
      // 2-stn with ME1
    case 6:
    case 7:
      //        c=(-.2999 * (eta*eta*eta) + 2.030 * (eta*eta) - 4.235 * eta + 3.024)+0.02;
      // Try just fixed value from fits to Pt turn-on curves
      c=0.3;
      break;
      // 2 or 3 stn without ME1
    case 5:
    case 8:
    case 9:
      c=(-2.484 * (eta*eta*eta) + 14.43 * (eta*eta) - 27.66 * eta + 18.47)*.72;
      break;
    case 10:
      c=2.0;
      break;
      // MB1 tracks
    case 11:
    case 12:
    case 13:
    case 14:
    case 15:
      c=0.3;
      break;
    default:
      c = 0.0;
    };
  // Calculation of 90% threshold from 50% one, assuming Gaussian resolution
  // 90% of Gaussian area extends to 1.28*sigma
  if (c>0.78) c = 0.78;
  //    return pt/(1.0 - 1.28*c);
  // Let's go back to old empirical method:
  return pt*(1.0 + 1.28*c);
}

---

Member Data Documentation

const float CSCTFPtMethods::AkHighEta_Fit1 [static]

Initial value:

 
  {
    
    {0.6275, -0.2133, 0.0, 0.0},                        
    {-1.124, 2.492, -1.347, 0.2204},            
    {0.5214, -0.7536, 0.3913, -.06739},         
    {1.394, -1.886, 0.8809, -0.1369},           
    {-0.5519E-02, 0.8496E-02, 0.0, 0.0},        
    {0.6644, -0.2312, 0.0, 0.0},                        
    {0.80, -0.54, 0.0, 0.0},                    
    {0.80, -0.54, 0.0, 0.0},                    
  }

Definition at line 31 of file CSCTFPtMethods.h.

Referenced by Pt3Stn().

const float CSCTFPtMethods::AkHighEta_Fit2 [static]

Initial value:

 
  {
    {0.6425, -0.2308, 0.0, 0.0},                                
    {0.6923, -0.2336, 0.0, 0.0},                                
    {-0.7147, .01117, -0.5255E-04, 0.7991E-07}, 
    {-1.377, 2.150, -1.046, 0.1667},                    
    {-.02974, .04816, -.02105, .003640},                
    {0.6425, -0.2308, 0.0, 0.0},                                
    {0.80, -0.54, 0.0, 0.0},                            
    {0.80, -0.54, 0.0, 0.0},                            
  }

Definition at line 23 of file CSCTFPtMethods.h.

Referenced by Pt2Stn().

const float CSCTFPtMethods::AkLowEta_Fit1 [static]

Initial value:

 
  {
    
    
    {3.900, -8.391, 6.062, -1.436},                             
    
    {8.03602, -17.5355, 12.7725, -3.0616 },     
    {0.8372, -2.220, 1.908, -0.5233},                   
    {-0.04389, 0.05174, 0.0, 0.0},                              
    {-0.5519E-02, 0.8496E-02, 0.0, 0.0},                
    {-4.228, 8.546, -5.617, 1.235},                             
    {0.80, -0.54, 0.0, 0.0},                            
    {0.80, -0.54, 0.0, 0.0},                            
  }

The three station pt measument only needs a one constant fit, but the dependence on eta is still there.

dphi = A/pt

Definition at line 30 of file CSCTFPtMethods.h.

Referenced by Pt3Stn().

const float CSCTFPtMethods::AkLowEta_Fit2 [static]

Initial value:

 
  {
    
    {-0.6131, 0.6870, 0.0285, -0.1147},                 
    
    {8.03602, -17.5355, 12.7725, -3.0616 },     
    {5.765, -12.85, 9.543, -2.345},                     
    {-4.068, 8.254, -5.491, 1.214},                     
    {-.02974, .04816, -.02105, .003640},        
    {-3.358, 7.038, -4.766, 1.079},                     
    {0.80, -0.54, 0.0, 0.0},                    
    {0.80, -0.54, 0.0, 0.0},                    
  }

First is the parameterizations of Acosta/McDonald.

The two station pt measument needs a two constant fit with a break in the detector depending on what detectors it hit (eta) dphi = A/pt + B/(pt^2)

Definition at line 22 of file CSCTFPtMethods.h.

Referenced by Pt2Stn().

const float CSCTFPtMethods::BkHighEta_Fit2 [static]

Initial value:

 
  {
    {14.79, -21.10, 10.16, -1.623},             
    {70.67, -105.7, 52.49, -8.631},             
    {15.70, -23.94, 11.89, -1.939},             
    {17.18, -26.38, 13.13, -2.139},             
    {0.0, 0.0, 0.0, 0.0},                               
    {14.79, -21.10, 10.16, -1.623},             
    {0.0, 0.0, 0.0, 0.0},                               
    {0.0, 0.0, 0.0, 0.0}                                
  }

Definition at line 26 of file CSCTFPtMethods.h.

Referenced by Pt2Stn().

const float CSCTFPtMethods::BkLowEta_Fit2 [static]

Initial value:

 
  {
    
    {7.0509, -7.3282, -0.1465, 1.308},  
    
    {0.0, 0.0, 0.0, 0.0},       
    {-38.55, 82.95, -59.66, 14.29},             
    {20.71, -45.34, 32.05, -7.415},             
    {0.0, 0.0, 0.0, 0.0},                               
    {-0.6167, 0.08421, 0.2454, 0.0},    
    {0.0, 0.0, 0.0, 0.0},                               
    {0.0, 0.0, 0.0, 0.0}                                
  }

Definition at line 25 of file CSCTFPtMethods.h.

Referenced by Pt2Stn().

const float CSCTFPtMethods::dphifr0 [static]

Definition at line 48 of file CSCTFPtMethods.h.

Referenced by Pt2StnChiSq(), and Pt3StnChiSq().

const float CSCTFPtMethods::dphifr1 [static]

Definition at line 49 of file CSCTFPtMethods.h.

Referenced by Pt2StnChiSq(), and Pt3StnChiSq().

const float CSCTFPtMethods::etabins [static]

Initial value:

 
  {0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 
   1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4}

Definition at line 47 of file CSCTFPtMethods.h.

Referenced by Pt2StnChiSq(), and Pt3StnChiSq().

const float CSCTFPtMethods::FRCorrHighEta [static]

Initial value:

 
  {


    {0.966, 1.176}, 
    {0.966, 1.176},   
    {1.0, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0}    
  }

Definition at line 36 of file CSCTFPtMethods.h.

Referenced by Pt2Stn(), and Pt3Stn().

const float CSCTFPtMethods::FRCorrLowEta [static]

Initial value:

 
  {


    {1.30, 1.0}, 
    {1.30, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0},   
    {1.0, 1.0}    
  }

Corrections for ME1 F/R bit.

Definition at line 35 of file CSCTFPtMethods.h.

Referenced by Pt2Stn(), and Pt3Stn().

const float CSCTFPtMethods::kGlobalScaleFactor = 1.36 [static]

Definition at line 33 of file CSCTFPtMethods.h.

Referenced by Pt2Stn(), and Pt3Stn().

const float CSCTFPtMethods::ptbins [static]

Initial value:

 
  {2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0, 7.0, 8.0, 
   10.0, 12.0, 14.0, 16.0, 18.0, 20.0, 25.0, 30.0, 35.0, 40.0, 
   45.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 120.0, 140.0}

Second are the parameterizations of Acosta/Yeh.

Definition at line 46 of file CSCTFPtMethods.h.

Referenced by Pt2StnChiSq(), and Pt3StnChiSq().

const float CSCTFPtMethods::sigmafr0 [static]

Definition at line 50 of file CSCTFPtMethods.h.

Referenced by Pt2StnChiSq(), and Pt3StnChiSq().

const float CSCTFPtMethods::sigmafr1 [static]

Definition at line 51 of file CSCTFPtMethods.h.

Referenced by Pt2StnChiSq(), and Pt3StnChiSq().

const L1MuTriggerPtScale* CSCTFPtMethods::trigger_scale [private]

Definition at line 77 of file CSCTFPtMethods.h.

Referenced by Pt2Stn(), and Pt3Stn().

---

The documentation for this class was generated from the following files:

• L1Trigger/CSCTrackFinder/interface/CSCTFPtMethods.h
• L1Trigger/CSCTrackFinder/src/CSCTFPtMethods.cc

---