Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

memory reduction most visible for HI processing, done mostly by dropping muon isodeposits in HI #12764

Merged
merged 7 commits into from Dec 14, 2015
Merged

memory reduction most visible for HI processing, done mostly by dropping muon isodeposits in HI #12764

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
2f75d96
more compact and possibly less memory-consuming valueMap filling
slava77devel Dec 9, 2015
aa62913
unique_ptr
slava77devel Dec 9, 2015
040b9e3
don't fill the vector of isodeposits if there is no plan to write it
slava77devel Dec 11, 2015
fb70567
don't save detector-based isodeposits from muons
slava77devel Dec 11, 2015
bee9370
better memory management
slava77devel Dec 12, 2015
1f193df
unused ID variables; now clearly unused after refactoring
slava77devel Dec 12, 2015
e66f643
unused ID variables; now clearly unused after refactoring
slava77devel Dec 12, 2015
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
12 changes: 6 additions & 6 deletions DQM/SiPixelMonitorTrack/src/SiPixelTrackResidualSource.cc
View file
Open in desktop
Expand Up @@ -875,26 +875,26 @@ void SiPixelTrackResidualSource::analyze(const edm::Event& iEvent, const edm::Ev
edm::Handle<std::vector<Trajectory> > trajCollectionHandle;
//iEvent.getByLabel(tracksrc_,trajCollectionHandle);
iEvent.getByToken ( tracksrcToken_, trajCollectionHandle );
const std::vector<Trajectory> trajColl = *(trajCollectionHandle.product());
const std::vector<Trajectory>& trajColl = *(trajCollectionHandle.product());

//get tracks
edm::Handle<std::vector<reco::Track> > trackCollectionHandle;
//iEvent.getByLabel(tracksrc_,trackCollectionHandle);
iEvent.getByToken( trackToken_, trackCollectionHandle );

const std::vector<reco::Track> trackColl = *(trackCollectionHandle.product());
const std::vector<reco::Track>& trackColl = *(trackCollectionHandle.product());

//get the map
edm::Handle<TrajTrackAssociationCollection> match;
//iEvent.getByLabel(tracksrc_,match);
iEvent.getByToken( trackAssociationToken_, match);
const TrajTrackAssociationCollection ttac = *(match.product());
const TrajTrackAssociationCollection& ttac = *(match.product());

// get clusters
edm::Handle< edmNew::DetSetVector<SiPixelCluster> > clusterColl;
//iEvent.getByLabel( clustersrc_, clusterColl );
iEvent.getByToken( clustersrcToken_, clusterColl );
const edmNew::DetSetVector<SiPixelCluster> clustColl = *(clusterColl.product());
const edmNew::DetSetVector<SiPixelCluster>& clustColl = *(clusterColl.product());

if(debug_){
std::cout << "Trajectories\t : " << trajColl.size() << std::endl;
Expand Down Expand Up @@ -923,7 +923,7 @@ void SiPixelTrackResidualSource::analyze(const edm::Event& iEvent, const edm::Ev

if(abs(d0)<15 && abs(dz)<50) crossesPixVol = true;

std::vector<TrajectoryMeasurement> tmeasColl =traj_iterator->measurements();
const std::vector<TrajectoryMeasurement>& tmeasColl =traj_iterator->measurements();
std::vector<TrajectoryMeasurement>::const_iterator tmeasIt;
//loop on measurements to find out whether there are bpix and/or fpix hits
for(tmeasIt = tmeasColl.begin();tmeasIt!=tmeasColl.end();tmeasIt++){
Expand All @@ -947,7 +947,7 @@ void SiPixelTrackResidualSource::analyze(const edm::Event& iEvent, const edm::Ev

if(crossesPixVol) meNofTracksInPixVol_->Fill(0,1);

std::vector<TrajectoryMeasurement> tmeasColl = traj_iterator->measurements();
const std::vector<TrajectoryMeasurement>& tmeasColl = traj_iterator->measurements();
for(std::vector<TrajectoryMeasurement>::const_iterator tmeasIt = tmeasColl.begin(); tmeasIt!=tmeasColl.end(); tmeasIt++){
if(! tmeasIt->updatedState().isValid()) continue;

Expand Down
30 changes: 15 additions & 15 deletions DQMOffline/EGamma/plugins/PhotonAnalyzer.cc
View file
Open in desktop
Expand Up @@ -615,45 +615,45 @@ void PhotonAnalyzer::analyze( const edm::Event& e, const edm::EventSetup& esup )
// Get the trigger results
bool validTriggerEvent=true;
edm::Handle<trigger::TriggerEvent> triggerEventHandle;
trigger::TriggerEvent triggerEvent;
const trigger::TriggerEvent dummyTE;
e.getByToken(triggerEvent_token_,triggerEventHandle);
if(!triggerEventHandle.isValid()) {
edm::LogInfo(fName_) << "Error! Can't get the product: triggerEvent_" << endl;
validTriggerEvent=false;
}
if(validTriggerEvent) triggerEvent = *(triggerEventHandle.product());
const trigger::TriggerEvent& triggerEvent(validTriggerEvent? *(triggerEventHandle.product()) : dummyTE);

// Get the reconstructed photons
// bool validPhotons=true;
Handle<reco::PhotonCollection> photonHandle;
reco::PhotonCollection photonCollection;
e.getByToken(photon_token_ , photonHandle);
if ( !photonHandle.isValid()) {
edm::LogInfo(fName_) << "Error! Can't get the product: photon_token_" << endl;
// validPhotons=false;
}
// if(validPhotons) photonCollection = *(photonHandle.product());
const reco::PhotonCollection& photonCollection(*(photonHandle.product()));


// Get the PhotonId objects
bool validloosePhotonID=true;
// bool validloosePhotonID=true;
Handle<edm::ValueMap<bool> > loosePhotonFlag;
edm::ValueMap<bool> loosePhotonID;
e.getByToken(PhotonIDLoose_token_, loosePhotonFlag);
if ( !loosePhotonFlag.isValid()) {
edm::LogInfo(fName_) << "Error! Can't get the product: PhotonIDLoose_token_" << endl;
validloosePhotonID=false;
// validloosePhotonID=false;
}
if (validloosePhotonID) loosePhotonID = *(loosePhotonFlag.product());
// edm::ValueMap<bool> dummyLPID;
// const edm::ValueMap<bool>& loosePhotonID(validloosePhotonID? *(loosePhotonFlag.product()) : dummyLPID);

bool validtightPhotonID=true;
// bool validtightPhotonID=true;
Handle<edm::ValueMap<bool> > tightPhotonFlag;
edm::ValueMap<bool> tightPhotonID;
e.getByToken(PhotonIDTight_token_, tightPhotonFlag);
if ( !tightPhotonFlag.isValid()) {
edm::LogInfo(fName_) << "Error! Can't get the product: PhotonIDTight_token_" << endl;
validtightPhotonID=false;
// validtightPhotonID=false;
}
if (validtightPhotonID) tightPhotonID = *(tightPhotonFlag.product());
// edm::ValueMap<bool> dummyTPI;
// const edm::ValueMap<bool>& tightPhotonID(validtightPhotonID ? *(tightPhotonFlag.product()) : dummyTPI);


edm::Handle<reco::VertexCollection> vtxH;
Expand Down Expand Up @@ -714,7 +714,7 @@ void PhotonAnalyzer::analyze( const edm::Event& e, const edm::EventSetup& esup )

/////////////////////////BEGIN LOOP OVER THE COLLECTION OF PHOTONS IN THE EVENT/////////////////////////
for(unsigned int iPho=0; iPho < photonHandle->size(); iPho++) {
reco::PhotonRef aPho(reco::PhotonRef(photonHandle, iPho));
const reco::Photon* aPho = &photonCollection[iPho];
// for( reco::PhotonCollection::const_iterator iPho = photonCollection.begin(); iPho != photonCollection.end(); iPho++) {

//for HLT efficiency plots
Expand Down Expand Up @@ -1083,7 +1083,7 @@ void PhotonAnalyzer::analyze( const edm::Event& e, const edm::EventSetup& esup )

if (isIsolated && aPho->et()>=invMassEtCut_){
for(unsigned int iPho2=iPho+1; iPho2 < photonHandle->size(); iPho2++) {
reco::PhotonRef aPho2(reco::PhotonRef(photonHandle, iPho2));
const reco::Photon* aPho2 = &photonCollection[iPho2];

// for (reco::PhotonCollection::const_iterator iPho2=iPho+1; iPho2!=photonCollection.end(); iPho2++){

Expand Down Expand Up @@ -1178,7 +1178,7 @@ void PhotonAnalyzer::fill3DHistoVector(vector<vector<vector<MonitorElement*> > >
histoVector[cut][type][part]->Fill(x,y);
}

bool PhotonAnalyzer::photonSelection(const reco::PhotonRef & pho)
bool PhotonAnalyzer::photonSelection(const reco::Photon* pho)
{
bool result=true;
if ( pho->pt() < minPhoEtCut_ ) result=false;
Expand Down
2 changes: 1 addition & 1 deletion DQMOffline/EGamma/plugins/PhotonAnalyzer.h
View file
Open in desktop
Expand Up @@ -136,7 +136,7 @@ class PhotonAnalyzer : public DQMEDAnalyzer
void fill3DHistoVector(std::vector<std::vector<std::vector<MonitorElement*> > >& histoVector,double x, int cut, int type, int part);
void fill3DHistoVector(std::vector<std::vector<std::vector<MonitorElement*> > >& histoVector,double x, double y, int cut, int type, int part);

bool photonSelection(const reco::PhotonRef & p);
bool photonSelection(const reco::Photon* p);
float phiNormalization(float& a);

//////////
Expand Down
54 changes: 27 additions & 27 deletions DQMOffline/JetMET/src/JetAnalyzer_HeavyIons_matching.cc
View file
Open in desktop
Expand Up @@ -183,9 +183,9 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm
// Get the Jet collection
//----------------------------------------------------------------------------

std::vector<Jet> recoJet1;
std::vector<const Jet*> recoJet1;
recoJet1.clear();
std::vector<Jet> recoJet2;
std::vector<const Jet*> recoJet2;
recoJet2.clear();

edm::Handle<CaloJetCollection> caloJet1;
Expand All @@ -196,36 +196,36 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm

if(std::string("VsCalo") == JetType1) {
mEvent.getByToken(caloJet1Token_, caloJet1);
for (unsigned ijet=0; ijet<caloJet1->size(); ++ijet) recoJet1.push_back((*caloJet1)[ijet]);
for (unsigned ijet=0; ijet<caloJet1->size(); ++ijet) recoJet1.push_back(&(*caloJet1)[ijet]);
}
if(std::string("PuCalo") == JetType1) {
mEvent.getByToken(caloJet2Token_, caloJet1);
for (unsigned ijet=0; ijet<caloJet1->size(); ++ijet) recoJet1.push_back((*caloJet1)[ijet]);
for (unsigned ijet=0; ijet<caloJet1->size(); ++ijet) recoJet1.push_back(&(*caloJet1)[ijet]);
}
if(std::string("VsPF") == JetType1) {
mEvent.getByToken(pfJetsToken_, pfJets);
for (unsigned ijet=0; ijet<pfJets->size(); ++ijet) recoJet1.push_back((*pfJets)[ijet]);
for (unsigned ijet=0; ijet<pfJets->size(); ++ijet) recoJet1.push_back(&(*pfJets)[ijet]);
}
if(std::string("PuPF") == JetType1) {
mEvent.getByToken(basicJetsToken_, basicJets);
for (unsigned ijet=0; ijet<basicJets->size(); ++ijet) recoJet1.push_back((*basicJets)[ijet]);
for (unsigned ijet=0; ijet<basicJets->size(); ++ijet) recoJet1.push_back(&(*basicJets)[ijet]);
}

if(std::string("VsCalo") == JetType2) {
mEvent.getByToken(caloJet1Token_, caloJet2);
for (unsigned ijet=0; ijet<caloJet2->size(); ++ijet) recoJet2.push_back((*caloJet2)[ijet]);
for (unsigned ijet=0; ijet<caloJet2->size(); ++ijet) recoJet2.push_back(&(*caloJet2)[ijet]);
}
if(std::string("PuCalo") == JetType2) {
mEvent.getByToken(caloJet2Token_, caloJet2);
for (unsigned ijet=0; ijet<caloJet2->size(); ++ijet) recoJet2.push_back((*caloJet2)[ijet]);
for (unsigned ijet=0; ijet<caloJet2->size(); ++ijet) recoJet2.push_back(&(*caloJet2)[ijet]);
}
if(std::string("VsPF") == JetType2) {
mEvent.getByToken(pfJetsToken_, pfJets);
for (unsigned ijet=0; ijet<pfJets->size(); ++ijet) recoJet2.push_back((*pfJets)[ijet]);
for (unsigned ijet=0; ijet<pfJets->size(); ++ijet) recoJet2.push_back(&(*pfJets)[ijet]);
}
if(std::string("PuPF") == JetType2) {
mEvent.getByToken(basicJetsToken_, basicJets);
for (unsigned ijet=0; ijet<basicJets->size(); ++ijet) recoJet2.push_back((*basicJets)[ijet]);
for (unsigned ijet=0; ijet<basicJets->size(); ++ijet) recoJet2.push_back(&(*basicJets)[ijet]);
}

// start to perform the matching - between recoJet1 and recoJet2.
Expand All @@ -244,13 +244,13 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm
for(unsigned ijet1 = 0; ijet1 < recoJet1.size(); ++ijet1){


if(recoJet1[ijet1].pt() < mRecoJetPtThreshold) continue;
if(fabs(recoJet1[ijet1].eta()) < mRecoJetEtaCut) continue;
if(recoJet1[ijet1]->pt() < mRecoJetPtThreshold) continue;
if(fabs(recoJet1[ijet1]->eta()) < mRecoJetEtaCut) continue;

MyJet JET1;
JET1.eta = recoJet1[ijet1].eta();
JET1.phi = recoJet1[ijet1].phi();
JET1.pt = recoJet1[ijet1].pt();
JET1.eta = recoJet1[ijet1]->eta();
JET1.phi = recoJet1[ijet1]->phi();
JET1.pt = recoJet1[ijet1]->pt();
JET1.id = ijet1;

vJet1.push_back(JET1);
Expand All @@ -260,13 +260,13 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm

for(unsigned ijet2 = 0; ijet2 < recoJet2.size(); ++ijet2){

if(recoJet2[ijet2].pt() < mRecoJetPtThreshold) continue;
if(fabs(recoJet2[ijet2].eta()) < mRecoJetEtaCut) continue;
if(recoJet2[ijet2]->pt() < mRecoJetPtThreshold) continue;
if(fabs(recoJet2[ijet2]->eta()) < mRecoJetEtaCut) continue;

MyJet JET2;
JET2.eta = recoJet2[ijet2].eta();
JET2.phi = recoJet2[ijet2].phi();
JET2.pt = recoJet2[ijet2].pt();
JET2.eta = recoJet2[ijet2]->eta();
JET2.phi = recoJet2[ijet2]->phi();
JET2.pt = recoJet2[ijet2]->pt();
JET2.id = ijet2;

vJet2.push_back(JET2);
Expand All @@ -290,7 +290,7 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm

int pj = (*iJet).id;

mpT_Jet1_unmatched->Fill(recoJet1[pj].pt());
mpT_Jet1_unmatched->Fill(recoJet1[pj]->pt());

}

Expand All @@ -300,7 +300,7 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm

int cj = (*iJet).id;

mpT_Jet2_unmatched->Fill(recoJet2[cj].pt());
mpT_Jet2_unmatched->Fill(recoJet2[cj]->pt());
}

}else if (bothJet1Jet2){
Expand All @@ -326,10 +326,10 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm

if( delr < mRecoDelRMatch && Jet1_ID[Aj.id] == 0){

mpT_ratio_Jet1Jet2->Fill((Float_t) recoJet2[Bj.id].pt()/recoJet1[Aj.id].pt());
mpT_ratio_Jet1Jet2->Fill((Float_t) recoJet2[Bj.id]->pt()/recoJet1[Aj.id]->pt());

mpT_Jet1_matched->Fill(recoJet1[Aj.id].pt());
mpT_Jet2_matched->Fill(recoJet2[Bj.id].pt());
mpT_Jet1_matched->Fill(recoJet1[Aj.id]->pt());
mpT_Jet2_matched->Fill(recoJet2[Bj.id]->pt());

Jet1_ID[Aj.id] = 1;
Jet2_ID[Bj.id] = 1;
Expand All @@ -351,7 +351,7 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm

if(Jet1_ID[Aj.id] == 0) {

mpT_Jet1_unmatched->Fill(recoJet1[Aj.id].pt());
mpT_Jet1_unmatched->Fill(recoJet1[Aj.id]->pt());
unmatchedJet1++;
Jet1_ID[Aj.id] = 1;

Expand Down Expand Up @@ -385,7 +385,7 @@ void JetAnalyzer_HeavyIons_matching::analyze(const edm::Event& mEvent, const edm

if(Jet2_ID[Bj.id] == 0) {

mpT_Jet2_unmatched->Fill(recoJet2[Bj.id].pt());
mpT_Jet2_unmatched->Fill(recoJet2[Bj.id]->pt());
unmatchedJet2++;
Jet2_ID[Bj.id] = 2;
if(std::string("VsCalo") == JetType2 || std::string("PuCalo") == JetType2){
Expand Down
44 changes: 22 additions & 22 deletions DQMOffline/Muon/src/MuonRecoOneHLT.cc
View file
Open in desktop
Expand Up @@ -150,7 +150,7 @@ void MuonRecoOneHLT::analyze(const edm::Event& iEvent, const edm::EventSetup& iS

edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
iEvent.getByToken(theBeamSpotLabel_,recoBeamSpotHandle);
reco::BeamSpot bs = *recoBeamSpotHandle;
const reco::BeamSpot& bs = *recoBeamSpotHandle;

posVtx = bs.position();
errVtx(0,0) = bs.BeamWidthX();
Expand All @@ -173,12 +173,12 @@ void MuonRecoOneHLT::analyze(const edm::Event& iEvent, const edm::EventSetup& iS
if(!muons.isValid()) return;

// Pick the leading lepton.
std::map<float,reco::Muon> muonMap;
std::map<float,const reco::Muon*> muonMap;
for (reco::MuonCollection::const_iterator recoMu = muons->begin(); recoMu!=muons->end(); ++recoMu){
muonMap[recoMu->pt()] = *recoMu;
muonMap[recoMu->pt()] = &*recoMu;
}
std::vector<reco::Muon> LeadingMuon;
for( std::map<float,reco::Muon>::reverse_iterator rit=muonMap.rbegin(); rit!=muonMap.rend(); ++rit){
std::vector<const reco::Muon*> LeadingMuon;
for( std::map<float,const reco::Muon*>::reverse_iterator rit=muonMap.rbegin(); rit!=muonMap.rend(); ++rit){
LeadingMuon.push_back( (*rit).second );
}

Expand All @@ -202,19 +202,19 @@ void MuonRecoOneHLT::analyze(const edm::Event& iEvent, const edm::EventSetup& iS
// if (_MuonEventFlag->on() && !(_MuonEventFlag->accept(iEvent,iSetup))) return;

// Check if Muon is Global
if(LeadingMuon[0].isGlobalMuon()) {
if((*LeadingMuon[0]).isGlobalMuon()) {
LogTrace(metname)<<"[MuonRecoOneHLT] The mu is global - filling the histos";
if(LeadingMuon[0].isTrackerMuon() && LeadingMuon[0].isStandAloneMuon()) muReco->Fill(1);
if(!(LeadingMuon[0].isTrackerMuon()) && LeadingMuon[0].isStandAloneMuon()) muReco->Fill(2);
if(!LeadingMuon[0].isStandAloneMuon())
if((*LeadingMuon[0]).isTrackerMuon() && (*LeadingMuon[0]).isStandAloneMuon()) muReco->Fill(1);
if(!((*LeadingMuon[0]).isTrackerMuon()) && (*LeadingMuon[0]).isStandAloneMuon()) muReco->Fill(2);
if(!(*LeadingMuon[0]).isStandAloneMuon())
LogTrace(metname)<<"[MuonRecoOneHLT] ERROR: the mu is global but not standalone!";

// get the track combinig the information from both the Tracker and the Spectrometer
reco::TrackRef recoCombinedGlbTrack = LeadingMuon[0].combinedMuon();
reco::TrackRef recoCombinedGlbTrack = (*LeadingMuon[0]).combinedMuon();
// get the track using only the tracker data
reco::TrackRef recoTkGlbTrack = LeadingMuon[0].track();
reco::TrackRef recoTkGlbTrack = (*LeadingMuon[0]).track();
// get the track using only the mu spectrometer data
reco::TrackRef recoStaGlbTrack = LeadingMuon[0].standAloneMuon();
reco::TrackRef recoStaGlbTrack = (*LeadingMuon[0]).standAloneMuon();

etaGlbTrack[0]->Fill(recoCombinedGlbTrack->eta());
etaGlbTrack[1]->Fill(recoTkGlbTrack->eta());
Expand All @@ -233,11 +233,11 @@ void MuonRecoOneHLT::analyze(const edm::Event& iEvent, const edm::EventSetup& iS
ptGlbTrack[2]->Fill(recoStaGlbTrack->pt());
}
// Check if Muon is Tight
if (muon::isTightMuon(LeadingMuon[0], vtx) ) {
if (muon::isTightMuon((*LeadingMuon[0]), vtx) ) {

LogTrace(metname)<<"[MuonRecoOneHLT] The mu is tracker only - filling the histos";

reco::TrackRef recoCombinedGlbTrack = LeadingMuon[0].combinedMuon();
reco::TrackRef recoCombinedGlbTrack = (*LeadingMuon[0]).combinedMuon();

etaTight->Fill(recoCombinedGlbTrack->eta());
phiTight->Fill(recoCombinedGlbTrack->phi());
Expand All @@ -246,13 +246,13 @@ void MuonRecoOneHLT::analyze(const edm::Event& iEvent, const edm::EventSetup& iS
}

// Check if Muon is Tracker but NOT Global
if(LeadingMuon[0].isTrackerMuon() && !(LeadingMuon[0].isGlobalMuon())) {
if((*LeadingMuon[0]).isTrackerMuon() && !((*LeadingMuon[0]).isGlobalMuon())) {
LogTrace(metname)<<"[MuonRecoOneHLT] The mu is tracker only - filling the histos";
if(LeadingMuon[0].isStandAloneMuon()) muReco->Fill(3);
if(!(LeadingMuon[0].isStandAloneMuon())) muReco->Fill(4);
if((*LeadingMuon[0]).isStandAloneMuon()) muReco->Fill(3);
if(!((*LeadingMuon[0]).isStandAloneMuon())) muReco->Fill(4);

// get the track using only the tracker data
reco::TrackRef recoTrack = LeadingMuon[0].track();
reco::TrackRef recoTrack = (*LeadingMuon[0]).track();

etaTrack->Fill(recoTrack->eta());
phiTrack->Fill(recoTrack->phi());
Expand All @@ -261,19 +261,19 @@ void MuonRecoOneHLT::analyze(const edm::Event& iEvent, const edm::EventSetup& iS
}

// Check if Muon is STA but NOT Global
if(LeadingMuon[0].isStandAloneMuon() && !(LeadingMuon[0].isGlobalMuon())) {
if((*LeadingMuon[0]).isStandAloneMuon() && !((*LeadingMuon[0]).isGlobalMuon())) {
LogTrace(metname)<<"[MuonRecoOneHLT] The mu is STA only - filling the histos";
if(!(LeadingMuon[0].isTrackerMuon())) muReco->Fill(5);
if(!((*LeadingMuon[0]).isTrackerMuon())) muReco->Fill(5);

// get the track using only the mu spectrometer data
reco::TrackRef recoStaTrack = LeadingMuon[0].standAloneMuon();
reco::TrackRef recoStaTrack = (*LeadingMuon[0]).standAloneMuon();

etaStaTrack->Fill(recoStaTrack->eta());
phiStaTrack->Fill(recoStaTrack->phi());
chi2OvDFStaTrack->Fill(recoStaTrack->normalizedChi2());
ptStaTrack->Fill(recoStaTrack->pt());
}
// Check if Muon is Only CaloMuon
if(LeadingMuon[0].isCaloMuon() && !(LeadingMuon[0].isGlobalMuon()) && !(LeadingMuon[0].isTrackerMuon()) && !(LeadingMuon[0].isStandAloneMuon()))
if((*LeadingMuon[0]).isCaloMuon() && !((*LeadingMuon[0]).isGlobalMuon()) && !((*LeadingMuon[0]).isTrackerMuon()) && !((*LeadingMuon[0]).isStandAloneMuon()))
muReco->Fill(6);
}