Upsilon production in Pb Pb Heavy ion collisions:

Comparing heavy ion muon reconstruction algorithms with standard pp muon reconstruction algorithms:

Try I. Using standard pp raw data and run heavy ion reconstruction algorithm on it.

1. Heavy Ion reconstruction is broken in CMSSW_3_1_2. All official production for Upsilon to Mu Mu is in CMSSW_3_1_2.

2. We Can make reco data in CMSSW_3_3_1 from raw data of CMSSW_3_1_2 but we have to change global tag. To make pp raw data MC_31X_V3 global tag is used and to run HI reconstruction on I use reconstruct_PbPb_cfg.py which have global tag MC_31X_v8::All. I used this global tag otherwise there were some errors in runing. I dont know what effect it will have.

3. PAT is broken on HI output file in CMSSW_3_3_1. So we compare efficiencies using reco:globalMuons.

4. Trying to calculate a rough idea for muon reconstruction efficiency with pp reconstruction and heavy ion reconstruction algorithms

i) Using official summer09 production Raw files and then reconstructing them with HI reconstruction algorithm

total no of events in Raw file on which HI algo is run = 50000 Now file reconstructed by HIAlgo have.

Genrated Dimuon pairs =51290

Reconstructed Dimuon pairs =16605

Overall dimuon reconstruction eff =32.17%

Total no of genrated muons =101067

Total no of reconstructed muon =57404

Overall muon reconstruction eff =56.7%

ii) using pp raw with pp Reco ie using summer09 official reconstructed data

total no of events run in =50000

total no of genrated dimuon pairs =51272

total no of reconstructed dimuon pairs =24218

dimuon reconstruction eff =47.23%

total no of genrated muons =101049

total no of reconstructed muons =72582

over all muon reconstruction eff =71.82%

I have also made some plots of efficiency vs dimuon pt,y and efficiency vs singlemuon pt and eta with pp and HI reconstruction algorithms. To see these plots one should

take care that at high pT statstics is very small to make any conclusion :).

muon reconstruction eff Vs single muon eta

muon reconstruction eff vs single muon pT

muon reconstruction eff vs dimuon pT

muon reconstruction eff vs Dimuon Y

To Be noted is here that there are some genfilters in the path which are

mumugenfilter = cms.EDFilter("MCParticlePairFilter",
Status = cms.untracked.vint32(1, 1),
MinPt = cms.untracked.vdouble(0.0, 0.0),
MaxEta = cms.untracked.vdouble(2.5, 2.5),
MinEta = cms.untracked.vdouble(-2.5, -2.5),
ParticleCharge = cms.untracked.int32(-1),
ParticleID1 = cms.untracked.vint32(13),
ParticleID2 = cms.untracked.vint32(13)
)

mugenfilter = cms.EDFilter("MCSmartSingleParticleFilter",
MaxDecayRadius = cms.untracked.vdouble(1500.0, 1500.0),
Status = cms.untracked.vint32(1, 1),
MinPt = cms.untracked.vdouble(2.5, 2.5),
ParticleID = cms.untracked.vint32(13, -13),
MaxEta = cms.untracked.vdouble(2.5, 2.5),
MinEta = cms.untracked.vdouble(-2.5, -2.5),
MaxDecayZ = cms.untracked.vdouble(3000.0, 3000.0),
MinDecayZ = cms.untracked.vdouble(-3000.0, -3000.0)
)

I dont know what exactly they are doing but what I guess second filter is selecting events which are having pt greater than 2.5 GeV. So if there is any increase for

low Pt muon reconstruction in HI Algo that we can not see.

Try II. Making a total offline reconstruction file (no HLT) for both pp and HI reconstruction:

Using pp Upsilon production cff file and cmsDriver:

To make reconstruction file I use Configuration/GenProduction/python/PYTHIA6_Upsilon1S_10TeV_cff.py and cmsDriver in CMSSW_3_3_1 as:

