CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Types | Private Member Functions | Private Attributes
l1t::L1TStage2CaloAnalyzer Class Reference
Inheritance diagram for l1t::L1TStage2CaloAnalyzer:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 L1TStage2CaloAnalyzer (const edm::ParameterSet &)
 
 ~L1TStage2CaloAnalyzer () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueue * globalLuminosityBlocksQueue ()
 
SerialTaskQueue * globalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 

Private Types

enum  ObjectType {
  Tower =1, Cluster =2, EG =3, Tau =4,
  Jet =5, SumET =6, SumETEm =7, SumHT =8,
  SumMET =9, SumMETHF =10, SumMHT =11, SumMHTHF =12,
  MPEG =13, MPTau =14, MPJet =15, MPSum =16,
  MPSumET =17, MPSumETHF =18, MPSumMETx =19, MPSumMETxHF =20,
  MPSumMETy =21, MPSumMETyHF =22, MPSumHT =23, MPSumHTHF =24,
  MPSumMHTx =25, MPSumMHTxHF =26, MPSumMHTy =27, MPSumMHTyHF =28,
  MPMinBiasHFP1 =29, MPMinBiasHFM1 =30, MPMinBiasHFP0 =31, MPMinBiasHFM0 =32,
  MinBiasHFP1 =33, MinBiasHFM1 =34, MinBiasHFP0 =35, MinBiasHFM0 =36,
  MPSumETEm = 37, MPSumHITowCount = 38, SumHITowCount = 39, SumCentrality = 40,
  SumAsymEt = 41, SumAsymEtHF = 42, SumAsymHt = 43, SumAsymHtHF = 44
}
 

Private Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void beginJob () override
 
void endJob () override
 

Private Attributes

std::map< ObjectType, TFileDirectory > dirs_
 
bool doHistos_
 
bool doText_
 
TFileDirectory evtDispDir_
 
std::map< ObjectType, TH1F * > hbx_
 
std::map< ObjectType, TH1F * > hem_
 
std::map< ObjectType, TH1F * > het_
 
std::map< ObjectType, TH1F * > heta_
 
std::map< ObjectType, TH2F * > hetaphi_
 
std::map< ObjectType, TH1F * > hhad_
 
std::map< ObjectType, TH1F * > hiso_
 
std::map< ObjectType, TH1F * > hphi_
 
std::map< ObjectType, TH1F * > hratio_
 
TH1F * hsort_
 
TH1F * hsortMP_
 
bool m_allBx = false
 
edm::EDGetToken m_clusterToken
 
int m_dmxBx = 0
 
bool m_doClusters
 
bool m_doEGs
 
bool m_doEvtDisp = false
 
bool m_doJets
 
bool m_doMPEGs
 
bool m_doMPJets
 
bool m_doMPSums
 
bool m_doMPTaus
 
bool m_doSums
 
bool m_doTaus
 
bool m_doTowers
 
edm::EDGetToken m_egToken
 
edm::EDGetToken m_jetToken
 
int m_mpBx = 0
 
edm::EDGetToken m_mpEGToken
 
edm::EDGetToken m_mpJetToken
 
edm::EDGetToken m_mpSumToken
 
edm::EDGetToken m_mpTauToken
 
edm::EDGetToken m_sumToken
 
edm::EDGetToken m_tauToken
 
edm::EDGetToken m_towerToken
 
std::vector< ObjectType > types_
 
std::vector< std::string > typeStr_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 29 of file L1TStage2CaloAnalyzer.cc.

Member Enumeration Documentation

enum l1t::L1TStage2CaloAnalyzer::ObjectType
private
Enumerator
Tower 
Cluster 
EG 
Tau 
Jet 
SumET 
SumETEm 
SumHT 
SumMET 
SumMETHF 
SumMHT 
SumMHTHF 
MPEG 
MPTau 
MPJet 
MPSum 
MPSumET 
MPSumETHF 
MPSumMETx 
MPSumMETxHF 
MPSumMETy 
MPSumMETyHF 
MPSumHT 
MPSumHTHF 
MPSumMHTx 
MPSumMHTxHF 
MPSumMHTy 
MPSumMHTyHF 
MPMinBiasHFP1 
MPMinBiasHFM1 
MPMinBiasHFP0 
MPMinBiasHFM0 
MinBiasHFP1 
MinBiasHFM1 
MinBiasHFP0 
MinBiasHFM0 
MPSumETEm 
MPSumHITowCount 
SumHITowCount 
SumCentrality 
SumAsymEt 
SumAsymEtHF 
SumAsymHt 
SumAsymHtHF 

Definition at line 73 of file L1TStage2CaloAnalyzer.cc.

73  {
74  Tower=1,
75  Cluster=2,
76  EG=3,
77  Tau=4,
78  Jet=5,
79  SumET=6,
80  SumETEm=7,
81  SumHT=8,
82  SumMET=9,
83  SumMETHF=10,
84  SumMHT=11,
85  SumMHTHF=12,
86  MPEG=13,
87  MPTau=14,
88  MPJet=15,
89  MPSum=16,
90  MPSumET=17,
91  MPSumETHF=18,
92  MPSumMETx=19,
93  MPSumMETxHF=20,
94  MPSumMETy=21,
95  MPSumMETyHF=22,
96  MPSumHT=23,
97  MPSumHTHF=24,
98  MPSumMHTx=25,
99  MPSumMHTxHF=26,
100  MPSumMHTy=27,
101  MPSumMHTyHF=28,
102  MPMinBiasHFP1=29,
103  MPMinBiasHFM1=30,
104  MPMinBiasHFP0=31,
105  MPMinBiasHFM0=32,
106  MinBiasHFP1=33,
107  MinBiasHFM1=34,
108  MinBiasHFP0=35,
109  MinBiasHFM0=36,
110  MPSumETEm = 37,
111  MPSumHITowCount = 38,
112  SumHITowCount = 39,
113  SumCentrality = 40,
114  SumAsymEt = 41,
115  SumAsymEtHF = 42,
116  SumAsymHt = 43,
117  SumAsymHtHF = 44
118  };
Jet
Definition: Jet.py:1
Tau
Definition: Tau.py:1

Constructor & Destructor Documentation

l1t::L1TStage2CaloAnalyzer::L1TStage2CaloAnalyzer ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 156 of file L1TStage2CaloAnalyzer.cc.

References Cluster, gather_cfg::cout, EG, edm::ParameterSet::getParameter(), mhtProducer_cfi::jetTag, edm::InputTag::label(), m_allBx, m_clusterToken, m_dmxBx, m_doClusters, m_doEGs, m_doEvtDisp, m_doJets, m_doMPEGs, m_doMPJets, m_doMPSums, m_doMPTaus, m_doSums, m_doTaus, m_doTowers, m_egToken, m_jetToken, m_mpBx, m_mpEGToken, m_mpJetToken, m_mpSumToken, m_mpTauToken, m_sumToken, m_tauToken, m_towerToken, MinBiasHFM0, MinBiasHFM1, MinBiasHFP0, MinBiasHFP1, MPEG, MPJet, MPMinBiasHFM0, MPMinBiasHFM1, MPMinBiasHFP0, MPMinBiasHFP1, MPSumET, MPSumETEm, MPSumETHF, MPSumHITowCount, MPSumHT, MPSumHTHF, MPSumMETx, MPSumMETxHF, MPSumMETy, MPSumMETyHF, MPSumMHTx, MPSumMHTxHF, MPSumMHTy, MPSumMHTyHF, MPTau, SumAsymEt, SumAsymEtHF, SumAsymHt, SumAsymHtHF, SumCentrality, SumET, SumETEm, SumHITowCount, SumHT, SumMET, SumMETHF, SumMHT, SumMHTHF, mhtProducer_cfi::tauTag, Tower, types_, and typeStr_.

156  :
157  doText_(iConfig.getUntrackedParameter<bool>("doText", true)),
158  doHistos_(iConfig.getUntrackedParameter<bool>("doHistos", true))
159  {
160  //now do what ever initialization is needed
161 
162  m_mpBx = iConfig.getParameter<int>("mpBx");
163  m_dmxBx = iConfig.getParameter<int>("dmxBx");
164  m_allBx = iConfig.getParameter<bool>("allBx");
165  m_doEvtDisp = iConfig.getParameter<bool>("doEvtDisp");
166 
167  // register what you consume and keep token for later access:
168  edm::InputTag nullTag("None");
169 
170  edm::InputTag towerTag = iConfig.getParameter<edm::InputTag>("towerToken");
171  m_towerToken = consumes<l1t::CaloTowerBxCollection>(towerTag);
172  m_doTowers = !(towerTag==nullTag);
173 
174  edm::InputTag clusterTag = iConfig.getParameter<edm::InputTag>("clusterToken");
175  m_clusterToken = consumes<l1t::CaloClusterBxCollection>(clusterTag);
176  m_doClusters = !(clusterTag==nullTag);
177 
178  edm::InputTag mpEGTag = iConfig.getParameter<edm::InputTag>("mpEGToken");
179  m_mpEGToken = consumes<l1t::EGammaBxCollection>(mpEGTag);
180  m_doMPEGs = !(mpEGTag==nullTag);
181 
182  edm::InputTag mpTauTag = iConfig.getParameter<edm::InputTag>("mpTauToken");
183  m_mpTauToken = consumes<l1t::TauBxCollection>(mpTauTag);
184  m_doMPTaus = !(mpTauTag==nullTag);
185 
186  edm::InputTag mpJetTag = iConfig.getParameter<edm::InputTag>("mpJetToken");
187  m_mpJetToken = consumes<l1t::JetBxCollection>(mpJetTag);
188  m_doMPJets = !(mpJetTag==nullTag);
189 
190  edm::InputTag mpSumTag = iConfig.getParameter<edm::InputTag>("mpEtSumToken");
191  m_mpSumToken = consumes<l1t::EtSumBxCollection>(mpSumTag);
192  m_doMPSums = !(mpSumTag==nullTag);
193 
194  edm::InputTag egTag = iConfig.getParameter<edm::InputTag>("egToken");
195  m_egToken = consumes<l1t::EGammaBxCollection>(egTag);
196  m_doEGs = !(egTag==nullTag);
197 
198  edm::InputTag tauTag = iConfig.getParameter<edm::InputTag>("tauToken");
199  m_tauToken = consumes<l1t::TauBxCollection>(tauTag);
200  m_doTaus = !(tauTag==nullTag);
201 
202  edm::InputTag jetTag = iConfig.getParameter<edm::InputTag>("jetToken");
203  m_jetToken = consumes<l1t::JetBxCollection>(jetTag);
204  m_doJets = !(jetTag==nullTag);
205 
206  edm::InputTag sumTag = iConfig.getParameter<edm::InputTag>("etSumToken");
207  m_sumToken = consumes<l1t::EtSumBxCollection>(sumTag);
208  m_doSums = !(sumTag==nullTag);
209 
210  std::cout << "Processing " << sumTag.label() << std::endl;
211 
212  types_.push_back( Tower );
213  types_.push_back( Cluster );
214  types_.push_back( MPEG );
215  types_.push_back( MPTau );
216  types_.push_back( MPJet );
217  types_.push_back( MPSumET );
218  types_.push_back( MPSumETHF );
219  types_.push_back( MPSumETEm );
220  types_.push_back( MPSumMETx );
221  types_.push_back( MPSumMETxHF );
222  types_.push_back( MPSumMETy );
223  types_.push_back( MPSumMETyHF );
224  types_.push_back( MPSumHT );
225  types_.push_back( MPSumHTHF );
226  types_.push_back( MPSumMHTx );
227  types_.push_back( MPSumMHTxHF );
228  types_.push_back( MPSumMHTy );
229  types_.push_back( MPSumMHTyHF );
230  types_.push_back( EG );
231  types_.push_back( Tau );
232  types_.push_back( Jet );
233  types_.push_back( SumET );
234  types_.push_back( SumETEm );
235  types_.push_back( SumHT );
236  types_.push_back( SumMET );
237  types_.push_back( SumMETHF );
238  types_.push_back( SumMHT );
239  types_.push_back( SumMHTHF );
240  types_.push_back( MPMinBiasHFP0 );
241  types_.push_back( MPMinBiasHFM0 );
242  types_.push_back( MPMinBiasHFP1 );
243  types_.push_back( MPMinBiasHFM1 );
244  types_.push_back( MinBiasHFP0 );
245  types_.push_back( MinBiasHFM0 );
246  types_.push_back( MinBiasHFP1 );
247  types_.push_back( MinBiasHFM1 );
248  types_.push_back( MPSumHITowCount );
249  types_.push_back( SumHITowCount );
250  types_.push_back( SumCentrality );
251  types_.push_back( SumAsymEt );
252  types_.push_back( SumAsymEtHF );
253  types_.push_back( SumAsymHt );
254  types_.push_back( SumAsymHtHF );
255 
256  typeStr_.push_back( "tower" );
257  typeStr_.push_back( "cluster" );
258  typeStr_.push_back( "mpeg" );
259  typeStr_.push_back( "mptau" );
260  typeStr_.push_back( "mpjet" );
261  typeStr_.push_back( "mpsumet" );
262  typeStr_.push_back( "mpsumethf" );
263  typeStr_.push_back( "mpsumetem" );
264  typeStr_.push_back( "mpsummetx" );
265  typeStr_.push_back( "mpsummetxhf" );
266  typeStr_.push_back( "mpsummety" );
267  typeStr_.push_back( "mpsummetyhf" );
268  typeStr_.push_back( "mpsumht" );
269  typeStr_.push_back( "mpsumhthf" );
270  typeStr_.push_back( "mpsummhtx" );
271  typeStr_.push_back( "mpsummhtxhf" );
272  typeStr_.push_back( "mpsummhty" );
273  typeStr_.push_back( "mpsummhtyhf" );
274  typeStr_.push_back( "eg" );
275  typeStr_.push_back( "tau" );
276  typeStr_.push_back( "jet" );
277  typeStr_.push_back( "sumet" );
278  typeStr_.push_back( "sumetem" );
279  typeStr_.push_back( "sumht" );
280  typeStr_.push_back( "summet" );
281  typeStr_.push_back( "summethf" );
282  typeStr_.push_back( "summht" );
283  typeStr_.push_back( "summhthf" );
284  typeStr_.push_back( "mpminbiashfp0" );
285  typeStr_.push_back( "mpminbiashfm0" );
286  typeStr_.push_back( "mpminbiashfp1" );
287  typeStr_.push_back( "mpminbiashfm1" );
288  typeStr_.push_back( "minbiashfp0" );
289  typeStr_.push_back( "minbiashfm0" );
290  typeStr_.push_back( "minbiashfp1" );
291  typeStr_.push_back( "minbiashfm1" );
292  typeStr_.push_back( "mpsumhitowercount");
293  typeStr_.push_back( "sumhitowercount");
294  typeStr_.push_back( "sumcentrality" );
295  typeStr_.push_back( "sumasymet" );
296  typeStr_.push_back( "sumasymethf" );
297  typeStr_.push_back( "sumasymht" );
298  typeStr_.push_back( "sumasymhthf" );
299  }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
mhtProducer_cfi.tauTag
tauTag
Definition: mhtProducer_cfi.py:10
mhtProducer_cfi.jetTag
jetTag
Definition: mhtProducer_cfi.py:7
Jet
Definition: Jet.py:1
l1t::L1TStage2CaloAnalyzer::typeStr_
std::vector< std::string > typeStr_
Definition: L1TStage2CaloAnalyzer.cc:121
l1t::L1TStage2CaloAnalyzer::types_
std::vector< ObjectType > types_
Definition: L1TStage2CaloAnalyzer.cc:120
Tau
Definition: Tau.py:1
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
gather_cfg.cout
cout
Definition: gather_cfg.py:144
l1t::L1TStage2CaloAnalyzer::~L1TStage2CaloAnalyzer ( )
override

