Flattening Event angle SW using track correction

Motivation

• Asymmetry of detector performance.

0. Prepare

##Make your working space before step##
cvs co UserCode/geonmo/HeavyIon
mv UserCode/geonmo/HeavyIon/* ./
rm -rf UserCode
cvs up -dP

1. Make track eta-phi map

• prepare

scram b
cd HITracksAna/HITrackAna01

• Make a table of eta-phi map.
  • Edit the hitrackana01_cfg.py
  • And Run. I tested this using 100k events.

process.CondDBCommon.connect = "sqlite_file:/net/hisrv0001/home/icpark/geonmo/UserCode/CmsHi/JulyExercise/data/CentralityTables.db" // Modify this line using your path.

process.PoolDBESSource = cms.ESSource("PoolDBESSource",
                                      process.CondDBCommon,
                                      toGet = cms.VPSet(cms.PSet(record = cms.string('HeavyIonRcd'),
                                      tag = cms.string('HFhits20_DataJulyExercise_AMPT2760GeV_MC_37Y_V5_ZS_v0')                  // Maybe only this line
                                      )))
process.TFileService = cms.Service("TFileService",
    fileName = cms.string('100730_trackana2.root')    // This is output file's name. 
)

• Using makeTable.C, you get the weight table of eta-phi map for correct.
  • default file name is "table.root"

2. Make resolution data / [or make flattened Evt angle]

• Edit hiresolutionplot_cfg.py

process.load("HIResolution.HIResolution.hiresolution_cfi")   // I think it need not modification
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10000) )
process.HIResolutionana = cms.EDAnalyzer('HIResolutionAna',
        centrality_bin = cms.int32(20)                             // It is very important value for histogram binning.
)

process.evtPlaneTrkCorr = cms.EDProducer('HICorrEvtPlane', inputfile = cms.string("table.root")      // Input file is weight table.
)
process.TFileService = cms.Service("TFileService",
                           fileName = cms.string('resolution_plot.root')                           // Define the name of output file.
)

• Then, cmsRun this file. We got the data about resolution profile. In addition, we got the data about Flattened Evt angle if we add output module.

3. Make V2 profile using RP method.

• Edit hiflow_cfg.py

process.evtPlaneTrkCorr = cms.EDProducer('HICorrEvtPlane', inputfile = cms.string("table.root"))     // it is a same on previous  step.
process.RPMethod = cms.EDAnalyzer('HIRPMethod',                                                                          
                                  centrality_bin = cms.int32(20),                                                                           // It is a same on previous step.
                                  resolutionfile= cms.string("resolution_plot.root")                                                // resolution plot file is needed. it is made on previous step.
)
process.TFileService = cms.Service("TFileService",
                           fileName = cms.string('V2_RP.root')                                                                          // Output histogram. It is a final data.

-- GunmoRyu - 10 Aug 2010