cmsDriver.py Configuration/GenProduction/python/PYTHIA6_Upsilon1S_10TeV_cff.py hiReco -s GEN,SIM,DIGI,L1,DIGI2RAW,RAW2DIGI,RECO --scenario HeavyIons --conditions FrontierConditions _GlobalTag,MC_31X_V8::All --fileout UpsilonMuMuReco.root --number 100 --mc --no_exec --datatier GEN-SIM-DIGI-RECO --eventcontent RECOSIM --processName RECO

and for pp reconstruction

cmsDriver.py Configuration/GenProduction/python/PYTHIA6_Upsilon1S_10TeV_cff.py -s GEN,SIM,DIGI,L1,DIGI2RAW,RAW2DIGI,RECO --conditions FrontierConditions _GlobalTag, MC_31X_V8::All --fileout UpsilonMuMuReco.root --number 100 --mc --no_exec --datatier GEN-SIM-DIGI-RECO --eventcontent RECOSIM --processName RECO

COMPARISON OF HI AND PP RECONSTRUCTION ALGORITHMS FOR PRIVATE UPSILON PRODUCTION :

reconstruction is run for 100000 events using both HI reco and pp reco: what are the main diff b/w pp and heavy ion reconstruction files

//---------------------------------------------------------------------------------------------------------------
23c23
< process.load('Configuration/StandardSequences/ReconstructionHeavyIons_cff')
---
> process.load('Configuration/StandardSequences/Reconstruction_cff')
26c26
< process.load('Configuration/EventContent/EventContentHeavyIons_cff')
---
> process.load('Configuration/EventContent/EventContent_cff')

< process.generation_step = cms.Path(process.pgen_hi)
>>> process.pgen_hi
cms.Sequence(randomEngineStateProducer+VtxSmeared+hiGenParticles)
---
> process.generation_step = cms.Path(process.pgen)
>>> process.pgen
cms.Sequence(randomEngineStateProducer+VtxSmeared+genParticles+genParticlesForJets*sisCone5GenJets+sisCone7GenJets+kt4GenJets+kt6GenJets+ite\
rativeCone5GenJets+ak5GenJets+ak7GenJets+genCandidatesForMET+genParticlesForMETAllVisible*genMetCalo+genMetCaloAndNonPrompt+genMetTrue*genMe\
tIC5GenJets)
162c164

< process.reconstruction_step = cms.Path(process.reconstructionHeavyIons)
>>> process.reconstructionHeavyIons
cms.Sequence(offlineBeamSpot*siPixelClusters*siPixelRecHits*siStripZeroSuppression*siStripClusters*siStripMatchedRecHits+dt1DRecHits*dt4DSegments

+csc2DRecHits*cscSegments+rpcRecHits+lumiProducer*ecalGlobalUncalibRecHit*ecalDetIdToBeRecovered*ecalRecHit+ecalPreshowerRecHit*hbherec
o+hfreco+horeco+zdcreco*hiPixel3ProtoTracks*hiPixelMedianVertex*hiSelectedProtoTracks*hiPixelAdaptiveVertex*hiSelectedVertex*hiPixel3PrimTra
cks*hiPixelTrackSeeds*hiPrimTrackCandidates*hiGlobalPrimTracks*hiSelectedTracks*islandBasicClusters*islandSuperClusters*correctedIslandBarre
lSuperClusters*correctedIslandEndcapSuperClusters*hybridSuperClusters*correctedHybridSuperClusters*correctedEndcapSuperClustersWithPreshower
*preshowerClusterShape*dynamicHybridSuperClusters*correctedDynamicHybridSuperClusters*multi5x5BasicClusters*multi5x5SuperClusters*multi5x5Su
perClustersWithPreshower*correctedMulti5x5SuperClustersWithPreshower*multi5x5PreshowerClusterShape*towerMaker+towerMakerWithHO*caloTowers*it
erativeConePu5CaloJets*ancientMuonSeed*standAloneMuons*globalMuons*tevMuons*glbTrackQual*muons*calomuons*muIsoDepositTk+muIsoDepositCalByAss
ociatorTowers+muIsoDepositJets*photonCore+photons*hiCentrality*hiEvtPlane)

