Monica Vazquez Acosta's working wiki
MAGIC
.tschrcsetenv 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 ./makeVersion recommended by J. Sitarek
wget http://root.cern.ch/download/root_v5.34.21.source.tar.gz ./configure --enable-builtin-zlib ./makeVersion 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
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.txtcoach
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: noOnline 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: 1Running 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
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 Topic revision: r27 - 2015-02-26 - MonicaVazquezAcosta
Home
Index
Search
Changes
Notifications
RSS Feed
Documentation
SupportWebs
Main
Main Archive
Plugins
Sandbox for tests
- ABATBEA
- ACPP
- ADCgroup
- AEGIS
- AfricaMap
- AgileInfrastructure
- ALICE
- AliceEbyE
- AliceSPD
- AliceSSD
- AliceTOF
- AliFemto
- ALPHA
- Altair
- ArdaGrid
- ASACUSA
- AthenaFCalTBAna
- Atlas
- AtlasLBNL
- AXIALPET
- CAE
- CALICE
- CDS
- CENF
- CERNSearch
- CLIC
- Cloud
- CloudServices
- CMS
- Controls
- CTA
- CvmFS
- DB
- DefaultWeb
- DESgroup
- DPHEP
- DM-LHC
- DSSGroup
- EGEE
- EgeePtf
- ELFms
- EMI
- ETICS
- FIOgroup
- FlukaTeam
- Frontier
- Gaudi
- GeneratorServices
- GuidesInfo
- HardwareLabs
- HCC
- HEPIX
- ILCBDSColl
- ILCTPC
- IMWG
- Inspire
- IPv6
- IT
- ItCommTeam
- ITCoord
- ITdeptTechForum
- ITDRP
- ITGT
- ITSDC
- LAr
- LCG
- LCGAAWorkbook
- Leade
- LHCAccess
- LHCAtHome
- LHCb
- LHCgas
- LHCONE
- LHCOPN
- LinuxSupport
- Main
- Medipix
- Messaging
- MPGD
- NA49
- NA61
- NA62
- NTOF
- Openlab
- PDBService
- Persistency
- PESgroup
- Plugins
- PSAccess
- PSBUpgrade
- R2Eproject
- RCTF
- RD42
- RFCond12
- RFLowLevel
- ROXIE
- Sandbox
- SocialActivities
- SPI
- SRMDev
- SSM
- Student
- SuperComputing
- Support
- SwfCatalogue
- TMVA
- TOTEM
- TWiki
- UNOSAT
- Virtualization
- VOBox
- WITCH
- XTCA
Welcome Guest Login or Register
 |
|
Copyright &© 2008-2024 by the contributing authors. All material on this collaboration platform is the property of the contributing authors. or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback