CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EcalSkim.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: EcalSkim
4 // Class: EcalSkim
5 //
6 //class EcalSkim EcalSkim.cc
7 //
8 // Original Author: Serena OGGERO
9 // Created: We May 14 10:10:52 CEST 2008
10 // Modified: Toyoko ORIMOTO
11 
12 #include <memory>
13 #include <vector>
14 #include <map>
15 #include <set>
16 
17 // user include files
19 
32 
33 using namespace edm;
34 using namespace std;
35 using namespace reco;
36 
38 {
39  BarrelClusterCollection = iConfig.getParameter<edm::InputTag>("barrelClusterCollection");
40  EndcapClusterCollection = iConfig.getParameter<edm::InputTag>("endcapClusterCollection");
41 
42  EnergyCutEB = iConfig.getUntrackedParameter<double>("energyCutEB");
43  EnergyCutEE = iConfig.getUntrackedParameter<double>("energyCutEE");
44 
45 }
46 
48 {
49 }
50 
52 {
53  int ievt = iEvent.id().event();
54 
56  edm::Handle<reco::SuperClusterCollection> eccHandle; // endcap
57 
58  iEvent.getByLabel("cosmicSuperClusters","CosmicBarrelSuperClusters", bccHandle);
59  if (!(bccHandle.isValid()))
60  {
61  LogWarning("EcalSkim") << BarrelClusterCollection << " not available in event " << ievt;
62  return false;
63  }
64  else
65  {
66  // edm::LogVerbatim("") << "I took the right barrel collection" ;
67  }
68  iEvent.getByLabel("cosmicSuperClusters","CosmicEndcapSuperClusters", eccHandle);
69 
70  if (!(eccHandle.isValid()))
71  {
72  LogWarning("EcalSkim") << EndcapClusterCollection << " not available";
73  //return false;
74  }
75  else
76  {
77  //edm::LogVerbatim("") << "I took the right endcap collection " ;
78  }
79 
80  bool accepted = false;
81  bool acceptedEB = false;
82  bool acceptedEE = false;
83 
84  // barrel
85  const reco::SuperClusterCollection *clusterCollectionEB = bccHandle.product();
86  for (reco::SuperClusterCollection::const_iterator clus = clusterCollectionEB->begin(); clus != clusterCollectionEB->end(); ++clus)
87  {
88  if ( clus->energy() >= EnergyCutEB )
89  {
90  acceptedEB = true;
91  break;
92  }
93  }
94 
95  // endcap
96  const reco::SuperClusterCollection *clusterCollectionEE = eccHandle.product();
97  for (reco::SuperClusterCollection::const_iterator clus = clusterCollectionEE->begin(); clus != clusterCollectionEE->end(); ++clus)
98  {
99  if ( clus->energy() >= EnergyCutEE )
100  {
101  acceptedEE = true;
102  break;
103  }
104  }
105 
106  // if there is at least one high energy cluster in EB OR EE, accept
107  if (acceptedEB||acceptedEE) accepted = true;
108 
109  return accepted;
110 
111 }
112 
113 //define this as a plug-in
T getParameter(std::string const &) const
EventNumber_t event() const
Definition: EventID.h:41
T getUntrackedParameter(std::string const &, T const &) const
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
int iEvent
Definition: GenABIO.cc:230
std::vector< SuperCluster > SuperClusterCollection
collection of SuperCluser objectr
virtual bool filter(edm::Event &, const edm::EventSetup &) override
Definition: EcalSkim.cc:51
EcalSkim(const edm::ParameterSet &)
Definition: EcalSkim.cc:37
bool isValid() const
Definition: HandleBase.h:75
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:420
T const * product() const
Definition: Handle.h:81
edm::EventID id() const
Definition: EventBase.h:60
~EcalSkim()
Definition: EcalSkim.cc:47