> process.reconstruction_step = cms.Path(process.reconstruction)

>>> process.reconstruction
cms.Sequence(siPixelClusters*siPixelRecHits*siStripZeroSuppression*siStripClusters*siStripMatchedRecHits+dt1DRecHits*dt4DSegments+csc2DRecHits*cscSegments+rp
cRecHits+ecalGlobalUncalibRecHit*ecalDetIdToBeRecovered*ecalRecHit+ecalPreshowerRecHit+hbhereco+hfreco+horeco+zdcreco+towerMakerPF*particleFlowRecHitECAL*par
ticleFlowClusterECAL*particleFlowRecHitHCAL*particleFlowClusterHCAL+particleFlowClusterHFHAD+particleFlowClusterHFEM*particleFlowRecHitPS*particleFlowCluster
PS+lumiProducer*offlineBeamSpot+pixelTracks*pixelVertices*newSeedFromTriplets*newTrackCandidateMaker*preFilterZeroStepTracks*zeroStepWithLooseQuality*zeroSte
pWithTightQuality*zeroStepTracksWithQuality*zeroStepFilter*newClusters*newPixelRecHits*newStripRecHits*newSeedFromPairs*stepOneTrackCandidateMaker*preFilterS
tepOneTracks*firstStepWithLooseQuality*firstStepWithTightQuality*preMergingFirstStepTracksWithQuality*firstStepTracksWithQuality*firstfilter*secClusters*secP
ixelRecHits*secStripRecHits*secTriplets*secTrackCandidates*secWithMaterialTracks*secStepVtxLoose*secStepTrkLoose*secStepLoose*secStepVtxTight*secStepTrkTight
*secStepTight*secStepVtx*secStepTrk*secStep*secfilter*thClusters*thPixelRecHits*thStripRecHits*thPLSeeds*thTrackCandidates*thWithMaterialTracks*thStepVtxLoos
e*thStepTrkLoose*thStepLoose*thStepVtxTight*thStepTrkTight*thStepTight*thStepVtx*thStepTrk*thStep*thfilter*fourthClusters*fourthPixelRecHits*fourthStripRecHi
ts*fourthPLSeeds*fourthTrackCandidates*fourthWithMaterialTracks*pixellessStepLoose*pixellessStepTight*pixellessStep*fourthfilter*fifthClusters*fifthPixelRecH
its*fifthStripRecHits*fifthSeeds*fifthTrackCandidates*fifthWithMaterialTracks*tobtecStepLoose*tobtecStepTight*tobtecStep*merge2nd3rdTracks*merge4th5thTracks*
iterTracks*generalTracks*newCombinedSeeds*dedxTruncated40+dedxMedian+dedxHarmonic2+hybridSuperClusters*correctedHybridSuperClusters*multi5x5BasicClusters*mul
ti5x5SuperClusters*multi5x5SuperClustersWithPreshower*correctedMulti5x5SuperClustersWithPreshower*multi5x5PreshowerClusterShape+towerMaker+towerMakerWithHO*o
fflinePrimaryVertices*offlinePrimaryVerticesWithBS*generalV0Candidates*sisCone5CaloJets+sisCone7CaloJets+kt4CaloJets+kt6CaloJets+iterativeCone5CaloJets+ak5Ca
loJets+ak7CaloJets*ak5JetID+ak7JetID+sc5JetID+sc7JetID+ic5JetID+kt4JetID+kt6JetID+ancientMuonSeed*standAloneMuons*globalMuons*tevMuons*glbTrackQual*muons*cal
omuons*SETMuonSeed*standAloneSETMuons*globalSETMuons*muIsoDepositTk+muIsoDepositCalByAssociatorTowers+muIsoDepositJets*muidTrackerMuonArbitrated+muidAllArbit
rated+muidGlobalMuonPromptTight+muidTMLastStationLoose+muidTMLastStationTight+muidTM2DCompatibilityLoose+muidTM2DCompatibilityTight+muidTMOneStationLoose+mui
dTMOneStationTight+muidTMLastStationOptimizedLowPtLoose+muidTMLastStationOptimizedLowPtTight+muidGMTkChiCompatibility+muidGMStaChiCompatibility+muidGMTkKinkT
ight+trackerDrivenElectronSeeds*ecalDrivenElectronSeeds*electronMergedSeeds*electronCkfTrackCandidates*electronGsfTracks*ic5JetTracksAssociatorAtVertex*itera
tiveCone5JetTracksAssociatorAtVertex*iterativeCone5JetTracksAssociatorAtCaloFace*iterativeCone5JetExtender+sisCone5JetTracksAssociatorAtVertex*sisCone5JetTra
cksAssociatorAtCaloFace*sisCone5JetExtender+kt4JetTracksAssociatorAtVertex*kt4JetTracksAssociatorAtCaloFace*kt4JetExtender+ak5JetTracksAssociatorAtVertex*ak5
JetTracksAssociatorAtCaloFace*ak5JetExtender+ak7JetTracksAssociatorAtVertex*ak7JetTracksAssociatorAtCaloFace*ak7JetExtender*caloRecoTauTagInfoProducer*caloRe
coTauProducer*caloRecoTauDiscriminationByLeadingTrackFinding*caloRecoTauDiscriminationByLeadingTrackPtCut*caloRecoTauDiscriminationByIsolation*caloRecoTauDis
criminationAgainstElectron*pfTrackElec*particleFlowBlock*particleFlow*pfElectronTranslator+pfElectronInterestingEcalDetIdEB+pfElectronInterestingEcalDetIdEE*
gsfElectronCores*gsfElectrons*conversionTrackCandidates*ckfOutInTracksFromConversions*ckfInOutTracksFromConversions*conversions*trackerOnlyConversions*photon
Core+photons*interestingEleIsoDetIdEB*interestingEleIsoDetIdEE*interestingGamIsoDetIdEB*interestingGamIsoDetIdEE*PhotonIDProd*eidRobustLoose+eidRobustTight+e
idRobustHighEnergy+eidLoose+eidTight*hfEMClusters+hfRecoEcalCandidate*met+metNoHF+metHO+metNoHFHO+calotoweroptmaker+metOpt+metOptNoHF+calotoweroptmakerWithHO
+metOptHO+metOptNoHFHO+htMetSC5+htMetSC7+htMetKT4+htMetKT6+htMetIC5+muonMETValueMapProducer+corMetGlobalMuons+muonTCMETValueMapProducer+tcMet+hcalnoise*inter
estingEcalDetIdEB*interestingEcalDetIdEE*reducedEcalRecHitsEB*reducedEcalRecHitsEE*impactParameterTagInfos*(trackCountingHighEffBJetTags+trackCountingHighPur
BJetTags+jetProbabilityBJetTags+jetBProbabilityBJetTags+secondaryVertexTagInfos*(simpleSecondaryVertexBJetTags+combinedSecondaryVertexBJetTags+combinedSecond
aryVertexMVABJetTags))+softPFElectrons*softElectronTagInfos*(softElectronByIP3dBJetTags+softElectronByPtBJetTags)+softMuonTagInfos*(softMuonBJetTags+softMuon
ByIP3dBJetTags+softMuonByPtBJetTags)*sisCone5PFJets+sisCone7PFJets+kt4PFJets+kt6PFJets+iterativeCone5PFJets+ak5PFJets+ak7PFJets*pfMet*ic5PFJetTracksAssociato
rAtVertex*pfRecoTauTagInfoProducer*shrinkingConePFTauProducer*shrinkingConePFTauDecayModeProducer*shrinkingConePFTauDecayModeIndexProducer*shrinkingConePFTau
DiscriminationByLeadingTrackFinding*shrinkingConePFTauDiscriminationByLeadingTrackPtCut*shrinkingConePFTauDiscriminationByLeadingPionPtCut*shrinkingConePFTau
DiscriminationByIsolation*shrinkingConePFTauDiscriminationByTrackIsolation*shrinkingConePFTauDiscriminationByECALIsolation*shrinkingConePFTauDiscriminationBy
IsolationUsingLeadingPion*shrinkingConePFTauDiscriminationByTrackIsolationUsingLeadingPion*shrinkingConePFTauDiscriminationByECALIsolationUsingLeadingPion*sh
rinkingConePFTauDiscriminationAgainstElectron*shrinkingConePFTauDiscriminationAgainstMuon+shrinkingConePFTauDiscriminationByTaNC*shrinkingConePFTauDiscrimina
tionByTaNCfrOnePercent*shrinkingConePFTauDiscriminationByTaNCfrHalfPercent*shrinkingConePFTauDiscriminationByTaNCfrQuarterPercent*shrinkingConePFTauDiscrimin
ationByTaNCfrTenthPercent+fixedConeHighEffPFTauProducer*fixedConeHighEffPFTauDiscriminationByLeadingTrackFinding*fixedConeHighEffPFTauDiscriminationByLeading
TrackPtCut*fixedConeHighEffPFTauDiscriminationByLeadingPionPtCut*fixedConeHighEffPFTauDiscriminationByIsolation*fixedConeHighEffPFTauDiscriminationByTrackIso
lation*fixedConeHighEffPFTauDiscriminationByECALIsolation*fixedConeHighEffPFTauDiscriminationByIsolationUsingLeadingPion*fixedConeHighEffPFTauDiscriminationB
yTrackIsolationUsingLeadingPion*fixedConeHighEffPFTauDiscriminationByECALIsolationUsingLeadingPion*fixedConeHighEffPFTauDiscriminationAgainstElectron*fixedCo
neHighEffPFTauDiscriminationAgainstMuon+fixedConePFTauProducer*fixedConePFTauDiscriminationByLeadingTrackFinding*fixedConePFTauDiscriminationByLeadingTrackPt
Cut*fixedConePFTauDiscriminationByLeadingPionPtCut*fixedConePFTauDiscriminationByIsolation*fixedConePFTauDiscriminationByTrackIsolation*fixedConePFTauDiscriminationByECALIsolation*

fixedConePFTauDiscriminationByIsolationUsingLeadingPion*fixedConePFTauDiscriminationByTrackIsolationUsingLeadingPion*fixedConePFTauDis
criminationByECALIsolationUsingLeadingPion*fixedConePFTauDiscriminationAgainstElectron*fixedConePFTauDiscriminationAgainstMuon*logErrorHarvester)


WAO


Now in this reconstruction We are not using any filter except onia filter

process.oniafilter = cms.EDFilter("PythiaFilter",
MaxEta = cms.untracked.double(1000.0),
Status = cms.untracked.int32(2),
MinEta = cms.untracked.double(-1000.0),
MinPt = cms.untracked.double(0.0),
ParticleID = cms.untracked.int32(553)
)

A rough idea about over all efficiencies:
Total no of events 100000

Data at CASTOR: /castor/cern.ch/user/k/kumarv/CMSSW330/UpsilonMuMu/UpsilonMuMu_HIReco_*.root

/castor/cern.ch/user/k/kumarv/CMSSW330/UpsilonMuMu/UpsilonMuMu_ppReco_*.root

And some primary plots are here ?

muon reconstruction eff Vs single muon eta

muon reconstruction eff vs single muon pT

muon reconstruction eff vs Dimuon pT

muon reconstruction eff vs Dimuon Y

One Query I Have what is the difference between pgen and hi_pgen and why all the modules which are coming in pgen are not coming in hi_pgen.

Puting genrated distributions and plots required by Catherine:

Few plots with upsilon inv mass cut


# reconstructed Upsilon / # generated Upsilons versus y

# reconstructed Upsilon / # generated Upsilon WITHIN |y|<2.4 versus pT

