Skip to content

Instantly share code, notes, and snippets.

@jdbrice
Last active March 6, 2017 20:27
Show Gist options
  • Save jdbrice/3e764e19684583fab82df4684d88bf57 to your computer and use it in GitHub Desktop.
Save jdbrice/3e764e19684583fab82df4684d88bf57 to your computer and use it in GitHub Desktop.
Accessing VPD Trigger Data from MuDst
#include "StEvent/StTriggerData2016.h"
// "event" is a class object storing the data into an nTuple
void StVpdPicoDstMaker::analyzeEvent(){
StMuEvent *muEvent = muDst->event();
const StTriggerData2016 *td = (StTriggerData2016*)muEvent->triggerData();
event.evtNumber = td->eventNumber();
event.npre = td->numberOfPreXing();
event.npost = td->numberOfPostXing();
event.nsum = event.npre + event.npost + 1;
int iPrePost = 0;
if ( event.npre != 0 || event.npost != 0){
cout << "Pre = " << event.npre << endl << endm;
cout << "Post = " << event.npost << endl << endm;
return;
}
for ( int iVPD = 0; iVPD < 16; iVPD ++ ){
event.bbqAdcEast[ iVPD ] = td->vpdADC( east, iVPD+1, iPrePost );
event.bbqAdcWest[ iVPD ] = td->vpdADC( west, iVPD+1, iPrePost );
event.bbqTdcEast[ iVPD ] = td->vpdTDC( east, iVPD+1, iPrePost );
event.bbqTdcWest[ iVPD ] = td->vpdTDC( west, iVPD+1, iPrePost );
// MXQ crate (for MTD)
event.mxqAdcEast[ iVPD ] = td->vpdADCHighThr( east, iVPD+1, iPrePost );
event.mxqAdcWest[ iVPD ] = td->vpdADCHighThr( west, iVPD+1, iPrePost );
event.mxqTdcEast[ iVPD ] = td->vpdTDCHighThr( east, iVPD+1, iPrePost );
event.mxqTdcWest[ iVPD ] = td->vpdTDCHighThr( west, iVPD+1, iPrePost );
}
event.bbqEarliestEast = td->vpdEarliestTDC( east, iPrePost );
event.bbqEarliestWest = td->vpdEarliestTDC( west, iPrePost );
event.grefmult = muEvent->grefmult();
event.refmult = muEvent->refMult();
// Run16 ONLY
event.bbcEastVpdNHits = (td->bbcVP101(4, iPrePost)>>12);
event.bbcEastVpdADCsum = (td->bbcVP101(4, iPrePost)&0xfff);
event.bbcEastVpdTACsum = td->bbcVP101(5, iPrePost);
event.bbcWestVpdNHits = (td->bbcVP101(6, iPrePost)>>12);
event.bbcWestVpdADCsum = (td->bbcVP101(6, iPrePost)&0xfff);
event.bbcWestVpdTACsum = td->bbcVP101(7, iPrePost);
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment