2 |
|
#define TRootSuperCluster_h |
3 |
|
|
4 |
|
/************************************************************************************************************************************ |
5 |
< |
SuperCluster Type: ijk |
6 |
< |
i = algo => 1=Island , 2=Hybrid, 3=multi5x5 |
7 |
< |
j = detector => 1=ECAL Barrel , 2=ECAL Endcap |
8 |
< |
k = corrections => 0=nopreshower / uncorrected , 1=no preshower / standard CMSSW corections |
9 |
< |
2=preshower / standard CMSSW corections , 3= preshower / uncorrected |
10 |
< |
|
11 |
< |
110 => Island Barrel SuperClusters |
12 |
< |
120 => Island Endcap SuperClusters |
13 |
< |
121 => Corrected Island Endcap SuperClusters (standard CMSSW corrections) |
14 |
< |
122 => Corrected Island Endcap SuperClusters with Preshower (standard CMSSW corrections) |
15 |
< |
210 => Hybrid Barrel SuperClusters |
16 |
< |
211 => Corrected Hybrid Barrel SuperClusters (standard CMSSW corrections) |
17 |
< |
320 => Multi5x5 Endcap SuperClusters |
18 |
< |
322 => Corrected Multi5x5 Endcap SuperClusters with Preshower (standard CMSSW corrections) |
19 |
< |
323 => Multi5x5 Endcap SuperClusters with Preshower |
20 |
< |
|
5 |
> |
SuperCluster Type: ijk |
6 |
> |
i = algo => 1=Island , 2=Hybrid, 3=multi5x5 |
7 |
> |
j = detector => 1=ECAL Barrel , 2=ECAL Endcap |
8 |
> |
k = corrections => 0=nopreshower / uncorrected , 1=no preshower / standard CMSSW corections |
9 |
> |
2=preshower / standard CMSSW corections , 3= preshower / uncorrected |
10 |
> |
|
11 |
> |
110 => Island Barrel SuperClusters |
12 |
> |
120 => Island Endcap SuperClusters |
13 |
> |
121 => Corrected Island Endcap SuperClusters (standard CMSSW corrections) |
14 |
> |
122 => Corrected Island Endcap SuperClusters with Preshower (standard CMSSW corrections) |
15 |
> |
210 => Hybrid Barrel SuperClusters |
16 |
> |
211 => Corrected Hybrid Barrel SuperClusters (standard CMSSW corrections) |
17 |
> |
320 => Multi5x5 Endcap SuperClusters |
18 |
> |
322 => Corrected Multi5x5 Endcap SuperClusters with Preshower (standard CMSSW corrections) |
19 |
> |
323 => Multi5x5 Endcap SuperClusters with Preshower |
20 |
> |
|
21 |
|
************************************************************************************************************************************/ |
22 |
|
|
23 |
|
#include <string> |
32 |
|
|
33 |
|
class TRootSuperCluster : public TRootCluster |
34 |
|
{ |
35 |
– |
|
36 |
– |
public: |
37 |
– |
|
38 |
– |
TRootSuperCluster() : |
39 |
– |
TRootCluster() |
40 |
– |
,nBasicClusters_(0) |
41 |
– |
,subClusterUID_() |
42 |
– |
,subClusterIndex_() |
43 |
– |
,seedClusterVectorIndex_(-1) |
44 |
– |
,preshowerEnergy_(0.) |
45 |
– |
,rawEnergy_(0.) |
46 |
– |
,photon_() |
47 |
– |
,electron_() |
48 |
– |
{;} |
49 |
– |
|
50 |
– |
TRootSuperCluster(const TRootSuperCluster& cluster) : |
51 |
– |
TRootCluster(cluster) |
52 |
– |
,nBasicClusters_(cluster.nBasicClusters_) |
53 |
– |
,subClusterUID_(cluster.subClusterUID_) |
54 |
– |
,subClusterIndex_(cluster.subClusterIndex_) |
55 |
– |
,seedClusterVectorIndex_(cluster.seedClusterVectorIndex_) |
56 |
– |
,preshowerEnergy_(cluster.preshowerEnergy_) |
57 |
– |
,rawEnergy_(cluster.rawEnergy_) |
58 |
– |
,photon_(cluster.photon_) |
59 |
– |
,electron_(cluster.photon_) |
60 |
– |
{;} |
61 |
– |
|
62 |
– |
TRootSuperCluster(Double_t energy, Double_t eta, Double_t phi) : |
63 |
– |
TRootCluster(energy, eta, phi) |
64 |
– |
,nBasicClusters_(0) |
65 |
– |
,subClusterUID_() |
66 |
– |
,subClusterIndex_() |
67 |
– |
,seedClusterVectorIndex_(-1) |
68 |
– |
,preshowerEnergy_(0.) |
69 |
– |
,rawEnergy_(0.) |
70 |
– |
,photon_() |
71 |
– |
,electron_() |
72 |
– |
{;} |
73 |
– |
|
74 |
– |
TRootSuperCluster(Double_t energy, Double_t eta, Double_t phi, Double_t x, Double_t y, Double_t z) : |
75 |
– |
TRootCluster(energy, eta, phi, x, y, z) |
76 |
– |
,nBasicClusters_(0) |
77 |
– |
,subClusterUID_() |
78 |
– |
,subClusterIndex_() |
79 |
– |
,seedClusterVectorIndex_(-1) |
80 |
– |
,preshowerEnergy_(0.) |
81 |
– |
,rawEnergy_(0.) |
82 |
– |
,photon_() |
83 |
– |
,electron_() |
84 |
– |
{;} |
85 |
– |
|
86 |
– |
TRootSuperCluster(Double_t energy, Double_t eta, Double_t phi, Double_t x, Double_t y, Double_t z, Int_t det) : |
87 |
– |
TRootCluster(energy, eta, phi, x, y, z, det) |
88 |
– |
,nBasicClusters_(0) |
89 |
– |
,subClusterUID_() |
90 |
– |
,subClusterIndex_() |
91 |
– |
,seedClusterVectorIndex_(-1) |
92 |
– |
,preshowerEnergy_(0.) |
93 |
– |
,rawEnergy_(0.) |
94 |
– |
,photon_() |
95 |
– |
,electron_() |
96 |
– |
{;} |
97 |
– |
|
98 |
– |
~TRootSuperCluster() {;} |
99 |
– |
|
100 |
– |
|
101 |
– |
Int_t nBasicClusters() const { return nBasicClusters_; } |
102 |
– |
std::vector<Int_t> subClusterUID() const { return subClusterUID_; } |
103 |
– |
std::vector<Int_t> subClusterIndex() const { return subClusterIndex_; } |
104 |
– |
Int_t seedIndex() const |
105 |
– |
{ |
106 |
– |
if( Int_t(subClusterIndex_.size())>seedClusterVectorIndex_) |
107 |
– |
{ |
108 |
– |
return subClusterIndex_.at(seedClusterVectorIndex_); |
109 |
– |
} |
110 |
– |
else |
111 |
– |
{ |
112 |
– |
cout << "No seed BasicCluster in this SuperCluster...." << endl; |
113 |
– |
return -1; |
114 |
– |
} |
115 |
– |
} |
116 |
– |
|
117 |
– |
Float_t preshowerEnergy() const { return preshowerEnergy_; } |
118 |
– |
Float_t rawEnergy() const { return rawEnergy_; } |
119 |
– |
TObject* photon() const { return photon_.GetObject();} |
120 |
– |
TObject* electron() const { return electron_.GetObject();} |
121 |
– |
Bool_t isPhoton() const { return ( photon_.GetObject() == 0 ? false : true); } |
122 |
– |
Bool_t isElectron() const { return ( electron_.GetObject() == 0 ? false : true); } |
123 |
– |
|
124 |
– |
void setNBasicClusters(Int_t nBasicClusters) { nBasicClusters_ = nBasicClusters; } |
125 |
– |
void addSubClusterUID(Int_t uid) { subClusterUID_.push_back(uid); } |
126 |
– |
void clearSubClusterUID() { subClusterUID_.clear(); } |
127 |
– |
void addSubClusterIndex(Int_t idx) { subClusterIndex_.push_back(idx); } |
128 |
– |
void clearSubClusterIndex() { subClusterIndex_.clear(); } |
129 |
– |
void setSeedClusterVectorIndex(Int_t seedClusterVectorIndex) { seedClusterVectorIndex_ = seedClusterVectorIndex; } |
130 |
– |
void setPreshowerEnergy(Float_t preshowerEnergy) { preshowerEnergy_ = preshowerEnergy; } |
131 |
– |
void setRawEnergy(Float_t rawEnergy) { rawEnergy_ = rawEnergy; } |
132 |
– |
void setPhoton(TObject* photon) { photon_ = photon; } |
133 |
– |
void setElectron(TObject* electron) { electron_ = electron; } |
134 |
– |
|
135 |
– |
|
136 |
– |
friend std::ostream& operator<< (std::ostream& stream, const TRootSuperCluster& clus) { |
137 |
– |
stream << "TRootSuperCluster - Type=" << clus.det_ << " (E,Et,eta,phi)=(" << clus.Mag() <<"," << clus.Pt() <<"," << clus.Eta() <<"," << clus.Phi() << ")" |
138 |
– |
<< " Calo position (x,y,z)=(" << clus.calX() << "," << clus.calY() << "," << clus.calZ() << ")" |
139 |
– |
<< " nBasicClusters=" << clus.nBasicClusters() << " e3x3=" << clus.e3x3() << " e5x5=" << clus.e5x5() << " eMax=" << clus.eMax()<< " e2nd=" << clus.e2nd() |
140 |
– |
<< " nXtals=" << clus.nXtals() << " preshowerEnergy=" << clus.preshowerEnergy() << " rawEnergy=" << clus.rawEnergy(); |
141 |
– |
return stream; |
142 |
– |
}; |
143 |
– |
|
144 |
– |
|
145 |
– |
|
146 |
– |
protected: |
147 |
– |
|
148 |
– |
Int_t nBasicClusters_; |
149 |
– |
std::vector<Int_t> subClusterUID_; |
150 |
– |
std::vector<Int_t> subClusterIndex_; |
151 |
– |
Int_t seedClusterVectorIndex_; // position in subClusterIndex_ vector containing the index of the seed basic cluster in the TRootBasicCluster array |
152 |
– |
Float_t preshowerEnergy_; |
153 |
– |
Float_t rawEnergy_; |
154 |
– |
TRef photon_; |
155 |
– |
TRef electron_; |
35 |
|
|
36 |
< |
|
37 |
< |
ClassDef (TRootSuperCluster,5); |
36 |
> |
public: |
37 |
> |
|
38 |
> |
TRootSuperCluster() : |
39 |
> |
TRootCluster() |
40 |
> |
,nBasicClusters_(0) |
41 |
> |
,subClusterUID_() |
42 |
> |
,subClusterIndex_() |
43 |
> |
,seedClusterVectorIndex_(-1) |
44 |
> |
,preshowerEnergy_(0.) |
45 |
> |
,rawEnergy_(0.) |
46 |
> |
,seedCluster_(0) |
47 |
> |
,photon_(0) |
48 |
> |
,electron_(0) |
49 |
> |
{;} |
50 |
> |
|
51 |
> |
TRootSuperCluster(const TRootSuperCluster& cluster) : |
52 |
> |
TRootCluster(cluster) |
53 |
> |
,nBasicClusters_(cluster.nBasicClusters_) |
54 |
> |
,subClusterUID_(cluster.subClusterUID_) |
55 |
> |
,subClusterIndex_(cluster.subClusterIndex_) |
56 |
> |
,seedClusterVectorIndex_(cluster.seedClusterVectorIndex_) |
57 |
> |
,preshowerEnergy_(cluster.preshowerEnergy_) |
58 |
> |
,rawEnergy_(cluster.rawEnergy_) |
59 |
> |
,seedCluster_(cluster.seedCluster_) |
60 |
> |
,photon_(cluster.photon_) |
61 |
> |
,electron_(cluster.photon_) |
62 |
> |
{;} |
63 |
> |
|
64 |
> |
TRootSuperCluster(Double_t energy, Double_t eta, Double_t phi) : |
65 |
> |
TRootCluster(energy, eta, phi) |
66 |
> |
,nBasicClusters_(0) |
67 |
> |
,subClusterUID_() |
68 |
> |
,subClusterIndex_() |
69 |
> |
,seedClusterVectorIndex_(-1) |
70 |
> |
,preshowerEnergy_(0.) |
71 |
> |
,rawEnergy_(0.) |
72 |
> |
,seedCluster_(0) |
73 |
> |
,photon_(0) |
74 |
> |
,electron_(0) |
75 |
> |
{;} |
76 |
> |
|
77 |
> |
TRootSuperCluster(Double_t energy, Double_t eta, Double_t phi, Double_t x, Double_t y, Double_t z) : |
78 |
> |
TRootCluster(energy, eta, phi, x, y, z) |
79 |
> |
,nBasicClusters_(0) |
80 |
> |
,subClusterUID_() |
81 |
> |
,subClusterIndex_() |
82 |
> |
,seedClusterVectorIndex_(-1) |
83 |
> |
,preshowerEnergy_(0.) |
84 |
> |
,rawEnergy_(0.) |
85 |
> |
,seedCluster_(0) |
86 |
> |
,photon_(0) |
87 |
> |
,electron_(0) |
88 |
> |
{;} |
89 |
> |
|
90 |
> |
TRootSuperCluster(Double_t energy, Double_t eta, Double_t phi, Double_t x, Double_t y, Double_t z, Int_t det) : |
91 |
> |
TRootCluster(energy, eta, phi, x, y, z, det) |
92 |
> |
,nBasicClusters_(0) |
93 |
> |
,subClusterUID_() |
94 |
> |
,subClusterIndex_() |
95 |
> |
,seedClusterVectorIndex_(-1) |
96 |
> |
,preshowerEnergy_(0.) |
97 |
> |
,rawEnergy_(0.) |
98 |
> |
,seedCluster_(0) |
99 |
> |
,photon_(0) |
100 |
> |
,electron_(0) |
101 |
> |
{;} |
102 |
> |
|
103 |
> |
~TRootSuperCluster() {;} |
104 |
> |
|
105 |
> |
|
106 |
> |
Int_t nBasicClusters() const { return nBasicClusters_; } |
107 |
> |
std::vector<Int_t> subClusterUID() const { return subClusterUID_; } |
108 |
> |
std::vector<Int_t> subClusterIndex() const { return subClusterIndex_; } |
109 |
> |
Int_t seedIndex() const |
110 |
> |
{ |
111 |
> |
if( Int_t(subClusterIndex_.size())>seedClusterVectorIndex_) |
112 |
> |
{ |
113 |
> |
return subClusterIndex_.at(seedClusterVectorIndex_); |
114 |
> |
} |
115 |
> |
else |
116 |
> |
{ |
117 |
> |
cout << "No seed BasicCluster in this SuperCluster...." << endl; |
118 |
> |
return -1; |
119 |
> |
} |
120 |
> |
} |
121 |
> |
|
122 |
> |
Float_t preshowerEnergy() const { return preshowerEnergy_; } |
123 |
> |
Float_t rawEnergy() const { return rawEnergy_; } |
124 |
> |
TObject* seedCluster() const { return seedCluster_.GetObject();} |
125 |
> |
TObject* photon() const { return photon_.GetObject();} |
126 |
> |
TObject* electron() const { return electron_.GetObject();} |
127 |
> |
Bool_t isPhoton() const { return ( photon_.GetObject() == 0 ? false : true); } |
128 |
> |
Bool_t isElectron() const { return ( electron_.GetObject() == 0 ? false : true); } |
129 |
> |
|
130 |
> |
void setNBasicClusters(Int_t nBasicClusters) { nBasicClusters_ = nBasicClusters; } |
131 |
> |
void addSubClusterUID(Int_t uid) { subClusterUID_.push_back(uid); } |
132 |
> |
void clearSubClusterUID() { subClusterUID_.clear(); } |
133 |
> |
void addSubClusterIndex(Int_t idx) { subClusterIndex_.push_back(idx); } |
134 |
> |
void clearSubClusterIndex() { subClusterIndex_.clear(); } |
135 |
> |
void setSeedClusterVectorIndex(Int_t seedClusterVectorIndex) { seedClusterVectorIndex_ = seedClusterVectorIndex; } |
136 |
> |
void setPreshowerEnergy(Float_t preshowerEnergy) { preshowerEnergy_ = preshowerEnergy; } |
137 |
> |
void setRawEnergy(Float_t rawEnergy) { rawEnergy_ = rawEnergy; } |
138 |
> |
void setSeedCluster(TObject* seedCluster) { seedCluster_ = seedCluster; } |
139 |
> |
void setPhoton(TObject* photon) { photon_ = photon; } |
140 |
> |
void setElectron(TObject* electron) { electron_ = electron; } |
141 |
> |
|
142 |
> |
|
143 |
> |
friend std::ostream& operator<< (std::ostream& stream, const TRootSuperCluster& clus) { |
144 |
> |
stream << "TRootSuperCluster - Type=" << clus.det_ << " (E,Et,eta,phi)=(" << clus.Mag() <<"," << clus.Pt() <<"," << clus.Eta() <<"," << clus.Phi() << ")" |
145 |
> |
<< " Calo position (x,y,z)=(" << clus.calX() << "," << clus.calY() << "," << clus.calZ() << ")" |
146 |
> |
<< " nBasicClusters=" << clus.nBasicClusters() << " e3x3=" << clus.e3x3() << " e5x5=" << clus.e5x5() << " eMax=" << clus.eMax()<< " e2nd=" << clus.e2nd() |
147 |
> |
<< " nXtals=" << clus.nXtals() << " preshowerEnergy=" << clus.preshowerEnergy() << " rawEnergy=" << clus.rawEnergy(); |
148 |
> |
return stream; |
149 |
> |
}; |
150 |
> |
|
151 |
> |
|
152 |
> |
protected: |
153 |
> |
|
154 |
> |
Int_t nBasicClusters_; |
155 |
> |
std::vector<Int_t> subClusterUID_; |
156 |
> |
std::vector<Int_t> subClusterIndex_; |
157 |
> |
Int_t seedClusterVectorIndex_; // position in subClusterIndex_ vector containing the index of the seed basic cluster in the TRootBasicCluster array |
158 |
> |
Float_t preshowerEnergy_; |
159 |
> |
Float_t rawEnergy_; |
160 |
> |
TRef seedCluster_; |
161 |
> |
TRef photon_; |
162 |
> |
TRef electron_; |
163 |
> |
|
164 |
> |
ClassDef (TRootSuperCluster,5); |
165 |
> |
|
166 |
|
}; |
167 |
|
|
168 |
|
#endif |