ATLFMiscMaker
class description - source file - inheritance tree
public:
ATLFMiscMaker ATLFMiscMaker()
ATLFMiscMaker ATLFMiscMaker(char* name, char* title)
ATLFMiscMaker ATLFMiscMaker(ATLFMiscMaker&)
void ~ATLFMiscMaker()
void Atanxy(Float_t xx, Float_t yy, Float_t& ang, Int_t& ierr)
void Atlut3(Int_t n, Float_t** p, Float_t& thr, Float_t& obl)
void Atplu3(Int_t indx, Float_t** plund, Float_t* pvec, char* option)
void Atrmx3(Float_t* vect, Float_t cp, Float_t sp, Float_t** rmx)
void Atrob3(Float_t** rmx, Float_t* vect, Float_t* rvec)
void Atrof3(Float_t** rmx, Float_t* vect, Float_t* rvec)
TH1F* CircEvent()
TH1F* CircJets()
void Circularity(Float_t& circjet, Float_t& circeve)
TClass* Class()
virtual void Clear(Option_t* option)
void EigenValue1(Int_t nm, Int_t n, Float_t** a, Float_t* d, Float_t* e, Float_t** z)
void EigenValue2(Int_t nm, Int_t n, Float_t* d, Float_t* e, Float_t** z, Int_t& ierr)
void EigenValues(Int_t nm, Int_t n, Float_t** ar, Float_t* wr, Float_t** zr, Int_t& ierr, Float_t* work)
virtual void Finish()
virtual void Init()
virtual TClass* IsA()
virtual Int_t Make()
virtual void MakeMisc()
virtual void MakeMissing()
Float_t MinET()
TH1F* Oblateness()
virtual void PrintInfo()
TH1F* PTmissPTnu()
TH1F* PTnu()
TH1F* PXmissPXnu()
TH1F* PXnu()
TH1F* PYmissPYnu()
TH1F* PYnu()
TH1F* RecPT()
TH1F* RecPTcells()
TH1F* RecPX()
TH1F* RecPXcells()
TH1F* RecPY()
TH1F* RecPYcells()
TH1F* ResolPX()
TH1F* ResolPX0()
TH1F* ResolPX100()
TH1F* ResolPX200()
TH1F* ResolPX25()
TH1F* ResolPX300()
TH1F* ResolPX50()
void SetMinET(Float_t val = 0.)
virtual void ShowMembers(TMemberInspector& insp, char* parent)
virtual void Streamer(TBuffer& b)
TH1F* Thrust()
void ThrustOblateness(Float_t& thrust, Float_t& obl)
private:
Float_t m_MinET Minimum E_T for unused cells
TH1F* m_CircJets circularity of jets
TH1F* m_CircEvent circularity of event
TH1F* m_Thrust thrust of event
TH1F* m_Oblateness oblateness of event
TH1F* m_RecPT reconstructed p_T
TH1F* m_RecPTcells reconstructed p_T +cells
TH1F* m_PTnu p_T nu
TH1F* m_PXnu pX nu
TH1F* m_PYnu pY nu
TH1F* m_RecPX reconstructed pX
TH1F* m_RecPY reconstructed pY
TH1F* m_RecPXcells reconstructed pX + cells
TH1F* m_RecPYcells reconstructed pY + cells
TH1F* m_ResolPX resol missing pX
TH1F* m_ResolPX0 resol missing pX (sumET=0-25)
TH1F* m_ResolPX25 resol missing pX (sumET=25-50)
TH1F* m_ResolPX50 resol missing pX (sumET=50-100)
TH1F* m_ResolPX100 resol missing pX (sumET=100-200)
TH1F* m_ResolPX200 resol missing pX (sumET=200-300)
TH1F* m_ResolPX300 resol missing pX (sumET=300-400)
TH1F* m_PTmissPTnu PTmiss/PTnu
TH1F* m_PXmissPXnu PXmiss/PXnu
TH1F* m_PYmissPYnu PYmiss/PYnu
ATLFast MiscMaker class.
ATLFMiscMaker()
ATLFMiscMaker(const char *name, const char *title)
:ATLFMaker(name,title)
~ATLFMiscMaker()
dummy
void Clear(Option_t *option)
Reset Misc Maker
void Init()
Create histograms
void Finish()
Function called by ATLFast::Finish at the end of the job
Int_t Make()
Compute Miscelleneous info
void MakeMisc()
.............................................
void MakeMissing()
.............................................
This function calculates the total missing transverse energy in the
event. The transverse energies of cells not used for cluster
reconstruction are smeared using the function ATLFClusterMaker::HadronResolution.
The total measured transverse momenta p_xobserv and p_yobserv are obtained by
summing over:
jets, b-jets, c-jets, isolated leptons and photons, non-isolated muons which
were not added to clusters, unused clusters, unused cells. The total missing
transverse momenta pxmiss, pymiss are calculated as :
pxmiss = - p_xobserv
pymiss = - p_yobserv
E_Tmiss = sqrt(pxmiss^2 + pymiss^2).
The sum of the momenta of all particles escaping detection is
also calculated. Results are stored in the object ATLFMisc, where
in addition to pxmiss and pymiss above, pxnu and pynu
denote the sum of momenta components of particles escaping detection
(muons outside detector acceptance, neutrinos and SUSY LSP).
.............................................
void PrintInfo()
void Circularity(Float_t &circj, Float_t &circev)
........................................................
...Purpose: to calculate sphericity of the event
Original routine SPHERIC from LUJET (and E. Richter-Was Spring 1996)
J Soderqvist 09/12/96
Rewrite to calculate Circularity of all jets (CIRCJ)
and of event (CIRCEV).
Use gen routine eisrs1.f to calculate eigenvalues
........................................................
void ThrustOblateness(Float_t &thrust, Float_t &obl)
........................................................
J Soderqvist 10/12/96
Function to calculate thrust and oblateness of events
in ATLFAST. Take clusters from ATLFAST and use
code that M. Pearce originally extracted from OPAL
to calculate thrust and oblateness.
(ATLUT3 is using JETSET thrust finding algorithm,
but is not dependent on the JETSET commons)
Modified 30/12/96 J.S.
to give values -777 if NCLU=0
The parameters to pass to ATLUT3 are:
function ATLUT3 (N,NRLUDM,P,THR,OBL)
INPUT : n Number of tracks
INPUT : nrludm First dimension of PLUND
INPUT : p 4-momenta in Jetset format (Px,Py,Pz,E)
OUTPUT : thr Thrust value
OUTPUT : obl Oblateness value
Use clusters identified by ATLFAST
........................................................
void EigenValues(Int_t nm, Int_t n, Float_t ar[2][2], Float_t *wr, Float_t zr[2][2],
Int_t &ierr, Float_t *work)
Compute all eigenvalues and corresponding eigenvectors of a real symmetric matrix.
translated to C++ from CERNLIB routine EISRS1
void EigenValue1(Int_t , Int_t n, Float_t a[2][2], Float_t *d, Float_t *e, Float_t z[2][2])
translated to C++ from CERNLIB routine TRED2
void EigenValue2(Int_t , Int_t n, Float_t *d, Float_t *e, Float_t z[2][2], Int_t &ierr)
translated to C++ from CERNLIB routine TQL2
void Atlut3(Int_t n,Float_t p[][5], Float_t &thr, Float_t &obl)
........................................................
. Extracted by M. Pearce from OPAL physics utilities (PX library)
. Originally authored by W. J. Gary
. This version converted for ATLAS use by
. M. Pearce and J. Soderqvist, KTH Stockholm, December 1996
.
. ------
. ATLUT3
. ------
. An "in-house" version of the Jetset thrust finding algorithm
. which works entirely through an argument list rather than
. with e.g. the Jetset common blocks. This routine calculates
. the standard linear thrust vectors and values. Its operation
. is entirely decoupled from any MST or MSTJ variables etc.
. which might be set by a user using Jetset.
. The main purpose of developing an in-house version of the
. Jetset thrust algorithm is so as to have a version
. which is compatible with both Jetset6.3 and Jetset7.1 etc.
. (because of the change in the Jetset common blocks between
. these two versions, the Jetset version of this thrust
. algorithm LUTHRU is version specific).
.
. The Jetset thrust algorithm implements an "approximate" method
. for thrust determination because not all particle combinations
. are tested. It is therefore logically possible that the thrust
. axes and values determined by this routine will correspond
. to local rather than to absolute maxima of the thrust function.
. However in practice this is unlikely because several starting
. points are used and the algorithm iterated to cross check one
. convergence vs. another. Thus this routine offers a very fast
. and in practice quite accurate algorithm for thrust (much faster
. than so-called "exact" algorithms).
. Usage :
.
.
. INPUT : ntrak Number of tracks
. INPUT : nrludm First dimension of PLUND
. INPUT : plund 4-momenta in Jetset format
. OUTPUT : thrust Thrust value
. OUTPUT : oblate Oblateness value
.
. CALLS : Atanxy, Atplu3, Atrmx3, Atrof3, Atrob3
. CALLED : Atlth4
.
. AUTHOR : Modified from LUTHRU (T.Sjostrand) by J.W.Gary (OPAL)
. CREATED : 31-Jan-89
. LAST MOD : Dec-96
.
. OPAL Modification Log.
. 04-Feb-89 In-house version for PX library J.W.Gary
. 12-Mar-89 Get rid of calls to RLU J.W.Gary
. 27-Nov-95 M.Schroder Clear part of the array P above tracks
. ATLAS Modification Log
. Dec-96 ATLAS naming of routines M. Pearce and J. Soderqvist
........................................................
void Atanxy(Float_t xx, Float_t yy, Float_t &ang, Int_t &ierr)
........................................................
. SOURCE: Jetset7.1 (T. Sjostrand)
. Reconstruct the azimuthal angle of a vector,
. given the X and Y components of a vector
. Usage :
.
.
. INPUT : xx The X component of a vector
. INPUT : yy The Y component of a vector
. OUTPUT : Ang The azimuthal angle
. OUTPUT : ierr = 0 if all is OK ; = -1 otherwise
........................................................
void Atplu3(Int_t indx, Float_t plund[][5], Float_t *pvec, char *option)
........................................................
. Extracted by M. Pearce from OPAL physics utilities (PX library)
. Originally authored by W. J. Gary
. This version converted for ATLAS use by
. M. Pearce and J. Soderqvist, KTH Stockholm, December 1996
.
. ------
. ATPLU3
. ------
. A utility routine to repack a Jetset 3-momentum as a
. standard ordered array or vice-versa. This routine
. is used to translate between the array conventions
. of the Jetset thrust algorithm and of the other routines
. in this library
. Usage :
.
.
. INPUT : indx The Jetset vector number
. IN/OUT : nrludm First argument of PLUND
. IN/OUT : plund The Jetset 3-momentum array
. IN/OUT : pvec The input or output array
. CONTROL : option ='U' means unpack Jetset array
. = anything else means pack Jetset array
.
. AUTHOR : J.W.Gary
. CREATED : 04-Feb-89
. LAST MOD : 04-Feb-89
.
. ATLAS Modification Log.
. Dec-96 ATLAS naming of routines M. Pearce and J. Soderqvist
........................................................
void Atrob3(Float_t rmx[3][3], Float_t *vect, Float_t *rvec)
........................................................
. SOURCE: HERWIG (B.Webber,G.Marchesini)
. Rotate 3-vector VECT by inverse of rotation matrix RMX,
. RVEC = (RMX)-1 * VECT
. Usage :
.
. Float_t vect (3.or.more),
. + rvec (3.or.more),
. + rmx (3,3.or.more)
.
.
. INPUT : rmx The rotation matrix
. INPUT : vect The vector to be rotated
. OUTPUT : rvec The rotated vector
........................................................
void Atrof3(Float_t rmx[3][3], Float_t *vect, Float_t *rvec)
........................................................
. SOURCE: HERWIG (B.Webber,G.Marchesini)
. Rotate 3-vector VECT by rotation matrix RMX,
. RVEC = RMX * VECT
. Usage :
.
. Float_t vect (3.or.more),
. + rvec (3.or.more),
. + rmx (3,3.or.more)
.
. CALL ATROF3 (RMX,VECT,RVEC)
.
. INPUT : rmx The rotation matrix
. INPUT : vect The vector to be rotated
. OUTPUT : rvec The rotated vector
........................................................
void Atrmx3(Float_t *vect, Float_t cp, Float_t sp, Float_t rmx[3][3])
........................................................
. SOURCE: HERWIG (B.Webber,G.Marchesini)
. Calculate the rotation matrix to get from vector VECT
. to the Z axis, followed by a rotation PHI about the Z axis,
. where CP, SP are the cosine and sine of PHI, respectively.
. Usage :
.
. INPUT : vect The vector for which the rotation matrix
. is to be calculated
. INPUT : cp Cosine of phi
. INPUT : sp Sine of phi
. OUTPUT : rmx The rotation matrix
........................................................
Inline Functions
Float_t MinET()
void SetMinET(Float_t val = 0.)
TH1F* CircJets()
TH1F* CircEvent()
TH1F* Thrust()
TH1F* Oblateness()
TH1F* RecPT()
TH1F* RecPTcells()
TH1F* PTnu()
TH1F* PXnu()
TH1F* PYnu()
TH1F* RecPX()
TH1F* RecPY()
TH1F* RecPXcells()
TH1F* RecPYcells()
TH1F* ResolPX()
TH1F* ResolPX0()
TH1F* ResolPX25()
TH1F* ResolPX50()
TH1F* ResolPX100()
TH1F* ResolPX200()
TH1F* ResolPX300()
TH1F* PTmissPTnu()
TH1F* PXmissPXnu()
TH1F* PYmissPYnu()
TClass* Class()
TClass* IsA()
void ShowMembers(TMemberInspector& insp, char* parent)
void Streamer(TBuffer& b)
ATLFMiscMaker ATLFMiscMaker(ATLFMiscMaker&)
ROOT page - Class index - Top of the page
This page has been automatically generated. If you have any comments or suggestions about the page layout send a mail to ROOT support, or contact the developers with any questions or problems regarding ROOT.