Monica Vazquez Acosta's working wiki

MAGIC

.tschrc

setenv ROOTSYS:<Path of ROOT directory>
setenv MARSSYS:<Path of MARS directory>
setenv LD_LIBRARY_PATH:$ROOTSYS/lib:$MARSSYS:$LD_LIBRARY_PATH
setenv DYLD_LIBRARY_PATH $LD_LIBRARY_PATH
setenv PATH "$PATH":$ROOTSYS/bin:$MARSSYS:$PATH

### setting up cvs in pic
setenv CVSROOT ":ext:USERNAME@cvs.magic.pic.es:/CVS"
setenv CVS_RSH "ssh"

.rootrc

Unix.*.Root.MacroPath: .:$(HOME)/macros:$(MARSSYS)/macros
Unix.*.Root.DynamicPath: .:$(MARSSYS)/lib
Unix.*.Gui.IconPath: $(MARSSYS)
Rint.Logon: rootlogon.C

Root

Version used online and recommended on the wiki:

wget http://root.cern.ch/download/root_v5.34.04.source.tar.gz
./configure --enable-builtin-zlib
./make

Version recommended by J. Sitarek

wget http://root.cern.ch/download/root_v5.34.21.source.tar.gz
./configure --enable-builtin-zlib
./make

Version used by Marina

wget http://root.cern.ch/download/root_v5.34.18.source.tar.gz
./configure --enable-builtin-zlib
./make

MARS

Version used in 4th Magic school: Mars_V2.14.3.
Latest recommended version: V2.14.2
cvs co -r Mars_V2-14-2 Mars
./make

MARS Commands

quate (star)
quate --useM1 --events --sizelow=50

  • Get config file:
          cp $MARSSYS/mrcfiles/quate.rc .
          

  • Edit config file to reduce zenith angle:
            emacs -nw  quate.rc 
            #ZdMax: 100
            ZdMax: 35
         

quate --config=./quate.rc --useM1 --applycuts

superstar
cp $MARSSYS/mrcfiles/superstar.rc .

superstar -f -q -b --config=./superstar.rc --ind1="../star/all/good/20*M1*root" --ind2="../star/all/good/20*M2*root" --out=./ --log="./superstar.log"

quate (superstar)
quate -b -f --stereo --events --log="./quate_1.log"
quate -b -f --config=./quate.rc --stereo --applycuts --log="./quate_cuts.log"

selectmc
OBSOLETE

selectmc -joinmc --split=0.3 --pathMC='./*.root' --out=./ --outname=LowZenithMC.root --log=log.txt

coach
cp $MARSSYS/mrcfiles/coach.rc .

RF.zdmin: 5.
#RF.zdmax: 50.
RF.zdmax: 35.

###GA_*1.root (MC training sample)
RF.mcdata: /scratch/monicava/MAGICMC/LOWZENITH/GA_*_1.root  
RF.data:  /scratch/monicava/MAGICDATA/1H0323+342/superstar/2013*.root
RF.outpath: /scratch/monicava/MAGICRF/

$MARSSYS/coach -f -b --config=./coach.rc -LUTs >& lut.txt &
$MARSSYS/coach -f -b --config=./coach.rc -RFgh >& RFgh.txt &
$MARSSYS/coach -f -b --config=./coach.rc -RFdisp >& RFdisp.txt &

$MARSSYS/coach -f -b --config=./coach_mod.rc -LUTs >& ST.03.05_MOD/lut.txt &
$MARSSYS/coach -f -b --config=./coach_mod.rc -RFgh >& ST.03.05_MOD/RFgh.txt &
$MARSSYS/coach -f -b --config=./coach_mod.rc -RFdisp >& ST.03.05_MOD/RFdisp.txt &


melibea
cp $MARSSYS/mrcfiles/melibea_stereo.rc .

# MJMelibea.CalcTrans: yes
 MJMelibea.CalcTrans: no

Online Data

$MARSSYS/melibea -f -q -b --stereo --rf --calcstereodisp --calc-disp-rf --calc-disp2-rf --disp-rf-sstrained -erec --rftree=/scratch/monicava/MAGICRF/RF.root --rfdisptree=/scratch/monicava/MAGICRF/disp1/DispRF.root --rfdisp2tree=/scratch/monicava/MAGICRF/disp2/DispRF.root --etab=/scratch/monicava/MAGICRF/Energy_Table.root --ind="/scratch/monicava/MAGICDATA/crab_clean/superstar/2013*root" --config=./melibea_stereo.rc --out=./ --log="./Melibea_OnData.log"

