#include <BmnUpstreamTracking.h>

Inheritance diagram for BmnUpstreamTracking:

Collaboration diagram for BmnUpstreamTracking:

Classes
struct	MC_points
struct	twotracks

Public Member Functions
	BmnUpstreamTracking ()
	BmnUpstreamTracking (Bool_t, Int_t)
virtual	~BmnUpstreamTracking ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()

Static Public Member Functions
static bool	compareSegments (const Smatch &a, const Smatch &b)

Private Member Functions
void	PrepareArraysToProcessEvent ()
void	ReadSiliconTracks (Double_t *, Double_t , Int_t &, vector< MC_points > &)
void	ReadSiliconHits (Double_t **, Int_t )
void	ReadMWPCSegments (Double_t *, Double_t , Int_t *, vector< MC_points > &)
void	ReadMWPCTracks (Double_t *, Double_t , Int_t *)
void	PCTracksSiTracksMatching (Double_t *, Double_t , Int_t &, Double_t *, Double_t , Int_t *, vector< Smatch > &, vector< Smatch > &)
void	PCSegmentsSiTracksMatching (Double_t *, Double_t , Int_t &, Double_t *, Double_t , Int_t *, vector< Smatch > &, vector< Smatch > &)
void	RecordingBest (Int_t &, vector< Smatch > &, vector< Smatch > &)
void	RecordCandidateSeg (Int_t &, Int_t &, Int_t &, Double_t , Double_t , Double_t , Double_t **, Double_t &, Double_t &, vector< Smatch > &)
void	GetAddSiHits (vector< Smatch > &, Double_t **, Int_t )
void	GetHitsToPoints (vector< Smatch > &, Double_t **, Int_t , Double_t **, Int_t , Double_t ***, Int_t &)
void	GlobalFit (Double_t , Double_t , Double_t *)
bool	InvertMatrix (Double_t , Double_t )
void	CalculationOfChi2 (Double_t ***, Int_t &, vector< UpTracks > &)
void	PrintAllTracks (vector< UpTracks > &)
void	TrackRecording (vector< UpTracks > &)
void	MCefficiencyCalculation (vector< MC_points > &, vector< UpTracks > &)

Private Attributes
Bool_t	fDebug = 0
UInt_t	fEventNo = 0
Int_t	fRunNumber
Bool_t	expData
TList	fList
TClonesArray *	fSiTracks
TClonesArray *	fSilHits
TClonesArray *	fMWPCTracks
TClonesArray *	fMWPCSegments
TString	fInputBranchName1
TString	fInputBranchName2
TString	fInputBranchName3
TString	fInputBranchHits
TClonesArray *	fSiTracksSim
TString	fInputBranchNameSimTrue
TClonesArray *	fBmnHitsArray
TString	fInputBranchName
TClonesArray *	fBmnUpstreamTracksArray
	Output array of Silicon Hits.
TString	fOutputTracksBranchName
TString	fOutputFileName
Double_t ***	par_ab_Ch
Double_t ***	par_ab_tr
Double_t ***	SiXYhits
Double_t ***	Points
Double_t **	par_ab_SiTr
Double_t **	par_Seg_z
Double_t **	par_ab_trz
Double_t **	Xforglfit
Double_t **	Amatr
Double_t *	par_SiTr_z
Double_t *	X_shift_seg
Double_t *	X_shiftUp
Double_t *	Y_shiftUp
Double_t *	Z_pair
Double_t *	Z_Chamber
Double_t *	AmatrInverted
Double_t *	rhs
Double_t *	AmatrArray
Double_t *	line
Int_t *	Nseg_Ch
Int_t *	NPCTracks
Int_t *	NSiXYhits
Int_t	NSiTracks
Int_t	NumPoints
Double_t	Shift_toCenterOfMagnetX
Double_t	Shift_toCenterOfMagnetY
Double_t	Shift_toCenterOfMagnetAX
Double_t	Shift_toCenterOfMagnetAY
Double_t	X_shift
Double_t	Y_shift
vector< Smatch >	vtmpSeg
vector< Smatch >	OutVector
vector< UpTracks >	vecUpTracks
vector< MC_points >	vec_points
vector< twotracks >	vec_twotracks
const Int_t	kBig = 200
const Int_t	kPoints = 5
const Int_t	kPoin_Par = 5
const Int_t	kNumPars = 10
const Int_t	kNumPairs = 2
const Int_t	kNumChambers = 4
const Int_t	kNumStSi = 4
const Int_t	kMaxMC = 100
const Int_t	fNstations = 4
const Int_t	kNPlanes = 6
const Double_t	kChi2_Max = 50.
const Double_t	dw = 0.25
const Double_t	SigmXmwpc = sqrt(2) * dw / 6.
const Double_t	SigmYmwpc = sqrt(2) * dw / 3.
const Double_t	cutAx = 0.02
const Double_t	cutAy = 0.02
const Double_t	cutX = 1.5
const Double_t	cutY = 2.0
const Double_t	Zcentr = -350.
const Double_t	kZ_target = -645.191
const Double_t	kZSi_cent = -392.5
const Int_t	PDG_Be7 = 1000040070
const Int_t	PDG_Li6 = 1000030060
const Int_t	PDG_Li7 = 1000030070
const Int_t	PDG_He3 = 1000020030
const Int_t	PDG_He4 = 1000020040
const Int_t	PDG_H2 = 1000010020
const Int_t	PDG_p = 2212
TH1D *	hAx_fitUp
TH1D *	hAy_fitUp
TH1D *	hx_fitUp
TH1D *	hy_fitUp
TH1D *	hchi2_fitUp
TH1D *	hNhitsUp
TH1D *	hdAx_tr_mc_comb
TH1D *	hdX_tr_mc_comb
TH1D *	hdAy_tr_mc_comb
TH1D *	hdY_tr_mc_comb
TH1D *	hdAx_uptr_mc
TH1D *	hdX_uptr_mc
TH1D *	hdAy_uptr_mc
TH1D *	hdY_uptr_mc
TH1D *	hDen_mcuptr
TH1D *	hNum_mcuptr
TH1D *	hEff_mcuptr
TH1D *	hAx_upmc
TH1D *	hAy_upmc
TH1D *	hX_upmc
TH1D *	hY_upmc
TH1D *	hNtr_reco
TH1D *	hNtr_mc
TH1D *	hNrecoTrif2mc
TH1D *	hAngle_reco
TH1D *	hAngle_recoifNmc2
TH1D *	hAngle_recoifNmc2Cases
TH2D *	hNtr_mc_vs_reco
TH2D *	hy_vs_x_Up
TH2D *	hY_vs_Xmctrue
TH2D *	hvertexXYUp
TH2D *	hTyTx_Up
vector< TH1D * >	hResXst
vector< TH1D * >	hResYst

Detailed Description

Definition at line 81 of file BmnUpstreamTracking.h.

Constructor & Destructor Documentation

◆ BmnUpstreamTracking() [1/2]

BmnUpstreamTracking::BmnUpstreamTracking ( )

Definition at line 1903 of file BmnUpstreamTracking.cxx.

References hAngle_reco, hAngle_recoifNmc2, hAngle_recoifNmc2Cases, hAx_fitUp, hAx_upmc, hAy_fitUp, hAy_upmc, hchi2_fitUp, hdAx_tr_mc_comb, hdAx_uptr_mc, hdAy_tr_mc_comb, hdAy_uptr_mc, hDen_mcuptr, hdX_tr_mc_comb, hdX_uptr_mc, hdY_tr_mc_comb, hdY_uptr_mc, hEff_mcuptr, hNhitsUp, hNrecoTrif2mc, hNtr_mc, hNtr_mc_vs_reco, hNtr_reco, hNum_mcuptr, hTyTx_Up, hvertexXYUp, hx_fitUp, hX_upmc, hy_fitUp, hY_upmc, hy_vs_x_Up, and hY_vs_Xmctrue.

◆ BmnUpstreamTracking() [2/2]

BmnUpstreamTracking::BmnUpstreamTracking	(	Bool_t	isExp,
		Int_t	runNumber )

Definition at line 1939 of file BmnUpstreamTracking.cxx.

References expData, fInputBranchHits, fInputBranchName, fInputBranchName1, fInputBranchName2, fInputBranchName3, fInputBranchNameSimTrue, fOutputTracksBranchName, and fRunNumber.

◆ ~BmnUpstreamTracking()

BmnUpstreamTracking::~BmnUpstreamTracking ( )

virtual

Definition at line 1957 of file BmnUpstreamTracking.cxx.

References Amatr, AmatrArray, AmatrInverted, i, kBig, kNumChambers, kNumPairs, kNumPars, kNumStSi, kPoints, line, NPCTracks, Nseg_Ch, NSiXYhits, par_ab_Ch, par_ab_SiTr, par_ab_tr, par_ab_trz, par_Seg_z, par_SiTr_z, Points, rhs, SiXYhits, X_shift_seg, X_shiftUp, Xforglfit, Y_shiftUp, Z_Chamber, and Z_pair.

Member Function Documentation

◆ CalculationOfChi2()

void BmnUpstreamTracking::CalculationOfChi2	(	Double_t ***	Points_,
		Int_t &	NumPoints_,
		vector< UpTracks > &	vec )

private

Definition at line 1214 of file BmnUpstreamTracking.cxx.

References Amatr, AmatrArray, AmatrInverted, UpTracks::Chi2, UpTracks::CoordX, UpTracks::CoordY, UpTracks::CoordZ, fDebug, GlobalFit(), i, InvertMatrix(), kPoin_Par, kPoints, UpTracks::Nhits, nhits, UpTracks::param, rhs, X_shiftUp, Xforglfit, Y_shiftUp, and Zcentr.

Referenced by Exec().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ compareSegments()

bool BmnUpstreamTracking::compareSegments	(	const Smatch &	a,
		const Smatch &	b )

inlinestatic

Definition at line 92 of file BmnUpstreamTracking.h.

Referenced by RecordingBest().

Here is the caller graph for this function:

◆ Exec()

void BmnUpstreamTracking::Exec ( Option_t * opt )

virtual

Definition at line 30 of file BmnUpstreamTracking.cxx.

References CalculationOfChi2(), expData, fDebug, fEventNo, GetAddSiHits(), GetHitsToPoints(), MCefficiencyCalculation(), NPCTracks, Nseg_Ch, NSiTracks, NSiXYhits, NumPoints, OutVector, par_ab_Ch, par_ab_SiTr, par_ab_tr, par_ab_trz, par_Seg_z, par_SiTr_z, PCSegmentsSiTracksMatching(), PCTracksSiTracksMatching(), Points, PrepareArraysToProcessEvent(), PrintAllTracks(), ReadMWPCSegments(), ReadMWPCTracks(), ReadSiliconHits(), ReadSiliconTracks(), SiXYhits, TrackRecording(), vec_points, vecUpTracks, vtmpSeg, and workTime.

Here is the call graph for this function:

◆ Finish()

void BmnUpstreamTracking::Finish ( )

virtual

Definition at line 2021 of file BmnUpstreamTracking.cxx.

References TFile::Close(), fDebug, fList, fOutputFileName, fRunNumber, hDen_mcuptr, hEff_mcuptr, hNum_mcuptr, and workTime.

Here is the call graph for this function:

◆ GetAddSiHits()

void BmnUpstreamTracking::GetAddSiHits	(	vector< Smatch > &	Out,
		Double_t ***	hits,
		Int_t *	Nh )

private

Definition at line 1061 of file BmnUpstreamTracking.cxx.

References fDebug, i, and kNumStSi.

Referenced by Exec().

Here is the caller graph for this function:

◆ GetHitsToPoints()

void BmnUpstreamTracking::GetHitsToPoints	(	vector< Smatch > &	Out,
		Double_t ***	ab_tr,
		Int_t *	NPCTracks_,
		Double_t ***	ab_Ch,
		Int_t *	Nseg_ch_,
		Double_t ***	Points_,
		Int_t &	NumPoints_ )

private

Definition at line 1109 of file BmnUpstreamTracking.cxx.

References fDebug, kNumStSi, SigmXmwpc, SigmYmwpc, Z_Chamber, and Z_pair.

Referenced by Exec().

Here is the caller graph for this function:

◆ GlobalFit()

void BmnUpstreamTracking::GlobalFit	(	Double_t **	XYZsigXY,
		Double_t **	Amatr_,
		Double_t *	rhs_ )

private

Definition at line 1527 of file BmnUpstreamTracking.cxx.

References kPoints.

Referenced by CalculationOfChi2().

Here is the caller graph for this function:

◆ Init()

InitStatus BmnUpstreamTracking::Init ( )

virtual

Definition at line 1577 of file BmnUpstreamTracking.cxx.

References Amatr, AmatrArray, AmatrInverted, expData, fBmnHitsArray, fBmnUpstreamTracksArray, fDebug, fInputBranchHits, fInputBranchName, fInputBranchName1, fInputBranchName2, fInputBranchName3, fInputBranchNameSimTrue, fList, fMWPCSegments, fMWPCTracks, fOutputTracksBranchName, fSilHits, fSiTracks, fSiTracksSim, hAngle_reco, hAngle_recoifNmc2, hAngle_recoifNmc2Cases, hAx_fitUp, hAx_upmc, hAy_fitUp, hAy_upmc, hchi2_fitUp, hdAx_tr_mc_comb, hdAx_uptr_mc, hdAy_tr_mc_comb, hdAy_uptr_mc, hDen_mcuptr, hdX_tr_mc_comb, hdX_uptr_mc, hdY_tr_mc_comb, hdY_uptr_mc, hEff_mcuptr, hNhitsUp, hNrecoTrif2mc, hNtr_mc, hNtr_mc_vs_reco, hNtr_reco, hNum_mcuptr, hResXst, hResYst, hTyTx_Up, hvertexXYUp, hx_fitUp, hX_upmc, hy_fitUp, hY_upmc, hy_vs_x_Up, hY_vs_Xmctrue, i, kBig, kNumChambers, kNumPairs, kNumPars, kNumStSi, kPoin_Par, kPoints, line, NPCTracks, Nseg_Ch, NSiXYhits, par_ab_Ch, par_ab_SiTr, par_ab_tr, par_ab_trz, par_Seg_z, par_SiTr_z, Points, rhs, Shift_toCenterOfMagnetAX, Shift_toCenterOfMagnetAY, Shift_toCenterOfMagnetX, Shift_toCenterOfMagnetY, SiXYhits, X_shift, X_shift_seg, X_shiftUp, Xforglfit, Y_shift, Y_shiftUp, Z_Chamber, and Z_pair.

◆ InvertMatrix()

bool BmnUpstreamTracking::InvertMatrix	(	Double_t *	m,
		Double_t *	invOut )

private

Definition at line 1478 of file BmnUpstreamTracking.cxx.

References i, and m.

Referenced by CalculationOfChi2().

Here is the caller graph for this function:

◆ MCefficiencyCalculation()

void BmnUpstreamTracking::MCefficiencyCalculation	(	vector< MC_points > &	points,
		vector< UpTracks > &	vecUp )

private

Definition at line 91 of file BmnUpstreamTracking.cxx.

References TMath::ACos(), fDebug, hAngle_reco, hAngle_recoifNmc2, hAngle_recoifNmc2Cases, hdAx_tr_mc_comb, hdAx_uptr_mc, hdAy_tr_mc_comb, hdAy_uptr_mc, hDen_mcuptr, hdX_tr_mc_comb, hdX_uptr_mc, hdY_tr_mc_comb, hdY_uptr_mc, hNrecoTrif2mc, hNtr_mc, hNtr_mc_vs_reco, hNtr_reco, hNum_mcuptr, hY_vs_Xmctrue, kMaxMC, PDG_Be7, PDG_H2, PDG_He3, PDG_He4, PDG_Li6, PDG_Li7, PDG_p, TMath::RadToDeg(), BmnUpstreamTracking::twotracks::track1, BmnUpstreamTracking::twotracks::track1_pdg, BmnUpstreamTracking::twotracks::track2, BmnUpstreamTracking::twotracks::track2_pdg, and vec_twotracks.

Referenced by Exec().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ PCSegmentsSiTracksMatching()

void BmnUpstreamTracking::PCSegmentsSiTracksMatching	(	Double_t **	par_Si,
		Double_t *	par_Siz,
		Int_t &	NTracksSi,
		Double_t ***	par_PC,
		Double_t **	par_PCz,
		Int_t *	NSegPC,
		vector< Smatch > &	vec_m,
		vector< Smatch > &	Out )

private

Definition at line 819 of file BmnUpstreamTracking.cxx.

References cutX, cutY, fabs(), fDebug, NSiTracks, RecordCandidateSeg(), and RecordingBest().

Referenced by Exec().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ PCTracksSiTracksMatching()

void BmnUpstreamTracking::PCTracksSiTracksMatching	(	Double_t **	par_Si,
		Double_t *	par_Siz,
		Int_t &	NTracksSi,
		Double_t ***	par_PC,
		Double_t **	par_PCz,
		Int_t *	NTracksPC,
		vector< Smatch > &	vec_m,
		vector< Smatch > &	Out )

private

Definition at line 732 of file BmnUpstreamTracking.cxx.

References Smatch::Chi2m, cutAx, cutAy, cutX, cutY, Smatch::distX, Smatch::distY, fabs(), fDebug, Smatch::Ind1, Smatch::Ind2, Smatch::Inds2, Smatch::Inds3, Smatch::Nhits1, Smatch::Nhits2, NPCTracks, Smatch::param1, Smatch::param2, RecordingBest(), vtmpSeg, Smatch::Z1, and Smatch::Z2.

Referenced by Exec().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ PrepareArraysToProcessEvent()

void BmnUpstreamTracking::PrepareArraysToProcessEvent ( )

private

Definition at line 1851 of file BmnUpstreamTracking.cxx.

References fBmnUpstreamTracksArray, i, kBig, kNumChambers, kNumPairs, kNumPars, kNumStSi, kPoin_Par, kPoints, NPCTracks, Nseg_Ch, NSiTracks, NSiXYhits, OutVector, par_ab_Ch, par_ab_SiTr, par_ab_tr, par_Seg_z, par_SiTr_z, Points, SiXYhits, vec_points, vec_twotracks, and vecUpTracks.

Referenced by Exec().

Here is the caller graph for this function:

◆ PrintAllTracks()

void BmnUpstreamTracking::PrintAllTracks ( vector< UpTracks > & vecUp )

private

Definition at line 1392 of file BmnUpstreamTracking.cxx.

References hAx_fitUp, hAy_fitUp, hchi2_fitUp, hNhitsUp, hResXst, hResYst, hTyTx_Up, hvertexXYUp, hx_fitUp, hy_fitUp, hy_vs_x_Up, kPoints, kZ_target, and Zcentr.

Referenced by Exec().

Here is the caller graph for this function:

◆ ReadMWPCSegments()

void BmnUpstreamTracking::ReadMWPCSegments	(	Double_t ***	par_ab,
		Double_t **	par_z,
		Int_t *	Nseg,
		vector< MC_points > &	vec )

private

Definition at line 476 of file BmnUpstreamTracking.cxx.

References expData, fBmnHitsArray, fDebug, fMWPCSegments, BmnMwpcHit::GetHitId(), BmnMwpcHit::GetMwpcId(), BmnTrack::GetNHits(), BmnTrack::GetParamFirst(), BmnMwpcHit::GetPdgId(), BmnMwpcHit::GetPlaneId(), hAx_upmc, hAy_upmc, hX_upmc, hY_upmc, i, kMaxMC, kNPlanes, kZ_target, X_shift, X_shift_seg, Y_shift, and Zcentr.

Referenced by Exec().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ReadMWPCTracks()

void BmnUpstreamTracking::ReadMWPCTracks	(	Double_t ***	par_ab,
		Double_t **	par_z,
		Int_t *	NTracks )

private

Definition at line 684 of file BmnUpstreamTracking.cxx.

References fDebug, fMWPCTracks, BmnTrack::GetNHits(), BmnTrack::GetParamFirst(), X_shift, and Y_shift.

Referenced by Exec().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ReadSiliconHits()

void BmnUpstreamTracking::ReadSiliconHits	(	Double_t ***	hits,
		Int_t *	NSihits )

private

Definition at line 449 of file BmnUpstreamTracking.cxx.

References fDebug, fSilHits, BmnHit::GetModule(), and BmnHit::GetStation().

Referenced by Exec().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ReadSiliconTracks()

void BmnUpstreamTracking::ReadSiliconTracks	(	Double_t **	par_ab,
		Double_t *	par_z,
		Int_t &	NTracks,
		vector< MC_points > &	vec )

private

Definition at line 334 of file BmnUpstreamTracking.cxx.

References expData, fDebug, fNstations, fSiTracks, fSiTracksSim, BmnHit::GetIndex(), BmnTrack::GetNHits(), BmnTrack::GetParamLast(), i, BmnUpstreamTracking::MC_points::Id, kMaxMC, BmnUpstreamTracking::MC_points::np_3si, BmnUpstreamTracking::MC_points::Np_Si, BmnUpstreamTracking::MC_points::x_si, BmnUpstreamTracking::MC_points::y_si, and BmnUpstreamTracking::MC_points::z_si.

Referenced by Exec().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ RecordCandidateSeg()

void BmnUpstreamTracking::RecordCandidateSeg	(	Int_t &	iSi,
		Int_t &	is2,
		Int_t &	is3,
		Double_t **	par_Si,
		Double_t ***	par_PC,
		Double_t *	par_Siz,
		Double_t **	ZPC,
		Double_t &	dx,
		Double_t &	dy,
		vector< Smatch > &	vtmp )