How to run CaloTrackCorr (Cath's receipee)

more details here : PhotonHadronCorrelations#Analysis_code

On the grid

  • in /alice/cern.ch/user/c/csilvest/PartIsoCorr/macros/JJ/
  • task AddTaskJJ.C
  • ana.C

  • MC Jet-Jet: LHC12a15e : pp, Jet-Jet Pythia6 (un-quenched), 0.5T, 2760 GeV, LHC11h anchors, 11 Pt bins, ID #300 link
    • Il faut que je creer les collection xml de tous les runs ici /alice/sim/2012/LHC12a15e ?
      • testing :=find -x collection /alice/sim/2012/LHC12a15e/169838/0/001 AliAOD.root > 0.xml=
cp -n file:0.xml .

      • a la main "find -x collection path AliESDs.root > nom.xml", ou tu utilises un script comme celle ci
/alice/cern.ch/user/g/gconesab/collection/scripts/makeCollections.sh, makeCollectionsMC.sh
    • example : /alice/cern.ch/user/g/gconesab/collection/
    • run makeCollectionsMC.sh 0  10 simu_pp2.76_JJ_EMCAL /alice/sim/2012/LHC12a15e/170040/ AOD
  • $2 =simu_pp2.76_JJ_EMCAL/, $1 =169838, 170040
  • > /alice/cern.ch/user/c/csilvest/PartIsoCorr/jdl/ >submit ana.jdl 0 simu_pp2.76_JJ_EMCAL/169838

at LPSC

  • runs with trunk
  • look at my working directory /cebaf/csilvestre/work/CaloTrackCorrelations/macros/ChargedNeutral/Corr
  • copy anaAllAOD.pl, ana.C, AddTaskCaloTrackCorr.C to your working directory
  • create a direcotry testing (or change the directory name in anaAllAOD.pl). The output root files will be copied in this testing directory.
    • anaAllAOD.pl : execute ana.C for each good run number. To run it do perl anaAllAOD.pl pp 0 300 (on pp, from the first runnumber to the 300th runnumber (or the last one if there are less runs than 300)
    • ana.C : runs on ESD the clusteriser and the correlation analysis
    • AddTaskCaloTrackCorr.C : defines the correlation analysis (cf maker) and the isolation criteria in ConfigureIsolationAnalysis()
      • if you want to change cone size or ptThreshold or the method change one of those parametters
    ic->SetPtThreshold(0.5);
    ic->SetConeSize(0.4);
   ic->SetPtFraction(0.1);
   ic->SetSumPtThreshold(1.0) ;
   ic->SetICMethod(thresType);
      • to set the type of particles used to isolate change SetParticleTypeInCone
  • It takes the night to run on all pp data. Then one needs to merge (using hadd for example).

-- CatherineSilvestre - 28-Apr-2012

Edit | Attach | Watch | Print version | History: r11 | r8 < r7 < r6 < r5 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r6 - 2012-06-04 - CatherineSilvestre
 
    • 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-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback