Flattening Event angle SW using track correction

Complete: 1


  • 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",
                                      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

Edit | Attach | Watch | Print version | History: r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r3 - 2010-08-11 - GunmoRyu
This site is powered by the TWiki collaboration platform Powered by Perl This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding KoreaCmsWiki? Send feedback