|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
#ifndef DATAFORMATS_METRECO_ECALHALODATA_H
#define DATAFORMATS_METRECO_ECALHALODATA_H
/*
[class]: EcalHaloData
[authors]: R. Remington, The University of Florida
[description]: Container class to store beam halo data specific to the Ecal subdetector
[date]: October 15, 2009
*/
#include "DataFormats/METReco/interface/PhiWedge.h"
#include <cmath>
#include "DataFormats/Common/interface/RefVector.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/Common/interface/RefProd.h"
#include "DataFormats/EgammaReco/interface/SuperCluster.h"
#include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
#include "DataFormats/EgammaReco/interface/BasicCluster.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/METReco/interface/HaloClusterCandidateECAL.h"
namespace reco {
class EcalHaloData {
public:
//Constructor
EcalHaloData();
//Destructor
~EcalHaloData() {}
// Number of Halo-like superclusters as a function of the roundness and angular cuts
//int NumberOfHaloSuperClusters(int iRoundness=2, int iAngle=2) const;
int NumberOfHaloSuperClusters(float roundness = 100., float angle = 4.0) const;
// Return collection of 5-degree Phi Wedges built from Ecal RecHits
const std::vector<PhiWedge>& GetPhiWedges() const { return PhiWedgeCollection; }
std::vector<PhiWedge>& GetPhiWedges() { return PhiWedgeCollection; }
// Get Reference to the SuperClusters
edm::RefVector<reco::SuperClusterCollection>& GetSuperClusters() { return TheSuperClusterRefs; }
const edm::RefVector<reco::SuperClusterCollection>& GetSuperClusters() const { return TheSuperClusterRefs; }
// Store Shower Shape variables for SuperClusters
edm::ValueMap<float>& GetShowerShapesRoundness() { return ShowerShapes_Roundness; }
const edm::ValueMap<float>& GetShowerShapesRoundness() const { return ShowerShapes_Roundness; }
edm::ValueMap<float>& GetShowerShapesAngle() { return ShowerShapes_Angle; }
const edm::ValueMap<float>& GetShowerShapesAngle() const { return ShowerShapes_Angle; }
const std::vector<HaloClusterCandidateECAL>& getHaloClusterCandidatesEB() const { return thehaloclustercands_eb; }
std::vector<HaloClusterCandidateECAL>& getHaloClusterCandidatesEB() { return thehaloclustercands_eb; }
const std::vector<HaloClusterCandidateECAL>& getHaloClusterCandidatesEE() const { return thehaloclustercands_ee; }
std::vector<HaloClusterCandidateECAL>& getHaloClusterCandidatesEE() { return thehaloclustercands_ee; }
void setHaloClusterCandidatesEB(const std::vector<HaloClusterCandidateECAL>& x) { thehaloclustercands_eb = x; }
void setHaloClusterCandidatesEE(const std::vector<HaloClusterCandidateECAL>& x) { thehaloclustercands_ee = x; }
private:
std::vector<PhiWedge> PhiWedgeCollection;
edm::RefVector<reco::SuperClusterCollection> TheSuperClusterRefs;
edm::ValueMap<float> ShowerShapes_Roundness;
edm::ValueMap<float> ShowerShapes_Angle;
std::vector<HaloClusterCandidateECAL> thehaloclustercands_eb;
std::vector<HaloClusterCandidateECAL> thehaloclustercands_ee;
};
} // namespace reco
#endif
|