MC: GA_*2.root (MC test sample)

$MARSSYS/melibea -f -q -b -mc --stereo --rf --calcstereodisp --calc-disp-rf --calc-disp2-rf --disp-rf-sstrained -erec --rftree=/scratch/monicava/MAGICRF/RF.root --rfdisptree=/scratch/monicava/MAGICRF/disp1/DispRF.root --rfdisp2tree=/scratch/monicava/MAGICRF/disp2/DispRF.root --etab=/scratch/monicava/MAGICRF/Energy_Table.root  --ind="/scratch/monicava/MAGICMC/LOWZENITH2/GA*_2.root" --config=./melibea_stereo.rc --out=./ --log="./Melibea_MC.log"

$MARSSYS/melibea -f -q -b --stereo --rf --calcstereodisp --calc-disp-rf --calc-disp2-rf --disp-rf-sstrained -erec --rftree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05/RF.root --rfdisptree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05/disp1/DispRF.root --rfdisp2tree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05/disp2/DispRF.root --etab=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05/Energy_Table.root --ind="/scratch/monicava/MAGICDATA/CRAB/superstar/all/2014*root" --config=./melibea_stereo.rc --out=./ST0305MC_training --log="ST0305MC_training/Melibea_OnData.log" >& ST0305MC_training/melibea_ondata.txt &

$MARSSYS/melibea -f -q -b -mc --stereo --rf --calcstereodisp --calc-disp-rf --calc-disp2-rf --disp-rf-sstrained -erec --rftree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05/RF.root --rfdisptree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05/disp1/DispRF.root --rfdisp2tree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05/disp2/DispRF.root --etab=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05/Energy_Table.root --ind="/scratch/monicava/MAGICMC/ST.03.05/GA*_2.root" --config=./melibea_stereo.rc --out=./ST0305MC_training --log="ST0305MC_training/Melibea_MC.log" >& ST0305MC_training/melibea_mc.txt &

$MARSSYS/melibea -f -q -b --stereo --rf --calcstereodisp --calc-disp-rf --calc-disp2-rf --disp-rf-sstrained -erec --rftree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05_MOD/RF.root --rfdisptree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05_MOD/disp1/DispRF.root --rfdisp2tree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05_MOD/disp2/DispRF.root --etab=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05_MOD/Energy_Table.root --ind="/scratch/monicava/MAGICDATA/CRAB/superstar/all/2014*root" --config=./melibea_stereo.rc --out=./ST0305MC_MOD_training --log="ST0305MC_MOD_training/Melibea_OnData.log" >& ST0305MC_MOD_training/melibea_ondata.txt &

$MARSSYS/melibea -f -q -b -mc --stereo --rf --calcstereodisp --calc-disp-rf --calc-disp2-rf --disp-rf-sstrained -erec --rftree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05_MOD/RF.root --rfdisptree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05_MOD/disp1/DispRF.root --rfdisp2tree=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05_MOD/disp2/DispRF.root --etab=/scratch/monicava/MAGICANALYSIS/RF/ST.03.05_MOD/Energy_Table.root --ind="/scratch/monicava/MAGICMC/ST.03.05/GA*_2.root" --config=./melibea_stereo.rc --out=./ST0305MC_MOD_training --log="ST0305MC_MOD_training/Melibea_MC.log" >& ST0305MC_MOD_training/melibea_mc.txt &

odie
cp $MARSSYS/mrcfiles/odie.rc .

Odie.dataName: /scratch/monicava/MELIBEA/20*W0.40*root
# For low energy (LE) analysis
# Sensitivity ~ 1.3% Crab 
#Odie.outFileName: OdieOutput_LE.root
#Odie.analysisEpoch: Jul13
#Odie.eRange: LE
#Odie.signalCut: 0.02
#Odie.psf40: 0.07

# For full range (FR) analysis
# Sensitivity ~ 0.7% Crab 
#Odie.outFileName: OdieOutput_FR.root
#Odie.analysisEpoch: Jul13
#Odie.eRange: FR
#Odie.signalCut: 0.009
#Odie.psf40: 0.056

# For high energy (HE) analysis
# Sensitivity ~ 1.3% Crab
# The signal cut is likely to be enlarged, see program output for recommendation
Odie.outFileName: OdieOutput_HE.root
Odie.analysisEpoch: Jul13
Odie.eRange: HE
Odie.signalCut: 0.007
Odie.psf40: 0.043

$MARSSYS/odie --config=./odie.rc --log="odie_LE.log"

flute
cp $MARSSYS/mrcfiles/flute.rc .
flute.mcdata:  /scratch/monicava/MELIBEA/GA_*_Q_*2.root
flute.data:    /scratch/monicava/MELIBEA/2013*_Q_*.root

$MARSSYS/flute --config=./flute.rc --log="flute.log"

combunfold

cp $MARSSYS/mrcfiles/combunfold.rc .

MCallUnfold.FlagUnfold:     1

# -----------------------
#    FlagUnfold    = 1   Schmelling (minimization by Gauss-Newton method) 
#                    2   Tikhonov   (minimization by MINUIT)
#                    3   Bertero (appropriate for under-constrained case)
#                    4   Forward unfolding
#                    5   Schmelling (minimization by MINUIT)
#                    6   BerteroW (approproate for over-constrained case)
MCallUnfold.FlagUnfold:     1

Running the macro:

root
.x /scratch/monicava/MARS/Mars/macros/CombUnfold.C("./combunfold.rc")

caspar and skymap
cp $MARSSYS/mrcfiles/caspar.rc .
Caspar.dataName: /scratch/monicava/MAGICANALYSIS/crab_sensitivity/melibea/20*W0.40*root
$MARSSYS/caspar --config=./caspar.rc 
cp $MARSSYS/mrcfiles/datacubes.rc .
Datacubes.dataNames0: /scratch/monicava/MAGICANALYSIS/crab_sensitivity/melibea/2*W0.4*035.root 
Datacubes.dataNames1: /scratch/monicava/MAGICANALYSIS/crab_sensitivity/melibea/2*W0.4*215.root 
$MARSSYS/datacubes --config=./datacubes.rc 

cp $MARSSYS/mrcfiles/skymapLRT.rc .
root -q -b '$MARSSYS/macros/skymapLRT.C+("skymapLRT.rc", "Output_datacubes.root", 0)' 

Getting MAGIC data

Find out runs to use for different sources/conditions in DB

Find the path here: data magic

Recommended MC: Recommended_Stereo_MC_productions

Example: Narrow-Line Seyfert 1 galaxy (NLSy1) 1H0323+342

Subdirectory: http://data.magic.pic.es/Data/Star/v1/1H0323+342-T/2013_09_03/
wget --user=MAGIC --password=XXX http://data.magic.pic.es/Data/Star/v1/1H0323+342-T/2013_09_03/star05029007.001_M2.root

Direct copy from web:

wget -r --no-parent --user=MAGIC --password=XXX -A.root http://data.magic.pic.es/Data/Star/v1/1H0323+342-T/2013_09_03/

Copy using lcg-tools:

ssh -Y mic.magic.pic.es
voms-proxy-init --voms magic
lcg-cp --verbose -n 5 "srm://srmcms.pic.es:8443/srm/managerv2?SFN=/pnfs/pic.es/data/magic/Data/Star/v1/1H0323+342-T/2013_09_03/star05029009.001_M2.root" "file:star05029009.001_M2.root"

Jobs at PIC

glite-wms-job-submit -a -o jobid.list superstar.jdl
glite-wms-job-status -i jobid.list 
glite-wms-job-output --dir . -i jobid.list

Parameters

// ----------
// fLength   [mm]       major axis of ellipse
// fWidth    [mm]       minor axis
// fDelta    [rad]      angle of major axis with x-axis
//                      by definition the major axis is pointing into
//                      the hemisphere x>0, thus -pi/2 < delta < pi/2
// fSize     [#CerPhot] total sum of pixels
// fMeanX    [mm]       x of center of ellipse
// fMeanY    [mm]       y of center of ellipse

// fNumCorePixels  number of pixels called core
// fNumUsedPixels  number of pixels which survived the cleaning

// fNumCorePixels  moved to MNewImagePar
// fNumUsedPixels  moved to MNewImagePar

Other

Monica Vazquez Acosta's IAC website

Monica Vazquez Acosta's Cern website

-- MonicaVazquezAcosta - 10 Aug 2006

Edit | Attach | Watch | Print version | History: r27 < r26 < r25 < r24 < r23 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r27 - 2015-02-26 - MonicaVazquezAcosta
 
    • 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