#include <MuonDTLocalMillepedeAlgorithm.h>

Inheritance diagram for MuonDTLocalMillepedeAlgorithm:

Classes
struct	Info1D

Public Member Functions
void	initialize (const edm::EventSetup &setup, AlignableTracker tracker, AlignableMuon muon, AlignmentParameterStore *store)
	Call at beginning of job. More...

	MuonDTLocalMillepedeAlgorithm (const edm::ParameterSet &cfg)
	Constructor. More...

void	run (const edm::EventSetup &setup, const EventInfo &eventInfo)
	Run the algorithm on trajectories and tracks. More...

void	terminate (void)
	Call at end of job. More...

	~MuonDTLocalMillepedeAlgorithm ()
	Destructor. More...

Public Member Functions inherited from AlignmentAlgorithmBase
virtual bool	addCalibrations (const Calibrations &)

bool	addCalibrations (const CalibrationsOwner &cals)

	AlignmentAlgorithmBase (const edm::ParameterSet &)
	Constructor. More...

virtual void	beginLuminosityBlock (const edm::EventSetup &setup)
	called at begin of luminosity block (no lumi block info passed yet) More...

virtual void	beginRun (const edm::Run &, const edm::EventSetup &, bool changed)
	called at begin of run More...

virtual void	endLuminosityBlock (const edm::EventSetup &setup)
	called at end of luminosity block (no lumi block info passed yet) More...

virtual void	endRun (const EndRunInfo &runInfo, const edm::EventSetup &setup)
	called at end of run - order of arguments like in EDProducer etc. More...

virtual void	initialize (const edm::EventSetup &setup, AlignableTracker tracker, AlignableMuon muon, AlignableExtras extras, AlignmentParameterStore store)=0
	Call at beginning of job (must be implemented in derived class) More...

virtual bool	processesEvents ()
	Returns whether algorithm proccesses events in current configuration. More...

virtual bool	setParametersForRunRange (const RunRange &rr)

virtual void	startNewLoop ()

virtual bool	storeAlignments ()
	Returns whether algorithm produced results to be stored. More...

virtual bool	supportsCalibrations ()

virtual void	terminate (const edm::EventSetup &iSetup)=0
	Call at end of each loop (must be implemented in derived class) More...

virtual	~AlignmentAlgorithmBase ()
	Destructor. More...

Private Member Functions
bool	build4DSegments ()

void	setBranchTrees ()

Private Attributes
float	charge

edm::InputTag	consTraj

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

edm::Service< TFileService >	fs

edm::InputTag	globalTracks

int	la [5][14]

float	meandxdz [5][4][14]

float	meandydz [5][4][14]

float	meanx [5][4][14]

float	meany [5][4][14]

Info1D	myTrack1D

int	nhits [5]

int	nMtxSection

int	nphihits [5]

float	nPhihits

int	nseg

int	nthetahits [5]

float	nThetahits

std::string	ntuplePath

int	numberOfRootFiles

float	numberOfSigmasDXDZ

float	numberOfSigmasDYDZ

float	numberOfSigmasX

float	numberOfSigmasY

float	p

float	phi

float	pt

float	ptMax

float	ptMin

float	sigmadxdz [5][4][14]

float	sigmadydz [5][4][14]

float	sigmax [5][4][14]

float	sigmay [5][4][14]

int	sl [5][14]

int	sr [5]

int	st [5]

TChain *	tali

AlignableNavigator *	theAlignableDetAccessor

align::Alignables	theAlignables

AlignmentParameterStore *	theAlignmentParameterStore

TTree *	ttreeOutput

int	wh [5]

int	workingmode

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]

Additional Inherited Members
Public Types inherited from AlignmentAlgorithmBase
typedef std::pair< const Trajectory , const reco::Track >	ConstTrajTrackPair

typedef std::vector< ConstTrajTrackPair >	ConstTrajTrackPairCollection

using	RunNumber = align::RunNumber

using	RunRange = align::RunRange

Detailed Description

Definition at line 30 of file MuonDTLocalMillepedeAlgorithm.h.

Constructor & Destructor Documentation

◆ MuonDTLocalMillepedeAlgorithm()

MuonDTLocalMillepedeAlgorithm::MuonDTLocalMillepedeAlgorithm ( const edm::ParameterSet & cfg )

Constructor.

Definition at line 41 of file MuonDTLocalMillepedeAlgorithm.cc.

                                                                                       :
   AlignmentAlgorithmBase( cfg )
 {
   
   edm::LogInfo("Alignment") << "[MuonDTLocalMillepedeAlgorithm] constructed.";
  
   //Parse parameters. In the future this section should be completed with more options
   ntuplePath = cfg.getParameter<std::string>( "ntuplePath" );
   numberOfRootFiles = cfg.getParameter<int>( "numberOfRootFiles" ); 
   ptMax = cfg.getParameter<double>( "ptMax" );
   ptMin = cfg.getParameter<double>( "ptMin" );
   numberOfSigmasX = cfg.getParameter<double>( "numberOfSigmasX" );
   numberOfSigmasDXDZ = cfg.getParameter<double>( "numberOfSigmasDXDZ" );
   numberOfSigmasY = cfg.getParameter<double>( "numberOfSigmasY" );
   numberOfSigmasDYDZ = cfg.getParameter<double>( "numberOfSigmasDYDZ" );
   nPhihits = cfg.getParameter<double>( "nPhihits" );
   nThetahits = cfg.getParameter<double>( "nThetahits" );
   workingmode = cfg.getParameter<int>( "workingMode" );
   nMtxSection = cfg.getParameter<int>( "nMtxSection" );
  
  
   //The algorithm has three working modes:
   //0.- aligment information is extracted from the events and stored in root files
   //1.- The SLtoSl algorithm
   //2.- The local MuonMillepede algorithm
   if(workingmode == 0) {
     edm::LogInfo("Alignment") << "[MuonDTLocalMillepedeAlgorithm] Running on production mode."; 
     char nameOfFile[200];
     snprintf(nameOfFile, sizeof(nameOfFile), "%s/MyNtupleResidual.root", ntuplePath.c_str());
     f = new TFile(nameOfFile, "RECREATE");
     f->cd();
     setBranchTrees();
   } else if (workingmode == 1) {
     edm::LogInfo("Alignment") << "[MuonDTLocalMillepedeAlgorithm] Running SLToSL algorithm."; 
   } else {
     edm::LogInfo("Alignment") << "[MuonDTLocalMillepedeAlgorithm] Running Local Millepede algorithm."; 
   }
  
 }

References looper::cfg, f, nMtxSection, nPhihits, nThetahits, ntuplePath, numberOfRootFiles, numberOfSigmasDXDZ, numberOfSigmasDYDZ, numberOfSigmasX, numberOfSigmasY, ptMax, ptMin, setBranchTrees(), AlCaHLTBitMon_QueryRunRegistry::string, and workingmode.

◆ ~MuonDTLocalMillepedeAlgorithm()

MuonDTLocalMillepedeAlgorithm::~MuonDTLocalMillepedeAlgorithm ( )

inline

Destructor.

Definition at line 39 of file MuonDTLocalMillepedeAlgorithm.h.

39 {};

Member Function Documentation

◆ build4DSegments()

bool MuonDTLocalMillepedeAlgorithm::build4DSegments ( )

private

Definition at line 214 of file MuonDTLocalMillepedeAlgorithm.cc.

                                                     {
  
   bool saveThis = false;
   
   //Set to 0
   int id[20][5];
   int numlayer[20][12];
   for(int s = 0; s < 20; ++s) {
     for(int k = 0; k < 5; ++k) id[s][k] = 0;
     for(int k = 0; k < 12; ++k) numlayer[s][k] = 0;
   }
  
   
   int nChambers = 0;
   for(int counter = 0; counter < myTrack1D.nhits; ++counter) {
     bool isNew = true;
     for(int counterCham = 0; counterCham < nChambers; counterCham++) {
       if(myTrack1D.wh[counter] == id[counterCham][0] &&
      myTrack1D.st[counter] == id[counterCham][1] &&
          myTrack1D.sr[counter] == id[counterCham][2]) {
     if(myTrack1D.sl[counter] == 2) { 
       id[counterCham][4]++;
     } else {
       id[counterCham][3]++;
     }
     for (int ila = 1; ila<=4; ila++)
       if (myTrack1D.la[counter]==ila) {
         int jla = (myTrack1D.sl[counter]-1)*4 + ila -1;
         numlayer[counterCham][jla]++;
       }
     isNew = false;
       } 
     }
     if(isNew) {
       id[nChambers][0] = myTrack1D.wh[counter];
       id[nChambers][1] = myTrack1D.st[counter];
       id[nChambers][2] = myTrack1D.sr[counter];
       if(myTrack1D.sl[counter] == 2) {
         id[nChambers][4]++;
       } else {
         id[nChambers][3]++;
       }
       for (int ila = 1; ila<=4; ila++)
     if (myTrack1D.la[counter]==ila) {
       int jla = (myTrack1D.sl[counter]-1)*4 + ila -1;
       numlayer[nChambers][jla]++;
     }
       nChambers++;
     }
   }
   
   for (int iseg = 0; iseg<MAX_SEGMENT; iseg++) 
     for (int ihit = 0; ihit<MAX_HIT_CHAM; ihit++) {
       xc[iseg][ihit] = -250.;
       yc[iseg][ihit] = -250.;
       zc[iseg][ihit] = -250.;
       ex[iseg][ihit] = -250.;
       xcp[iseg][ihit] = -250.;
       ycp[iseg][ihit] = -250.;
       excp[iseg][ihit] = -250.;
       eycp[iseg][ihit] = -250.;
       sl[iseg][ihit] = 0;
       la[iseg][ihit] = 0;
     }
  
   nseg = 0;
   for(int counter = 0; counter < nChambers; ++counter) {
     
     bool GoodPhiChamber = true, GoodThetaChamber = true;
     for (int ila = 1; ila<=12; ila++) {
       if (numlayer[counter][ila-1]!=1 && (ila<5 || ila>8)) GoodPhiChamber = false;
       if (numlayer[counter][ila-1]!=1 && (ila<9 || ila>4) && id[counter][1]!=4) GoodThetaChamber = false;
     }
  
     if(id[counter][3] >= nPhihits && (id[counter][4] >= nThetahits || id[counter][1] == 4) &&
        GoodPhiChamber && GoodThetaChamber) {
  
       TMatrixD phiProjection(2,2);
       TMatrixD thetaProjection(2,2);
       TMatrixD bphiProjection(2,1);
       TMatrixD bthetaProjection(2,1);
       
       TMatrixD phiProjectionSL1(2,2);
       TMatrixD bphiProjectionSL1(2,1);
       TMatrixD phiProjectionSL3(2,2);
       TMatrixD bphiProjectionSL3(2,1);
      
       float SL1_z_ave = 0;
       float SL3_z_ave = 0;
  
       int numh1 = 0, numh2 = 0, numh3 = 0;
       for(int counterH = 0; counterH < myTrack1D.nhits; ++counterH) {
         if(myTrack1D.wh[counterH] == id[counter][0] && myTrack1D.st[counterH] == id[counter][1] &&
            myTrack1D.sr[counterH] == id[counter][2]) {
       if(myTrack1D.sl[counterH] == 2) {
         numh2++;
             thetaProjection(0,0) += 1.0/myTrack1D.erx[counterH];
             thetaProjection(0,1) += myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
             thetaProjection(1,0) += myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
             thetaProjection(1,1) += myTrack1D.zc[counterH]*myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
             bthetaProjection(0,0) += myTrack1D.yc[counterH]/myTrack1D.erx[counterH];
             bthetaProjection(1,0) += myTrack1D.yc[counterH]*myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
           } else {
             phiProjection(0,0) += 1.0/myTrack1D.erx[counterH];
             phiProjection(0,1) += myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
             phiProjection(1,0) += myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
             phiProjection(1,1) += myTrack1D.zc[counterH]*myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
             bphiProjection(0,0) += myTrack1D.xc[counterH]/myTrack1D.erx[counterH];
             bphiProjection(1,0) += myTrack1D.xc[counterH]*myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
         if(myTrack1D.sl[counterH] == 1) {
           numh1++;
               phiProjectionSL1(0,0) += 1.0/myTrack1D.erx[counterH];
               phiProjectionSL1(0,1) += myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
               phiProjectionSL1(1,0) += myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
               phiProjectionSL1(1,1) += myTrack1D.zc[counterH]*myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
               bphiProjectionSL1(0,0) += myTrack1D.xc[counterH]/myTrack1D.erx[counterH];
               bphiProjectionSL1(1,0) += myTrack1D.xc[counterH]*myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
               SL1_z_ave += myTrack1D.zc[counterH];
             } else {
           numh3++;
               phiProjectionSL3(0,0) += 1.0/myTrack1D.erx[counterH];
               phiProjectionSL3(0,1) += myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
               phiProjectionSL3(1,0) += myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
               phiProjectionSL3(1,1) += myTrack1D.zc[counterH]*myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
               bphiProjectionSL3(0,0) += myTrack1D.xc[counterH]/myTrack1D.erx[counterH];
               bphiProjectionSL3(1,0) += myTrack1D.xc[counterH]*myTrack1D.zc[counterH]/myTrack1D.erx[counterH];
               SL3_z_ave += myTrack1D.zc[counterH];
             } 
           }
     }
       }
  
       SL1_z_ave /= 4.0;
       SL3_z_ave /= 4.0;
       
       if (phiProjection(0,0) != 0 && phiProjectionSL1(0,0) != 0 && phiProjectionSL3(0,0) != 0 && 
       (thetaProjection(0,0) != 0 || id[counter][1] == 4)) {
     
     wh[nseg] = id[counter][0];
     st[nseg] = id[counter][1];
     sr[nseg] = id[counter][2];
     
     if(thetaProjection(0,0) != 0 && id[counter][1] != 4) { // Already asked (almost)
       thetaProjection.Invert();
       TMatrixD solution = thetaProjection*bthetaProjection;
       ySl[nseg] = solution(0,0);
       dydzSl[nseg] = solution(1,0);
       eySl[nseg] = thetaProjection(0,0);
       edydzSl[nseg] = thetaProjection(1,1);
       eydydzSl[nseg] = thetaProjection(0,1);
     }
     phiProjection.Invert();
     phiProjectionSL1.Invert();
     phiProjectionSL3.Invert();
     TMatrixD solution = phiProjection*bphiProjection;
     TMatrixD solutionSL1 = phiProjectionSL1*bphiProjectionSL1;
     TMatrixD solutionSL3 = phiProjectionSL3*bphiProjectionSL3;
     xSl[nseg] = solution(0,0);
     dxdzSl[nseg] = solution(1,0);
     exSl[nseg] = phiProjection(0,0);
     edxdzSl[nseg] = phiProjection(1,1);
     exdxdzSl[nseg] = phiProjection(0,0);
     xSlSL1[nseg] = solutionSL1(0,0);
     dxdzSlSL1[nseg] = solutionSL1(1,0);
     exSlSL1[nseg] = phiProjectionSL1(0,0);
     edxdzSlSL1[nseg] = phiProjectionSL1(1,1);
     exdxdzSlSL1[nseg] = phiProjectionSL1(0,0);
         xSL1SL3[nseg] = solutionSL1(0,0) + SL3_z_ave * solutionSL1(1,0);
     xSlSL3[nseg] = solutionSL3(0,0);
     dxdzSlSL3[nseg] = solutionSL3(1,0);
     exSlSL3[nseg] = phiProjectionSL3(0,0);
     edxdzSlSL3[nseg] = phiProjectionSL3(1,1);
     exdxdzSlSL3[nseg] = phiProjectionSL3(0,0);
         xSL3SL1[nseg] = solutionSL3(0,0) + SL1_z_ave * solutionSL3(1,0);
         int hitcounter = 0;
         for(int counterH = 0; counterH < myTrack1D.nhits; ++counterH) {
           if(myTrack1D.wh[counterH] == wh[nseg] && myTrack1D.st[counterH] == st[nseg] &&
              myTrack1D.sr[counterH] == sr[nseg]) {
             xc[nseg][hitcounter] = myTrack1D.xc[counterH];
             yc[nseg][hitcounter] = myTrack1D.yc[counterH];
             zc[nseg][hitcounter] = myTrack1D.zc[counterH];
             ex[nseg][hitcounter] = myTrack1D.erx[counterH];
             xcp[nseg][hitcounter] = xSl[nseg]+dxdzSl[nseg]*myTrack1D.zc[counterH];
             ycp[nseg][hitcounter] = ySl[nseg]+dydzSl[nseg]*myTrack1D.zc[counterH];
             excp[nseg][hitcounter] = exSl[nseg]*exSl[nseg]+ (edxdzSl[nseg]*edxdzSl[nseg])*myTrack1D.zc[counterH];
             eycp[nseg][hitcounter] = eySl[nseg]*eySl[nseg]+ (edydzSl[nseg]*edydzSl[nseg])*myTrack1D.zc[counterH];
             sl[nseg][hitcounter] = myTrack1D.sl[counterH];
             la[nseg][hitcounter] = myTrack1D.la[counterH];
             saveThis = true;
             hitcounter++;
           }
         }
     nphihits[nseg] = id[counter][3];
     nthetahits[nseg] = id[counter][4];
     nhits[nseg] = hitcounter;
     nseg++;
       }
     }
   }
   return saveThis;
 }

Referenced by run().

◆ initialize()

void MuonDTLocalMillepedeAlgorithm::initialize	(	const edm::EventSetup &	setup,
		AlignableTracker *	tracker,
		AlignableMuon *	muon,
		AlignmentParameterStore *	store
	)

Call at beginning of job.

Definition at line 86 of file MuonDTLocalMillepedeAlgorithm.cc.

 {
   
   // accessor Det->AlignableDet
   if ( !muon )
     theAlignableDetAccessor = new AlignableNavigator(tracker);
   else if ( !tracker )
     theAlignableDetAccessor = new AlignableNavigator(muon);
   else 
     theAlignableDetAccessor = new AlignableNavigator(tracker,muon);
   
   // set alignmentParameterStore
   theAlignmentParameterStore=store;
   
   // get alignables
   theAlignables = theAlignmentParameterStore->alignables();
   
 }

References AlignmentParameterStore::alignables(), theAlignableDetAccessor, theAlignables, theAlignmentParameterStore, and PbPb_ZMuSkimMuonDPG_cff::tracker.

◆ run()

void MuonDTLocalMillepedeAlgorithm::run	(	const edm::EventSetup &	setup,
		const EventInfo &	eventInfo
	)

virtual

Run the algorithm on trajectories and tracks.

Implements AlignmentAlgorithmBase.

Definition at line 132 of file MuonDTLocalMillepedeAlgorithm.cc.

 {
  
   //Only important in the production mode
   if(workingmode != 0) return;
  
   const ConstTrajTrackPairCollection &tracks = eventInfo.trajTrackPairs();
   for( ConstTrajTrackPairCollection::const_iterator it=tracks.begin();
        it!=tracks.end();it++) {
  
     const Trajectory *traj = (*it).first;
     const reco::Track *track = (*it).second;
     
     p    = track->p();
     pt    = track->pt();
     eta   = track->eta();
     phi   = track->phi();
     charge = track->charge();
  
  
     if(pt < ptMin || pt > ptMax) continue;
     
     vector<const TransientTrackingRecHit*> hitvec;
     vector<TrajectoryMeasurement> measurements = traj->measurements();
  
     //In this loop the measurements and hits are extracted and put into two vectors 
     int ch_muons = 0;
     for (vector<TrajectoryMeasurement>::iterator im=measurements.begin();
      im!=measurements.end(); im++) {
       TrajectoryMeasurement meas = *im;
       const TransientTrackingRecHit* hit = &(*meas.recHit());
       //We are not very strict at this point
       if (hit->isValid()) {
         if(hit->det()->geographicalId().det() == 2 && hit->det()->geographicalId().subdetId() == 1) {
           hitvec.push_back(hit);
           ch_muons++;
     } 
       }
     }
  
  
     vector<const TransientTrackingRecHit*>::const_iterator ihit=hitvec.begin();
     //Information is stored temporally in the myTrack1D object, which is analyzed
     //in the build4DSegments method in order to associate hits to segments. 
     int ch_counter = 0;
     while (ihit != hitvec.end()) 
       {
     const GeomDet* det=(*ihit)->det();
         if(det->geographicalId().det() == 2) {
       if(det->geographicalId().subdetId() == 1) {
             DTLayerId mLayer(det->geographicalId().rawId());
             DTChamberId mChamber(mLayer.wheel(), mLayer.station(), mLayer.sector());
             AlignableDet *aliDet = theAlignableDetAccessor->alignableDetFromDetId(mChamber);
         myTrack1D.wh[ch_counter] = mLayer.wheel();
             myTrack1D.st[ch_counter] = mLayer.station();
             myTrack1D.sr[ch_counter] = mLayer.sector();
             myTrack1D.sl[ch_counter] = mLayer.superlayer();
             myTrack1D.la[ch_counter] = mLayer.layer();
             myTrack1D.erx[ch_counter] = (*ihit)->localPositionError().xx();
             align::GlobalPoint globhit = det->toGlobal((*ihit)->localPosition());
             align::LocalPoint seghit = aliDet->surface().toLocal(globhit);
             myTrack1D.xc[ch_counter] = seghit.x();
         myTrack1D.yc[ch_counter] = seghit.y();
         myTrack1D.zc[ch_counter] = seghit.z();
         ch_counter++;
       }
         }
         ihit++;
       }
     myTrack1D.nhits = ch_counter;
     if(build4DSegments()) ttreeOutput->Fill();
   }
   
 }

Referenced by DTWorkflow.DTWorkflow::all(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), and o2olib.O2OTool::execute().

◆ setBranchTrees()

void MuonDTLocalMillepedeAlgorithm::setBranchTrees ( )

private

Definition at line 419 of file MuonDTLocalMillepedeAlgorithm.cc.

                                                    {
  
   ttreeOutput = new TTree("InfoTuple", "InfoTuple");
   
   ttreeOutput->Branch("p", &p, "p/F");
   ttreeOutput->Branch("pt", &pt, "pt/F");
   ttreeOutput->Branch("eta", &eta, "eta/F");
   ttreeOutput->Branch("phi", &phi, "phi/F");
   ttreeOutput->Branch("charge", &charge, "charge/F");
   ttreeOutput->Branch("nseg", &nseg, "nseg/I");
   ttreeOutput->Branch("nphihits", nphihits, "nphihits[nseg]/I");
   ttreeOutput->Branch("nthetahits", nthetahits, "nthetahits[nseg]/I");
   ttreeOutput->Branch("nhits", nhits, "nhits[nseg]/I");
   ttreeOutput->Branch("xSl", xSl, "xSl[nseg]/F");
   ttreeOutput->Branch("dxdzSl", dxdzSl, "dxdzSl[nseg]/F");
   ttreeOutput->Branch("exSl", exSl, "exSl[nseg]/F");
   ttreeOutput->Branch("edxdzSl", edxdzSl, "edxdzSl[nseg]/F");
   ttreeOutput->Branch("exdxdzSl", edxdzSl, "exdxdzSl[nseg]/F");
   ttreeOutput->Branch("ySl", ySl, "ySl[nseg]/F");
   ttreeOutput->Branch("dydzSl", dydzSl, "dydzSl[nseg]/F");
   ttreeOutput->Branch("eySl", eySl, "eySl[nseg]/F");
   ttreeOutput->Branch("edydzSl", edydzSl, "edydzSl[nseg]/F");
   ttreeOutput->Branch("eydydzSl", eydydzSl, "eydydzSl[nseg]/F");
   ttreeOutput->Branch("xSlSL1", xSlSL1, "xSlSL1[nseg]/F");
   ttreeOutput->Branch("dxdzSlSL1", dxdzSlSL1, "dxdzSlSL1[nseg]/F");
   ttreeOutput->Branch("exSlSL1", exSlSL1, "exSlSL1[nseg]/F");
   ttreeOutput->Branch("edxdzSlSL1", edxdzSlSL1, "edxdzSlSL1[nseg]/F");
   ttreeOutput->Branch("xSL1SL3", xSL1SL3, "xSL1SL3[nseg]/F");
   ttreeOutput->Branch("xSlSL3", xSlSL3, "xSlSL3[nseg]/F");
   ttreeOutput->Branch("dxdzSlSL3", dxdzSlSL3, "dxdzSlSL3[nseg]/F");
   ttreeOutput->Branch("exSlSL3", exSlSL3, "exSlSL3[nseg]/F");
   ttreeOutput->Branch("edxdzSlSL3", edxdzSlSL3, "edxdzSlSL3[nseg]/F");
   ttreeOutput->Branch("xSL3SL1", xSL3SL1, "xSL3SL1[nseg]/F");
   ttreeOutput->Branch("xc", xc, "xc[nseg][14]/F");
   ttreeOutput->Branch("yc", yc, "yc[nseg][14]/F");
   ttreeOutput->Branch("zc", zc, "zc[nseg][14]/F");
   ttreeOutput->Branch("ex", ex, "ex[nseg][14]/F");
   ttreeOutput->Branch("xcp", xcp, "xcp[nseg][14]/F");
   ttreeOutput->Branch("ycp", ycp, "ycp[nseg][14]/F");
   ttreeOutput->Branch("excp", excp, "excp[nseg][14]/F");
   ttreeOutput->Branch("eycp", eycp, "eycp[nseg][14]/F");
   ttreeOutput->Branch("wh", wh, "wh[nseg]/I");
   ttreeOutput->Branch("st", st, "st[nseg]/I");
   ttreeOutput->Branch("sr", sr, "sr[nseg]/I");
   ttreeOutput->Branch("sl", sl, "sl[nseg][14]/I");
   ttreeOutput->Branch("la", la, "la[nseg][14]/I");
  
 }

References charge, dxdzSl, dxdzSlSL1, dxdzSlSL3, dydzSl, edxdzSl, edxdzSlSL1, edxdzSlSL3, edydzSl, eta, ex, excp, exSl, exSlSL1, exSlSL3, eycp, eydydzSl, eySl, la, nhits, nphihits, nseg, nthetahits, p, phi, pt, sl, sr, st, ttreeOutput, wh, xc, xcp, xSl, xSL1SL3, xSL3SL1, xSlSL1, xSlSL3, yc, ycp, ySl, and zc.

Referenced by MuonDTLocalMillepedeAlgorithm().

◆ terminate()

void MuonDTLocalMillepedeAlgorithm::terminate ( void )

virtual

Call at end of job.

Reimplemented from AlignmentAlgorithmBase.

Definition at line 110 of file MuonDTLocalMillepedeAlgorithm.cc.

 {
  
   //If workingmode equals 1 or 2, the algorithms are run before saving.   
   if(workingmode == 1) {
     edm::LogInfo("Alignment") << "[MuonDTLocalMillepedeAlgorithm] Starting SLToSL algorithm";
     DTMuonSLToSL mySLToSL(ntuplePath, numberOfRootFiles, ptMax, ptMin, f);
   } else if(workingmode >= 2) {
     edm::LogInfo("Alignment") << "[MuonDTLocalMillepedeAlgorithm] Starting local MuonMillepede algorithm";
     DTMuonMillepede myMillepede(ntuplePath, numberOfRootFiles, ptMax, ptMin, nPhihits, nThetahits, workingmode, nMtxSection);
   } 
  
   if (workingmode==0) {
     f->Write();
     f->Close();
   }
  
 }

References f, nMtxSection, nPhihits, nThetahits, ntuplePath, numberOfRootFiles, ptMax, ptMin, and workingmode.

Member Data Documentation

◆ charge

float MuonDTLocalMillepedeAlgorithm::charge

private

Definition at line 78 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by run(), and setBranchTrees().

◆ consTraj

edm::InputTag MuonDTLocalMillepedeAlgorithm::consTraj

private

Definition at line 135 of file MuonDTLocalMillepedeAlgorithm.h.

◆ dxdzSl

float MuonDTLocalMillepedeAlgorithm::dxdzSl[5]

private

Definition at line 84 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ dxdzSlSL1

float MuonDTLocalMillepedeAlgorithm::dxdzSlSL1[5]

private

Definition at line 94 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ dxdzSlSL3

float MuonDTLocalMillepedeAlgorithm::dxdzSlSL3[5]

private

Definition at line 100 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ dydzSl

float MuonDTLocalMillepedeAlgorithm::dydzSl[5]

private

Definition at line 89 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ edxdzSl

float MuonDTLocalMillepedeAlgorithm::edxdzSl[5]

private

Definition at line 86 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ edxdzSlSL1

float MuonDTLocalMillepedeAlgorithm::edxdzSlSL1[5]

private

Definition at line 96 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ edxdzSlSL3

float MuonDTLocalMillepedeAlgorithm::edxdzSlSL3[5]

private

Definition at line 102 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ edydzSl

float MuonDTLocalMillepedeAlgorithm::edydzSl[5]

private

Definition at line 91 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ eta

float MuonDTLocalMillepedeAlgorithm::eta

private

Definition at line 78 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), Jet.Jet::puJetId(), run(), and setBranchTrees().

◆ ex

float MuonDTLocalMillepedeAlgorithm::ex[5][14]

private

Definition at line 108 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ excp

float MuonDTLocalMillepedeAlgorithm::excp[5][14]

private

Definition at line 111 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ exdxdzSl

float MuonDTLocalMillepedeAlgorithm::exdxdzSl[5]

private

Definition at line 87 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments().

◆ exdxdzSlSL1

float MuonDTLocalMillepedeAlgorithm::exdxdzSlSL1[5]

private

Definition at line 97 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments().

◆ exdxdzSlSL3

float MuonDTLocalMillepedeAlgorithm::exdxdzSlSL3[5]

private

Definition at line 103 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments().

◆ exSl

float MuonDTLocalMillepedeAlgorithm::exSl[5]

private

Definition at line 85 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ exSlSL1

float MuonDTLocalMillepedeAlgorithm::exSlSL1[5]

private

Definition at line 95 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ exSlSL3

float MuonDTLocalMillepedeAlgorithm::exSlSL3[5]

private

Definition at line 101 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ eycp

float MuonDTLocalMillepedeAlgorithm::eycp[5][14]

private

Definition at line 112 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ eydydzSl

float MuonDTLocalMillepedeAlgorithm::eydydzSl[5]

private

Definition at line 92 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ eySl

float MuonDTLocalMillepedeAlgorithm::eySl[5]

private

Definition at line 90 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ f

TFile* MuonDTLocalMillepedeAlgorithm::f

private

Definition at line 122 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by ztail.Decoder::follow(), ztail.Decoder::initial_synchronize(), MuonDTLocalMillepedeAlgorithm(), DQMIO2histo.DQMIO::print_index(), terminate(), electronCompare.flushfile::write(), electronStore.flushfile::write(), and DQMIO2histo.DQMIO::write_to_file().

◆ fs

edm::Service<TFileService> MuonDTLocalMillepedeAlgorithm::fs

private

Definition at line 132 of file MuonDTLocalMillepedeAlgorithm.h.

◆ globalTracks

edm::InputTag MuonDTLocalMillepedeAlgorithm::globalTracks

private

Definition at line 134 of file MuonDTLocalMillepedeAlgorithm.h.

◆ la

int MuonDTLocalMillepedeAlgorithm::la[5][14]

private

Definition at line 115 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ meandxdz

float MuonDTLocalMillepedeAlgorithm::meandxdz[5][4][14]

private

Definition at line 154 of file MuonDTLocalMillepedeAlgorithm.h.

◆ meandydz

float MuonDTLocalMillepedeAlgorithm::meandydz[5][4][14]

private

Definition at line 158 of file MuonDTLocalMillepedeAlgorithm.h.

◆ meanx

float MuonDTLocalMillepedeAlgorithm::meanx[5][4][14]

private

Definition at line 152 of file MuonDTLocalMillepedeAlgorithm.h.

◆ meany

float MuonDTLocalMillepedeAlgorithm::meany[5][4][14]

private

Definition at line 156 of file MuonDTLocalMillepedeAlgorithm.h.

◆ myTrack1D

Info1D MuonDTLocalMillepedeAlgorithm::myTrack1D

private

Definition at line 73 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and run().

◆ nhits

int MuonDTLocalMillepedeAlgorithm::nhits[5]

private

Definition at line 82 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ nMtxSection

int MuonDTLocalMillepedeAlgorithm::nMtxSection

private

Definition at line 143 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm(), and terminate().

◆ nphihits

int MuonDTLocalMillepedeAlgorithm::nphihits[5]

private

Definition at line 80 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ nPhihits

float MuonDTLocalMillepedeAlgorithm::nPhihits

private

Definition at line 139 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), MuonDTLocalMillepedeAlgorithm(), and terminate().

◆ nseg

int MuonDTLocalMillepedeAlgorithm::nseg

private

Definition at line 79 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ nthetahits

int MuonDTLocalMillepedeAlgorithm::nthetahits[5]

private

Definition at line 81 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ nThetahits

float MuonDTLocalMillepedeAlgorithm::nThetahits

private

Definition at line 140 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), MuonDTLocalMillepedeAlgorithm(), and terminate().

◆ ntuplePath

std::string MuonDTLocalMillepedeAlgorithm::ntuplePath

private

Definition at line 136 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm(), and terminate().

◆ numberOfRootFiles

int MuonDTLocalMillepedeAlgorithm::numberOfRootFiles

private

Definition at line 142 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm(), and terminate().

◆ numberOfSigmasDXDZ

float MuonDTLocalMillepedeAlgorithm::numberOfSigmasDXDZ

private

Definition at line 148 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm().

◆ numberOfSigmasDYDZ

float MuonDTLocalMillepedeAlgorithm::numberOfSigmasDYDZ

private

Definition at line 150 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm().

◆ numberOfSigmasX

float MuonDTLocalMillepedeAlgorithm::numberOfSigmasX

private

Definition at line 147 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm().

◆ numberOfSigmasY

float MuonDTLocalMillepedeAlgorithm::numberOfSigmasY

private

Definition at line 149 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm().

◆ p

float MuonDTLocalMillepedeAlgorithm::p

private

Definition at line 78 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by Electron.Electron::ptErr(), run(), and setBranchTrees().

◆ phi

float MuonDTLocalMillepedeAlgorithm::phi

private

Definition at line 78 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by Particle.Particle::__str__(), ntupleDataFormat.Track::phiPull(), run(), and setBranchTrees().

◆ pt

float MuonDTLocalMillepedeAlgorithm::pt

private

Definition at line 78 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by Particle.Particle::__str__(), ZMuMuRochCorAnalyzer.DiMuon::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), run(), setBranchTrees(), and Jet.Jet::setCorrP4().

◆ ptMax

float MuonDTLocalMillepedeAlgorithm::ptMax

private

Definition at line 137 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm(), run(), and terminate().

◆ ptMin

float MuonDTLocalMillepedeAlgorithm::ptMin

private

Definition at line 138 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm(), and terminate().

◆ sigmadxdz

float MuonDTLocalMillepedeAlgorithm::sigmadxdz[5][4][14]

private

Definition at line 155 of file MuonDTLocalMillepedeAlgorithm.h.

◆ sigmadydz

float MuonDTLocalMillepedeAlgorithm::sigmadydz[5][4][14]

private

Definition at line 159 of file MuonDTLocalMillepedeAlgorithm.h.

◆ sigmax

float MuonDTLocalMillepedeAlgorithm::sigmax[5][4][14]

private

Definition at line 153 of file MuonDTLocalMillepedeAlgorithm.h.

◆ sigmay

float MuonDTLocalMillepedeAlgorithm::sigmay[5][4][14]

private

Definition at line 157 of file MuonDTLocalMillepedeAlgorithm.h.

◆ sl

int MuonDTLocalMillepedeAlgorithm::sl[5][14]

private

Definition at line 114 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ sr

int MuonDTLocalMillepedeAlgorithm::sr[5]

private

Definition at line 113 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ st

int MuonDTLocalMillepedeAlgorithm::st[5]

private

Definition at line 113 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ tali

TChain* MuonDTLocalMillepedeAlgorithm::tali

private

Definition at line 124 of file MuonDTLocalMillepedeAlgorithm.h.

◆ theAlignableDetAccessor

AlignableNavigator* MuonDTLocalMillepedeAlgorithm::theAlignableDetAccessor

private

Definition at line 129 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by initialize(), and run().

◆ theAlignables

align::Alignables MuonDTLocalMillepedeAlgorithm::theAlignables

private

Definition at line 128 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by initialize().

◆ theAlignmentParameterStore

AlignmentParameterStore* MuonDTLocalMillepedeAlgorithm::theAlignmentParameterStore

private

Definition at line 127 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by initialize().

◆ ttreeOutput

TTree* MuonDTLocalMillepedeAlgorithm::ttreeOutput

private

Definition at line 123 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by run(), and setBranchTrees().

◆ wh

int MuonDTLocalMillepedeAlgorithm::wh[5]

private

Definition at line 113 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ workingmode

int MuonDTLocalMillepedeAlgorithm::workingmode

private

Definition at line 141 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by MuonDTLocalMillepedeAlgorithm(), run(), and terminate().

◆ xc

float MuonDTLocalMillepedeAlgorithm::xc[5][14]

private

Definition at line 105 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), geometryXMLparser.Alignable::covariance(), and setBranchTrees().

◆ xcp

float MuonDTLocalMillepedeAlgorithm::xcp[5][14]

private

Definition at line 109 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ xSl

float MuonDTLocalMillepedeAlgorithm::xSl[5]

private

Definition at line 83 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ xSL1SL3

float MuonDTLocalMillepedeAlgorithm::xSL1SL3[5]

private

Definition at line 98 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ xSL3SL1

float MuonDTLocalMillepedeAlgorithm::xSL3SL1[5]

private

Definition at line 104 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ xSlSL1

float MuonDTLocalMillepedeAlgorithm::xSlSL1[5]

private

Definition at line 93 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ xSlSL3

float MuonDTLocalMillepedeAlgorithm::xSlSL3[5]

private

Definition at line 99 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ yc

float MuonDTLocalMillepedeAlgorithm::yc[5][14]

private

Definition at line 106 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), geometryXMLparser.Alignable::covariance(), and setBranchTrees().

◆ ycp

float MuonDTLocalMillepedeAlgorithm::ycp[5][14]

private

Definition at line 110 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ ySl

float MuonDTLocalMillepedeAlgorithm::ySl[5]

private

Definition at line 88 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), and setBranchTrees().

◆ zc

float MuonDTLocalMillepedeAlgorithm::zc[5][14]

private

Definition at line 107 of file MuonDTLocalMillepedeAlgorithm.h.

Referenced by build4DSegments(), geometryXMLparser.Alignable::covariance(), and setBranchTrees().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MuonDTLocalMillepedeAlgorithm()

◆ ~MuonDTLocalMillepedeAlgorithm()

Member Function Documentation

◆ build4DSegments()

◆ initialize()

◆ run()

◆ setBranchTrees()

◆ terminate()

Member Data Documentation

◆ charge

◆ consTraj

◆ dxdzSl

◆ dxdzSlSL1

◆ dxdzSlSL3

◆ dydzSl

◆ edxdzSl

◆ edxdzSlSL1

◆ edxdzSlSL3

◆ edydzSl

◆ eta

◆ ex

◆ excp

◆ exdxdzSl

◆ exdxdzSlSL1

◆ exdxdzSlSL3

◆ exSl

◆ exSlSL1

◆ exSlSL3

◆ eycp

◆ eydydzSl

◆ eySl

◆ f

◆ fs

◆ globalTracks

◆ la

◆ meandxdz

◆ meandydz

◆ meanx

◆ meany

◆ myTrack1D

◆ nhits

◆ nMtxSection

◆ nphihits

◆ nPhihits

◆ nseg

◆ nthetahits

◆ nThetahits

◆ ntuplePath

◆ numberOfRootFiles

◆ numberOfSigmasDXDZ

◆ numberOfSigmasDYDZ

◆ numberOfSigmasX

◆ numberOfSigmasY

◆ p

◆ phi

◆ pt

◆ ptMax

◆ ptMin

◆ sigmadxdz

◆ sigmadydz

◆ sigmax

◆ sigmay

◆ sl

◆ sr

◆ st

◆ tali

◆ theAlignableDetAccessor

◆ theAlignables

◆ theAlignmentParameterStore

◆ ttreeOutput

◆ wh

◆ workingmode

◆ xc