Conv Class Reference

#include <CSCxTalk.h>

Public Member Functions
	Conv ()
void	convolution (float *xleft_a, float *xleft_b, float *min_left, float *xright_a, float *xright_b, float *min_right, float *pTime)
float	elec (float t, float vs)
void	mkbins (float vs)
	~Conv ()

Public Attributes
float	conv [3][120]
float	convd [3][120]
float	conve [120]
float	nconvd [3][120]

Detailed Description

Definition at line 5 of file CSCxTalk.h.

Constructor & Destructor Documentation

Conv::Conv ( )

inline

Definition at line 9 of file CSCxTalk.h.

{}

Conv::~Conv ( )

inline

Definition at line 149 of file CSCxTalk.h.

{}

Member Function Documentation

void Conv::convolution ( float *  xleft_a, float *  xleft_b, float *  min_left, float *  xright_a, float *  xright_b, float *  min_right, float *  pTime  )

inline

Definition at line 39 of file CSCxTalk.h.

References conv, convd, conve, cross2(), i, j, relval_steps::k, prof2calltree::l, bookConverter::max, and makeMuonMisalignmentScenario::sum_x.

                                                                                                                                     {
    //void(convolution){  
    
    float max, cross0,cross2,min_l,min_r,sum_x=0.0,sumx2=0.;
    float sum_y_left=0.0,sum_y_right=0.0,sum_xy_left=0.0,sum_xy_right=0.0;
    float a_left=0.0,a_right=0.0,b_left=0.0,b_right=0.0,chi2_left=0.0,chi2_right=0.0,chi_left=0.0,chi_right=0.0;
    float aleft=0.0,aright=0.0,bleft=0.0,bright=0.0;
    int i,j,k,l,imax=0;
    
    for(l=0;l<3;l++){
      for(i=0;i<119;i++)conv[l][i]=0.0;
      for(j=0;j<119;j++){
    for(k=0;k<119;k++){
      if(j+k<119)conv[l][j+k]=conv[l][j+k]+convd[l][j]*conve[k];
    }
      }
    }
    max=0;
    min_l=9999999999999999.0;
    min_r=9999999999999999.0;
    for(i=0;i<119;i++){
      if(conv[1][i]>max){ 
    max=conv[1][i];
    imax=i;
      }
    }
    
    //find the max peak time from 3 timebins when line intersects x axis a+b*x=0 -> x=-a/b 
    float time1=-999.0, time2=-999.0, time3=-999.0;
    float data1=-999.0, data2=-999.0, data3=-999.0;
    float peakTime=0.0;
    
    time1=imax-1;
    time2=imax;
    time3=imax+1;
    
    data1=conv[1][imax-1];
    data2=conv[1][imax];
    data3=conv[1][imax+1];
    
    peakTime=(0.5)*((time1*time1*(data3-data2)+time2*time2*(data1-data3)+time3*time3*(data2-data1))/(time1*(data3-data2)+time2*(data1-data3)+time3*(data2-data1)))*6.25;
    
    for(l=0;l<3;l++){
      for(i=0;i<119;i++)conv[l][i]=conv[l][i]/max;
    }
    
    int nobs = 0;
    for (int j=0; j<119; j++){
      if (conv[1][j]>0.6) nobs++;
    }
    
    for(i=0;i<119;i++){
      cross0=0.0;
      cross2=0.0;
      
      if(conv[1][i]>0.6){
    cross0=conv[0][i]/(conv[0][i]+conv[1][i]+conv[2][i]);
    cross2=conv[2][i]/(conv[0][i]+conv[1][i]+conv[2][i]);
    
    sum_x += i;
    sum_y_left += cross0;
    sum_y_right += cross2;
    sumx2 += i*i;
    sum_xy_left += i*cross0;
    sum_xy_right += i*cross2;
      }
    }  
    
    //LMS fitting straight line y=a+b*x
    
    bleft  = ((nobs*sum_xy_left) - (sum_x * sum_y_left))/((nobs*sumx2) - (sum_x*sum_x));
    bright = ((nobs*sum_xy_right) - (sum_x * sum_y_right))/((nobs*sumx2) - (sum_x*sum_x));
    
    aleft  = ((sum_y_left*sumx2)-(sum_x*sum_xy_left))/((nobs*sumx2)-(sum_x*sum_x));
    aright = ((sum_y_right*sumx2)-(sum_x*sum_xy_right))/((nobs*sumx2)-(sum_x*sum_x));
    
    for(i=0;i<119;i++ ){
      chi2_left  += (cross0 -(aleft+(bleft*i)))*(cross0 -(aleft+(bleft*i)));
      chi2_right += (cross2 -(aright+(bright*i)))*(cross2 -(aright+(bright*i)));
    }   
    
    if(chi_left<min_l){
      min_l=chi_left;
      bleft=bleft;
      aleft=aleft;
    }
    if(chi_right<min_r){
      min_r=chi_right;
      bright=bright;
      aright=aright;
    }
    
    
    //Now calculating parameters in ns to compensate for drift in peak time  
    b_left  = bleft/6.25;
    b_right = bright/6.25;
    
    a_left  = aleft  + (bleft*peakTime/6.25);
    a_right = aright + (bright*peakTime/6.25);
    
    *xleft_a   = a_left; 
    *xleft_b   = b_left;
    *min_left  = chi2_left;
    *xright_a  = a_right;
    *xright_b  = b_right;
    *min_right = chi2_right;
    *pTime     = peakTime;
    
  } //end CONVOLUTION  

float Conv::elec ( float  t, float  vs  )

inline

Definition at line 12 of file CSCxTalk.h.

References f, and create_public_lumi_plots::log.

Referenced by mkbins().

                              {
    float f;
    if (t<=vs){
      f=log(t+2.5)-log(2.5);
    }
    else{
      f=log(t+2.5)-log(t-50+2.5);
    }
    return f;
  }//end elec

void Conv::mkbins ( float  vs )

inline

Definition at line 25 of file CSCxTalk.h.

References conve, elec(), i, relval_steps::k, and tree::t.

                       {
    int i,k;
    float t;
    for(i=0;i<120;i++) conve[i]=0.0;
    for(i=0;i<120;i++){
      for(k=0;k<16;k++){
    t=(6.25*i)+(k*0.625);
    conve[i]=conve[i]+elec(t,vs);
      }
    }
  } //end mkbins

Member Data Documentation

float Conv::conv[3][120]

Definition at line 155 of file CSCxTalk.h.

Referenced by convolution().

float Conv::convd[3][120]

Definition at line 152 of file CSCxTalk.h.

Referenced by convolution().

float Conv::conve[120]

Definition at line 154 of file CSCxTalk.h.

Referenced by convolution(), and mkbins().

float Conv::nconvd[3][120]

Definition at line 153 of file CSCxTalk.h.