BOSS 7.0.4
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
EmcRecShowerEnergy Class Reference

#include <EmcRecShowerEnergy.h>

Public Member Functions

void Energy (RecEmcShower &aShower)
 
void E3x3 (RecEmcShower &aShower)
 
void E5x5 (RecEmcShower &aShower)
 
void ELepton (RecEmcShower &aShower)
 
void EAll (RecEmcShower &aShower)
 
RecEmcEnergy ECorrection (const RecEmcEnergy eIn)
 
RecEmcEnergy ECorrTheta (const RecEmcEnergy eIn, const RecEmcID &id)
 
void Energy (RecEmcShower &aShower)
 
void E3x3 (RecEmcShower &aShower)
 
void E5x5 (RecEmcShower &aShower)
 
void ELepton (RecEmcShower &aShower)
 
void EAll (RecEmcShower &aShower)
 
RecEmcEnergy ECorrection (const RecEmcEnergy eIn)
 
RecEmcEnergy ECorrTheta (const RecEmcEnergy eIn, const RecEmcID &id)
 

Detailed Description

Definition at line 11 of file InstallArea/include/EmcRec/EmcRec/EmcRecShowerEnergy.h.

Member Function Documentation

◆ E3x3() [1/2]

void EmcRecShowerEnergy::E3x3 ( RecEmcShower &  aShower)

◆ E3x3() [2/2]

void EmcRecShowerEnergy::E3x3 ( RecEmcShower &  aShower)

◆ E5x5() [1/2]

void EmcRecShowerEnergy::E5x5 ( RecEmcShower &  aShower)

◆ E5x5() [2/2]

void EmcRecShowerEnergy::E5x5 ( RecEmcShower &  aShower)

◆ EAll() [1/2]

void EmcRecShowerEnergy::EAll ( RecEmcShower &  aShower)

◆ EAll() [2/2]

void EmcRecShowerEnergy::EAll ( RecEmcShower &  aShower)

◆ ECorrection() [1/2]

RecEmcEnergy EmcRecShowerEnergy::ECorrection ( const RecEmcEnergy  eIn)

Definition at line 125 of file EmcRecShowerEnergy.cxx.

126{
127 if(eIn>3.) return eIn;
128
129 EmcRecParameter& Para=EmcRecParameter::GetInstance();
130
131 RecEmcEnergy eOut=0;
132 double par[4];
133 for(int i=0;i<4;i++) {
134 par[i]=Para.ECorr(i);
135 }
136
137 eOut = eIn/(par[0]+par[1]*eIn+par[2]*eIn*eIn+par[3]*eIn*eIn*eIn);
138 return eOut;
139}
EmcRecParameter::ECorr
double ECorr(int n) const
Definition: EmcRecParameter.cxx:619
EmcRecParameter::GetInstance
static EmcRecParameter & GetInstance()
Definition: EmcRecParameter.cxx:537

◆ ECorrection() [2/2]

RecEmcEnergy EmcRecShowerEnergy::ECorrection ( const RecEmcEnergy  eIn)

◆ ECorrTheta() [1/2]

RecEmcEnergy EmcRecShowerEnergy::ECorrTheta ( const RecEmcEnergy  eIn,
const RecEmcID &  id 
)

Definition at line 141 of file EmcRecShowerEnergy.cxx.

142{
143 EmcRecParameter& Para=EmcRecParameter::GetInstance();
144 RecEmcEnergy eOut=eIn;
145
146 unsigned int npart = EmcID::barrel_ec(id);
147 unsigned int ntheta = EmcID::theta_module(id);
148
149 if(npart==1) {
150 eOut *= 1.843/Para.Peak(ntheta);
151 } else if(npart==0) {
152 eOut *= 1.843/Para.Peak(ntheta+44);
153 } else if(npart==2) {
154 eOut *= 1.843/Para.Peak(ntheta+50);
155 }
156
157 return eOut;
158}
EmcID::barrel_ec
static unsigned int barrel_ec(const Identifier &id)
Values of different levels (failure returns 0)
Definition: EmcID.cxx:38
EmcID::theta_module
static unsigned int theta_module(const Identifier &id)
Definition: EmcID.cxx:43
EmcRecParameter::Peak
double Peak(int n) const
Definition: EmcRecParameter.cxx:654

◆ ECorrTheta() [2/2]

RecEmcEnergy EmcRecShowerEnergy::ECorrTheta ( const RecEmcEnergy  eIn,
const RecEmcID &  id 
)

◆ ELepton() [1/2]

void EmcRecShowerEnergy::ELepton ( RecEmcShower &  aShower)

◆ ELepton() [2/2]

void EmcRecShowerEnergy::ELepton ( RecEmcShower &  aShower)

◆ Energy() [1/2]

void EmcRecShowerEnergy::Energy ( RecEmcShower &  aShower)

Definition at line 6 of file EmcRecShowerEnergy.cxx.

7{
8 RecEmcFractionMap::const_iterator cit;
9 RecEmcEnergy e1=0;
10 RecEmcEnergy e9=0;
11 RecEmcEnergy e25=0;
12 RecEmcEnergy elepton=0;
13 RecEmcEnergy eall=0;
14
15 EmcRecNeighbor nhb;
16
17 RecEmcID CellId=aShower.getShowerId();
18 int module=EmcID::barrel_ec(CellId);
19 RecEmcIDVector NearCell=nhb.GetNeighbors(CellId);
20 RecEmcIDVector NextNearCell=nhb.GetNextNeighbors(CellId);
21 RecEmcIDVector tmpNearCell;
22 RecEmcIDVector tmpNextNearCell;
23 i_RecEmcIDVector pNearCell;
24 i_RecEmcIDVector pNextNearCell;
25 vector<RecEmcEnergy> eVec;
26 vector<RecEmcEnergy>::const_iterator ciVec;
27
28 tmpNearCell.push_back(CellId);
29 tmpNextNearCell.push_back(CellId);
30
31 cit=aShower.Find(CellId);
32 //int time_seed = cit->second.getTime();
33 e1=(cit->second.getEnergy())*(cit->second.getFraction());
34 e9+=(cit->second.getEnergy())*(cit->second.getFraction());
35 e25+=(cit->second.getEnergy())*(cit->second.getFraction());
36
37 //e3x3
38 for(pNearCell=NearCell.begin();
39 pNearCell!=NearCell.end();
40 pNearCell++) {
41 cit=aShower.Find(*pNearCell);
42 if(cit!=aShower.End()) {
43 tmpNearCell.push_back(*pNearCell);
44 tmpNextNearCell.push_back(*pNearCell);
45 e9+=cit->second.getEnergy()*cit->second.getFraction();
46 e25+=cit->second.getEnergy()*cit->second.getFraction();
47 }
48 }
49
50 //e5x5
51 for(pNextNearCell=NextNearCell.begin();
52 pNextNearCell!=NextNearCell.end();
53 pNextNearCell++) {
54 cit=aShower.Find(*pNextNearCell);
55 if(cit!=aShower.End()) {
56 tmpNextNearCell.push_back(*pNextNearCell);
57 e25+=cit->second.getEnergy()*cit->second.getFraction();
58 }
59 }
60
61 //eall
62 for(cit=aShower.Begin();cit!=aShower.End();++cit) {
63 eall+=(cit->second.getEnergy())*(cit->second.getFraction());
64 eVec.push_back(cit->second.getEnergy()*cit->second.getFraction());
65 }
66
67 //calculate number of hits from MC
68 int nHit,n;
69 EmcRecParameter& Para=EmcRecParameter::GetInstance();
70 nHit=(int)(Para.HitNb(0)*log(Para.HitNb(1)*e9+Para.HitNb(2)));
71 n=0;
72
73 //sort by energy
74 sort(eVec.begin(), eVec.end(), greater<RecEmcEnergy>());
75
76 for(ciVec=eVec.begin();ciVec!=eVec.end();ciVec++) {
77 if(n<nHit) {
78 elepton+=*ciVec;
79 n++;
80 }
81 }
82
83 //energy correction
84 //RecEmcEnergy eCorr=ECorrTheta(e25,CellId);
85 //RecEmcEnergy eCorr=ECorrection(e25);
86 int getthetaid = EmcID::theta_module(CellId);
87 int getmodule = EmcID::barrel_ec(CellId);
88 if(getthetaid>21)getthetaid=43-getthetaid;
89 if(getmodule==1)getthetaid=getthetaid+6;
90 double dthetaid=double(getthetaid);
91 double eCorr = Para.ECorrMC(e25,dthetaid);
92
93 //energy error
94 RecEmcEnergy de,de1,de2,de3;
95 de1 = Para.SigE(0)/eCorr;
96 de2 = Para.SigE(1)/pow(eCorr,0.25);
97 de3 = Para.SigE(2);
98 de = sqrt(de1*de1+de2*de2+de3*de3)*perCent*eCorr;
99
100 double err = Para.ErrMC(e25,dthetaid);
101 if(err>0) de = err*e25;
102
103 aShower.setTrackId(-1);
104 aShower.setCellId(CellId);
105 aShower.setModule(module);
106 aShower.setNumHits(aShower.getSize());
107 aShower.setDE(de);
108 aShower.CellId3x3(tmpNearCell);
109 aShower.CellId5x5(tmpNextNearCell);
110 aShower.setESeed(e1);
111 aShower.setE3x3(e9);
112 aShower.setE5x5(e25);
113 aShower.ELepton(elepton);
114 aShower.EAll(eall);
115 aShower.setEnergy(eCorr);
116
117 //cout<<"e1="<<aShower.eSeed()
118 // <<"\te9="<<aShower.e3x3()
119 // <<"\te25="<<aShower.e5x5()
120 // <<"\telepton="<<aShower.getELepton()
121 // <<"\teall="<<aShower.getEAll()
122 // <<"\tenergy="<<aShower.energy()<<endl;
123}
n
const Int_t n
Definition: DataBase/tau_mode.c:65
e1
Double_t e1
Definition: DataBase/tau_mode.c:7
EmcRecNeighbor::GetNeighbors
RecEmcIDVector GetNeighbors(const Identifier &id)
Definition: EmcRecNeighbor.cxx:54
EmcRecNeighbor::GetNextNeighbors
RecEmcIDVector GetNextNeighbors(const Identifier &id)
Definition: EmcRecNeighbor.cxx:319
EmcRecParameter::SigE
double SigE(int n) const
Definition: EmcRecParameter.cxx:624
EmcRecParameter::ECorrMC
double ECorrMC(double eg, double theid) const
Definition: EmcRecParameter.cxx:839
EmcRecParameter::HitNb
double HitNb(int n) const
Definition: EmcRecParameter.cxx:639
EmcRecParameter::ErrMC
double ErrMC(double eg, double theid) const
Definition: EmcRecParameter.cxx:856
RecEmcShower::End
RecEmcFractionMap::const_iterator End() const
Definition: RecEmcShower.cxx:136
RecEmcShower::CellId3x3
void CellId3x3(RecEmcIDVector &id3x3)
Definition: RecEmcShower.cxx:119
RecEmcShower::EAll
RecEmcEnergy EAll(RecEmcEnergy e)
Definition: RecEmcShower.cxx:278
RecEmcShower::ELepton
RecEmcEnergy ELepton(RecEmcEnergy e)
Definition: RecEmcShower.cxx:282
RecEmcShower::Begin
RecEmcFractionMap::const_iterator Begin() const
Definition: RecEmcShower.cxx:131
RecEmcShower::Find
RecEmcFractionMap::const_iterator Find(const RecEmcID &CellId) const
Definition: RecEmcShower.cxx:141
RecEmcShower::getSize
unsigned int getSize() const
Definition: RecEmcShower.cxx:209
RecEmcShower::CellId5x5
void CellId5x5(RecEmcIDVector &id5x5)
Definition: RecEmcShower.cxx:125

Referenced by EmcRecSplitWeighted::Split().

◆ Energy() [2/2]

void EmcRecShowerEnergy::Energy ( RecEmcShower &  aShower)
The documentation for this class was generated from the following files: