DTMuonSLToSL Class Reference

#include <DTMuonSLToSL.h>

Inheritance diagram for DTMuonSLToSL:

Public Member Functions
void	calculationSLToSL ()
	DTMuonSLToSL (std::string, int, float, float, TFile *)
TMatrixD	returnbSLMatrix (float, float, float)
TMatrixD	returnCSLMatrix (float, float, float)
void	setBranchTree ()
	~DTMuonSLToSL ()
Public Member Functions inherited from DTMuonLocalAlignment
	DTMuonLocalAlignment ()
void	initNTuples (int)
void	setBranchAddressTree ()
	~DTMuonLocalAlignment ()

Private Attributes
float	cov [3][3]
float	dx
float	dz
float	phiy
float	ptMax
float	ptMin
int	srC
int	stC
TTree *	ttreeOutput
int	whC

Additional Inherited Members
Public Attributes inherited from DTMuonLocalAlignment
float	charge
float	dxdzSl [5]
float	dxdzSlSL1 [5]
float	dxdzSlSL3 [5]
float	dydzSl [5]
float	edxdzSl [5]
float	edxdzSlSL1 [5]
float	edxdzSlSL3 [5]
float	edydzSl [5]
float	eta
float	ex [5][14]
float	excp [5][14]
float	exdxdzSl [5]
float	exdxdzSlSL1 [5]
float	exdxdzSlSL3 [5]
float	exSl [5]
float	exSlSL1 [5]
float	exSlSL3 [5]
float	eycp [5][14]
float	eydydzSl [5]
float	eySl [5]
TFile *	f
int	la [5][14]
int	nhits [5]
int	nphihits [5]
int	nseg
int	nthetahits [5]
std::string	ntuplePath
int	numberOfRootFiles
float	p
float	phi
float	pt
int	sl [5][14]
int	sr [5]
int	st [5]
TChain *	tali
int	wh [5]
float	xc [5][14]
float	xcp [5][14]
float	xSl [5]
float	xSL1SL3 [5]
float	xSL3SL1 [5]
float	xSlSL1 [5]
float	xSlSL3 [5]
float	yc [5][14]
float	ycp [5][14]
float	ySl [5]
float	zc [5][14]

Detailed Description

$Date$

Revision

1.3

Author

Luca Scodellaro Luca..nosp@m.Scod.nosp@m.ellar.nosp@m.o@ce.nosp@m.rn.ch

Definition at line 18 of file DTMuonSLToSL.h.

Constructor & Destructor Documentation

DTMuonSLToSL()

DTMuonSLToSL::DTMuonSLToSL	(	std::string	path,
		int	n_files,
		float	MaxPt,
		float	MinPt,
		TFile *	f_
	)

Definition at line 3 of file DTMuonSLToSL.cc.

References calculationSLToSL(), DTMuonLocalAlignment::f, DTMuonLocalAlignment::initNTuples(), genfragment_ptgun_cfg::MaxPt, electronAnalyzer_cfi::MinPt, DTMuonLocalAlignment::ntuplePath, DTMuonLocalAlignment::numberOfRootFiles, EnsembleCalibrationLA_cfg::path, ptMax, ptMin, and setBranchTree().

                                                                                           {
  ntuplePath = path;
  numberOfRootFiles = n_files;

  f = f_;

  ptMax = MaxPt;
  ptMin = MinPt;

  setBranchTree();

  initNTuples(0);

  calculationSLToSL();
}

~DTMuonSLToSL()

DTMuonSLToSL::~DTMuonSLToSL

(

)

Definition at line 19 of file DTMuonSLToSL.cc.

{}

Member Function Documentation

calculationSLToSL()

void DTMuonSLToSL::calculationSLToSL

(

)

Definition at line 21 of file DTMuonSLToSL.cc.

References DummyCfis::c, counter, cov, dx, DTMuonLocalAlignment::dxdzSlSL1, DTMuonLocalAlignment::dxdzSlSL3, dz, DTMuonLocalAlignment::f, mps_fire::i, DTMuonLocalAlignment::nhits, DTMuonLocalAlignment::nseg, phiy, DTMuonLocalAlignment::pt, ptMax, ptMin, returnbSLMatrix(), returnCSLMatrix(), alignCSCRings::s, nano_mu_digi_cff::sector, DTMuonLocalAlignment::sr, srC, DTMuonLocalAlignment::st, relativeConstraints::station, stC, DTMuonLocalAlignment::tali, ttreeOutput, DTMuonLocalAlignment::wh, whC, makeMuonMisalignmentScenario::wheel, DTMuonLocalAlignment::xSL1SL3, DTMuonLocalAlignment::xSL3SL1, DTMuonLocalAlignment::xSlSL1, DTMuonLocalAlignment::xSlSL3, and DTMuonLocalAlignment::zc.

Referenced by DTMuonSLToSL().

                                     {
  TMatrixD ****C13 = new TMatrixD ***[5];
  TMatrixD ****b13 = new TMatrixD ***[5];
  TMatrixD ****C31 = new TMatrixD ***[5];
  TMatrixD ****b31 = new TMatrixD ***[5];

  for (int whI = -2; whI < 3; ++whI) {
    C13[whI + 2] = new TMatrixD **[4];
    b13[whI + 2] = new TMatrixD **[4];
    C31[whI + 2] = new TMatrixD **[4];
    b31[whI + 2] = new TMatrixD **[4];
    for (int stI = 1; stI < 5; ++stI) {
      C13[whI + 2][stI - 1] = new TMatrixD *[14];
      b13[whI + 2][stI - 1] = new TMatrixD *[14];
      C31[whI + 2][stI - 1] = new TMatrixD *[14];
      b31[whI + 2][stI - 1] = new TMatrixD *[14];
      for (int seI = 1; seI < 15; ++seI) {
        if (seI > 12 && stI != 4)
          continue;
        C13[whI + 2][stI - 1][seI - 1] = new TMatrixD(3, 3);
        b13[whI + 2][stI - 1][seI - 1] = new TMatrixD(3, 1);
        C31[whI + 2][stI - 1][seI - 1] = new TMatrixD(3, 3);
        b31[whI + 2][stI - 1][seI - 1] = new TMatrixD(3, 1);
      }
    }
  }

  //Run over the TTree
  Int_t nentries = (Int_t)tali->GetEntries();
  for (Int_t i = 0; i < nentries; i++) {
    tali->GetEntry(i);
    //Basic cuts
    if (pt > ptMax || pt < ptMin)
      continue;

    bool repeatedHits = false;
    for (int counter = 0; counter < nseg; ++counter) {
      //Make sure there are no repeated hits
      for (int counterHi = 0; counterHi < nhits[counter]; counterHi++) {
        for (int counterHj = 0; counterHj < nhits[counter]; counterHj++) {
          if (counterHi == counterHj)
            continue;
          if (zc[counter][counterHi] == zc[counter][counterHj]) {
            repeatedHits = true;
          }
        }
      }
      if (repeatedHits == true)
        continue;

      float x_13 = xSlSL3[counter];
      float xp_13 = xSL1SL3[counter];
      float x_31 = xSlSL1[counter];
      float xp_31 = xSL3SL1[counter];
      //float tanphi = dxdzSl[counter];
      float tanphi_13 = dxdzSlSL1[counter];
      float tanphi_31 = dxdzSlSL3[counter];
      int wheel = wh[counter];
      int station = st[counter];
      int sector = sr[counter];

      if (fabs(x_13 - xp_13) < 3 && fabs(x_31 - xp_31) && fabs(tanphi_13 - tanphi_31) < 0.06) {
        *(C13[wheel + 2][station - 1][sector - 1]) += returnCSLMatrix(x_13, xp_13, tanphi_13);
        *(b13[wheel + 2][station - 1][sector - 1]) += returnbSLMatrix(x_13, xp_13, tanphi_13);

        *(C31[wheel + 2][station - 1][sector - 1]) += returnCSLMatrix(x_31, xp_31, tanphi_31);
        *(b31[wheel + 2][station - 1][sector - 1]) += returnbSLMatrix(x_31, xp_31, tanphi_31);
      }
    }
  }

  for (int wheel = -2; wheel < 3; ++wheel) {
    for (int station = 1; station < 5; ++station) {
      for (int sector = 1; sector < 15; ++sector) {
        if (sector > 12 && station != 4)
          continue;
        TMatrixD solution13(3, 1);
        TMatrixD solution31(3, 1);
        TMatrixD C31_copy = *(C31[wheel + 2][station - 1][sector - 1]);
        TMatrixD C13_copy = *(C13[wheel + 2][station - 1][sector - 1]);
        TMatrixD b31_copy = *(b31[wheel + 2][station - 1][sector - 1]);
        TMatrixD b13_copy = *(b13[wheel + 2][station - 1][sector - 1]);

        C31_copy.Invert();
        C13_copy.Invert();
        solution13 = C13_copy * b13_copy;
        solution31 = C31_copy * b31_copy;
        whC = wheel;
        stC = station;
        srC = sector;
        dx = solution13(0, 0);
        dz = solution13(1, 0);
        phiy = solution13(2, 0);
        for (int c = 0; c < 3; ++c) {
          for (int s = 0; s < 3; ++s) {
            cov[c][s] = C13_copy(c, s);
          }
        }
        ttreeOutput->Fill();
      }
    }
  }
  f->Write();
}

returnbSLMatrix()

TMatrixD DTMuonSLToSL::returnbSLMatrix	(	float	x,
		float	xp,
		float	tanphi
	)

Definition at line 142 of file DTMuonSLToSL.cc.

References makeMuonMisalignmentScenario::matrix, and x.

Referenced by calculationSLToSL().

                                                                      {
  TMatrixD matrix(3, 1);

  matrix(0, 0) = -(x - xp);
  matrix(1, 0) = -(x - xp) * tanphi;
  matrix(2, 0) = -(x - xp) * tanphi * xp;

  return matrix;
}

returnCSLMatrix()

TMatrixD DTMuonSLToSL::returnCSLMatrix	(	float	x,
		float	xp,
		float	tanphi
	)

Definition at line 126 of file DTMuonSLToSL.cc.

References makeMuonMisalignmentScenario::matrix.

Referenced by calculationSLToSL().

                                                                      {
  TMatrixD matrix(3, 3);

  matrix(0, 0) = 1.0;
  matrix(1, 0) = tanphi;
  matrix(0, 1) = tanphi;
  matrix(1, 1) = tanphi * tanphi;
  matrix(0, 2) = tanphi * xp;
  matrix(2, 0) = tanphi * xp;
  matrix(2, 2) = tanphi * tanphi * xp * xp;
  matrix(2, 1) = tanphi * tanphi * xp;
  matrix(1, 2) = tanphi * tanphi * xp;

  return matrix;
}

setBranchTree()

void DTMuonSLToSL::setBranchTree

(

)

Definition at line 152 of file DTMuonSLToSL.cc.

References cov, dx, dz, phiy, srC, stC, ttreeOutput, and whC.

Referenced by DTMuonSLToSL().

                                 {
  ttreeOutput = new TTree("DTSLToSLResult", "DTSLToSLResult");

  ttreeOutput->Branch("wh", &whC, "wh/F");
  ttreeOutput->Branch("st", &stC, "st/F");
  ttreeOutput->Branch("sr", &srC, "sr/F");
  ttreeOutput->Branch("dx", &dx, "dx/F");
  ttreeOutput->Branch("dz", &dz, "dz/F");
  ttreeOutput->Branch("phiy", &phiy, "phiy/F");
  ttreeOutput->Branch("cov", cov, "cov[3][3]/F");
}

Member Data Documentation

cov

float DTMuonSLToSL::cov[3][3]

private

Definition at line 37 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and setBranchTree().

dx

float DTMuonSLToSL::dx

private

Definition at line 36 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and setBranchTree().

dz

float DTMuonSLToSL::dz

private

Definition at line 36 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), Electron.Electron::cutBasedId(), ntupleDataFormat.Track::dzPull(), and setBranchTree().

phiy

float DTMuonSLToSL::phiy

private

Definition at line 36 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and setBranchTree().

ptMax

float DTMuonSLToSL::ptMax

private

Definition at line 40 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and DTMuonSLToSL().

ptMin

float DTMuonSLToSL::ptMin

private

Definition at line 40 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and DTMuonSLToSL().

srC

int DTMuonSLToSL::srC

private

Definition at line 35 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and setBranchTree().

stC

int DTMuonSLToSL::stC

private

Definition at line 35 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and setBranchTree().

ttreeOutput

TTree* DTMuonSLToSL::ttreeOutput

private

Definition at line 42 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and setBranchTree().

whC

int DTMuonSLToSL::whC

private

Definition at line 35 of file DTMuonSLToSL.h.

Referenced by calculationSLToSL(), and setBranchTree().