Difference: CSCOfflineAnalysisTools (16 vs. 17)

Revision 172018-07-09 - StefanRaimundMaschek

Line: 1 to 1
 
META TOPICPARENT name="KalliopiIordanidou"
%CERTIFY%
Line: 334 to 334
 The package CscSegmentPerformance contains an executable csc_segment_performance.exe and csc_4d_segment_performance.exe that read a segment prediction tree and a segment measurement tree and, event by event, matches these and write the output to the cluster performance tree csc_segperf.root and csc_4d_segperf.root.

Creating Ntuples

Added:
>
>

21.0.73

fork athena and run the scipt in the appendix.

 kinit -5 ${USER}@CERN.CH
 aklog -c CERN.CH -noprdb

#~ cd `mktemp -d` # test in new dir



CSCDir=${PWD}
#setting up git and get kerberos ticket
 lsetup git
 mkdir build run
# fork athena before doing the next step
# clone the athena repository locally


#~ # using a full checkout
#~ git clone https://:@gitlab.cern.ch:8443/${USER}/athena.git
# or using the sparse checkout

 git atlas init-workdir https://:@gitlab.cern.ch:8443/${USER}/athena.git  ${CSCDir}/athena
 cd ${CSCDir}/athena
# Checkout packages that might need to be changed
 git atlas addpkg CscSegmentMakers
 git atlas addpkg  MuonRIO_OnTrack
 git atlas addpkg  MuonSegment
 git atlas addpkg  MuonRecHelperTools
 git atlas addpkg  MuonCalibITools
 git atlas addpkg  CscClusterization
#~  git atlas addpkg MuonPrepRawData #proposed by Fede

 git fetch upstream > ${CSCDir}/AvailableRelease.txt
 git checkout -b CSCntuple release/21.0.73 --no-track
# Clone the CSCntuple package. In the end this should be a athena package
cd  ${CSCDir}/
 git clone https://gitlab.cern.ch/smaschek/CSCntuple
mv CSCntuple/ ${CSCDir}/athena/
#~ rm -fr ${CSCDir}/athena/CSCntuple/.git/


# set enviroment and compile
cd ${CSCDir}/build/
 asetup 21.0.73,Athena,here
 cmake ${CSCDir}/athena/Projects/WorkDir/
make -j 64
 source ${CSCDir}/build/x86_64-slc6-gcc62-opt/setup.sh

# download a test sample from rucio and test
mkdir ${CSCDir}/data/
cd ${CSCDir}/data/
lsetup rucio
voms-proxy-init voms atlas
rucio get --nrandom 1 data17_13TeV:data17_13TeV.00339957.physics_Main.merge.DESDM_MCP.r10258_p3412
cd ${CSCDir}/run
testOption=${CSCDir}/run/runJobs.test.py
cp ../source/CSCntuple/run/runJobs.CSCspoiledCluster.grid.py ${testOption}
testfile=${CSCDir}/data/*root*
sed -i -e "s:InputFile\ \= \".*\":InputFile = \"${testfile}\":g"  -e "s:NumberOfEvents\ \= .*$:NumberOfEvents\ \= 1000:g" ${testOption}

#now you can run locally by using
athena ${testOption}

# or on grid by using
pathena ${testOption} --inDS data17_13TeV:data17_13TeV.00339957.physics_Main.merge.DESDM_MCP.r10258_p3412 --outDS user.${USER}.data.00339957.CSCntuple.1 --supStream GLOBAL --extOutFile ntuple.root --excludeFile='*.git','*.svn','*.pyc'

 

21.0.37

testdir=<some dir>
Line: 419 to 492
 
---+ Variables in the Ntuples | *csc_calib* | stored calibration values from DB | | *Variable* | *Meaning* | | zsec | z-sector: -1 for z<0, +1 for z>0 | | istation | 1=CSS, 2=CSL | | phisec | Phi-sector: 1-8 | | sector | zsec*(2*phisec-istation+1) | | wlay | Wire layer 1, 2, 3, 4 increasing with abs(z) | | measphi | 0 for eta-measurement, 1 for phi-measurement | | maxstrip | Number of strips in the layer, 48 or 192 | | nchan | number of channels per layer per measphi | | istrip[200] | strip Id (1-192 or 1-48) | | ped[200] |pedestal from DB in unit of charge | | sigma[200]| sigma from DB in unit of charge | | rms[200] | rms from DB in unit of charge| | f001[200] | f001 from DB in unit of charge (NOTE: it is not pedestal subtracted!!)| | pedADC[200] | same but in unit of ADC| |sigmaADC[200]| " | |rmsADC[200]| " | |f001ADC[200]| " | |thr[200]| threshold applied based on noise option and max(T0, k*noise)| |stat[200]| status bit from DB (0th: hot channel, 1st: dead channel, 2nd: stuck bit channel) | |address[200]| equivalent to the word from Bytestream data|  | *csc_strip* || | *Variable* | *Meaning* | | run | Run number | | evt | Event number | | bcid | Bunch Crossind id | | zsec | z-sector: -1 for z<0, +1 for z>0 | | istation | 1=CSS, 2=CSL | | phisec | Phi-sector: 1-8 | | sector | zsec*(2*phisec-istation+1) | | wlay | Wire layer 1, 2, 3, 4 increasing with abs(z) | | measphi | 0 for eta-measurement, 1 for phi-measurement | | maxstrip | Number of strips in the layer, 48 or 192 | | nstrip | Number of strips with charge above threshold | | istrip[200] | Indices for these strips | | qstrip[200] | Charges on these strips | | phase[200] | Sampling phase, 0 or 1 | | astrip[200] | true if calculated charge by parabola interpolation is above threshold Max(T0,k0*f001) | | tstrip[200] | Time at the maximum on these strips | | sfit[200] | strip fit return code | | stat[200] | status bit from DB (0th: hot channel, 1st: dead channel, 2nd: stuck bit channel) | | samples[200][4] | four sample charges in the unit of electrons (e-) |  | *csc_cluster* || | *Variable* | *Meaning* | | run | Run number | | evt | Event number | | lumiblock | Luminosity block | | bcid | Bunch Crossing id | | nentry | Number of entries (clusters) for this event | | phase[MAXENT] | Sampling phase per Segment | | pos[MAXENT] | Cluster position in local coordinates (mm) | | error[MAXENT] | Error assigned to the cluster position (mm) | | x[MAXENT] | Cluster global x-position from MuonGeoModel | | y[MAXENT] | Cluster global y-position from MuonGeoModel | | z[MAXENT] | Cluster global z-position from MuonGeoModel | | nstrip[MAXENT] | number of strips in the cluster | | qstrip[MAXENT][MAXSTRIP] | Charge on each strip in the cluster | | tstrip[MAXENT][MAXSTRIP] | Time on each strip in the cluster | | qsum[MAXENT] | Sum of three strip charges | | qpeak[MAXENT] | Peak strip charge | | speak[MAXENT] | Peak strip status | | qleft[MAXENT] | Charge on the strip left of the peak | | sleft[MAXENT] | Status of the strip left of the peak strip | | qright[MAXENT] | Charge on the strip right of the peak strip | | sright[MAXENT] | Status of the strip right of the peak strip | | tpeak[MAXENT] | Cluster time (peak strip) | | phase[MAXENT] | Sampling phase | | sfit[MAXENT] | Cluster status | | posrefit[MAXENT][MAXPOS] | Refit position (mm) | | errrefit[MAXENT][MAXPOS] | Refit position error (mm) | | srefit[MAXENT][MAXPOS] | Cluster refit status | | measphi[MAXENT] | 0 for eta-measurement, 1 for phi-measurement | | wlay[MAXENT] | Wire layer 1, 2, 3, 4 increasing with abs(z) | | istation[MAXENT] | 1=CSS, 2=CSL | | zsec[MAXENT] | z-sector: -1 for z<0, +1 for z>0 | | phisec[MAXENT] | Phi-sector: 1-8 | | sector[MAXENT] | zsec*(2*phisec-istation+1) | | pitch[MAXENT] | Strip pitch | | maxstrip[MAXENT] | Number of strips in this cathode plane | | strip0[MAXENT] | Channel number for the first strip in the cluster | | pstrip[MAXENT] | Offset of the peak strip from the first | | data[MAXDATA][]MAXENT] | Extra data from cluster fit |  | *csc_segment* || | *Variable* | *Meaning* | | run | Run number | | evt | Event number | | lumiblock | Luminosity block (not added yet) | | bcid | Bunch Crossing id | | nchamber | Number of chambers with segments in this event | | nseg [MAXCHAMBER] | Number of segments in each chamber | | nentry | Number of entries (segments) in this event | | iseg[MAXENT] | index of segment in each chamber | | nspoil[MAXENT]| # spoiled measurements on segment. | | nunspoil[MAXENT]| # not spoiled measurements on segment. | | x[MAXENT]| Segment global coordinate x | | y[MAXENT]| Segment global coordinate y | | z[MAXENT]| Segment global coordinate z | | s0[MAXENT]| Segment position | | s1[MAXENT]| Segment slope dy/dx | | sa[MAXENT]| Segment angle atan(dx/dy) degree (normal to chamber 0 deg)| | d0[MAXENT]| Position error. | | d1[MAXENT]| Slope error. | | da[MAXENT]| Angle error. | | d01[MAXENT]| Position-slope error correlation.| | d0a[MAXENT]| Position-angle error correlation.| | chsq[MAXENT]| Segment chi-square. | | phase[MAXENT]| sampling phase (in principle dependent on chamber (ROD)) | | qsum[MAXENT][4]| sum of qleft+qpeak+qright associated 2dsegment | | qpeak[MAXENT][4]| qpeak associated 2dsegment | | qleft[MAXENT][4]| qleft associated 2dsegment | | qright[MAXENT][4]| qright associated 2dsegment | | pos[MAXENT][4] | hit position associated 2dsegment | | dpos[MAXENT][4]| error of hit position associated 2dsegment | | sfit[MAXENT][4]| sfit (cluster quality) associated 2dsegment | | str0[MAXENT][4]| First strip id of cluster (strip#) | | nstr[MAXENT][4]| Cluster size (strip#) | | pstr[MAXENT][4]| Peak strip id (strip#) | | tstr[MAXENT][4]| Time of peak strip (ns) | | res[MAXENT][4]| Residual. | | dres[MAXENT][4]| Residual error. | | posAsSid[MAXENT][4]| position in unit of stripId (1-192/48) converted from position (mm) | | refAsSid[MAXENT][4]| position in unit of stripId (1-192/48) converted from position (mm) | | posInStr[MAXENT][4]| position in unit of strip UNIT (-0.5 - 0.5 in case of precision fit)| | refInStr[MAXENT][4]| Reference position from the other three hits in unit of strip UNIT | | - | (-0.5 - 0.5 in case of precision fit) | | rs2[MAXENT]| residual from the first 3 hits | | rs3[MAXENT]| residual from the second 3 hits | | drs2[MAXENT]| error of rs2 | | drs3[MAXENT]| error of rs3 | | measphi[MAXENT]|0 for r-measurement, 1 for phi-measurement | | sector[MAXENT]| zsec*(2*phisec-istation+1) | | hitmap[MAXENT]| hit3*2^3 + hit2*2^2 + hit1*2^1 + hit0*2^0 | | spoilmap[MAXENT]| hit3*2^3 + hit2*2^2 + hit1*2^1 + hit0*2^0 | | istation[MAXENT]| 1=CSS, 2=CSL | | zsec[MAXENT]| Z-sector: -1 for z<0, +1 for z>0 | | phisec[MAXENT]| Phi-sector: 1-8 |  | *csc_4d_segment* || | *Variable* | *Meaning* | | run| run# | | evt| event # | | bcid| beam crossing id | | nentry| # segments for this event | | nchamber| # chambers with segments for this event | | nseg[MAXCHAMBER]| # segments in the chamber | | gx[MAXENT]| Segment global coordinate x (mm) | | gy[MAXENT]| Segment global coordinate y (mm) | | gz[MAXENT]| Segment global coordinate z (mm) | | eta[MAXENT]| Segment global pseudorapidity | | theta[MAXENT]| Segment global azimuthal angle degree | | phi[MAXENT]| Segment global polar angle degree | | y[MAXENT]|Segment local position phi (mm) | | z[MAXENT]| Segment local position eta (mm) | | ay[MAXENT]|Segment angle degree (normal to chamber 0 deg) | | az[MAXENT]| | | dy[MAXENT]| Position error. | | dz[MAXENT]| | | day[MAXENT]| Angle error. degree | | daz[MAXENT]| | | eyz[MAXENT]| Error correlations. | | eyay[MAXENT]| | | eyaz[MAXENT]| | | ezay[MAXENT]| | | ezaz[MAXENT]| | | eayaz[MAXENT]| | | chsq[MAXENT]| Segment chi-square. | | nhit[MAXENT]| Number of hits either 6,7,8 | | ihit[MAXENT][8]|index of hit 0-7 | | str0[MAXENT][8]|First strip id of cluster (strip#) | | nstr[MAXENT][8]|Cluster size (strip#) | | pstr[MAXENT][8]|Peak strip id (strip#) | | sfit[MAXENT][8]|Cluster Status tag (strip#) | | qsum[MAXENT][8]| Sum over strip charges | | qleft[MAXENT][8]| Qleft | | qright[MAXENT][8]| Qright | | qpeak[MAXENT][8]|Qpeak | | tpeak[MAXENT][8]|Qpeak | | dqleft[MAXENT][8]| Qleft | | dqright[MAXENT][8]| Qright | | dqpeak[MAXENT][8]| Qpeak | | pos[MAXENT][8]| Strip position (mm) | | dpos[MAXENT][8]|Position error (mm) | | res[MAXENT][8]| Residual from 4 hits (mm)| | dres[MAXENT][8]|Residual error (mm) | | rs[MAXENT][8]| Residual from 3 hits (mm) | | drs[MAXENT][8]| Residual error (mm) | | posAsSid[MAXENT][8]| position in unit of stripId (1-192/48) converted from position (mm) | | refAsSid[MAXENT][8]| position in unit of stripId (1-192/48) converted from position (mm) | | posInStr[MAXENT][8]| position in unit of strip UNIT (-0.5 - 0.5 in case of precision fit)| | refInStr[MAXENT][8]| Reference position from the other three hits | | |-|in unit of strip UNIT (-0.5 - 0.5 in case of precision fit) | | istation[MAXENT]|1=CSS, 2=CSL | | zsec[MAXENT]| Z-sector: -1 for z<0, +1 for z>0 | | phisec[MAXENT]| Phi-sector: 1-8 | | sector[MAXENT]| zsec*(2*phisec-istation+1) | | hitmap[MAXENT]| hit3*2^3 + hit2*2^2 + hit1*2^1 + hit0*2^0 (X) | |-| hit7*2^7 + hit6*2^6 + hit5*2^5 + hit4*2^4 | | phase[MAXENT]| sampling phase (in principle dependent on chamber (ROD)) | ==================================================================== ----- *Major updates*:
<!--Future editors should add their signatures beneath yours!-->
-- Main.dadams - 29 Jan 2007 -- Main.StevenGoldfarb - 03 Jul 2007 (reviewed & certified) --Kalliopi Iordanidou - 6 Feb 2013
<!--Please add the name of someone who is responsible for this page so that he/she can be contacted if changes are needed.
The creator's name will be added by default, but this can be replaced if appropriate.-->
%RESPONSIBLE% Main.dadams
<!--Once this page has been reviewed, please add the name and the date e.g. Main.StephenHaywood - 31 Oct 2006 -->
%REVIEW% Main.StevenGoldfarb - 03 Jul 2007
Added:
>
>
  • setupCSC.sh: a setup script to get the CSCntuple package

META FILEATTACHMENT attachment="setupCSC.sh" attr="" comment="a setup script to get the CSCntuple package" date="1531168711" name="setupCSC.sh" path="setupCSC.sh" size="2444" user="smaschek" version="1"
 
META TOPICMOVED by="kalliopi" date="1360241407" from="Main.OfflineAnalysisTools" to="Main.CSCOfflineAnalysisTools"
 
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