- dolv719's home page
- Posts
- 2017
- 2016
- December (1)
- November (3)
- October (2)
- September (5)
- August (3)
- July (1)
- June (3)
- April (2)
- March (2)
- February (2)
- January (5)
- 2015
- December (1)
- November (1)
- October (2)
- September (2)
- August (1)
- July (1)
- June (1)
- May (1)
- April (1)
- March (3)
- February (1)
- 2014
- 2013
- My blog
- Post new blog entry
- All blogs
ADC Plots
The program generates plots by running over a daq file for single run. A tree is created and stores info such as maxadc, strip info, disk, layer, etc.
There are no cuts or any requirements, they idea is to accept as much data as possible.
All the information is gathered from the StFgtStrip class, using StGeomDecode function in the .ccx. This all follows the formula of StFgtClusterPlotter code.
Creating the plots is very versatile, I am able to plot data over all 24 quadrants for both the R and Phi layer, or for a specific quadrant and layer.
The plots below are for all of disk 1 for R and Phi layer.
Questions
1. Why is there a gap in the R plot from ~285-400, but not in the Phi plot?
->Answer by Xuan, see link: http://drupal.star.bnl.gov/STAR/blog/balewski/2012/jan/28/master-mapping-fgt-2012
2. Should the maxadc values be higher, in the jplots adc hits ~3000?
The full code can be found at
/star/u/dolv719/FgtOffline/StRoot/StFgtPool/StFgtClusterTools/StFgtClusterSaver.cxx
The "main" loop is below
for(int iDx=0;iDx<kFgtNumDiscs;iDx++)
{
StFgtStripCollection& strips=*(mFgtCollectionPtr->getStripCollection( iDx ));
fgtcol=mFgtCollectionPtr;
for( StSPtrVecFgtStripIterator it=strips.getStripVec().begin();it!=strips.getStripVec().end(); ++it)
{
Short_t quad, disc, strip;
Char_t layer;
Int_t layernum;
StFgtGeom::decodeGeoId((*it)->getGeoId(),disc, quad, layer, strip);
Double_t maxadc = (*it)->getMaxAdc();
Double_t geoid = (*it)->getGeoId();
if(layer=='R')
layernum=0;
else layernum=1;
StFgtHitTrivia* hit = new StFgtHitTrivia();
hit->setDisk(disc);
hit->setStrip(strip);
hit->setMaxAdc(maxadc);
hit->setGeoId(geoid);
hit->setLayer(layernum);
hit->setQuad(quad);
vec->push_back(hit);
}
- dolv719's blog
- Login or register to post comments