Definition at line 302 of file L1TStage2CaloAnalyzer.cc.

303  {
304 
305  // do anything here that needs to be done at desctruction time
306  // (e.g. close files, deallocate resources etc.)
307 
308  }

Member Function Documentation

void l1t::L1TStage2CaloAnalyzer::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivate

Definition at line 317 of file L1TStage2CaloAnalyzer.cc.

References Cluster, fastPrimaryVertexProducer_cfi::clusters, gather_cfg::cout, dir, doText_, EG, edm::EventID::event(), evtDispDir_, edm::Event::getByToken(), hbx_, hem_, het_, heta_, hetaphi_, hhad_, hiso_, hphi_, hratio_, hsort_, hsortMP_, mps_fire::i, edm::EventBase::id(), fwrapper::jets, L1Analysis::kAsymEt, L1Analysis::kAsymEtHF, L1Analysis::kAsymHt, L1Analysis::kAsymHtHF, L1Analysis::kCentrality, L1Analysis::kMinBiasHFM0, L1Analysis::kMinBiasHFM1, L1Analysis::kMinBiasHFP0, L1Analysis::kMinBiasHFP1, L1Analysis::kMissingEt, L1Analysis::kMissingEtHF, L1Analysis::kMissingHt, L1Analysis::kMissingHtHF, L1Analysis::kTotalEt, L1Analysis::kTotalEtEm, L1Analysis::kTotalEtHF, L1Analysis::kTotalEtx, L1Analysis::kTotalEtxHF, L1Analysis::kTotalEty, L1Analysis::kTotalEtyHF, L1Analysis::kTotalHt, L1Analysis::kTotalHtHF, L1Analysis::kTotalHtx, L1Analysis::kTotalHtxHF, L1Analysis::kTotalHty, L1Analysis::kTotalHtyHF, L1Analysis::kTowerCount, m_allBx, m_clusterToken, m_dmxBx, m_doClusters, m_doEGs, m_doEvtDisp, m_doJets, m_doMPEGs, m_doMPJets, m_doMPSums, m_doMPTaus, m_doSums, m_doTaus, m_doTowers, m_egToken, m_jetToken, m_mpBx, m_mpEGToken, m_mpJetToken, m_mpSumToken, m_mpTauToken, m_sumToken, m_tauToken, m_towerToken, TFileDirectory::make(), MinBiasHFM0, MinBiasHFM1, MinBiasHFP0, MinBiasHFP1, TFileDirectory::mkdir(), MPEG, MPJet, MPMinBiasHFM0, MPMinBiasHFM1, MPMinBiasHFP0, MPMinBiasHFP1, MPSumET, MPSumETEm, MPSumETHF, MPSumHITowCount, MPSumHT, MPSumHTHF, MPSumMETx, MPSumMETxHF, MPSumMETy, MPSumMETyHF, MPSumMHTx, MPSumMHTxHF, MPSumMHTy, MPSumMHTyHF, MPTau, edm::Event::run(), SumAsymEt, SumAsymEtHF, SumAsymHt, SumAsymHtHF, SumCentrality, SumET, SumETEm, SumHITowCount, SumHT, SumMET, SumMETHF, SumMHT, SumMHTHF, nano_cff::taus, RecoTauValidation_cfi::text, and Tower.

318  {
319  using namespace edm;
320 
321  std::stringstream text;
322 
323  TH2F* hEvtTow = new TH2F();
324  TH2F* hEvtMPEG = new TH2F();
325  TH2F* hEvtMPTau = new TH2F();
326  TH2F* hEvtMPJet = new TH2F();
327  TH2F* hEvtDemuxEG = new TH2F();
328  TH2F* hEvtDemuxTau = new TH2F();
329  TH2F* hEvtDemuxJet = new TH2F();
330 
331 
332  if (m_doEvtDisp) {
333  std::stringstream ss;
334  ss << iEvent.run() << "-" << iEvent.id().event();
335  TFileDirectory dir = evtDispDir_.mkdir(ss.str());
336  hEvtTow = dir.make<TH2F>("Tower", "", 83, -41.5, 41.5, 72, .5, 72.5);
337  hEvtMPEG = dir.make<TH2F>("MPEG", "", 83, -41.5, 41.5, 72, .5, 72.5);
338  hEvtMPTau = dir.make<TH2F>("MPTau", "", 83, -41.5, 41.5, 72, .5, 72.5);
339  hEvtMPJet = dir.make<TH2F>("MPJet", "", 83, -41.5, 41.5, 72, .5, 72.5);
340  hEvtDemuxEG = dir.make<TH2F>("DemuxEG", "", 227, -113.5, 113.5, 144, -0.5, 143.5);
341  hEvtDemuxTau = dir.make<TH2F>("DemuxTau", "", 227, -113.5, 113.5, 144, -0.5, 143.5);
342  hEvtDemuxJet = dir.make<TH2F>("DemuxJet", "", 227, -113.5, 113.5, 144, -0.5, 143.5);
343  }
344 
345  // get TPs ?
346  // get regions ?
347  // get RCT clusters ?
348 
349  //check mpbx and dmxbx
350  if(m_mpBx < -2 || m_mpBx > 2 || m_dmxBx < -2 || m_dmxBx > 2)
351  edm::LogError("L1T") << "Selected MP Bx or Demux Bx to fill histograms is outside of range -2,2. Histos will be empty!";
352 
353 
354  // get towers
355  if (m_doTowers) {
356  Handle< BXVector<l1t::CaloTower> > towers;
357  iEvent.getByToken(m_towerToken,towers);
358 
359  for ( int ibx=towers->getFirstBX(); ibx<=towers->getLastBX(); ++ibx) {
360 
361  if ( !m_allBx && ibx != m_mpBx ) continue;
362 
363  for ( auto itr = towers->begin(ibx); itr !=towers->end(ibx); ++itr ) {
364 
365  if (itr->hwPt()<=0) continue;
366 
367  hbx_.at(Tower)->Fill( ibx );
368  het_.at(Tower)->Fill( itr->hwPt() );
369  heta_.at(Tower)->Fill( itr->hwEta() );
370  hphi_.at(Tower)->Fill( itr->hwPhi() );
371  hem_.at(Tower)->Fill( itr->hwEtEm() );
372  hhad_.at(Tower)->Fill( itr->hwEtHad() );
373  hratio_.at(Tower)->Fill( itr->hwEtRatio() );
374  hetaphi_.at(Tower)->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
375 
376  text << "Tower : " << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << " iem=" << itr->hwEtEm() << " ihad=" << itr->hwEtHad() << " iratio=" << itr->hwEtRatio() << std::endl;
377 
378  if (m_doEvtDisp) hEvtTow->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
379 
380  }
381 
382  }
383 
384  }
385  // get cluster
386  if (m_doClusters) {
387  Handle< BXVector<l1t::CaloCluster> > clusters;
388  iEvent.getByToken(m_clusterToken,clusters);
389 
390  for ( int ibx=clusters->getFirstBX(); ibx<=clusters->getLastBX(); ++ibx) {
391 
392  if ( !m_allBx && ibx != m_mpBx ) continue;
393 
394  for ( auto itr = clusters->begin(ibx); itr !=clusters->end(ibx); ++itr ) {
395  hbx_.at(Cluster)->Fill( ibx );
396  het_.at(Cluster)->Fill( itr->hwPt() );
397  heta_.at(Cluster)->Fill( itr->hwEta() );
398  hphi_.at(Cluster)->Fill( itr->hwPhi() );
399  hetaphi_.at(Cluster)->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
400  text << "Cluster : " << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
401  }
402 
403  }
404  }
405 
406  // get EG
407  if (m_doMPEGs) {
408  Handle< BXVector<l1t::EGamma> > mpegs;
409  iEvent.getByToken(m_mpEGToken,mpegs);
410 
411  for ( int ibx=mpegs->getFirstBX(); ibx<=mpegs->getLastBX(); ++ibx) {
412 
413  if ( !m_allBx && ibx != m_mpBx ) continue;
414 
415  for ( auto itr = mpegs->begin(ibx); itr != mpegs->end(ibx); ++itr ) {
416  hbx_.at(MPEG)->Fill( ibx );
417  het_.at(MPEG)->Fill( itr->hwPt() );
418  heta_.at(MPEG)->Fill( itr->hwEta() );
419  hphi_.at(MPEG)->Fill( itr->hwPhi() );
420  hetaphi_.at(MPEG)->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
421  hiso_.at(MPEG)->Fill( itr->hwIso() );
422 
423  text << "MP EG : " << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
424 
425  if (m_doEvtDisp) hEvtMPEG->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
426 
427  }
428  }
429 
430  }
431 
432  // get tau
433  if (m_doMPTaus) {
434  Handle< BXVector<l1t::Tau> > mptaus;
435  iEvent.getByToken(m_mpTauToken,mptaus);
436 
437  for ( int ibx=mptaus->getFirstBX(); ibx<=mptaus->getLastBX(); ++ibx) {
438 
439  if ( !m_allBx && ibx != m_mpBx ) continue;
440 
441  for ( auto itr = mptaus->begin(ibx); itr != mptaus->end(ibx); ++itr ) {
442  hbx_.at(MPTau)->Fill( ibx );
443  het_.at(MPTau)->Fill( itr->hwPt() );
444  heta_.at(MPTau)->Fill( itr->hwEta() );
445  hphi_.at(MPTau)->Fill( itr->hwPhi() );
446  hetaphi_.at(MPTau)->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
447 
448  text << "MP Tau : " << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
449 
450  if (m_doEvtDisp) hEvtMPTau->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
451  }
452 
453  }
454 
455  }
456 
457  // get jet
458  int njetmp=0;
459  std::vector<l1t::Jet> thejets_poseta;
460  std::vector<l1t::Jet> thejets_negeta;
461 
462  if (m_doMPJets) {
463  Handle< BXVector<l1t::Jet> > mpjets;
464  iEvent.getByToken(m_mpJetToken,mpjets);
465 
466  //Handle<BXVector<l1t::Jet>> jets;
467  //iEvent.getByToken(m_jetToken,jets);
468  //if (mpjets->size(0) == jets->size(0)) std::cout<<"******notequal"<<std::endl;
469  for ( int ibx=mpjets->getFirstBX(); ibx<=mpjets->getLastBX(); ++ibx) {
470 
471  if ( !m_allBx && ibx != m_mpBx ) continue;
472 
473  for ( auto itr = mpjets->begin(ibx); itr != mpjets->end(ibx); ++itr ) {
474  njetmp+=1;
475  hbx_.at(MPJet)->Fill( ibx );
476  het_.at(MPJet)->Fill( itr->hwPt() );
477  heta_.at(MPJet)->Fill( itr->hwEta() );
478  hphi_.at(MPJet)->Fill( itr->hwPhi() );
479  hetaphi_.at(MPJet)->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
480 
481  text << "MP Jet : " << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
482 
483  if (m_doEvtDisp) hEvtMPJet->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
484 
485  itr->hwEta()>0 ? thejets_poseta.push_back(*itr) : thejets_negeta.push_back(*itr);
486  }
487 
488  }
489 
490  }
491 
492  if (!thejets_poseta.empty()) {
493  for (unsigned int i=0; i<thejets_poseta.size()-1; i++) {
494  for (unsigned int j=i+1; j<thejets_poseta.size(); j++) {
495  hsortMP_->Fill(thejets_poseta.at(i).hwPt()-thejets_poseta.at(j).hwPt());
496  }
497  }
498  }
499 
500  if (!thejets_negeta.empty()) {
501  for (unsigned int i=0; i<thejets_negeta.size()-1; i++) {
502  for (unsigned int j=i+1; j<thejets_negeta.size(); j++) {
503  hsortMP_->Fill(thejets_negeta.at(i).hwPt()-thejets_negeta.at(j).hwPt());
504  }
505  }
506  }
507 
508  //std::cout<<"njetmp "<<njetmp<<std::endl;
509 
510  // get sums
511  if (m_doMPSums) {
512  Handle< BXVector<l1t::EtSum> > mpsums;
513  iEvent.getByToken(m_mpSumToken,mpsums);
514 
515  for ( int ibx=mpsums->getFirstBX(); ibx<=mpsums->getLastBX(); ++ibx) {
516 
517  if ( !m_allBx && ibx != m_mpBx ) continue;
518 
519  for ( auto itr = mpsums->begin(ibx); itr != mpsums->end(ibx); ++itr ) {
520 
521  switch(itr->getType()){
522  case l1t::EtSum::EtSumType::kTotalEt: het_.at(MPSumET) ->Fill( itr->hwPt() ); break;
523  case l1t::EtSum::EtSumType::kTotalEtHF: het_.at(MPSumETHF) ->Fill( itr->hwPt() ); break;
524  case l1t::EtSum::EtSumType::kTotalEtEm: het_.at(MPSumETEm) ->Fill( itr->hwPt() ); break;
525  case l1t::EtSum::EtSumType::kTotalEtx: het_.at(MPSumMETx) ->Fill( itr->hwPt() ); break;
526  case l1t::EtSum::EtSumType::kTotalEtxHF: het_.at(MPSumMETxHF) ->Fill( itr->hwPt() ); break;
527  case l1t::EtSum::EtSumType::kTotalEty: het_.at(MPSumMETy) ->Fill( itr->hwPt() ); break;
528  case l1t::EtSum::EtSumType::kTotalEtyHF: het_.at(MPSumMETyHF) ->Fill( itr->hwPt() ); break;
529  case l1t::EtSum::EtSumType::kTotalHt: het_.at(MPSumHT) ->Fill( itr->hwPt() ); break;
530  case l1t::EtSum::EtSumType::kTotalHtHF: het_.at(MPSumHTHF) ->Fill( itr->hwPt() ); break;
531  case l1t::EtSum::EtSumType::kTotalHtx: het_.at(MPSumMHTx) ->Fill( itr->hwPt() ); break;
532  case l1t::EtSum::EtSumType::kTotalHtxHF: het_.at(MPSumMHTxHF) ->Fill( itr->hwPt() ); break;
533  case l1t::EtSum::EtSumType::kTotalHty: het_.at(MPSumMHTy) ->Fill( itr->hwPt() ); break;
534  case l1t::EtSum::EtSumType::kTotalHtyHF: het_.at(MPSumMHTyHF) ->Fill( itr->hwPt() ); break;
535  case l1t::EtSum::EtSumType::kMinBiasHFP0: het_.at(MPMinBiasHFP0) ->Fill( itr->hwPt() ); break;
536  case l1t::EtSum::EtSumType::kMinBiasHFM0: het_.at(MPMinBiasHFM0) ->Fill( itr->hwPt() ); break;
537  case l1t::EtSum::EtSumType::kMinBiasHFP1: het_.at(MPMinBiasHFP1) ->Fill( itr->hwPt() ); break;
538  case l1t::EtSum::EtSumType::kMinBiasHFM1: het_.at(MPMinBiasHFM1) ->Fill( itr->hwPt() ); break;
539  case l1t::EtSum::EtSumType::kTowerCount: het_.at(MPSumHITowCount) ->Fill( itr->hwPt() ); break;
540 
541  default: std::cout<<"wrong type of MP sum"<<std::endl;
542  }
543  text << "MP Sum : " << " type=" << itr->getType() << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
544  }
545 
546  }
547 
548  }
549 
550  // get EG
551  if (m_doEGs) {
552  Handle< BXVector<l1t::EGamma> > egs;
553  iEvent.getByToken(m_egToken,egs);
554 
555  for ( int ibx=egs->getFirstBX(); ibx<=egs->getLastBX(); ++ibx) {
556 
557  if ( !m_allBx && ibx != m_dmxBx ) continue;
558 
559  for ( auto itr = egs->begin(ibx); itr != egs->end(ibx); ++itr ) {
560  hbx_.at(EG)->Fill( ibx );
561  het_.at(EG)->Fill( itr->hwPt() );
562  heta_.at(EG)->Fill( itr->hwEta() );
563  hphi_.at(EG)->Fill( itr->hwPhi() );
564  hetaphi_.at(EG)->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
565  hiso_.at(EG)->Fill( itr->hwIso() );
566 
567  text << "EG : " << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
568 
569  if (m_doEvtDisp) hEvtDemuxEG->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
570  }
571 
572  }
573 
574  }
575 
576  // get tau
577  if (m_doTaus) {
578  Handle< BXVector<l1t::Tau> > taus;
579  iEvent.getByToken(m_tauToken,taus);
580 
581  for ( int ibx=taus->getFirstBX(); ibx<=taus->getLastBX(); ++ibx) {
582 
583  if ( !m_allBx && ibx != m_dmxBx ) continue;
584 
585  for ( auto itr = taus->begin(ibx); itr != taus->end(ibx); ++itr ) {
586  hbx_.at(Tau)->Fill( ibx );
587  het_.at(Tau)->Fill( itr->hwPt() );
588  heta_.at(Tau)->Fill( itr->hwEta() );
589  hphi_.at(Tau)->Fill( itr->hwPhi() );
590  hetaphi_.at(Tau)->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
591 
592  text << "Tau : " << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
593 
594  if (m_doEvtDisp) hEvtDemuxTau->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
595  }
596 
597  }
598 
599  }
600 
601 
602  // get jet
603  int njetdem=0;
604  std::vector<l1t::Jet> thejets;
605 
606  if (m_doJets) {
607  Handle< BXVector<l1t::Jet> > jets;
608  iEvent.getByToken(m_jetToken,jets);
609 
610  for ( int ibx=jets->getFirstBX(); ibx<=jets->getLastBX(); ++ibx) {
611 
612  if ( !m_allBx && ibx != m_dmxBx ) continue;
613 
614  for ( auto itr = jets->begin(ibx); itr != jets->end(ibx); ++itr ) {
615  njetdem+=1;
616  hbx_.at(Jet)->Fill( ibx );
617  het_.at(Jet)->Fill( itr->hwPt() );
618  heta_.at(Jet)->Fill( itr->hwEta() );
619  hphi_.at(Jet)->Fill( itr->hwPhi() );
620  hetaphi_.at(Jet)->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
621 
622  text << "Jet : " << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
623 
624  if (m_doEvtDisp) hEvtDemuxJet->Fill( itr->hwEta(), itr->hwPhi(), itr->hwPt() );
625 
626  thejets.push_back(*itr);
627  }
628 
629  }
630 
631  }
632 
633  if (!thejets.empty()) {
634  for (unsigned int i=0; i<thejets.size()-1; i++) {
635  for (unsigned int j=i+1; j<thejets.size(); j++) {
636  hsort_->Fill(thejets.at(i).hwPt()-thejets.at(j).hwPt());
637  }
638  }
639  }
640 
641  // get sums
642  if (m_doSums) {
643  Handle< BXVector<l1t::EtSum> > sums;
644  iEvent.getByToken(m_sumToken,sums);
645 
646  for ( int ibx=sums->getFirstBX(); ibx<=sums->getLastBX(); ++ibx) {
647 
648  if ( !m_allBx && ibx != m_dmxBx ) continue;
649 
650  for ( auto itr = sums->begin(ibx); itr != sums->end(ibx); ++itr ) {
651 
652  switch(itr->getType()){
653  case l1t::EtSum::EtSumType::kTotalEt: het_.at(SumET) ->Fill( itr->hwPt() ); break;
654  case l1t::EtSum::EtSumType::kTotalEtHF: break;
655  case l1t::EtSum::EtSumType::kTotalEtEm: het_.at(SumETEm) ->Fill( itr->hwPt() ); break;
656  case l1t::EtSum::EtSumType::kTotalHt: het_.at(SumHT) ->Fill( itr->hwPt() ); break;
657  case l1t::EtSum::EtSumType::kTotalHtHF: break;
658  case l1t::EtSum::EtSumType::kMissingEt: het_.at(SumMET) ->Fill( itr->hwPt() ); hphi_.at(SumMET) ->Fill( itr->hwPhi() ); break;
659  case l1t::EtSum::EtSumType::kMissingEtHF: het_.at(SumMETHF) ->Fill( itr->hwPt() ); hphi_.at(SumMETHF) ->Fill( itr->hwPhi() ); break;
660  case l1t::EtSum::EtSumType::kMissingHt: het_.at(SumMHT) ->Fill( itr->hwPt() ); hphi_.at(SumMHT) ->Fill( itr->hwPhi() ); break;
661  case l1t::EtSum::EtSumType::kMissingHtHF: het_.at(SumMHTHF) ->Fill( itr->hwPt() ); hphi_.at(SumMHTHF) ->Fill( itr->hwPhi() ); break;
662  case l1t::EtSum::EtSumType::kMinBiasHFP0: het_.at(MinBiasHFP0) ->Fill( itr->hwPt() ); break;
663  case l1t::EtSum::EtSumType::kMinBiasHFM0: het_.at(MinBiasHFM0) ->Fill( itr->hwPt() ); break;
664  case l1t::EtSum::EtSumType::kMinBiasHFP1: het_.at(MinBiasHFP1) ->Fill( itr->hwPt() ); break;
665  case l1t::EtSum::EtSumType::kMinBiasHFM1: het_.at(MinBiasHFM1) ->Fill( itr->hwPt() ); break;
666  case l1t::EtSum::EtSumType::kTowerCount: het_.at(SumHITowCount) ->Fill( itr->hwPt() ); break;
667  case l1t::EtSum::EtSumType::kCentrality: het_.at(SumCentrality) ->Fill( itr->hwPt() ); break;
668  case l1t::EtSum::EtSumType::kAsymEt: het_.at(SumAsymEt) ->Fill( itr->hwPt() ); break;
669  case l1t::EtSum::EtSumType::kAsymHt: het_.at(SumAsymHt) ->Fill( itr->hwPt() ); break;
670  case l1t::EtSum::EtSumType::kAsymEtHF: het_.at(SumAsymEtHF) ->Fill( itr->hwPt() ); break;
671  case l1t::EtSum::EtSumType::kAsymHtHF: het_.at(SumAsymHtHF) ->Fill( itr->hwPt() ); break;
672 
673  default: std::cout<<"wrong type of demux sum: "<< itr->getType() << std::endl;
674  }
675 
676  text << "Sum : " << " type=" << itr->getType() << " BX=" << ibx << " ipt=" << itr->hwPt() << " ieta=" << itr->hwEta() << " iphi=" << itr->hwPhi() << std::endl;
677  }
678  }
679  }
680 
681  if (doText_) edm::LogVerbatim("L1TCaloEvents") << text.str();
682 
683  }
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
l1t::L1TStage2CaloAnalyzer::het_
std::map< ObjectType, TH1F * > het_
Definition: L1TStage2CaloAnalyzer.cc:124
l1t::L1TStage2CaloAnalyzer::hhad_
std::map< ObjectType, TH1F * > hhad_
Definition: L1TStage2CaloAnalyzer.cc:129
mps_fire.i
i
Definition: mps_fire.py:338
l1t::L1TStage2CaloAnalyzer::heta_
std::map< ObjectType, TH1F * > heta_
Definition: L1TStage2CaloAnalyzer.cc:125
l1t::L1TStage2CaloAnalyzer::hbx_
std::map< ObjectType, TH1F * > hbx_
Definition: L1TStage2CaloAnalyzer.cc:127
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
l1t::L1TStage2CaloAnalyzer::hetaphi_
std::map< ObjectType, TH2F * > hetaphi_
Definition: L1TStage2CaloAnalyzer.cc:131
RecoTauValidation_cfi.text
text
Definition: RecoTauValidation_cfi.py:338
Jet
Definition: Jet.py:1
l1t::L1TStage2CaloAnalyzer::hem_
std::map< ObjectType, TH1F * > hem_
Definition: L1TStage2CaloAnalyzer.cc:128
fwrapper::jets
vector< PseudoJet > jets
Definition: fastjetfortran_madfks.cc:42
edm::Event::run
RunNumber_t run() const
Definition: Event.h:101
l1t::L1TStage2CaloAnalyzer::hratio_
std::map< ObjectType, TH1F * > hratio_
Definition: L1TStage2CaloAnalyzer.cc:130
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:59
l1t::L1TStage2CaloAnalyzer::hphi_
std::map< ObjectType, TH1F * > hphi_
Definition: L1TStage2CaloAnalyzer.cc:126
l1t::L1TStage2CaloAnalyzer::hiso_
std::map< ObjectType, TH1F * > hiso_
Definition: L1TStage2CaloAnalyzer.cc:132
Tau
Definition: Tau.py:1
fastPrimaryVertexProducer_cfi.clusters
clusters
Definition: fastPrimaryVertexProducer_cfi.py:5
TFileDirectory::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileDirectory.cc:162
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:59
edm
HLT enums.
Definition: AlignableModifier.h:17
gather_cfg.cout
cout
Definition: gather_cfg.py:144
dir
dbl *** dir
Definition: mlp_gen.cc:35
nano_cff.taus
taus
Definition: nano_cff.py:36
void l1t::L1TStage2CaloAnalyzer::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 688 of file L1TStage2CaloAnalyzer.cc.

