|
|
|
#include <HSCPValidation/HSCPValidator/src/HSCPValidator.cc>
Description: [one line class summary]
Implementation: [Notes on implementation]
Definition at line 41 of file HSCPValidator.h.
| HSCPValidator::HSCPValidator | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 98 of file HSCPValidator.cc.
References cluSizeDistribution_, dedxVsp, digiOccupancyMapEB_, digiOccupancyMapEEM_, digiOccupancyMapEEP_, digisEcalNumHistEB_, digisEcalNumHistEE_, efficiencyRPCRecHitSimDigis_, h_genhscp_met, h_genhscp_met_nohscp, h_genhscp_scaloret, h_genhscp_scaloret_nohscp, hltjet, hltmet, hltmu, particleBetaHist_, particleBetaInverseHist_, particleEtaHist_, particleMassHist_, particlePhiHist_, particlePHist_, particlePtHist_, particleStatusHist_, RecoHSCPPtVsGenPt, residualsRPCRecHitSimDigis_, rpcBXBarrel_, rpcBXEndCap_, rpcTimeOfFlightBarrel_, rpcTimeOfFlightEndCap_, simHitsEcalDigiMatchEnergyHistEB_, simHitsEcalDigiMatchEnergyHistEE_, simHitsEcalDigiMatchEnergyVsTimeHistEB_, simHitsEcalDigiMatchEnergyVsTimeHistEE_, simHitsEcalDigiMatchIEtaHist_, simHitsEcalDigiMatchIPhiHist_, simHitsEcalDigiMatchTimeHistEB_, simHitsEcalDigiMatchTimeHistEE_, simHitsEcalEnergyHistEB_, simHitsEcalEnergyHistEE_, simHitsEcalEnergyVsTimeHistEB_, simHitsEcalEnergyVsTimeHistEE_, simHitsEcalNumHistEB_, simHitsEcalNumHistEE_, simHitsEcalTimeHistEB_, simHitsEcalTimeHistEE_, simTrackParticleBetaHist_, simTrackParticleEtaHist_, simTrackParticlePhiHist_, simTrackParticlePHist_, and simTrackParticlePtHist_.
: doGenPlots_ (iConfig.getParameter<bool>("MakeGenPlots")), doHLTPlots_ (iConfig.getParameter<bool>("MakeHLTPlots")), doSimTrackPlots_ (iConfig.getParameter<bool>("MakeSimTrackPlots")), doSimDigiPlots_ (iConfig.getParameter<bool>("MakeSimDigiPlots")), doRecoPlots_ (iConfig.getParameter<bool>("MakeRecoPlots")), label_ (iConfig.getParameter<edm::InputTag>("generatorLabel")), particleIds_ (iConfig.getParameter< std::vector<int> >("particleIds")), particleStatus_ (iConfig.getUntrackedParameter<int>("particleStatus",1)), ebSimHitTag_ (iConfig.getParameter<edm::InputTag>("EBSimHitCollection")), eeSimHitTag_ (iConfig.getParameter<edm::InputTag>("EESimHitCollection")), simTrackTag_ (iConfig.getParameter<edm::InputTag>("SimTrackCollection")), EBDigiCollection_ (iConfig.getParameter<edm::InputTag>("EBDigiCollection")), EEDigiCollection_ (iConfig.getParameter<edm::InputTag>("EEDigiCollection")), RPCRecHitTag_ (iConfig.getParameter<edm::InputTag>("RPCRecHitTag")) { //now do what ever initialization is needed // GEN particleEtaHist_ = fileService->make<TH1F>("particleEta","Eta of gen particle",100,-5,5); particlePhiHist_ = fileService->make<TH1F>("particlePhi","Phi of gen particle",180,-3.15,3.15); particlePHist_ = fileService->make<TH1F>("particleP","Momentum of gen particle",500,0,2000); particlePtHist_ = fileService->make<TH1F>("particlePt","P_{T} of gen particle",500,0,2000); particleMassHist_ = fileService->make<TH1F>("particleMass","Mass of gen particle",1000,0,2000); particleStatusHist_ = fileService->make<TH1F>("particleStatus","Status of gen particle",10,0,10); particleBetaHist_ = fileService->make<TH1F>("particleBeta","Beta of gen particle",100,0,1); particleBetaInverseHist_ = fileService->make<TH1F>("particleBetaInverse","1/#beta of gen particle",100,0,5); h_genhscp_met = fileService->make<TH1F>( "hscp_met" , "missing E_{T} hscp" , 100, 0., 1500. ); h_genhscp_met_nohscp = fileService->make<TH1F>( "hscp_met_nohscp" , "missing E_{T} w/o hscp" , 100, 0., 1500. ); h_genhscp_scaloret = fileService->make<TH1F>( "hscp_scaloret" , "scalor E_{T} sum" , 100, 0., 1500. ); h_genhscp_scaloret_nohscp = fileService->make<TH1F>( "hscp_scaloret_nohscp" , "scalor E_{T} sum w/o hscp" , 100, 0., 1500. ); //SIM track Info simTrackParticleEtaHist_ = fileService->make<TH1F>("simTrackParticleEta","Eta of simTrackParticle",100,-5,5); simTrackParticlePhiHist_ = fileService->make<TH1F>("simTrackParticlePhi","Phi of simTrackParticle",180,-3.15,3.15); simTrackParticlePHist_ = fileService->make<TH1F>("simTrackParticleP","Momentum of simTrackParticle",500,0,2000); simTrackParticlePtHist_ = fileService->make<TH1F>("simTrackParticlePt","P_{T} of simTrackParticle",500,0,2000); simTrackParticleBetaHist_ = fileService->make<TH1F>("simTrackParticleBeta","Beta of simTrackParticle",100,0,1); //reco track Info RecoHSCPPtVsGenPt= fileService->make<TH2F>("Recovsgenpt","RecovsGen",100,0,1000,100,0,1000); dedxVsp = fileService->make<TH2F>("dedxvsp","dedxvsp",100,0,1000,100,0,10); //HLT Info hltmet = fileService->make<TH1F>("HLT_MET","MET",3,-1,2); hltjet = fileService->make<TH1F>("HLT_JET","JET",3,-1,2); hltmu = fileService->make<TH1F>("HLT_Mu","Mu",3,-1,2); // SIM-DIGI: ECAL simHitsEcalEnergyHistEB_ = fileService->make<TH1F>("ecalEnergyOfSimHitsEB","HSCP SimTrack-matching SimHit energy EB [GeV]",125,-1,4); simHitsEcalEnergyHistEE_ = fileService->make<TH1F>("ecalEnergyOfSimHitsEE","HSCP SimTrack-matching SimHit energy EE [GeV]",125,-1,4); simHitsEcalTimeHistEB_ = fileService->make<TH1F>("ecalTimingOfSimHitsEB","HSCP SimTrack-matching SimHit time EB [ns]",115,-15,100); simHitsEcalTimeHistEE_ = fileService->make<TH1F>("ecalTimingOfSimHitsEE","HSCP SimTrack-matching SimHit time EE [ns]",115,-15,100); simHitsEcalNumHistEB_ = fileService->make<TH1F>("ecalNumberOfSimHitsEB","Number of HSCP SimTrack-matching EB sim hits in event",100,0,200); simHitsEcalNumHistEE_ = fileService->make<TH1F>("ecalNumberOfSimHitsEE","Number of HSCP SimTrack-matching EE sim hits in event",100,0,200); simHitsEcalEnergyVsTimeHistEB_ = fileService->make<TH2F>("ecalEnergyVsTimeOfSimHitsEB","Energy vs. time of HSCP SimTrack-matching EB sim hits in event",115,-15,100,125,-1,4); simHitsEcalEnergyVsTimeHistEE_ = fileService->make<TH2F>("ecalEnergyVsTimeOfSimHitsEE","Energy vs. time of HSCP SimTrack-matching EE sim hits in event",115,-15,100,125,-1,4); simHitsEcalDigiMatchEnergyHistEB_ = fileService->make<TH1F>("ecalEnergyOfDigiMatSimHitsEB","HSCP digi-matching SimHit energy EB [GeV]",125,-1,4); simHitsEcalDigiMatchEnergyHistEE_ = fileService->make<TH1F>("ecalEnergyOfDigiMatSimHitsEE","HSCP digi-matching SimHit energy EE [GeV]",125,-1,4); simHitsEcalDigiMatchTimeHistEB_ = fileService->make<TH1F>("ecalTimingOfDigiMatSimHitsEB","HSCP digi-matching SimHit time EB [ns]",115,-15,100); simHitsEcalDigiMatchTimeHistEE_ = fileService->make<TH1F>("ecalTimingOfDigiMatSimHitsEE","HSCP digi-matching SimHit time EE [ns]",115,-15,100); simHitsEcalDigiMatchEnergyVsTimeHistEB_ = fileService->make<TH2F>("ecalEnergyVsTimeOfDigiMatSimHitsEB","HSCP digi-matching EB SimHit energy vs. time",115,-15,100,125,-1,4); simHitsEcalDigiMatchEnergyVsTimeHistEE_ = fileService->make<TH2F>("ecalEnergyVsTimeOfDigiMatSimHitsEE","HSCP digi-matching EE SimHit energy vs. time",115,-15,100,125,-1,4); simHitsEcalDigiMatchIEtaHist_ = fileService->make<TH1F>("ecalIEtaOfDigiMatchSimHits","iEta of digi-matching Ecal simHits (EB)",171,-85,86); simHitsEcalDigiMatchIPhiHist_ = fileService->make<TH1F>("ecalIPhiOfDigiMatchSimHits","iPhi of digi-matching Ecal simHits (EB)",360,1,361); digisEcalNumHistEB_ = fileService->make<TH1F>("ecalDigisNumberEB","Number of EB digis matching simhits in event",200,0,1000); digisEcalNumHistEE_ = fileService->make<TH1F>("ecalDigisNumberEE","Number of EE digis matching simhits in event",200,0,1000); digiOccupancyMapEB_ = fileService->make<TH2F>("ecalDigiOccupancyMapEB","Occupancy of simhit-matching digis EB;i#phi;i#eta",360,1,361,171,-85,86); digiOccupancyMapEEP_ = fileService->make<TH2F>("ecalDigiOccupancyMapEEM","Occupancy of simhit-matching digis EEM;ix;iy",100,1,100,100,1,100); digiOccupancyMapEEM_ = fileService->make<TH2F>("ecalDigiOccupancyMapEEP","Occupancy of simhit-matching digis EEP;ix;iy",100,1,100,100,1,100); // SIM-DIGI: RPC residualsRPCRecHitSimDigis_ = fileService->make<TH1F>("residualsRPCRecHitSimDigis","HSCP SimHit - Clossest RPC RecHit",100,-5,5); efficiencyRPCRecHitSimDigis_ = fileService->make<TH1F>("efficiencyRPCRecHitSimDigis","HSCP SimHits RecHits Efficiency",2,-0.5,1.5); cluSizeDistribution_ = fileService->make<TH1F>("RPCCluSizeDistro","RPC HSCP CluSize Distribution",11,-0.5,10.5); rpcTimeOfFlightBarrel_[0] = fileService->make<TH1F>("RPCToFLayer1","RPC HSCP Time Of Flight Layer 1",50,5,100); rpcTimeOfFlightBarrel_[1] = fileService->make<TH1F>("RPCToFLayer2","RPC HSCP Time Of Flight Layer 2",50,5,100); rpcTimeOfFlightBarrel_[2] = fileService->make<TH1F>("RPCToFLayer3","RPC HSCP Time Of Flight Layer 3",50,5,100); rpcTimeOfFlightBarrel_[3] = fileService->make<TH1F>("RPCToFLayer4","RPC HSCP Time Of Flight Layer 4",50,5,100); rpcTimeOfFlightBarrel_[4] = fileService->make<TH1F>("RPCToFLayer5","RPC HSCP Time Of Flight Layer 5",50,5,100); rpcTimeOfFlightBarrel_[5] = fileService->make<TH1F>("RPCToFLayer6","RPC HSCP Time Of Flight Layer 6",50,5,100); rpcBXBarrel_[0] = fileService->make<TH1F>("RPCBXLayer1","RPC HSCP BX Layer 1",5,-0.5,4.5); rpcBXBarrel_[1] = fileService->make<TH1F>("RPCBXLayer2","RPC HSCP BX Layer 2",5,-0.5,4.5); rpcBXBarrel_[2] = fileService->make<TH1F>("RPCBXLayer3","RPC HSCP BX Layer 3",5,-0.5,4.5); rpcBXBarrel_[3] = fileService->make<TH1F>("RPCBXLayer4","RPC HSCP BX Layer 4",5,-0.5,4.5); rpcBXBarrel_[4] = fileService->make<TH1F>("RPCBXLayer5","RPC HSCP BX Layer 5",5,-0.5,4.5); rpcBXBarrel_[5] = fileService->make<TH1F>("RPCBXLayer6","RPC HSCP BX Layer 6",5,-0.5,4.5); rpcTimeOfFlightEndCap_[0]= fileService->make<TH1F>("RPCToFDisk1","RPC HSCP Time Of Flight Disk 1",50,5,100); rpcTimeOfFlightEndCap_[1]= fileService->make<TH1F>("RPCToFDisk2","RPC HSCP Time Of Flight Disk 2",50,5,100); rpcTimeOfFlightEndCap_[2]= fileService->make<TH1F>("RPCToFDisk3","RPC HSCP Time Of Flight Disk 3",50,5,100); rpcBXEndCap_[0] = fileService->make<TH1F>("RPCBXDisk1","RPC HSCP BX Disk 1",5,-0.5,4.5); rpcBXEndCap_[1] = fileService->make<TH1F>("RPCBXDisk2","RPC HSCP BX Disk 2",5,-0.5,4.5); rpcBXEndCap_[2] = fileService->make<TH1F>("RPCBXDisk3","RPC HSCP BX Disk 3",5,-0.5,4.5); }
| HSCPValidator::~HSCPValidator | ( | ) |
Definition at line 197 of file HSCPValidator.cc.
{
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
// particleEtaHist_ = fileService->make<TH1F>("particleEta","Eta of gen particle",100,-5,5);
// particlePhiHist_ = fileService->make<TH1F>("particlePhi","Phi of gen particle",180,-3.15,3.15);
// particlePHist_ = fileService->make<TH1F>("particleP","Momentum of gen particle",500,0,2000);
// particlePtHist_ = fileService->make<TH1F>("particlePt","P_{T} of gen particle",500,0,2000);
// particleMassHist_ = fileService->make<TH1F>("particleMass","Mass of gen particle",1000,0,2000);
// particleStatusHist_ = fileService->make<TH1F>("particleStatus","Status of gen particle",10,0,10);
// particleBetaHist_ = fileService->make<TH1F>("particleBeta","Beta of gen particle",100,0,1);
// particleBetaInverseHist_ = fileService->make<TH1F>("particleBetaInverse","1/#beta of gen particle",100,0,5);
}
| void HSCPValidator::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 220 of file HSCPValidator.cc.
References doGenPlots_, doHLTPlots_, doRecoPlots_, doSimDigiPlots_, doSimTrackPlots_, edm::EventSetup::get(), makeGenPlots(), makeHLTPlots(), makeRecoPlots(), makeSimDigiPlotsECAL(), makeSimDigiPlotsRPC(), makeSimTrackPlots(), and rpcGeo.
{
using namespace edm;
iSetup.get<MuonGeometryRecord>().get(rpcGeo);
if(doGenPlots_)
makeGenPlots(iEvent);
if(doSimTrackPlots_)
makeSimTrackPlots(iEvent);
if(doSimDigiPlots_){
makeSimDigiPlotsECAL(iEvent);
makeSimDigiPlotsRPC(iEvent);
}
if(doHLTPlots_){
makeHLTPlots(iEvent);
}
if(doRecoPlots_){
makeRecoPlots(iEvent);
}
}
| void HSCPValidator::beginJob | ( | void | ) | [private, virtual] |
| void HSCPValidator::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 251 of file HSCPValidator.cc.
References gather_cfg::cout, intToString(), and particleIdsFoundMap_.
{
std::string frequencies = "";
for(std::map<int,int>::const_iterator itr = particleIdsFoundMap_.begin();
itr != particleIdsFoundMap_.end(); ++itr)
{
frequencies+="PDG ID: ";
frequencies+=intToString(itr->first);
frequencies+=" Frequency: ";
frequencies+=intToString(itr->second);
frequencies+="\n";
}
std::cout << "Found PDGIds: " << "\n\n" << frequencies << std::endl;
}
| bool HSCPValidator::IncreasedTreshold | ( | const trigger::TriggerEvent & | trEv, |
| const edm::InputTag & | InputPath, | ||
| double | NewThreshold, | ||
| double | etaCut, | ||
| int | NObjectAboveThreshold, | ||
| bool | averageThreshold | ||
| ) | [private] |
Definition at line 837 of file HSCPValidator.cc.
References eta(), trigger::TriggerEvent::filterIds(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), i, max(), n, trigger::TriggerEvent::sizeFilters(), and python::multivaluedict::sort().
Referenced by makeHLTPlots().
{
unsigned int filterIndex = trEv.filterIndex(InputPath);
//if(filterIndex<trEv.sizeFilters())printf("SELECTED INDEX =%i --> %s XXX %s\n",filterIndex,trEv.filterTag(filterIndex).label().c_str(), trEv.filterTag(filterIndex).process().c_str());
if (filterIndex<trEv.sizeFilters()){
const trigger::Vids& VIDS(trEv.filterIds(filterIndex));
const trigger::Keys& KEYS(trEv.filterKeys(filterIndex));
const int nI(VIDS.size());
const int nK(KEYS.size());
assert(nI==nK);
const int n(std::max(nI,nK));
const trigger::TriggerObjectCollection& TOC(trEv.getObjects());
if(!averageThreshold){
int NObjectAboveThresholdObserved = 0;
for (int i=0; i!=n; ++i) {
if(TOC[KEYS[i]].pt()> NewThreshold && fabs(TOC[KEYS[i]].eta())<etaCut) NObjectAboveThresholdObserved++;
//cout << " " << i << " " << VIDS[i] << "/" << KEYS[i] << ": "<< TOC[KEYS[i]].id() << " " << TOC[KEYS[i]].pt() << " " << TOC[KEYS[i]].eta() << " " << TOC[KEYS[i]].phi() << " " << TOC[KEYS[i]].mass()<< endl;
}
if(NObjectAboveThresholdObserved>=NObjectAboveThreshold)return true;
}else{
std::vector<double> ObjPt;
for (int i=0; i!=n; ++i) {
ObjPt.push_back(TOC[KEYS[i]].pt());
//cout << " " << i << " " << VIDS[i] << "/" << KEYS[i] << ": "<< TOC[KEYS[i]].id() << " " << TOC[KEYS[i]].pt() << " " << TOC[KEYS[i]].eta() << " " << TOC[KEYS[i]].phi() << " " << TOC[KEYS[i]].mass()<< endl;
}
if((int)(ObjPt.size())<NObjectAboveThreshold)return false;
std::sort(ObjPt.begin(), ObjPt.end());
double Average = 0;
for(int i=0; i<NObjectAboveThreshold;i++){
Average+= ObjPt[ObjPt.size()-1-i];
}Average/=NObjectAboveThreshold;
//cout << "AVERAGE = " << Average << endl;
if(Average>NewThreshold)return true;
}
}
return false;
}
| std::string HSCPValidator::intToString | ( | int | num | ) | [private] |
Definition at line 824 of file HSCPValidator.cc.
Referenced by endJob().
{
using namespace std;
ostringstream myStream;
myStream << num << flush;
return(myStream.str()); //returns the string form of the stringstream object
}
| void HSCPValidator::makeGenPlots | ( | const edm::Event & | iEvent | ) | [private] |
Definition at line 268 of file HSCPValidator.cc.
References abs, spr::find(), edm::Event::getByLabel(), h_genhscp_met, h_genhscp_met_nohscp, h_genhscp_scaloret, h_genhscp_scaloret_nohscp, label_, mag(), AlCaHLTBitMon_ParallelJobs::p, particleBetaHist_, particleBetaInverseHist_, particleEtaHist_, particleIds_, particleIdsFoundMap_, particleMassHist_, particlePhiHist_, particlePHist_, particlePtHist_, particleStatus_, particleStatusHist_, funct::pow(), and mathSSE::sqrt().
Referenced by analyze().
{
using namespace edm;
double missingpx=0;
double missingpy=0;
double missingpx_nohscp=0;
double missingpy_nohscp=0;
double scalorEt=0;
double scalorEt_nohscp=0;
Handle<HepMCProduct> evt;
iEvent.getByLabel(label_, evt);
HepMC::GenEvent * myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
for(HepMC::GenEvent::particle_iterator p = myGenEvent->particles_begin();
p != myGenEvent->particles_end(); ++p )
{
if((*p)->status() != particleStatus_)
continue;
//calculate MET(neutrino as MET)
if(abs((*p)->pdg_id())!=12 && abs((*p)->pdg_id())!=14 && abs((*p)->pdg_id())!=16){ //for non-neutrino particles.
missingpx-=(*p)->momentum().px();
missingpy-=(*p)->momentum().py();
scalorEt+=(*p)->momentum().perp();
}
// Check if the particleId is in our R-hadron list
std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(),particleIds_.end(),(*p)->pdg_id());
if(partIdItr==particleIds_.end()){
//calculate MET(neutrino+ HSCP as MET)
if(abs((*p)->pdg_id())!=12 && abs((*p)->pdg_id())!=14 && abs((*p)->pdg_id())!=16){ //for non-neutrino particles.
missingpx_nohscp-=(*p)->momentum().px();
missingpy_nohscp-=(*p)->momentum().py();
scalorEt_nohscp+=(*p)->momentum().perp();
}
}
else{
particleStatusHist_->Fill((*p)->status());
std::pair<std::map<int,int>::iterator,bool> pair = particleIdsFoundMap_.insert(std::make_pair<int,int>((*p)->pdg_id(),1));
if(!pair.second)
{
++(pair.first->second);
}
double mag = sqrt(pow((*p)->momentum().px(),2) + pow((*p)->momentum().py(),2) + pow((*p)->momentum().pz(),2) );
particleEtaHist_->Fill((*p)->momentum().eta());
particlePhiHist_->Fill((*p)->momentum().phi());
particlePHist_->Fill(mag);
particlePtHist_->Fill((*p)->momentum().perp());
particleMassHist_->Fill((*p)->generated_mass());
float particleP = mag;
float particleM = (*p)->generated_mass();
particleBetaHist_->Fill(particleP/sqrt(particleP*particleP+particleM*particleM));
particleBetaInverseHist_->Fill(sqrt(particleP*particleP+particleM*particleM)/particleP);
}
}
h_genhscp_met->Fill(sqrt(missingpx*missingpx+missingpy*missingpy));
h_genhscp_met_nohscp->Fill(sqrt(missingpx_nohscp*missingpx_nohscp+missingpy_nohscp*missingpy_nohscp));
h_genhscp_scaloret->Fill(scalorEt);
h_genhscp_scaloret_nohscp->Fill(scalorEt_nohscp);
delete myGenEvent;
}
| void HSCPValidator::makeHLTPlots | ( | const edm::Event & | iEvent | ) | [private] |
Definition at line 372 of file HSCPValidator.cc.
References gather_cfg::cout, edm::Event::getByLabel(), hltjet, hltmet, hltmu, i, IncreasedTreshold(), and edm::Event::triggerResultsByName().
Referenced by analyze().
{
using namespace edm;
//get HLT infos
edm::TriggerResultsByName tr = iEvent.triggerResultsByName("HLT");
if(!tr.isValid()){
std::cout<<"Tirgger Results not available"<<std::endl;
}
edm::Handle< trigger::TriggerEvent > trEvHandle;
iEvent.getByLabel("hltTriggerSummaryAOD", trEvHandle);
trigger::TriggerEvent trEv = *trEvHandle;
unsigned int TrIndex_Unknown = tr.size();
// HLT TRIGGER BASED ON 1 MUON!
if(TrIndex_Unknown != tr.triggerIndex("HLT_Mu40_v1")){
if(tr.accept(tr.triggerIndex("HLT_Mu40_v1"))) hltmu->Fill(1);
else {hltmu->Fill(0);}
}
else{
if(TrIndex_Unknown != tr.triggerIndex("HLT_Mu30_v1")){
if(IncreasedTreshold(trEv, InputTag("hltSingleMu30L3Filtered30","","HLT"), 40,2.1, 1, false)) hltmu->Fill(1);
else hltmu->Fill(0);
}else{
printf("BUG with HLT_Mu\n");
std::cout<<"trigger names are : ";
for(unsigned int i=0;i<tr.size();i++){
std::cout<<" "<<tr.triggerName(i);
}
std::cout<<std::endl;
}
}
// HLT TRIGGER BASED ON MET!
if(TrIndex_Unknown != tr.triggerIndex("HLT_PFMHT150_v3")){
if(tr.accept(tr.triggerIndex("HLT_PFMHT150_v3")))hltmet->Fill(1);
else hltmet->Fill(0);
}else{
if(TrIndex_Unknown != tr.triggerIndex("HLT_PFMHT150_v2")){
if(tr.accept(tr.triggerIndex("HLT_PFMHT150_v2"))) hltmet->Fill(1);
else hltmet->Fill(0);
}
else{
printf("BUG with HLT_MET\n");
}
}
// HLT TRIGGER BASED ON 1 JET!
if(TrIndex_Unknown != tr.triggerIndex("HLT_Jet370_v1")){
if(tr.accept(tr.triggerIndex("HLT_Jet370_v1")))hltjet->Fill(1);
else hltjet->Fill(0);
}else{
if(TrIndex_Unknown != tr.triggerIndex("HLT_Jet100U")){
if(IncreasedTreshold(trEv, InputTag("hlt1jet100U","","HLT"), 140, 5.,1, false))hltjet->Fill(1);
else hltjet->Fill(0);
}else{
if(TrIndex_Unknown != tr.triggerIndex("HLT_Jet70U")){
if(IncreasedTreshold(trEv, InputTag("hlt1jet70U","","HLT"), 140, 5.,1, false))hltjet->Fill(1);
else hltjet->Fill(0);
}else{
if(TrIndex_Unknown != tr.triggerIndex("HLT_Jet50U")){
if(IncreasedTreshold(trEv, InputTag("hlt1jet50U","","HLT"), 140,2.5, 1, false))hltjet->Fill(1);
else hltjet->Fill(0);
}else{
printf("BUG with HLT_Jet\n");
}
}
}
}
}
| void HSCPValidator::makeRecoPlots | ( | const edm::Event & | iEvent | ) | [private] |
Definition at line 671 of file HSCPValidator.cc.
References dedxVsp, spr::find(), edm::Event::getByLabel(), i, label_, AlCaHLTBitMon_ParallelJobs::p, particleIds_, particleStatus_, funct::pow(), dt_dqm_sourceclient_common_cff::reco, RecoHSCPPtVsGenPt, and mathSSE::sqrt().
Referenced by analyze().
{
using namespace edm;
using namespace reco;
Handle<HepMCProduct> evt;
iEvent.getByLabel(label_, evt);
Handle<TrackCollection> tkTracks;
iEvent.getByLabel("generalTracks",tkTracks);
const reco::TrackCollection tkTC = *(tkTracks.product());
Handle<ValueMap<DeDxData> > dEdxTrackHandle;
iEvent.getByLabel("dedxHarmonic2", dEdxTrackHandle);
const ValueMap<DeDxData> dEdxTrack = *dEdxTrackHandle.product();
for(size_t i=0; i<tkTracks->size(); i++){
reco::TrackRef trkRef = reco::TrackRef(tkTracks, i);
if(trkRef->pt()<5 || trkRef->normalizedChi2()>10) continue;
double minR= 999;
double hscpgenPt =-1;
HepMC::GenEvent * myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
for(HepMC::GenEvent::particle_iterator p = myGenEvent->particles_begin();
p != myGenEvent->particles_end(); ++p )
{
if((*p)->status() != particleStatus_)
continue;
// Check if the particleId is in our R-hadron list
std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(),particleIds_.end(),(*p)->pdg_id());
if(partIdItr!=particleIds_.end()){
//calculate DeltaR
double distance =pow((*p)->momentum().eta()-trkRef->eta(),2)+pow((*p)->momentum().phi()-trkRef->phi(),2);
distance =sqrt(distance);
if(distance <minR ){
minR = distance;
hscpgenPt= (*p)->momentum().perp();
}
}
}
RecoHSCPPtVsGenPt->Fill(trkRef->pt(),hscpgenPt);
delete myGenEvent;
double dedx = dEdxTrack[trkRef].dEdx();
dedxVsp->Fill( trkRef->p(),dedx);
}
}
| void HSCPValidator::makeSimDigiPlotsECAL | ( | const edm::Event & | iEvent | ) | [private] |
Definition at line 459 of file HSCPValidator.cc.
References EcalMGPASample::adc(), gather_cfg::cout, digiOccupancyMapEB_, digiOccupancyMapEEM_, digiOccupancyMapEEP_, digisEcalNumHistEB_, digisEcalNumHistEE_, EBDigiCollection_, ebSimHitTag_, EEDigiCollection_, eeSimHitTag_, spr::find(), edm::Event::getByLabel(), i, particleIds_, EcalDataFrame::sample(), simHitsEcalDigiMatchEnergyHistEB_, simHitsEcalDigiMatchEnergyHistEE_, simHitsEcalDigiMatchEnergyVsTimeHistEB_, simHitsEcalDigiMatchEnergyVsTimeHistEE_, simHitsEcalDigiMatchIEtaHist_, simHitsEcalDigiMatchIPhiHist_, simHitsEcalDigiMatchTimeHistEB_, simHitsEcalDigiMatchTimeHistEE_, simHitsEcalEnergyHistEB_, simHitsEcalEnergyHistEE_, simHitsEcalEnergyVsTimeHistEB_, simHitsEcalEnergyVsTimeHistEE_, simHitsEcalNumHistEB_, simHitsEcalNumHistEE_, simHitsEcalTimeHistEB_, simHitsEcalTimeHistEE_, simTrackTag_, and EcalDataFrame::size().
Referenced by analyze().
{
using namespace edm;
// EB SimHits
Handle<PCaloHitContainer> ebSimHits;
iEvent.getByLabel(ebSimHitTag_, ebSimHits);
if(!ebSimHits.isValid())
{
std::cout << "Cannot get EBSimHits from event!" << std::endl;
return;
}
// EE SimHits
Handle<PCaloHitContainer> eeSimHits;
iEvent.getByLabel(eeSimHitTag_, eeSimHits);
if(!eeSimHits.isValid())
{
std::cout << "Cannot get EESimHits from event!" << std::endl;
return;
}
// SimTracks
Handle<SimTrackContainer> simTracks;
iEvent.getByLabel(simTrackTag_,simTracks);
if(!simTracks.isValid())
{
std::cout << "Cannot get SimTracks from event!" << std::endl;
return;
}
// EB Digis
Handle<EBDigiCollection> ebDigis;
iEvent.getByLabel(EBDigiCollection_,ebDigis);
if(!ebDigis.isValid())
{
std::cout << "Cannot get EBDigis from event!" << std::endl;
return;
}
// EE Digis
Handle<EEDigiCollection> eeDigis;
iEvent.getByLabel(EEDigiCollection_,eeDigis);
if(!eeDigis.isValid())
{
std::cout << "Cannot get EEDigis from event!" << std::endl;
return;
}
// EB first
// 1) Look at SimTracks, getting only the HSCP tracks
// 2) Match to PCaloHits
// 3) Match to digis
int numMatchedSimHitsEventEB = 0;
int numMatchedDigisEventEB = 0;
const PCaloHitContainer* phitsEB=0;
phitsEB = ebSimHits.product();
for(SimTrackContainer::const_iterator simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack)
{
// Check if the particleId is in our list
std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(),particleIds_.end(),simTrack->type());
if(partIdItr==particleIds_.end())
continue;
PCaloHitContainer mySimHitsEB;
std::vector<EBDataFrame> myDigisEB;
//int particleId = simTrack->type();
int trackId = simTrack->trackId();
PCaloHitContainer::const_iterator simHitItr = phitsEB->begin();
while(simHitItr != phitsEB->end())
{
if(simHitItr->geantTrackId()==trackId)
mySimHitsEB.push_back(*simHitItr);
++simHitItr;
}
if(mySimHitsEB.size()==0)
{
std::cout << "Could not find matching EB PCaloHits for SimTrack id: " << trackId << ". Skipping this SimTrack" << std::endl;
continue;
}
// Loop over matching PCaloHits
for(simHitItr = mySimHitsEB.begin(); simHitItr != mySimHitsEB.end(); ++simHitItr)
{
simHitsEcalEnergyHistEB_->Fill(simHitItr->energy());
simHitsEcalTimeHistEB_->Fill(simHitItr->time());
simHitsEcalEnergyVsTimeHistEB_->Fill(simHitItr->time(),simHitItr->energy());
EBDetId simHitId = EBDetId(simHitItr->id());
std::cout << "SimHit DetId found: " << simHitId << " for PDGid: " << simTrack->type() << std::endl;
//std::cout << "SimHit hashedIndex: " << simHitId.hashedIndex() << std::endl;
std::cout << "SimHit energy: " << simHitItr->energy() << " time: " << simHitItr->time() << std::endl;
++numMatchedSimHitsEventEB;
EBDigiCollection::const_iterator digiItr = ebDigis->begin();
while(digiItr != ebDigis->end() && (digiItr->id() != simHitId))
++digiItr;
if(digiItr==ebDigis->end())
{
// Commented out for debugging ease, Aug 3 2009
std::cout << "Could not find simHit detId: " << simHitId << "in EBDigiCollection!" << std::endl;
continue;
}
std::vector<EBDataFrame>::const_iterator myDigiItr = myDigisEB.begin();
while(myDigiItr != myDigisEB.end() && (digiItr->id() != myDigiItr->id()))
++myDigiItr;
if(myDigiItr!=myDigisEB.end())
continue; // if this digi is already in the list, skip it
++numMatchedDigisEventEB;
EBDataFrame df = *digiItr;
myDigisEB.push_back(df);
std::cout << "SAMPLE ADCs: " << "\t";
for(int i=0; i<10;++i)
std::cout << i << "\t";
std::cout << std::endl << "\t\t";
for(int i=0; i < df.size(); ++i)
{
std::cout << df.sample(i).adc() << "\t";
}
std::cout << std::endl << std::endl;
simHitsEcalDigiMatchEnergyHistEB_->Fill(simHitItr->energy());
simHitsEcalDigiMatchTimeHistEB_->Fill(simHitItr->time());
simHitsEcalDigiMatchEnergyVsTimeHistEB_->Fill(simHitItr->time(),simHitItr->energy());
simHitsEcalDigiMatchIEtaHist_->Fill(((EBDetId)digiItr->id()).ieta());
simHitsEcalDigiMatchIPhiHist_->Fill(((EBDetId)digiItr->id()).iphi());
digiOccupancyMapEB_->Fill(((EBDetId)digiItr->id()).iphi(),((EBDetId)digiItr->id()).ieta());
}
}
simHitsEcalNumHistEB_->Fill(numMatchedSimHitsEventEB);
digisEcalNumHistEB_->Fill(numMatchedDigisEventEB);
// EE next
int numMatchedSimHitsEventEE = 0;
int numMatchedDigisEventEE = 0;
const PCaloHitContainer* phitsEE=0;
phitsEE = eeSimHits.product();
for(SimTrackContainer::const_iterator simTrack = simTracks->begin(); simTrack != simTracks->end(); ++simTrack)
{
// Check if the particleId is in our list
std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(),particleIds_.end(),simTrack->type());
if(partIdItr==particleIds_.end())
continue;
PCaloHitContainer mySimHitsEE;
std::vector<EEDataFrame> myDigisEE;
//int particleId = simTrack->type();
int trackId = simTrack->trackId();
PCaloHitContainer::const_iterator simHitItr = phitsEE->begin();
while(simHitItr != phitsEE->end())
{
if(simHitItr->geantTrackId()==trackId)
mySimHitsEE.push_back(*simHitItr);
++simHitItr;
}
if(mySimHitsEE.size()==0)
{
std::cout << "Could not find matching EE PCaloHits for SimTrack id: " << trackId << ". Skipping this SimTrack" << std::endl;
continue;
}
// Loop over matching PCaloHits
for(simHitItr = mySimHitsEE.begin(); simHitItr != mySimHitsEE.end(); ++simHitItr)
{
simHitsEcalEnergyHistEE_->Fill(simHitItr->energy());
simHitsEcalTimeHistEE_->Fill(simHitItr->time());
simHitsEcalEnergyVsTimeHistEE_->Fill(simHitItr->time(),simHitItr->energy());
EEDetId simHitId = EEDetId(simHitItr->id());
std::cout << "SimHit DetId found: " << simHitId << " for PDGid: " << simTrack->type() << std::endl;
//std::cout << "SimHit hashedIndex: " << simHitId.hashedIndex() << std::endl;
std::cout << "SimHit energy: " << simHitItr->energy() << " time: " << simHitItr->time() << std::endl;
++numMatchedSimHitsEventEE;
EEDigiCollection::const_iterator digiItr = eeDigis->begin();
while(digiItr != eeDigis->end() && (digiItr->id() != simHitId))
++digiItr;
if(digiItr==eeDigis->end())
{
// Commented out for debugging ease, Aug 3 2009
std::cout << "Could not find simHit detId: " << simHitId << "in EEDigiCollection!" << std::endl;
continue;
}
std::vector<EEDataFrame>::const_iterator myDigiItr = myDigisEE.begin();
while(myDigiItr != myDigisEE.end() && (digiItr->id() != myDigiItr->id()))
++myDigiItr;
if(myDigiItr!=myDigisEE.end())
continue; // if this digi is already in the list, skip it
++numMatchedDigisEventEE;
EEDataFrame df = *digiItr;
myDigisEE.push_back(df);
std::cout << "SAMPLE ADCs: " << "\t";
for(int i=0; i<10;++i)
std::cout << i << "\t";
std::cout << std::endl << "\t\t";
for(int i=0; i < df.size(); ++i)
{
std::cout << df.sample(i).adc() << "\t";
}
std::cout << std::endl << std::endl;
simHitsEcalDigiMatchEnergyHistEE_->Fill(simHitItr->energy());
simHitsEcalDigiMatchTimeHistEE_->Fill(simHitItr->time());
simHitsEcalDigiMatchEnergyVsTimeHistEE_->Fill(simHitItr->time(),simHitItr->energy());
if(((EEDetId)digiItr->id()).zside() > 0)
digiOccupancyMapEEP_->Fill(((EEDetId)digiItr->id()).ix(),((EEDetId)digiItr->id()).iy());
else if(((EEDetId)digiItr->id()).zside() < 0)
digiOccupancyMapEEM_->Fill(((EEDetId)digiItr->id()).ix(),((EEDetId)digiItr->id()).iy());
}
}
simHitsEcalNumHistEE_->Fill(numMatchedSimHitsEventEE);
digisEcalNumHistEE_->Fill(numMatchedDigisEventEE);
}
| void HSCPValidator::makeSimDigiPlotsRPC | ( | const edm::Event & | iEvent | ) | [private] |
Definition at line 727 of file HSCPValidator.cc.
References cluSizeDistribution_, gather_cfg::cout, efficiencyRPCRecHitSimDigis_, PV3DBase< T, PVType, FrameType >::eta(), HcalObjRepresent::Fill(), spr::find(), edm::Event::getByLabel(), edm::Event::getManyByType(), i, RPCDetId::layer(), PV3DBase< T, PVType, FrameType >::mag(), DetId::Muon, RPCGeomServ::name(), particleIds_, PV3DBase< T, PVType, FrameType >::phi(), RPCDetId::region(), residualsRPCRecHitSimDigis_, dedefs::RPC, rpcBXBarrel_, rpcBXEndCap_, rpcGeo, python::rpcRecHits_cfi::rpcRecHits, rpcTimeOfFlightBarrel_, rpcTimeOfFlightEndCap_, RPCDetId::station(), GeomDet::surface(), Surface::toGlobal(), and PV3DBase< T, PVType, FrameType >::x().
Referenced by analyze().
{
using namespace edm;
//std::cout << " Getting the SimHits " <<std::endl;
std::vector<Handle<edm::PSimHitContainer> > theSimHitContainers;
iEvent.getManyByType(theSimHitContainers);
//std::cout << " The Number of sim Hits is " << theSimHitContainers.size() <<std::endl;
Handle<RPCRecHitCollection> rpcRecHits;
iEvent.getByLabel("rpcRecHits","",rpcRecHits);
//SimTrack Stuff
std::vector<PSimHit> theSimHits;
for (int i = 0; i < int(theSimHitContainers.size()); i++){
theSimHits.insert(theSimHits.end(),theSimHitContainers.at(i)->begin(),theSimHitContainers.at(i)->end());
}
for (std::vector<PSimHit>::const_iterator iHit = theSimHits.begin(); iHit != theSimHits.end(); iHit++){
std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(),particleIds_.end(),(*iHit).particleType());
if(partIdItr==particleIds_.end())
continue;
DetId theDetUnitId((*iHit).detUnitId());
DetId simdetid= DetId((*iHit).detUnitId());
if(simdetid.det()==DetId::Muon && simdetid.subdetId()== MuonSubdetId::RPC){//Only RPCs
RPCDetId rollId(theDetUnitId);
RPCGeomServ rpcsrv(rollId);
//std::cout << " Reading the Roll"<<std::endl;
const RPCRoll* rollasociated = rpcGeo->roll(rollId);
//std::cout << " Getting the Surface"<<std::endl;
const BoundPlane & RPCSurface = rollasociated->surface();
GlobalPoint SimHitInGlobal = RPCSurface.toGlobal((*iHit).localPosition());
std::cout<<"\t\t We have an RPC Sim Hit! in t="<<(*iHit).timeOfFlight()<<"ns "<<rpcsrv.name()<<" Global postition="<<SimHitInGlobal<<std::endl;
int layer = 0;
if(rollId.station()==1&&rollId.layer()==1) layer = 1;
else if(rollId.station()==1&&rollId.layer()==2) layer = 2;
else if(rollId.station()==2&&rollId.layer()==1) layer = 3;
else if(rollId.station()==2&&rollId.layer()==2) layer = 4;
else if(rollId.station()==3) layer = 5;
else if(rollId.station()==4) layer = 6;
if(rollId.region()==0){
rpcTimeOfFlightBarrel_[layer-1]->Fill((*iHit).timeOfFlight());
}else{
rpcTimeOfFlightEndCap_[rollId.station()-1]->Fill((*iHit).timeOfFlight());
}
std::cout<<"\t\t r="<<SimHitInGlobal.mag()<<" phi="<<SimHitInGlobal.phi()<<" eta="<<SimHitInGlobal.eta()<<std::endl;
int cluSize = 0;
int bx = 100;
float minres = 3000.;
std::cout<<"\t \t \t \t Getting RecHits in Roll Asociated"<<std::endl;
typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
rangeRecHits recHitCollection = rpcRecHits->get(rollId);
RPCRecHitCollection::const_iterator recHit;
efficiencyRPCRecHitSimDigis_->Fill(0);
for(recHit = recHitCollection.first; recHit != recHitCollection.second ; recHit++){
LocalPoint recHitPos=recHit->localPosition();
float res=(*iHit).localPosition().x()- recHitPos.x();
if(fabs(res)<fabs(minres)){
minres=res;
cluSize = recHit->clusterSize();
bx=recHit->BunchX();
std::cout<<"\t New Min Res "<<res<<"cm."<<std::endl;
}
}
if(minres<3000.){
residualsRPCRecHitSimDigis_->Fill(minres);
efficiencyRPCRecHitSimDigis_->Fill(1);
cluSizeDistribution_->Fill(cluSize);
if(rollId.region()==0) rpcBXBarrel_[layer-1]->Fill(bx);
else rpcBXEndCap_[rollId.station()-1]->Fill(bx);
}
}
}
}
| void HSCPValidator::makeSimTrackPlots | ( | const edm::Event & | iEvent | ) | [private] |
Definition at line 345 of file HSCPValidator.cc.
References spr::find(), edm::Event::getByLabel(), particleIds_, simTrackParticleBetaHist_, simTrackParticleEtaHist_, simTrackParticlePhiHist_, simTrackParticlePHist_, and simTrackParticlePtHist_.
Referenced by analyze().
{ using namespace edm;
//get sim track infos
Handle<edm::SimTrackContainer> simTracksHandle;
iEvent.getByLabel("g4SimHits",simTracksHandle);
const SimTrackContainer simTracks = *(simTracksHandle.product());
SimTrackContainer::const_iterator simTrack;
for (simTrack = simTracks.begin(); simTrack != simTracks.end(); ++simTrack){
// Check if the particleId is in our list
std::vector<int>::const_iterator partIdItr = find(particleIds_.begin(),particleIds_.end(),simTrack->type());
if(partIdItr==particleIds_.end()) continue;
simTrackParticleEtaHist_->Fill((*simTrack).momentum().eta());
simTrackParticlePhiHist_->Fill((*simTrack).momentum().phi());
simTrackParticlePHist_->Fill((*simTrack).momentum().P());
simTrackParticlePtHist_->Fill((*simTrack).momentum().pt());
simTrackParticleBetaHist_->Fill((*simTrack).momentum().P()/(*simTrack).momentum().e());
// std::cout<<"Particle:"<<simTrack->type()<<" Charge:"<<simTrack->charge()<<std::endl;
}
}
TH1F* HSCPValidator::cluSizeDistribution_ [private] |
Definition at line 125 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().
TH2F* HSCPValidator::dedxVsp [private] |
Definition at line 136 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeRecoPlots().
TH2F* HSCPValidator::digiOccupancyMapEB_ [private] |
Definition at line 119 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH2F* HSCPValidator::digiOccupancyMapEEM_ [private] |
Definition at line 121 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH2F* HSCPValidator::digiOccupancyMapEEP_ [private] |
Definition at line 120 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::digisEcalNumHistEB_ [private] |
Definition at line 117 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::digisEcalNumHistEE_ [private] |
Definition at line 118 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
bool HSCPValidator::doGenPlots_ [private] |
Definition at line 60 of file HSCPValidator.h.
Referenced by analyze().
bool HSCPValidator::doHLTPlots_ [private] |
Definition at line 61 of file HSCPValidator.h.
Referenced by analyze().
bool HSCPValidator::doRecoPlots_ [private] |
Definition at line 64 of file HSCPValidator.h.
Referenced by analyze().
bool HSCPValidator::doSimDigiPlots_ [private] |
Definition at line 63 of file HSCPValidator.h.
Referenced by analyze().
bool HSCPValidator::doSimTrackPlots_ [private] |
Definition at line 62 of file HSCPValidator.h.
Referenced by analyze().
Definition at line 96 of file HSCPValidator.h.
Referenced by makeSimDigiPlotsECAL().
edm::InputTag HSCPValidator::ebSimHitTag_ [private] |
Definition at line 93 of file HSCPValidator.h.
Referenced by makeSimDigiPlotsECAL().
Definition at line 97 of file HSCPValidator.h.
Referenced by makeSimDigiPlotsECAL().
edm::InputTag HSCPValidator::eeSimHitTag_ [private] |
Definition at line 94 of file HSCPValidator.h.
Referenced by makeSimDigiPlotsECAL().
TH1F* HSCPValidator::efficiencyRPCRecHitSimDigis_ [private] |
Definition at line 124 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().
TH1F* HSCPValidator::h_genhscp_met [private] |
Definition at line 79 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::h_genhscp_met_nohscp [private] |
Definition at line 80 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::h_genhscp_scaloret [private] |
Definition at line 81 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::h_genhscp_scaloret_nohscp [private] |
Definition at line 82 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::hltjet [private] |
Definition at line 132 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeHLTPlots().
TH1F* HSCPValidator::hltmet [private] |
Definition at line 131 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeHLTPlots().
TH1F* HSCPValidator::hltmu [private] |
Definition at line 133 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeHLTPlots().
edm::InputTag HSCPValidator::label_ [private] |
Definition at line 67 of file HSCPValidator.h.
Referenced by makeGenPlots(), and makeRecoPlots().
TH1F* HSCPValidator::particleBetaHist_ [private] |
Definition at line 77 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::particleBetaInverseHist_ [private] |
Definition at line 78 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::particleEtaHist_ [private] |
Definition at line 71 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
std::vector<int> HSCPValidator::particleIds_ [private] |
Definition at line 68 of file HSCPValidator.h.
Referenced by makeGenPlots(), makeRecoPlots(), makeSimDigiPlotsECAL(), makeSimDigiPlotsRPC(), and makeSimTrackPlots().
std::map<int,int> HSCPValidator::particleIdsFoundMap_ [private] |
Definition at line 70 of file HSCPValidator.h.
Referenced by endJob(), and makeGenPlots().
TH1F* HSCPValidator::particleMassHist_ [private] |
Definition at line 75 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::particlePhiHist_ [private] |
Definition at line 72 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::particlePHist_ [private] |
Definition at line 73 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH1F* HSCPValidator::particlePtHist_ [private] |
Definition at line 74 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
int HSCPValidator::particleStatus_ [private] |
Definition at line 69 of file HSCPValidator.h.
Referenced by makeGenPlots(), and makeRecoPlots().
TH1F* HSCPValidator::particleStatusHist_ [private] |
Definition at line 76 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeGenPlots().
TH2F* HSCPValidator::RecoHSCPPtVsGenPt [private] |
Definition at line 135 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeRecoPlots().
TH1F* HSCPValidator::residualsRPCRecHitSimDigis_ [private] |
Definition at line 123 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().
TH1F* HSCPValidator::rpcBXBarrel_[6] [private] |
Definition at line 127 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().
TH1F* HSCPValidator::rpcBXEndCap_[3] [private] |
Definition at line 129 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().
edm::ESHandle<RPCGeometry> HSCPValidator::rpcGeo [private] |
Definition at line 99 of file HSCPValidator.h.
Referenced by analyze(), and makeSimDigiPlotsRPC().
edm::InputTag HSCPValidator::RPCRecHitTag_ [private] |
Definition at line 98 of file HSCPValidator.h.
TH1F* HSCPValidator::rpcTimeOfFlightBarrel_[6] [private] |
Definition at line 126 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().
TH1F* HSCPValidator::rpcTimeOfFlightEndCap_[3] [private] |
Definition at line 128 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsRPC().
TH1F* HSCPValidator::simHitsEcalDigiMatchEnergyHistEB_ [private] |
Definition at line 109 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalDigiMatchEnergyHistEE_ [private] |
Definition at line 112 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH2F* HSCPValidator::simHitsEcalDigiMatchEnergyVsTimeHistEB_ [private] |
Definition at line 111 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH2F* HSCPValidator::simHitsEcalDigiMatchEnergyVsTimeHistEE_ [private] |
Definition at line 114 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalDigiMatchIEtaHist_ [private] |
Definition at line 115 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalDigiMatchIPhiHist_ [private] |
Definition at line 116 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalDigiMatchTimeHistEB_ [private] |
Definition at line 110 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalDigiMatchTimeHistEE_ [private] |
Definition at line 113 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalEnergyHistEB_ [private] |
Definition at line 101 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalEnergyHistEE_ [private] |
Definition at line 105 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH2F* HSCPValidator::simHitsEcalEnergyVsTimeHistEB_ [private] |
Definition at line 104 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH2F* HSCPValidator::simHitsEcalEnergyVsTimeHistEE_ [private] |
Definition at line 108 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalNumHistEB_ [private] |
Definition at line 103 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalNumHistEE_ [private] |
Definition at line 107 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalTimeHistEB_ [private] |
Definition at line 102 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simHitsEcalTimeHistEE_ [private] |
Definition at line 106 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimDigiPlotsECAL().
TH1F* HSCPValidator::simTrackParticleBetaHist_ [private] |
Definition at line 90 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimTrackPlots().
TH1F* HSCPValidator::simTrackParticleEtaHist_ [private] |
Definition at line 86 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimTrackPlots().
TH1F* HSCPValidator::simTrackParticlePhiHist_ [private] |
Definition at line 87 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimTrackPlots().
TH1F* HSCPValidator::simTrackParticlePHist_ [private] |
Definition at line 88 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimTrackPlots().
TH1F* HSCPValidator::simTrackParticlePtHist_ [private] |
Definition at line 89 of file HSCPValidator.h.
Referenced by HSCPValidator(), and makeSimTrackPlots().
edm::InputTag HSCPValidator::simTrackTag_ [private] |
Definition at line 95 of file HSCPValidator.h.
Referenced by makeSimDigiPlotsECAL().
1.7.3