|
|
|
#include <src/L1RCTTestAnalyzer/src/L1RCTTestAnalyzer.cc>
Public Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
| L1RCTTestAnalyzer (const edm::ParameterSet &) | |
| ~L1RCTTestAnalyzer () | |
Static Private Member Functions | |
| static bool | compareEmCands (const L1CaloEmCand &cand1, const L1CaloEmCand &cand2) |
Private Attributes | |
| edm::InputTag | ecalDigisLabel |
| std::vector< int > | emIeta |
| std::vector< int > | emIphi |
| std::vector< int > | emIso |
| std::vector< int > | emRank |
| TTree * | emTree |
| TH1F * | h_ecalTimeSample |
| TH1F * | h_emCandTimeSample |
| TH1F * | h_emIeta |
| TH1F * | h_emIphi |
| TH1F * | h_emIso |
| TH2F * | h_emIsoInIetaIphi |
| TH2F * | h_emNonIsoInIetaIphi |
| TH1F * | h_emRank |
| TH2F * | h_emRankInIetaIphi |
| TH1F * | h_emRankOutOfTime |
| TH1F * | h_hcalTimeSample |
| TH2F * | h_regionFGInIetaIphi |
| TH1F * | h_regionIeta |
| TH1F * | h_regionIphi |
| TH1F * | h_regionMip |
| TH1F * | h_regionSum |
| TH2F * | h_regionSumInIetaIphi |
| TH1F * | h_towerMip |
| edm::InputTag | hcalDigisLabel |
| std::string | m_HepMCProduct |
| edm::InputTag | rctDigisLabel |
| bool | showEmCands |
| bool | showRegionSums |
Description: <one line="" class="" summary>="">
Implementation: <Notes on="" implementation>="">
Definition at line 56 of file L1RCTTestAnalyzer.h.
| L1RCTTestAnalyzer::L1RCTTestAnalyzer | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 27 of file L1RCTTestAnalyzer.cc.
References emIeta, emIphi, emIso, emRank, emTree, h_ecalTimeSample, h_emCandTimeSample, h_emIeta, h_emIphi, h_emIso, h_emIsoInIetaIphi, h_emNonIsoInIetaIphi, h_emRank, h_emRankInIetaIphi, h_emRankOutOfTime, h_hcalTimeSample, h_regionFGInIetaIphi, h_regionIeta, h_regionIphi, h_regionMip, h_regionSum, h_regionSumInIetaIphi, and h_towerMip.
: showEmCands(iConfig.getUntrackedParameter<bool>("showEmCands")), showRegionSums(iConfig.getUntrackedParameter<bool>("showRegionSums")), ecalDigisLabel(iConfig.getParameter<edm::InputTag>("ecalDigisLabel")), hcalDigisLabel(iConfig.getParameter<edm::InputTag>("hcalDigisLabel")), rctDigisLabel(iConfig.getParameter<edm::InputTag>("rctDigisLabel")) { //now do what ever initialization is needed edm::Service<TFileService> fs; emTree = fs->make<TTree>("emTree","L1 RCT EM tree"); // emTree->Branch("emRank",emRank,"emRank/I"); // emTree->Branch("emIeta",emIeta,"emIeta/I"); // emTree->Branch("emIphi",emIphi,"emIphi/I"); // emTree->Branch("emIso" ,emIso ,"emIso/I"); emTree->Branch("emRank",&emRank); emTree->Branch("emIeta",&emIeta); emTree->Branch("emIphi",&emIphi); emTree->Branch("emIso" ,&emIso); h_emRank = fs->make<TH1F>( "emRank", "emRank", 64, 0., 64. ); h_emRankOutOfTime = fs->make<TH1F>( "emRankOutOfTime", "emRankOutOfTime", 64, 0., 64. ); h_emIeta = fs->make<TH1F>( "emIeta", "emIeta", 22, 0., 22. ); h_emIphi = fs->make<TH1F>( "emIphi", "emIphi", 18, 0., 18. ); h_emIso = fs->make<TH1F>( "emIso", "emIso", 2, 0., 2. ); h_emRankInIetaIphi = fs->make<TH2F>( "emRank2D", "emRank2D", 22, 0., 22., 18, 0., 18. ); h_emIsoInIetaIphi = fs->make<TH2F>( "emIso2D", "emIso2D", 22, 0., 22., 18, 0., 18. ); h_emNonIsoInIetaIphi = fs->make<TH2F>( "emNonIso2D", "emNonIso2D", 22, 0., 22., 18, 0., 18. ); h_emCandTimeSample = fs->make<TH1F>( "emCandTimeSample", "emCandTimeSample", 5, -2., 2.); h_regionSum = fs->make<TH1F>( "regionSum", "regionSum", 100, 0., 100. ); h_regionIeta = fs->make<TH1F>( "regionIeta", "regionIeta", 22, 0., 22. ); h_regionIphi = fs->make<TH1F>( "regionIphi", "regionIphi", 18, 0., 18. ); h_regionMip = fs->make<TH1F>( "regionMip", "regionMipBit", 2, 0., 2. ); h_regionSumInIetaIphi = fs->make<TH2F>( "regionSum2D", "regionSum2D", 22, 0., 22., 18, 0., 18. ); h_regionFGInIetaIphi = fs->make<TH2F>( "regionFG2D", "regionFG2D", 22, 0., 22., 18, 0., 18. ); h_towerMip = fs->make<TH1F>( "towerMip", "towerMipBit", 2, 0., 2. ); h_ecalTimeSample = fs->make<TH1F>( "ecalTimeSample", "ecalTimeSample", 10, 0., 10. ); h_hcalTimeSample = fs->make<TH1F>( "hcalTimeSample", "hcalTimeSample", 10, 0., 10. ); // get names of modules, producing object collections }
| L1RCTTestAnalyzer::~L1RCTTestAnalyzer | ( | ) |
Definition at line 83 of file L1RCTTestAnalyzer.cc.
{
// do anything here that needs to be done at destruction time
// (e.g. close files, deallocate resources etc.)
}
| void L1RCTTestAnalyzer::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 98 of file L1RCTTestAnalyzer.cc.
References compareEmCands(), gather_cfg::cout, patCandidatesForDimuonsSequences_cff::ecal, ecalDigisLabel, emIeta, emIphi, emIso, emRank, emTree, edm::EventSetup::get(), edm::Event::getByLabel(), h_ecalTimeSample, h_emCandTimeSample, h_emIeta, h_emIphi, h_emIso, h_emIsoInIetaIphi, h_emNonIsoInIetaIphi, h_emRank, h_emRankInIetaIphi, h_emRankOutOfTime, h_hcalTimeSample, h_regionFGInIetaIphi, h_regionIeta, h_regionIphi, h_regionMip, h_regionSum, h_regionSumInIetaIphi, h_towerMip, patCandidatesForDimuonsSequences_cff::hcal, hcalDigisLabel, rctDigisLabel, compare_using_db::sample, showEmCands, showRegionSums, and python::multivaluedict::sort().
{
using namespace edm;
#ifdef THIS_IS_AN_EVENT_EXAMPLE
Handle<ExampleData> pIn;
iEvent.getByLabel("example",pIn);
#endif
#ifdef THIS_IS_AN_EVENTSETUP_EXAMPLE
ESHandle<SetupData> pSetup;
iSetup.get<SetupRecord>().get(pSetup);
#endif
// as in L1GctTestAnalyzer.cc
Handle<L1CaloEmCollection> rctEmCands;
Handle<L1CaloRegionCollection> rctRegions;
Handle<EcalTrigPrimDigiCollection> ecalColl;
Handle<HcalTrigPrimDigiCollection> hcalColl;
L1CaloEmCollection::const_iterator em;
L1CaloRegionCollection::const_iterator rgn;
EcalTrigPrimDigiCollection::const_iterator ecal;
HcalTrigPrimDigiCollection::const_iterator hcal;
iEvent.getByLabel(rctDigisLabel, rctEmCands);
iEvent.getByLabel(rctDigisLabel, rctRegions);
iEvent.getByLabel(ecalDigisLabel, ecalColl);
iEvent.getByLabel(hcalDigisLabel, hcalColl);
// for sorting later
L1CaloEmCollection * myL1EmColl = new L1CaloEmCollection;
for (ecal=ecalColl->begin(); ecal!=ecalColl->end(); ecal++)
{
for (unsigned short sample = 0; sample < (*ecal).size(); sample++)
{
h_ecalTimeSample->Fill(sample);
}
}
for (hcal=hcalColl->begin(); hcal!=hcalColl->end(); hcal++)
{
h_towerMip->Fill( (*hcal).SOI_fineGrain() );
for (unsigned short sample = 0; sample < (*hcal).size(); sample++)
{
h_hcalTimeSample->Fill(sample);
}
}
if(showEmCands)
{
// std::cout << std::endl << "L1 RCT EmCand objects" << std::endl;
}
for (em=rctEmCands->begin(); em!=rctEmCands->end(); em++){
// std::cout << "(Analyzer)\n" << (*em) << std::endl;
L1CaloEmCand * myL1EmCand = new L1CaloEmCand(*em);
(*myL1EmColl).push_back(*myL1EmCand);
delete myL1EmCand;
h_emCandTimeSample->Fill((*em).bx());
if ((*em).bx() == 0)
{
unsigned short n_emcands = 0;
//std::cout << std::endl << "rank: " << (*em).rank() ;
if ((*em).rank() > 0)
{
h_emRank->Fill( (*em).rank() );
h_emIeta->Fill( (*em).regionId().ieta() );
h_emIphi->Fill( (*em).regionId().iphi() );
h_emIso->Fill( (*em).isolated() );
h_emRankInIetaIphi->Fill( (*em).regionId().ieta(),
(*em).regionId().iphi(),
(*em).rank() );
if ((*em).isolated())
{
h_emIsoInIetaIphi->Fill( (*em).regionId().ieta(),
(*em).regionId().iphi() );
}
else
{
h_emNonIsoInIetaIphi->Fill( (*em).regionId().ieta(),
(*em).regionId().iphi() );
}
}
if (showEmCands)
{
if ((*em).rank() > 0)
{
// std::cout << std::endl << "rank: " << (*em).rank();
unsigned short rgnPhi = 999;
unsigned short rgn = (unsigned short) (*em).rctRegion();
unsigned short card = (unsigned short) (*em).rctCard();
unsigned short crate = (unsigned short) (*em).rctCrate();
if (card == 6)
{
rgnPhi = rgn;
}
else if (card < 6)
{
rgnPhi = (card % 2);
}
else
{
std::cout << "rgnPhi not assigned (still " << rgnPhi << ") -- Weird card number! " << card ;
}
//unsigned short phi_bin = ((crate % 9) * 2) + rgnPhi;
short eta_bin = (card/2) * 2 + 1;
if (card < 6)
{
eta_bin = eta_bin + rgn;
}
if (crate < 9)
{
eta_bin = -eta_bin;
}
n_emcands++;
// std::cout << /* "rank: " << (*em).rank() << */ " eta_bin: " << eta_bin << " phi_bin: " << phi_bin << ". crate: " << crate << " card: " << card << " region: " << rgn << ". isolated: " << (*em).isolated();
}
}
}
else
{
h_emRankOutOfTime->Fill( (*em).rank() );
}
}
if(showEmCands)
{
// std::cout << std::endl;
}
// next: SORT THESE GUYS so they're entered into the tree highest first
// std::sort(rctEmCands->begin(),rctEmCands->end(),compareEmCands);
// for (em=rctEmCands->begin(); em!=rctEmCands->end(); em++)
std::sort(myL1EmColl->begin(),myL1EmColl->end(),compareEmCands);
std::reverse(myL1EmColl->begin(),myL1EmColl->end()); // whoops!
for (em=myL1EmColl->begin(); em!=myL1EmColl->end(); em++)
{
emRank.push_back( (*em).rank() );
emIeta.push_back( (*em).regionId().ieta() );
emIphi.push_back( (*em).regionId().iphi() );
emIso.push_back( (*em).isolated() );
}
emTree->Fill();
if(showRegionSums)
{
std::cout << "Regions" << std::endl;
}
for (rgn=rctRegions->begin(); rgn!=rctRegions->end(); rgn++)
{
if ((*rgn).bx() == 0)
{
if(showRegionSums&&(*rgn).et()>0)
{
std::cout << /* "(Analyzer)\n" << */ (*rgn) << std::endl;
}
if ( (*rgn).et() > 0 )
{
h_regionSum->Fill( (*rgn).et() );
h_regionIeta->Fill( (*rgn).gctEta() );
h_regionIphi->Fill( (*rgn).gctPhi() );
h_regionSumInIetaIphi->Fill( (*rgn).gctEta(), (*rgn).gctPhi(),
(*rgn).et() );
h_regionFGInIetaIphi->Fill( (*rgn).gctEta(), (*rgn).gctPhi(),
(*rgn).fineGrain() );
}
h_regionMip->Fill( (*rgn).mip() );
}
}
if(showRegionSums)
{
std::cout << std::endl;
}
emRank.clear();
emIeta.clear();
emIphi.clear();
emIso.clear();
delete myL1EmColl;
}
| bool L1RCTTestAnalyzer::compareEmCands | ( | const L1CaloEmCand & | cand1, |
| const L1CaloEmCand & | cand2 | ||
| ) | [static, private] |
Definition at line 288 of file L1RCTTestAnalyzer.cc.
References L1CaloEmCand::rank().
Referenced by analyze().
Definition at line 70 of file L1RCTTestAnalyzer.h.
Referenced by analyze().
std::vector<int> L1RCTTestAnalyzer::emIeta [private] |
Definition at line 80 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
std::vector<int> L1RCTTestAnalyzer::emIphi [private] |
Definition at line 81 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
std::vector<int> L1RCTTestAnalyzer::emIso [private] |
Definition at line 82 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
std::vector<int> L1RCTTestAnalyzer::emRank [private] |
Definition at line 79 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TTree* L1RCTTestAnalyzer::emTree [private] |
Definition at line 74 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_ecalTimeSample [private] |
Definition at line 106 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_emCandTimeSample [private] |
Definition at line 93 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_emIeta [private] |
Definition at line 86 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_emIphi [private] |
Definition at line 87 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_emIso [private] |
Definition at line 88 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH2F* L1RCTTestAnalyzer::h_emIsoInIetaIphi [private] |
Definition at line 91 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH2F* L1RCTTestAnalyzer::h_emNonIsoInIetaIphi [private] |
Definition at line 92 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_emRank [private] |
Definition at line 84 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH2F* L1RCTTestAnalyzer::h_emRankInIetaIphi [private] |
Definition at line 89 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_emRankOutOfTime [private] |
Definition at line 85 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_hcalTimeSample [private] |
Definition at line 107 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH2F* L1RCTTestAnalyzer::h_regionFGInIetaIphi [private] |
Definition at line 102 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_regionIeta [private] |
Definition at line 96 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_regionIphi [private] |
Definition at line 97 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_regionMip [private] |
Definition at line 98 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_regionSum [private] |
Definition at line 95 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH2F* L1RCTTestAnalyzer::h_regionSumInIetaIphi [private] |
Definition at line 99 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
TH1F* L1RCTTestAnalyzer::h_towerMip [private] |
Definition at line 104 of file L1RCTTestAnalyzer.h.
Referenced by analyze(), and L1RCTTestAnalyzer().
Definition at line 71 of file L1RCTTestAnalyzer.h.
Referenced by analyze().
std::string L1RCTTestAnalyzer::m_HepMCProduct [private] |
Definition at line 67 of file L1RCTTestAnalyzer.h.
Definition at line 72 of file L1RCTTestAnalyzer.h.
Referenced by analyze().
bool L1RCTTestAnalyzer::showEmCands [private] |
Definition at line 68 of file L1RCTTestAnalyzer.h.
Referenced by analyze().
bool L1RCTTestAnalyzer::showRegionSums [private] |
Definition at line 69 of file L1RCTTestAnalyzer.h.
Referenced by analyze().
1.7.3