References Cluster, dirs_, EG, evtDispDir_, hbx_, hem_, het_, heta_, hetaphi_, hhad_, hiso_, hphi_, hratio_, hsort_, hsortMP_, m_doEvtDisp, TFileService::make(), MinBiasHFM0, MinBiasHFM1, MinBiasHFP0, MinBiasHFP1, TFileService::mkdir(), MPEG, MPJet, MPMinBiasHFM0, MPMinBiasHFM1, MPMinBiasHFP0, MPMinBiasHFP1, MPSumET, MPSumETEm, MPSumETHF, MPSumHITowCount, MPSumHT, MPSumHTHF, MPSumMETx, MPSumMETxHF, MPSumMETy, MPSumMETyHF, MPSumMHTx, MPSumMHTxHF, MPSumMHTy, MPSumMHTyHF, MPTau, str, SumET, SumETEm, SumHITowCount, SumHT, SumMET, SumMETHF, SumMHT, SumMHTHF, Tower, types_, and typeStr_.

689  {
690 
691  edm::Service<TFileService> fs;
692 
693  auto itr = types_.cbegin();
694  auto str = typeStr_.cbegin();
695 
696  for (; itr!=types_.end(); ++itr, ++str ) {
697 
698  dirs_.insert( std::pair< ObjectType, TFileDirectory >(*itr, fs->mkdir(*str) ) );
699 
700  if (*itr==MPSumMETx || *itr == MPSumMETxHF || *itr==MPSumMETy || *itr==MPSumMETyHF){
701  het_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("et", "", 2200000, -2199999500, 2200000500) ));
702  }
703  else if (*itr==MPSumMHTx || *itr==MPSumMHTxHF || *itr==MPSumMHTy || *itr==MPSumMHTyHF ) {
704  het_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("et", "", 2200000, -2199999500, 2200000500) ));
705  }
706  else if (*itr==SumET || *itr==SumETEm || *itr==MPSumETEm || *itr==MPSumET || *itr==MPSumETHF || *itr==SumHT || *itr==MPSumHT || *itr==MPSumHTHF ) {
707  het_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("et", "", 100000, -0.5, 99999.5) ));
708  }
709  else if (*itr==MPMinBiasHFP0 ||
710  *itr==MPMinBiasHFM0 ||
711  *itr==MPMinBiasHFP1 ||
712  *itr==MPMinBiasHFM1 ||
713  *itr==MinBiasHFP0 ||
714  *itr==MinBiasHFM0 ||
715  *itr==MinBiasHFP1 ||
716  *itr==MinBiasHFM1) {
717  het_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("et", "", 16, -0.5, 15.5) ));
718  }
719  else if (*itr==MPSumHITowCount || *itr==SumHITowCount){
720  het_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("et", "", 5904, -0.5, 5903.5)));
721  }
722  else if (*itr==EG || *itr==Tau || *itr==MPEG || *itr==MPTau){
723  het_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("et", "", 1500, -0.5, 1499.5) ));
724  }
725  else {
726  het_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("et", "", 100000, -0.5, 99999.5) ));
727  }
728 
729  hbx_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("bx", "", 11, -5.5, 5.5) ));
730 
731  if (*itr==EG || *itr==Jet || *itr==Tau) {
732  heta_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("eta", "", 227, -113.5, 113.5) ));
733  hphi_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("phi", "", 144, -0.5, 143.5) ));
734  hetaphi_.insert( std::pair< ObjectType, TH2F* >(*itr, dirs_.at(*itr).make<TH2F>("etaphi", "", 227, -113.5, 113.5, 144, -0.5, 143.5) ));
735  if(*itr==EG) hiso_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("iso", "", 4, -0.5, 3.5) ));
736  }
737  else if (*itr==Tower || *itr==Cluster || *itr==MPEG || *itr==MPJet || *itr==MPTau) {
738  heta_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("eta", "", 83, -41.5, 41.5) ));
739  hphi_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("phi", "", 72, 0.5, 72.5) ));
740  hetaphi_.insert( std::pair< ObjectType, TH2F* >(*itr, dirs_.at(*itr).make<TH2F>("etaphi", "", 83, -41.5, 41.5, 72, .5, 72.5) ));
741  if(*itr==MPEG) hiso_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("iso", "", 4, -0.5, 3.5) ));
742  }
743  else if (*itr==SumMET || *itr==SumMETHF || *itr==SumMHT || *itr==SumMHTHF) {
744  hphi_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("phi", "", 1008, -0.5, 1007.5) ));
745  }
746 
747  if (*itr==Tower) {
748  hem_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("em", "", 101, -0.5, 100.5) ));
749  hhad_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("had", "", 101, -0.5, 100.5) ));
750  hratio_.insert( std::pair< ObjectType, TH1F* >(*itr, dirs_.at(*itr).make<TH1F>("ratio", "", 11, -0.5, 10.5) ));
751  }
752 
753  }
754 
755  if (m_doEvtDisp) {
756  evtDispDir_ = fs->mkdir("Events");
757  }
758 
759  hsort_ = fs->make<TH1F>("sort","",201,-100.5,100.5);
760  hsortMP_ = fs->make<TH1F>("sortMP","",201,-100.5,100.5);
761 
762  }
l1t::L1TStage2CaloAnalyzer::het_
std::map< ObjectType, TH1F * > het_
Definition: L1TStage2CaloAnalyzer.cc:124
l1t::L1TStage2CaloAnalyzer::hhad_
std::map< ObjectType, TH1F * > hhad_
Definition: L1TStage2CaloAnalyzer.cc:129
l1t::L1TStage2CaloAnalyzer::heta_
std::map< ObjectType, TH1F * > heta_
Definition: L1TStage2CaloAnalyzer.cc:125
l1t::L1TStage2CaloAnalyzer::hbx_
std::map< ObjectType, TH1F * > hbx_
Definition: L1TStage2CaloAnalyzer.cc:127
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
l1t::L1TStage2CaloAnalyzer::hetaphi_
std::map< ObjectType, TH2F * > hetaphi_
Definition: L1TStage2CaloAnalyzer.cc:131
l1t::L1TStage2CaloAnalyzer::dirs_
std::map< ObjectType, TFileDirectory > dirs_
Definition: L1TStage2CaloAnalyzer.cc:123
Jet
Definition: Jet.py:1
l1t::L1TStage2CaloAnalyzer::hem_
std::map< ObjectType, TH1F * > hem_
Definition: L1TStage2CaloAnalyzer.cc:128
l1t::L1TStage2CaloAnalyzer::hratio_
std::map< ObjectType, TH1F * > hratio_
Definition: L1TStage2CaloAnalyzer.cc:130
l1t::L1TStage2CaloAnalyzer::typeStr_
std::vector< std::string > typeStr_
Definition: L1TStage2CaloAnalyzer.cc:121
l1t::L1TStage2CaloAnalyzer::hphi_
std::map< ObjectType, TH1F * > hphi_
Definition: L1TStage2CaloAnalyzer.cc:126
l1t::L1TStage2CaloAnalyzer::types_
std::vector< ObjectType > types_
Definition: L1TStage2CaloAnalyzer.cc:120
l1t::L1TStage2CaloAnalyzer::hiso_
std::map< ObjectType, TH1F * > hiso_
Definition: L1TStage2CaloAnalyzer.cc:132
Tau
Definition: Tau.py:1
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
str
#define str(s)
Definition: TestProcessor.cc:48
void l1t::L1TStage2CaloAnalyzer::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 766 of file L1TStage2CaloAnalyzer.cc.