# reconstructed Upsilon WITHIN |y|<0.8 / # generated Upsilon WITHIN |y|<0.8 versus pT

HI reco eff/pp reco eff vs Upsilon Eta

HI reco eff/ pp reco eff vs Upsilon Y

HI reco eff / pp reco eff vs Upsilon pT

Working on to put error bars on eff histograms.Done

PlotsWithErrorBars

ReconstructionOfUpsilonFlatInPt

HowToRunCondorInMit

Topic attachments
I Attachment History Action Size Date Who Comment
Unknown file formateps EffRatioVsEta.eps r1 manage 8.4 K 2009-11-12 - 20:13 VineetKumar HI and PP reco eff ratio vs Upsilon eta
Unknown file formateps EffRatioVsPt.eps r1 manage 8.1 K 2009-11-12 - 20:13 VineetKumar Eff ratio vs Upsilon Pt
Unknown file formateps EffRatioVsY.eps r1 manage 7.6 K 2009-11-12 - 20:14 VineetKumar Eff ratio vs upsilon rapidity
Unknown file formateps EffVsPtDimuonPr.eps r2 r1 manage 7.5 K 2009-11-11 - 19:00 VineetKumar  
Unknown file formateps EffVsPtDimuonSt.eps r2 r1 manage 7.7 K 2009-11-11 - 19:16 VineetKumar  
Unknown file formateps EffVsPtGenYLT2.4.eps r1 manage 8.7 K 2009-11-12 - 20:15 VineetKumar # reconstructed Upsilon / # generated Upsilon WITHIN |y|<2.4 versus pT
Unknown file formateps EffVsPtYLT08.eps r1 manage 8.6 K 2009-11-12 - 20:16 VineetKumar # reconstructed Upsilon WITHIN |y|<0.8 / # generated Upsilon WITHIN |y|<0.8 versus pT
PNGpng EffVsYDiMuon.png r1 manage 11.1 K 2009-11-12 - 20:17 VineetKumar - # reconstructed Upsilon / # generated Upsilons versus y
Unknown file formateps EffVsYDimuon.eps r1 manage 7.7 K 2009-11-11 - 13:03 VineetKumar  
Unknown file formateps EffVsYDimuonPr.eps r1 manage 7.7 K 2009-11-11 - 14:58 VineetKumar  
Unknown file formateps EffVsYDimuonSt.eps r1 manage 8.1 K 2009-11-11 - 14:58 VineetKumar  
Unknown file formateps EffVsetaSingleMuPr.eps r1 manage 8.5 K 2009-11-11 - 14:56 VineetKumar  
Unknown file formateps EffVsetaSingleMuSt.eps r1 manage 8.5 K 2009-11-11 - 14:56 VineetKumar  
Unknown file formateps EffVsptSingleMuPr.eps r2 r1 manage 8.9 K 2009-11-11 - 19:00 VineetKumar  
Unknown file formateps EffVsptSingleMuSt.eps r2 r1 manage 7.5 K 2009-11-11 - 19:10 VineetKumar  
PNGpng GenDiMuEta.png r2 r1 manage 7.4 K 2009-11-20 - 20:50 VineetKumar GenDiMuEta
PNGpng GenDiMuPt.png r2 r1 manage 6.9 K 2009-11-20 - 20:50 VineetKumar gen di mu pt
PNGpng GenDiMuY.png r2 r1 manage 8.0 K 2009-11-20 - 20:50 VineetKumar  
PNGpng GenMuEta.png r2 r1 manage 7.1 K 2009-11-20 - 20:51 VineetKumar  
PNGpng GenMuPt.png r2 r1 manage 7.1 K 2009-11-20 - 20:51 VineetKumar  
Unknown file formateps RRR.eps r1 manage 7.7 K 2009-11-11 - 19:16 VineetKumar  
Edit | Attach | Watch | Print version | History: r8 < r7 < r6 < r5 < r4 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r8 - 2009-12-17 - VineetKumar
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback