CMS 3D CMS Logo

PFECALSuperClusterProducer.cc
Go to the documentation of this file.
2 
3 #include <memory>
4 
6 
8 
10 
15 
19 
22 
25 
28 
30 
31 #include "TVector2.h"
33 
34 using namespace std;
35 using namespace edm;
36 
37 namespace {
38  const std::string ClusterType__BOX("Box");
39  const std::string ClusterType__Mustache("Mustache");
40 
41  const std::string EnergyWeight__Raw("Raw");
42  const std::string EnergyWeight__CalibratedNoPS("CalibratedNoPS");
43  const std::string EnergyWeight__CalibratedTotal("CalibratedTotal");
44 } // namespace
45 
47  verbose_ = iConfig.getUntrackedParameter<bool>("verbose", false);
48 
49  superClusterAlgo_.setUseRegression(iConfig.getParameter<bool>("useRegression"));
50 
51  isOOTCollection_ = iConfig.getParameter<bool>("isOOTCollection");
52  superClusterAlgo_.setIsOOTCollection(isOOTCollection_);
53 
54  std::string _typename = iConfig.getParameter<std::string>("ClusteringType");
55  if (_typename == ClusterType__BOX) {
56  _theclusteringtype = PFECALSuperClusterAlgo::kBOX;
57  } else if (_typename == ClusterType__Mustache) {
58  _theclusteringtype = PFECALSuperClusterAlgo::kMustache;
59  } else {
60  throw cms::Exception("InvalidClusteringType") << "You have not chosen a valid clustering type,"
61  << " please choose from \"Box\" or \"Mustache\"!";
62  }
63  superClusterAlgo_.setClusteringType(_theclusteringtype);
64  superClusterAlgo_.setUseDynamicDPhi(iConfig.getParameter<bool>("useDynamicDPhiWindow"));
65  // clusteringType and useDynamicDPhi need to be defined before setting the tokens in order to esConsume only the necessary records
66  superClusterAlgo_.setTokens(iConfig, consumesCollector());
67 
68  std::string _weightname = iConfig.getParameter<std::string>("EnergyWeight");
69  if (_weightname == EnergyWeight__Raw) {
70  _theenergyweight = PFECALSuperClusterAlgo::kRaw;
71  } else if (_weightname == EnergyWeight__CalibratedNoPS) {
72  _theenergyweight = PFECALSuperClusterAlgo::kCalibratedNoPS;
73  } else if (_weightname == EnergyWeight__CalibratedTotal) {
75  } else {
76  throw cms::Exception("InvalidClusteringType") << "You have not chosen a valid energy weighting scheme,"
77  << " please choose from \"Raw\", \"CalibratedNoPS\", or"
78  << " \"CalibratedTotal\"!";
79  }
80 
81  // parameters for clustering
82  bool seedThresholdIsET = iConfig.getParameter<bool>("seedThresholdIsET");
83 
84  double threshPFClusterSeedBarrel = iConfig.getParameter<double>("thresh_PFClusterSeedBarrel");
85  double threshPFClusterBarrel = iConfig.getParameter<double>("thresh_PFClusterBarrel");
86 
87  double threshPFClusterSeedEndcap = iConfig.getParameter<double>("thresh_PFClusterSeedEndcap");
88  double threshPFClusterEndcap = iConfig.getParameter<double>("thresh_PFClusterEndcap");
89 
90  double phiwidthSuperClusterBarrel = iConfig.getParameter<double>("phiwidth_SuperClusterBarrel");
91  double etawidthSuperClusterBarrel = iConfig.getParameter<double>("etawidth_SuperClusterBarrel");
92 
93  double phiwidthSuperClusterEndcap = iConfig.getParameter<double>("phiwidth_SuperClusterEndcap");
94  double etawidthSuperClusterEndcap = iConfig.getParameter<double>("etawidth_SuperClusterEndcap");
95 
96  //double threshPFClusterMustacheOutBarrel = iConfig.getParameter<double>("thresh_PFClusterMustacheOutBarrel");
97  //double threshPFClusterMustacheOutEndcap = iConfig.getParameter<double>("thresh_PFClusterMustacheOutEndcap");
98 
99  double doSatelliteClusterMerge = iConfig.getParameter<bool>("doSatelliteClusterMerge");
100  double satelliteClusterSeedThreshold = iConfig.getParameter<double>("satelliteClusterSeedThreshold");
101  double satelliteMajorityFraction = iConfig.getParameter<double>("satelliteMajorityFraction");
102  bool dropUnseedable = iConfig.getParameter<bool>("dropUnseedable");
103 
104  superClusterAlgo_.setVerbosityLevel(verbose_);
105  superClusterAlgo_.setEnergyWeighting(_theenergyweight);
106  superClusterAlgo_.setUseETForSeeding(seedThresholdIsET);
107 
108  superClusterAlgo_.setThreshSuperClusterEt(iConfig.getParameter<double>("thresh_SCEt"));
109 
110  superClusterAlgo_.setThreshPFClusterSeedBarrel(threshPFClusterSeedBarrel);
111  superClusterAlgo_.setThreshPFClusterBarrel(threshPFClusterBarrel);
112 
113  superClusterAlgo_.setThreshPFClusterSeedEndcap(threshPFClusterSeedEndcap);
114  superClusterAlgo_.setThreshPFClusterEndcap(threshPFClusterEndcap);
115 
116  superClusterAlgo_.setPhiwidthSuperClusterBarrel(phiwidthSuperClusterBarrel);
117  superClusterAlgo_.setEtawidthSuperClusterBarrel(etawidthSuperClusterBarrel);
118 
119  superClusterAlgo_.setPhiwidthSuperClusterEndcap(phiwidthSuperClusterEndcap);
120  superClusterAlgo_.setEtawidthSuperClusterEndcap(etawidthSuperClusterEndcap);
121 
122  superClusterAlgo_.setSatelliteMerging(doSatelliteClusterMerge);
123  superClusterAlgo_.setSatelliteThreshold(satelliteClusterSeedThreshold);
124  superClusterAlgo_.setMajorityFraction(satelliteMajorityFraction);
125  superClusterAlgo_.setDropUnseedable(dropUnseedable);
126  //superClusterAlgo_.setThreshPFClusterMustacheOutBarrel( threshPFClusterMustacheOutBarrel );
127  //superClusterAlgo_.setThreshPFClusterMustacheOutEndcap( threshPFClusterMustacheOutEndcap );
128 
129  //Load the ECAL energy calibration
130  thePFEnergyCalibration_ = std::make_shared<PFEnergyCalibration>();
131  superClusterAlgo_.setPFClusterCalibration(thePFEnergyCalibration_);
132 
133  bool applyCrackCorrections_ = iConfig.getParameter<bool>("applyCrackCorrections");
134  superClusterAlgo_.setCrackCorrections(applyCrackCorrections_);
135 
136  PFBasicClusterCollectionBarrel_ = iConfig.getParameter<string>("PFBasicClusterCollectionBarrel");
137  PFSuperClusterCollectionBarrel_ = iConfig.getParameter<string>("PFSuperClusterCollectionBarrel");
138 
139  PFBasicClusterCollectionEndcap_ = iConfig.getParameter<string>("PFBasicClusterCollectionEndcap");
140  PFSuperClusterCollectionEndcap_ = iConfig.getParameter<string>("PFSuperClusterCollectionEndcap");
141 
142  PFBasicClusterCollectionPreshower_ = iConfig.getParameter<string>("PFBasicClusterCollectionPreshower");
143  PFSuperClusterCollectionEndcapWithPreshower_ =
144  iConfig.getParameter<string>("PFSuperClusterCollectionEndcapWithPreshower");
145 
146  PFClusterAssociationEBEE_ = "PFClusterAssociationEBEE";
147  PFClusterAssociationES_ = "PFClusterAssociationES";
148 
149  produces<reco::SuperClusterCollection>(PFSuperClusterCollectionBarrel_);
150  produces<reco::SuperClusterCollection>(PFSuperClusterCollectionEndcapWithPreshower_);
151  produces<reco::CaloClusterCollection>(PFBasicClusterCollectionBarrel_);
152  produces<reco::CaloClusterCollection>(PFBasicClusterCollectionEndcap_);
153  produces<reco::CaloClusterCollection>(PFBasicClusterCollectionPreshower_);
154  produces<edm::ValueMap<reco::CaloClusterPtr>>(PFClusterAssociationEBEE_);
155  produces<edm::ValueMap<reco::CaloClusterPtr>>(PFClusterAssociationES_);
156 }
157 
159 
161  superClusterAlgo_.update(iE);
162 }
163 
165  // update SC parameters
166  superClusterAlgo_.updateSCParams(iSetup);
167  // do clustering
168  superClusterAlgo_.loadAndSortPFClusters(iEvent);
169  superClusterAlgo_.run();
170 
171  //build collections of output CaloClusters from the used PFClusters
172  auto caloClustersEB = std::make_unique<reco::CaloClusterCollection>();
173  auto caloClustersEE = std::make_unique<reco::CaloClusterCollection>();
174  auto caloClustersES = std::make_unique<reco::CaloClusterCollection>();
175 
176  std::map<reco::CaloClusterPtr, unsigned int> pfClusterMapEB; //maps of pfclusters to caloclusters
177  std::map<reco::CaloClusterPtr, unsigned int> pfClusterMapEE;
178  std::map<reco::CaloClusterPtr, unsigned int> pfClusterMapES;
179 
180  //fill calocluster collections and maps
181  for (const auto& ebsc : *(superClusterAlgo_.getEBOutputSCCollection())) {
182  for (reco::CaloCluster_iterator pfclus = ebsc.clustersBegin(); pfclus != ebsc.clustersEnd(); ++pfclus) {
183  if (!pfClusterMapEB.count(*pfclus)) {
184  reco::CaloCluster caloclus(**pfclus);
185  caloClustersEB->push_back(caloclus);
186  pfClusterMapEB[*pfclus] = caloClustersEB->size() - 1;
187  } else {
188  throw cms::Exception("PFECALSuperClusterProducer::produce")
189  << "Found an EB pfcluster matched to more than one EB supercluster!" << std::dec << std::endl;
190  }
191  }
192  }
193  for (const auto& eesc : *(superClusterAlgo_.getEEOutputSCCollection())) {
194  for (reco::CaloCluster_iterator pfclus = eesc.clustersBegin(); pfclus != eesc.clustersEnd(); ++pfclus) {
195  if (!pfClusterMapEE.count(*pfclus)) {
196  reco::CaloCluster caloclus(**pfclus);
197  caloClustersEE->push_back(caloclus);
198  pfClusterMapEE[*pfclus] = caloClustersEE->size() - 1;
199  } else {
200  throw cms::Exception("PFECALSuperClusterProducer::produce")
201  << "Found an EE pfcluster matched to more than one EE supercluster!" << std::dec << std::endl;
202  }
203  }
204  for (reco::CaloCluster_iterator pfclus = eesc.preshowerClustersBegin(); pfclus != eesc.preshowerClustersEnd();
205  ++pfclus) {
206  if (!pfClusterMapES.count(*pfclus)) {
207  reco::CaloCluster caloclus(**pfclus);
208  caloClustersES->push_back(caloclus);
209  pfClusterMapES[*pfclus] = caloClustersES->size() - 1;
210  } else {
211  throw cms::Exception("PFECALSuperClusterProducer::produce")
212  << "Found an ES pfcluster matched to more than one EE supercluster!" << std::dec << std::endl;
213  }
214  }
215  }
216 
217  //create ValueMaps from output CaloClusters back to original PFClusters
218  auto pfClusterAssociationEBEE = std::make_unique<edm::ValueMap<reco::CaloClusterPtr>>();
219  auto pfClusterAssociationES = std::make_unique<edm::ValueMap<reco::CaloClusterPtr>>();
220 
221  //vectors to fill ValueMaps
222  std::vector<reco::CaloClusterPtr> clusptrsEB(caloClustersEB->size());
223  std::vector<reco::CaloClusterPtr> clusptrsEE(caloClustersEE->size());
224  std::vector<reco::CaloClusterPtr> clusptrsES(caloClustersES->size());
225 
226  //put calocluster output collections in event and get orphan handles to create ptrs
227  const edm::OrphanHandle<reco::CaloClusterCollection>& caloClusHandleEB =
228  iEvent.put(std::move(caloClustersEB), PFBasicClusterCollectionBarrel_);
229  const edm::OrphanHandle<reco::CaloClusterCollection>& caloClusHandleEE =
230  iEvent.put(std::move(caloClustersEE), PFBasicClusterCollectionEndcap_);
231  const edm::OrphanHandle<reco::CaloClusterCollection>& caloClusHandleES =
232  iEvent.put(std::move(caloClustersES), PFBasicClusterCollectionPreshower_);
233 
234  //relink superclusters to output caloclusters and fill vectors for ValueMaps
235  for (auto& ebsc : *(superClusterAlgo_.getEBOutputSCCollection())) {
236  reco::CaloClusterPtr seedptr(caloClusHandleEB, pfClusterMapEB[ebsc.seed()]);
237  ebsc.setSeed(seedptr);
238 
240  for (reco::CaloCluster_iterator pfclus = ebsc.clustersBegin(); pfclus != ebsc.clustersEnd(); ++pfclus) {
241  int caloclusidx = pfClusterMapEB[*pfclus];
242  reco::CaloClusterPtr clusptr(caloClusHandleEB, caloclusidx);
243  clusters.push_back(clusptr);
244  clusptrsEB[caloclusidx] = *pfclus;
245  }
246  ebsc.setClusters(clusters);
247  }
248  for (auto& eesc : *(superClusterAlgo_.getEEOutputSCCollection())) {
249  reco::CaloClusterPtr seedptr(caloClusHandleEE, pfClusterMapEE[eesc.seed()]);
250  eesc.setSeed(seedptr);
251 
253  for (reco::CaloCluster_iterator pfclus = eesc.clustersBegin(); pfclus != eesc.clustersEnd(); ++pfclus) {
254  int caloclusidx = pfClusterMapEE[*pfclus];
255  reco::CaloClusterPtr clusptr(caloClusHandleEE, caloclusidx);
256  clusters.push_back(clusptr);
257  clusptrsEE[caloclusidx] = *pfclus;
258  }
259  eesc.setClusters(clusters);
260 
261  reco::CaloClusterPtrVector psclusters;
262  for (reco::CaloCluster_iterator pfclus = eesc.preshowerClustersBegin(); pfclus != eesc.preshowerClustersEnd();
263  ++pfclus) {
264  int caloclusidx = pfClusterMapES[*pfclus];
265  reco::CaloClusterPtr clusptr(caloClusHandleES, caloclusidx);
266  psclusters.push_back(clusptr);
267  clusptrsES[caloclusidx] = *pfclus;
268  }
269  eesc.setPreshowerClusters(psclusters);
270  }
271 
272  //fill association maps from output CaloClusters back to original PFClusters
273  edm::ValueMap<reco::CaloClusterPtr>::Filler fillerEBEE(*pfClusterAssociationEBEE);
274  fillerEBEE.insert(caloClusHandleEB, clusptrsEB.begin(), clusptrsEB.end());
275  fillerEBEE.insert(caloClusHandleEE, clusptrsEE.begin(), clusptrsEE.end());
276  fillerEBEE.fill();
277 
278  edm::ValueMap<reco::CaloClusterPtr>::Filler fillerES(*pfClusterAssociationES);
279  fillerES.insert(caloClusHandleES, clusptrsES.begin(), clusptrsES.end());
280  fillerES.fill();
281 
282  //store in the event
283  iEvent.put(std::move(pfClusterAssociationEBEE), PFClusterAssociationEBEE_);
284  iEvent.put(std::move(pfClusterAssociationES), PFClusterAssociationES_);
285  iEvent.put(std::move(superClusterAlgo_.getEBOutputSCCollection()), PFSuperClusterCollectionBarrel_);
286  iEvent.put(std::move(superClusterAlgo_.getEEOutputSCCollection()), PFSuperClusterCollectionEndcapWithPreshower_);
287 }
288 
291  desc.add<std::string>("PFSuperClusterCollectionEndcap", "particleFlowSuperClusterECALEndcap");
292  desc.add<bool>("doSatelliteClusterMerge", false);
293  desc.add<double>("thresh_PFClusterBarrel", 0.0);
294  desc.add<std::string>("PFBasicClusterCollectionBarrel", "particleFlowBasicClusterECALBarrel");
295  desc.add<bool>("useRegression", true);
296  desc.add<double>("satelliteMajorityFraction", 0.5);
297  desc.add<double>("thresh_PFClusterEndcap", 0.0);
298  desc.add<edm::InputTag>("ESAssociation", edm::InputTag("particleFlowClusterECAL"));
299  desc.add<std::string>("PFBasicClusterCollectionPreshower", "particleFlowBasicClusterECALPreshower");
300  desc.add<bool>("use_preshower", true);
301  desc.addUntracked<bool>("verbose", false);
302  desc.add<double>("thresh_SCEt", 4.0);
303  desc.add<double>("etawidth_SuperClusterEndcap", 0.04);
304  desc.add<double>("phiwidth_SuperClusterEndcap", 0.6);
305  desc.add<bool>("useDynamicDPhiWindow", true);
306  desc.add<std::string>("PFSuperClusterCollectionBarrel", "particleFlowSuperClusterECALBarrel");
307  {
309  psd0.add<bool>("isHLT", false);
310  psd0.add<bool>("applySigmaIetaIphiBug", false);
311  psd0.add<edm::InputTag>("ecalRecHitsEE", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
312  psd0.add<edm::InputTag>("ecalRecHitsEB", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
313  psd0.add<std::string>("regressionKeyEB", "pfscecal_EBCorrection_offline_v2");
314  psd0.add<std::string>("regressionKeyEE", "pfscecal_EECorrection_offline_v2");
315  psd0.add<std::string>("uncertaintyKeyEB", "pfscecal_EBUncertainty_offline_v2");
316  psd0.add<std::string>("uncertaintyKeyEE", "pfscecal_EEUncertainty_offline_v2");
317  psd0.add<edm::InputTag>("vertexCollection", edm::InputTag("offlinePrimaryVertices"));
318  psd0.add<double>("eRecHitThreshold", 1.);
319  desc.add<edm::ParameterSetDescription>("regressionConfig", psd0);
320  }
321  desc.add<bool>("applyCrackCorrections", false);
322  desc.add<double>("satelliteClusterSeedThreshold", 50.0);
323  desc.add<double>("etawidth_SuperClusterBarrel", 0.04);
324  desc.add<std::string>("PFBasicClusterCollectionEndcap", "particleFlowBasicClusterECALEndcap");
325  desc.add<edm::InputTag>("PFClusters", edm::InputTag("particleFlowClusterECAL"));
326  desc.add<double>("thresh_PFClusterSeedBarrel", 1.0);
327  desc.add<std::string>("ClusteringType", "Mustache");
328  desc.add<std::string>("EnergyWeight", "Raw");
329  desc.add<edm::InputTag>("BeamSpot", edm::InputTag("offlineBeamSpot"));
330  desc.add<double>("thresh_PFClusterSeedEndcap", 1.0);
331  desc.add<double>("phiwidth_SuperClusterBarrel", 0.6);
332  desc.add<double>("thresh_PFClusterES", 0.0);
333  desc.add<bool>("seedThresholdIsET", true);
334  desc.add<bool>("isOOTCollection", false);
335  desc.add<edm::InputTag>("barrelRecHits", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
336  desc.add<edm::InputTag>("endcapRecHits", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
337  desc.add<std::string>("PFSuperClusterCollectionEndcapWithPreshower",
338  "particleFlowSuperClusterECALEndcapWithPreshower");
339  desc.add<bool>("dropUnseedable", false);
340  descriptions.add("particleFlowSuperClusterECALMustache", desc);
341 }
ConfigurationDescriptions.h
PFECALSuperClusterProducer::beginLuminosityBlock
void beginLuminosityBlock(const edm::LuminosityBlock &, const edm::EventSetup &) override
Definition: PFECALSuperClusterProducer.cc:160
edm::helper::Filler::insert
void insert(const H &h, I begin, I end)
Definition: ValueMap.h:53
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
MessageLogger.h
ESHandle.h
edm::LuminosityBlock
Definition: LuminosityBlock.h:50
edm
HLT enums.
Definition: AlignableModifier.h:19
GBRWrapperRcd.h
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89287
edm::helper::Filler::fill
void fill()
Definition: ValueMap.h:65
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
edm::PtrVectorItr
Definition: PtrVector.h:51
HLT_FULL_cff.satelliteClusterSeedThreshold
satelliteClusterSeedThreshold
Definition: HLT_FULL_cff.py:14917
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
PFECALSuperClusterProducer::PFECALSuperClusterProducer
PFECALSuperClusterProducer(const edm::ParameterSet &)
Definition: PFECALSuperClusterProducer.cc:46
GBRForest.h
PFECALSuperClusterAlgo.h
PFECALSuperClusterAlgo::kMustache
Definition: PFECALSuperClusterAlgo.h:53
PFRecHit.h
HLT_FULL_cff.dropUnseedable
dropUnseedable
Definition: HLT_FULL_cff.py:14924
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
edm::PtrVector< CaloCluster >
reco::CaloCluster
Definition: CaloCluster.h:31
CaloClusterFwd.h
PFECALSuperClusterProducer::produce
void produce(edm::Event &, const edm::EventSetup &) override
Definition: PFECALSuperClusterProducer.cc:164
ParameterSetDescription.h
edm::ConfigurationDescriptions
Definition: ConfigurationDescriptions.h:28
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Vertex.h
PFECALSuperClusterAlgo::kRaw
Definition: PFECALSuperClusterAlgo.h:54
bsc_activity_cfg.clusters
clusters
Definition: bsc_activity_cfg.py:36
edm::PtrVector::push_back
void push_back(Ptr< T > const &iPtr)
Definition: PtrVector.h:149
PFECALSuperClusterProducer::~PFECALSuperClusterProducer
~PFECALSuperClusterProducer() override
Definition: PFECALSuperClusterProducer.cc:158
edm::ParameterSet
Definition: ParameterSet.h:47
PFECALSuperClusterProducer.h
PFECALSuperClusterProducer::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: PFECALSuperClusterProducer.cc:289
deltaR.h
iEvent
int iEvent
Definition: GenABIO.cc:224
PFClusterFwd.h
edm::EventSetup
Definition: EventSetup.h:57
EcalRecHit.h
edm::Ptr< CaloCluster >
ValueMap.h
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
eostools.move
def move(src, dest)
Definition: eostools.py:511
std
Definition: JetResolutionObject.h:76
PFECALSuperClusterAlgo::kCalibratedTotal
Definition: PFECALSuperClusterAlgo.h:54
edm::OrphanHandle
Definition: EDProductfwd.h:39
SuperClusterFwd.h
HLT_FULL_cff.seedThresholdIsET
seedThresholdIsET
Definition: HLT_FULL_cff.py:14930
SuperCluster.h
Exception
Definition: hltDiff.cc:246
HLT_FULL_cff.satelliteMajorityFraction
satelliteMajorityFraction
Definition: HLT_FULL_cff.py:14894
PFECALSuperClusterAlgo::kBOX
Definition: PFECALSuperClusterAlgo.h:53
EventSetup.h
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
edm::helper::Filler
Definition: ValueMap.h:22
HLT_FULL_cff.doSatelliteClusterMerge
doSatelliteClusterMerge
Definition: HLT_FULL_cff.py:14890
edm::Event
Definition: Event.h:73
PFECALSuperClusterAlgo::kCalibratedNoPS
Definition: PFECALSuperClusterAlgo.h:54
TauDecayModes.dec
dec
Definition: TauDecayModes.py:143
edm::InputTag
Definition: InputTag.h:15
CaloCluster.h