767  {
768  }
void l1t::L1TStage2CaloAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 804 of file L1TStage2CaloAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setUnknown().

804  {
805  //The following says we do not know what parameters are allowed so do no validation
806  // Please change this to state exactly what you do use, even if it is no parameters
807  edm::ParameterSetDescription desc;
808  desc.setUnknown();
809  descriptions.addDefault(desc);
810  }
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99

Member Data Documentation

std::map< ObjectType, TFileDirectory > l1t::L1TStage2CaloAnalyzer::dirs_
private

Definition at line 123 of file L1TStage2CaloAnalyzer.cc.

Referenced by beginJob().

bool l1t::L1TStage2CaloAnalyzer::doHistos_
private

Definition at line 71 of file L1TStage2CaloAnalyzer.cc.

bool l1t::L1TStage2CaloAnalyzer::doText_
private

Definition at line 70 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze().

TFileDirectory l1t::L1TStage2CaloAnalyzer::evtDispDir_
private

Definition at line 134 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH1F* > l1t::L1TStage2CaloAnalyzer::hbx_
private

Definition at line 127 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH1F* > l1t::L1TStage2CaloAnalyzer::hem_
private

Definition at line 128 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH1F* > l1t::L1TStage2CaloAnalyzer::het_
private

Definition at line 124 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH1F* > l1t::L1TStage2CaloAnalyzer::heta_
private

Definition at line 125 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH2F* > l1t::L1TStage2CaloAnalyzer::hetaphi_
private

Definition at line 131 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH1F* > l1t::L1TStage2CaloAnalyzer::hhad_
private

Definition at line 129 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH1F* > l1t::L1TStage2CaloAnalyzer::hiso_
private

Definition at line 132 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH1F* > l1t::L1TStage2CaloAnalyzer::hphi_
private

Definition at line 126 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::map< ObjectType, TH1F* > l1t::L1TStage2CaloAnalyzer::hratio_
private

Definition at line 130 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1F * l1t::L1TStage2CaloAnalyzer::hsort_
private

Definition at line 136 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1F* l1t::L1TStage2CaloAnalyzer::hsortMP_
private

Definition at line 136 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool l1t::L1TStage2CaloAnalyzer::m_allBx = false
private

Definition at line 140 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_clusterToken
private

Definition at line 49 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

int l1t::L1TStage2CaloAnalyzer::m_dmxBx = 0
private

Definition at line 139 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doClusters
private

Definition at line 60 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doEGs
private

Definition at line 65 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doEvtDisp = false
private

Definition at line 141 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), beginJob(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doJets
private

Definition at line 67 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doMPEGs
private

Definition at line 61 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doMPJets
private

Definition at line 63 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doMPSums
private

Definition at line 64 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doMPTaus
private

Definition at line 62 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doSums
private

Definition at line 68 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doTaus
private

Definition at line 66 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

bool l1t::L1TStage2CaloAnalyzer::m_doTowers
private

Definition at line 59 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_egToken
private

Definition at line 54 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_jetToken
private

Definition at line 56 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

int l1t::L1TStage2CaloAnalyzer::m_mpBx = 0
private

Definition at line 138 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_mpEGToken
private

Definition at line 50 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_mpJetToken
private

Definition at line 52 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_mpSumToken
private

Definition at line 53 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_mpTauToken
private

Definition at line 51 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_sumToken
private

Definition at line 57 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_tauToken
private

Definition at line 55 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

edm::EDGetToken l1t::L1TStage2CaloAnalyzer::m_towerToken
private

Definition at line 48 of file L1TStage2CaloAnalyzer.cc.

Referenced by analyze(), and L1TStage2CaloAnalyzer().

std::vector< ObjectType > l1t::L1TStage2CaloAnalyzer::types_
private

Definition at line 120 of file L1TStage2CaloAnalyzer.cc.

Referenced by beginJob(), and L1TStage2CaloAnalyzer().

std::vector< std::string > l1t::L1TStage2CaloAnalyzer::typeStr_
private

Definition at line 121 of file L1TStage2CaloAnalyzer.cc.

Referenced by beginJob(), and L1TStage2CaloAnalyzer().

Generated for CMSSW Reference Manual by doxygen 1.8.11