arich modules

Collaboration diagram for arich modules:

Classes
class	arichBtestModule
	The UserTutorial module. More...
class	ARICHChannelMaskModule
	Testing module for collection of calibration data. More...
class	ARICHDigitizerModule
	ARICH digitizer module. More...
class	ARICHDQMModule
	Make summary of data quality from reconstruction. More...
class	ARICHFillHitsModule
	Fill ARICHHit collection from ARICHDigits. More...
class	ARICHMCParticlesModule
	Module to match ARICH hits to MCParticles. More...
class	ARICHNtupleModule
	ARICH reconstruction efficiency test module. More...
class	ARICHPackerModule
	Raw data packer for ARICH. More...
class	ARICHRateCalModule
	Fill ARICHHit collection from ARICHDigits. More...
class	ARICHRawUnpackerModule
	Fill ARICHHit collection from ARICHDigits. More...
class	ARICHReconstruction
	Internal ARICH track reconstruction. More...
class	ARICHReconstructorModule
	ARICH subdetector main module. More...
class	ARICHRelateModule
	Creates relations between ARICHAeroHits and ExtHits. More...
class	arichToNtupleModule
	This module extends existing variablesToNtuple to append detailed arich information to candidates in the analysis output ntuple. More...
struct	ARICHRawHeader
class	ARICHUnpackerModule
	Raw data unpacker for ARICH. More...

Macros
#define	ARICH_BUFFER_NWORDS   252
#define	ARICHFEB_HEADER_SIZE   10
#define	ARICHRAW_HEADER_SIZE   12

Functions
	REG_MODULE (arichBtest)
	Register the Module.
	REG_MODULE (ARICHDigitizer)
	Register the Module.
TH2 *	moduleHitMap (TH1 *hitMap, int moduleID)
	Make hit map in HAPD view (12*12 channels) More...
TH2 *	moduleDeadMap (TH1 *hitMap, int moduleID)
	Make chip dead/alive map in HAPD view (2*2 chips) More...
TH1 *	mergerClusterHitMap1D (TH1 *hitMap, int mergerID)
	Make 1D hit map of specified Merger Board. More...
TCanvas *	mergerClusterHitMap2D (TH1 *hitMap, int mergerID)
	Make display of 6 HAPDs' 2D hit map of the Merger Board. More...
TCanvas *	sectorHitMap (TH1 *hitMap, int sector)
	Make display of 70 HAPDs' 2D hit map of the sector. More...
TCanvas *	sectorDeadMap (TH1 *hitMap, int sector)
	Make display of 70 HAPDs' 2D dead/alive map of the sector. More...
	REG_MODULE (ARICHDQM)
void	deadPalette ()
	REG_MODULE (ARICHFillHits)
	Register module.
	REG_MODULE (ARICHMCParticles)
	Register module.
	REG_MODULE (ARICHNtuple)
	Register module.
	REG_MODULE (ARICHPacker)
	Register module.
	REG_MODULE (ARICHRateCal)
	REG_MODULE (ARICHRawUnpacker)
	REG_MODULE (ARICHReconstructor)
	Register the Module.
	REG_MODULE (ARICHRelate)
	Register module.
	REG_MODULE (ARICHUnpacker)
	Register module.
	arichBtestModule ()
	Constructor. More...
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
int	skipdata (gzFile fp)
	Skip the data part of the record.
void	readmwpc (unsigned int *dbuf, unsigned int len, int print=0)
	Read the MWPC information from the data buffer.
int	readhapd (unsigned int len, unsigned int *data)
	Read the HAPD hits from the data buffer.
int	getTrack (int mask, TVector3 &r, TVector3 &dir)
	Beamtest Track reconstruction.
int	readdata (gzFile fp, int rec_id, int print)
	Read the data from the file (can be compressed)
virtual void	event () override
	Running over all events. More...
virtual void	endRun () override
	Is called after processing the last event of a run. More...
virtual void	terminate () override
	Is called at the end of your Module. More...
	ARICHDigitizerModule ()
	Constructor.
virtual	~ARICHDigitizerModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor. More...
void	magFieldDistorsion (TVector2 &hit, int copyno)
	Apply correction to hit position due to non-perpendicular component of magnetic field. More...
	ARICHDQMModule ()
	Constructor.
virtual	~ARICHDQMModule ()
	Destructor.
virtual void	defineHisto () override
	Function to define histograms. More...
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
	ARICHFillHitsModule ()
	Constructor.
virtual	~ARICHFillHitsModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
void	magFieldCorrection (TVector3 &hitpos)
	Corrects hit position for distorsion due to non-perpendicular magnetic field component.
	ARICHMCParticlesModule ()
	Constructor.
virtual	~ARICHMCParticlesModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
	ARICHNtupleModule ()
	Constructor.
virtual	~ARICHNtupleModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
virtual void	terminate () override
	Termination action. More...
	ARICHPackerModule ()
	Constructor.
virtual	~ARICHPackerModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
void	writeHeader (int *buffer, unsigned &ibyte, const ARICHRawHeader &head)
	TODO!
unsigned int	calbyte (const int *buf)
unsigned int	calword (const int *buf)
	ARICHRateCalModule ()
	Constructor.
virtual	~ARICHRateCalModule ()
	Destructor.
virtual void	defineHisto () override
	Function to define histograms. More...
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
unsigned int	calbyte (const int *buf)
unsigned int	calword (const int *buf)
	ARICHRawUnpackerModule ()
	Constructor.
virtual	~ARICHRawUnpackerModule ()
	Destructor.
virtual void	defineHisto () override
	Function to define histograms. More...
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
	ARICHReconstruction (int storePhotons=0)
	Constructor.
void	initialize ()
	read geomerty parameters from xml and initialize class memebers
int	InsideDetector (TVector3 a, int copyno)
	Returns 1 if vector "a" lies on "copyno"-th detector active surface of detector and 0 else.
int	smearTrack (ARICHTrack &arichTrack)
	Smeares track parameters ("simulate" the uncertainties of tracking).
TVector3	FastTracking (TVector3 dirf, TVector3 r, double *refind, double *z, int n, int opt)
	Calculates the intersection of the Cherenkov photon emitted from point "r" in "dirf" direction with the detector plane. More...
TVector3	HitVirtualPosition (const TVector3 &hitpos, int mirrorID)
	Returns the hit virtual position, assuming that it was reflected from mirror. More...
bool	HitsMirror (const TVector3 &pos, const TVector3 &dir, int mirrorID)
	returns true if photon at position pos with direction dir hits mirror plate with ID mirrorID
int	CherenkovPhoton (TVector3 r, TVector3 rh, TVector3 &rf, TVector3 &dirf, double *refind, double *z, int n, int mirrorID=0)
	Calculates the direction of photon emission. More...
int	likelihood2 (ARICHTrack &arichTrack, const StoreArray< ARICHHit > &arichHits, ARICHLikelihood &arichLikelihood)
	Computes the value of identity likelihood function for different particle hypotheses.
void	setTrackPositionResolution (double pRes)
	Sets track position resolution (from tracking)
void	setTrackAngleResolution (double aRes)
	Sets track direction resolution (from tracking)
TVector3	getTrackMeanEmissionPosition (const ARICHTrack &track, int iAero)
	Returns mean emission position of Cherenkov photons from i-th aerogel layer.
TVector3	getTrackPositionAtZ (const ARICHTrack &track, double zout)
	Returns track direction at point with z coordinate "zout" (assumes straight track).
void	transformTrackToLocal (ARICHTrack &arichTrack, bool align)
	Transforms track parameters from global Belle2 to ARICH local frame.
TVector3	getMirrorPoint (int mirrorID)
	Returns point on the mirror plate with id mirrorID.
TVector3	getMirrorNorm (int mirrorID)
	Returns normal vector of the mirror plate with id mirrorID.
void	correctEmissionPoint (int tileID, double r)
	correct mean emission point z position
	ARICHReconstructorModule ()
	Constructor.
virtual	~ARICHReconstructorModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module.
virtual void	beginRun () override
	Called when entering a new run.
virtual void	event () override
	Event processor.
void	printModuleParams ()
	Print module parameters.
	ARICHRelateModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
unsigned int	calbyte (const int *buf)
	calculate number of bytes in raw Unpacker
unsigned int	cal2byte (const int *buf)
	calculate number of lines (2 bytes) in raw Unpacker
unsigned int	calword (const int *buf)
	calculate number of words in raw Unpacker
	ARICHUnpackerModule ()
	Constructor.
virtual	~ARICHUnpackerModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor. More...
void	readHeader (const int *buffer, unsigned &ibyte, ARICHRawHeader &head)
void	readFEHeader (const int *buffer, unsigned &ibyte, ARICHRawHeader &head)
	read Merger header
void	printBits (const int *buffer, int bufferSize)
	Unpack raw data given in production format. More...

Variables
TNtuple *	m_tuple
	ntuple containing hapd hits
TH1F *	hapd [6]
	histogram of hits for each hapd
TH1F *	mwpc_tdc [4][5]
	tdc information from mwpcs
TH1F *	mwpc_diff [4][2]
	tdc difference from mwpcs
TH1F *	mwpc_sum [4][2]
	tdc sum from mwpcs
TH1F *	mwpc_sum_cut [4][2]
	tdc sum from mwpcs, with sum cut applied
TH1F *	mwpc_residuals [4][2]
	residuals from mwpcs
TH2F *	mwpc_xy [4]
	calculated x-y track positions
TH2F *	mwpc_residualsz [4][2]
	z-residuals from mwpcs

Detailed Description

Function Documentation

arichBtestModule()

arichBtestModule

(

)

Constructor.

Sets the module parameters.

Definition at line 62 of file arichBtestModule.cc.

                                     : Module(), m_end(0), m_events(0), m_file(NULL), m_timestart(0), m_mwpc(NULL)
  {
    //Set module properties
    setDescription("Module for the ARICH Beamtest data analysis. It creates track form the MWPC hits and reads the HAPD hits");
 
    //Parameter definition
    addParam("outputFileName", m_outFile, "Output Root Filename", std::string("output.root"));
    std::vector<std::string> defaultList;
    addParam("runList", m_runList, "Data Filenames.", defaultList);
    std::vector<int> defaultMask;
    addParam("mwpcTrackMask", m_MwpcTrackMask, "Create track from MWPC layers", defaultMask);
    m_fp = NULL;
    addParam("beamMomentum", m_beamMomentum, "Momentum of the beam [GeV]", 0.0);
 
    for (int i = 0; i < 32; i++) m_tdc[i] = 0;
 
  }

beginRun() [1/4]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

At the beginning of each run, the function gives you the chance to change run dependent constants like alignment parameters, etc.

Reimplemented from Module.

Definition at line 148 of file arichBtestModule.cc.

beginRun() [2/4]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

Set run dependent things like run header parameters, alignment, etc.

Reimplemented from Module.

Definition at line 89 of file ARICHDigitizerModule.cc.

  {
    if (m_simPar->getNBkgHits() > 0)  m_bkgLevel = m_simPar->getNBkgHits();
  }

beginRun() [3/4]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

Set run dependent things like run header parameters, alignment, etc.

Reimplemented from HistoModule.

Definition at line 176 of file ARICHDQMModule.cc.

  {
 
    h_chStat->Reset();
    h_aeroStat->Reset();
    h_chDigit->Reset();
    h_chipDigit->Reset();
    h_hapdDigit->Reset();
 
    h_chHit->Reset();
    h_chipHit->Reset();
    h_hapdHit->Reset();
    h_mergerHit->Reset();
    h_bitsPerMergerNorm->Reset();
    h_bitsPerHapdMerger->Reset();
 
    h_aerogelHit->Reset();
    h_bits->Reset();
    h_bitsPerChannel->Reset();
    h_hitsPerTrack2D->Reset();
    h_tracks2D->Reset();
    h_aerogelHits3D->Reset();
 
    h_hitsPerEvent->Reset();
    h_theta->Reset();
    h_hitsPerTrack->Reset();
    h_trackPerEvent->Reset();
    h_flashPerAPD->Reset();
    h_mirrorThetaPhi->Reset();
    h_thetaPhi->Reset();
 
    for (int i = 0; i < 6; i++) {
      h_secTheta[i]->Reset();
      h_secHitsPerTrack[i]->Reset();
    }
 
    h_ARICHOccAfterInjLer->Reset();
    h_ARICHEOccAfterInjLer->Reset();
    h_ARICHOccAfterInjHer->Reset();
    h_ARICHEOccAfterInjHer->Reset();
  }

beginRun() [4/4]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

Set run dependent things like run header parameters, alignment, etc.

Reimplemented from HistoModule.

Definition at line 104 of file ARICHRateCalModule.cc.

  {
    StoreObjPtr<EventMetaData> evtmetadata;
    //int expno = evtmetadata->getExperiment();
    int runno = evtmetadata->getRun();
    if (runno == 100 && !m_internalmode) {
      terminate();
    }
  }

CherenkovPhoton()

int CherenkovPhoton ( TVector3  r, TVector3  rh, TVector3 &  rf, TVector3 &  dirf, double *  refind, double *  z, int  n, int  mirrorID = 0  )

private

Calculates the direction of photon emission.

Giving the Cherenkov photon emission point "r" and its position on detector plane "rh" (hit position) this methods calculates the direction "dirf" in which photon was emitted, under the assumption that it was reflected from "nmir"-th mirror plate (nmir=-1 for no reflection).

Parameters

r	vector of photon emission point
rh	photon hit position
dirf	vector of photon emission direction (this is output of method)
rf	vector of photon position on aerogel exit
refind	array of layers refractive indices
n	number of aerogel layers through which photon passes
z	array of z coordinates of borders between layers
nmir	id of mirror from which the foton was reflected (assuming).

Definition at line 250 of file ARICHReconstruction.cc.

  {
    //
    // Description:
    // The method calculates the direction of the cherenkov photon
    // and intersection with the aerogel exit point
    //
    // Called by: CerenkovAngle
 
 
    // Arguments:
    // Input:
    // dir, r track position
    // rh photon hit
 
 
    // Output:
    // rf aerogel exit from which the  photon was emitted
    // dirf photon direction in aerogel
    static TVector3 norm(0, 0, 1); // detector plane normal vector
 
    double dwin    = m_arichgp->getHAPDGeometry().getWinThickness();
    double refractiveInd0 = m_arichgp->getHAPDGeometry().getWinRefIndex();
 
    // iteration is stoped when the difference of photon positions on first aerogel exit
    // between two iterations is smaller than this value.
    const double dmin  = 0.0000001;
    const int    niter = 100; // maximal number of iterations
    TVector3  dirw;
    TVector3 rh1 = rh - dwin * norm;
 
    std::vector <TVector3> rf0(n + 2);
    //  rf0[0] .. track point
    //  rf0[1] 1. 1st aerogel exit
    //  rf0[n] n.  aerogel exit ...
    std::vector <TVector3> dirf0(n + 2);
    //  dirf0[0] .. 1st aerogel direction
    //  dirf0[1] .. 2nd aerogel direction
    //  dirf0[n] .. direction after aerogels
 
    //  z[0] .. 1st aerogel exit
    //  z[n-1] .. 2nd aerogel exit
    //  z[n]    .. detector position
    //  z[n+1]  .. detector + window
 
    rf0[0] = r;
    rf0[1] = rf;
 
    for (int iter = 0; iter < niter; iter++) {
 
      // direction in the space between aerogels and detector
      // *************************************
      if (iter == 0) dirf0[n] = (rh1 - rf0[1]).Unit();
      else  dirf0[n] = (rh1 - rf0[n]).Unit();
 
      // *************************************
      // n-layers of aerogel // refractiveInd relative refractive index
      for (int a = n - 1; a >= 0 ; a--) {
        double rind = refractiveInd[a] / refractiveInd[a + 1];
        dirf0[a] = Refraction(dirf0[a + 1], rind);
      }
 
      double path = (z[0] - r.Z()) / dirf0[0].Z();
      double x1 = rf0[1].X();
      double y1 = rf0[1].Y();
      for (int a = 0; a < n ; a++) {
        rf0[a + 1] = rf0[a] + dirf0[a] * path;
        path = (z[a + 1] - rf0[a + 1].Z()) / dirf0[a + 1].Z();
      }
 
      Refraction(dirf0[n], norm, refractiveInd0, dirw);
 
      // *************************************
 
      path = dwin / (dirw * norm);
      rh1 = rh - dirw * path;
 
      double d2 = (rf0[1].X() - x1) * (rf0[1].X() - x1) + (rf0[1].Y() - y1) * (rf0[1].Y() - y1);
 
      if ((d2 < dmin) && (iter)) {
 
        // if mirror hypothesis check if reflection point lies on mirror plate
        if (mirrorID) {
          if (!HitsMirror(rf0[n], dirf0[n], mirrorID)) return -1;
        }
 
        rf = rf0[1];
        dirf = dirf0[0];
        return iter;
      }
    }
    return -1;
  }

defineHisto() [1/3]

void defineHisto ( )

overridevirtual

Function to define histograms.

This function is hooked to HistoManager by calling RbTupleManager::Instance().register_module ( this ) or using a macro REG_HISTOGRAM. It is supposed to be done in initialize() function.

Reimplemented from HistoModule.

Definition at line 62 of file ARICHDQMModule.cc.

defineHisto() [2/3]

void defineHisto ( )

overridevirtual

Function to define histograms.

This function is hooked to HistoManager by calling RbTupleManager::Instance().register_module ( this ) or using a macro REG_HISTOGRAM. It is supposed to be done in initialize() function.

Reimplemented from HistoModule.

Definition at line 65 of file ARICHRateCalModule.cc.

defineHisto() [3/3]

void defineHisto ( )

overridevirtual

Function to define histograms.

This function is hooked to HistoManager by calling RbTupleManager::Instance().register_module ( this ) or using a macro REG_HISTOGRAM. It is supposed to be done in initialize() function.

Reimplemented from HistoModule.

Definition at line 46 of file ARICHRawUnpackerModule.cc.

  {
    h_rate_a_all = new TH1F("h_rate_a_all", ";Channel ID", 144 * 6, 0, 144 * 6); //yone
    h_rate_b_all = new TH1F("h_rate_b_all", ";Channel ID", 144 * 6, 0, 144 * 6); //yone
    h_rate_c_all = new TH1F("h_rate_c_all", ";Channel ID", 144 * 6, 0, 144 * 6); //yone
    h_rate_d_all = new TH1F("h_rate_d_all", ";Channel ID", 144 * 6, 0, 144 * 6); //yone
  }

endRun() [1/2]

void endRun ( void  )

overridevirtual

Is called after processing the last event of a run.

Good e.g. for storing stuff, that you want to aggregate over one run.

Reimplemented from Module.

Definition at line 554 of file arichBtestModule.cc.

endRun() [2/2]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from HistoModule.

Definition at line 398 of file ARICHDQMModule.cc.

event() [1/3]

void event ( void  )

overridevirtual

Running over all events.

Function is called for each evRunning over all events This means, this function is called very often, and good performance of the code is of strong interest.

Reimplemented from Module.

Definition at line 465 of file arichBtestModule.cc.

event() [2/3]

void event ( void  )

overridevirtual

Event processor.

Convert ARICHSimHits of the event to arichDigits.

Reimplemented from Module.

Definition at line 94 of file ARICHDigitizerModule.cc.

event() [3/3]

void event ( void  )

overridevirtual

Event processor.

temporary! FEB Slots on merger should be 1-6 (now 0-5). Remove when firmware is updated!

Reimplemented from Module.

Definition at line 90 of file ARICHUnpackerModule.cc.

  {
 
    StoreArray<RawARICH> rawData(m_inputRawDataName);
    StoreArray<ARICHDigit> digits(m_outputDigitsName);
    StoreArray<ARICHRawDigit> rawdigits(m_outputRawDigitsName);
    StoreObjPtr<ARICHInfo> arichinfo(m_outputarichinfoName);
    arichinfo.create();
    StoreObjPtr<EventMetaData> evtMetaData;
 
    digits.clear();
    bool m_pciedata = false;
    int trgtype = 16;
    double vth_thscan = 0.0;
 
    if (m_debug) {
      std::cout << std::endl << "------------------------" << std::endl;
      std::cout << "Run: " << evtMetaData->getRun() << " Event: " << evtMetaData->getEvent()  << std::endl;
      std::cout << "------------------------" << std::endl << std::endl;
    }
 
    unsigned thscan_mode = 0;
    // regular Unpacker mode, fill ARICHDigit
//    if (m_disable_unpacker == 0) {
 
    for (auto& raw : rawData) {
      // Check PCIe40 data or Copper data
      if (raw.GetMaxNumOfCh(0) == 48) { m_pciedata = true; } // Could be 36 or 48
      else if (raw.GetMaxNumOfCh(0) == 4) { m_pciedata = false; }
      else { B2FATAL("ARICHUnpackerModule: Invalid value of GetMaxNumOfCh from raw data: " << LogVar("Number of ch: ", raw.GetMaxNumOfCh(0))); }
 
      for (int finesse = 0; finesse < raw.GetMaxNumOfCh(0); finesse++) {
        const int* buffer = raw.GetDetectorBuffer(0, finesse);
        int bufferSize = raw.GetDetectorNwords(0, finesse);
 
        if (bufferSize < 1)
          continue;
 
        // record the trigger type from the B2L data
        trgtype = raw.GetTRGType(0);
 
        // read merger header
        unsigned ibyte = 0;
        ARICHRawHeader head;
 
        readHeader(buffer, ibyte, head);
 
        if (m_debug > 1) printBits(buffer, bufferSize);
 
        if (m_debug) {
          std::cout << "Merger header" << std::endl;
          head.print();
        }
        //-- RawDigit for Merger info
        int type = (int)head.type;
        int ver = (int)head.version;
        int boardid = (int)head.mergerID;
        int febno = (int)head.FEBSlot;
        unsigned int length_all = (unsigned int)head.length;
        unsigned int mrg_evtno = (unsigned int)head.trigger;
        ARICHRawDigit* rawdigit = rawdigits.appendNew(type, ver, boardid, febno, length_all, mrg_evtno);
 
        if (!m_pciedata) {
          rawdigit->setCopperId(raw.GetNodeID(0));
          rawdigit->setHslbId(finesse);
        } else {
          if (raw.GetNodeID(0) == 0x4000001) {
            rawdigit->setCopperId(raw.GetNodeID(0) + (int)(finesse / 4));
          } else if (raw.GetNodeID(0) == 0x4000002) {
            rawdigit->setCopperId(0x400000A + (int)(finesse / 4));
          } else {
            B2FATAL("ARICHUnpackerModule: Invalid Node ID from readout: " <<  LogVar("NodeID: ", raw.GetNodeID(0)));
          }
 
          rawdigit->setHslbId(finesse % 4);
          rawdigit->setPcieId(raw.GetNodeID(0));
          rawdigit->setPcieChId(finesse);
        }
 
        //-- end of RawDigit for Merger info
 
        // record the ibyte here
        unsigned begin = ibyte;
 
        while (ibyte < head.length) {
 
          // new feb
          ARICHRawHeader febHead;
          readFEHeader(buffer, ibyte, febHead);
          if (febHead.thscan_mode) {thscan_mode++;}
          if (m_debug) febHead.print();
 
          if (/*febHead.type != head.type ||*/ febHead.version != head.version || febHead.mergerID != head.mergerID
                                               || febHead.trigger != head.trigger) {
            B2ERROR("ARICHUnpackerModule: data in FEB header not consistent with data in merger HEADER " << LogVar("FEB ID",
                    (unsigned)febHead.FEBSlot) <<
                    LogVar("merger ID", (unsigned)head.mergerID)); break;
          }
 
          // feb header shift
          ibyte += ARICHFEB_HEADER_SIZE;
          int dataLen = febHead.length - ARICHFEB_HEADER_SIZE;
 
          febHead.FEBSlot += 1; 
 
          unsigned mergID = m_mergerMap->getMergerIDfromSN((unsigned)head.mergerID);
 
          if (mergID == 99) { B2ERROR("ARICHUnpackerModule: unknown merger number. Merger data will be skipped.  " << LogVar("merger ID", mergID) << LogVar("Serial Number", (unsigned)head.mergerID)); break;}
 
          unsigned moduleID = m_mergerMap->getModuleID(mergID, (unsigned)febHead.FEBSlot);
 
          if (!moduleID) { B2ERROR("ARICHUnpackerModule: no merger to FEB mapping. Merger data will be skipped. " << LogVar("merger ID", mergID) << LogVar("Serial Number", (unsigned)head.mergerID) << LogVar("FEB slot", (unsigned)febHead.FEBSlot)); break;}
 
          // read data
          if (m_debug) std::cout << "Hit channels: " << std::endl;
          if (febHead.type == 1) {
            for (int i = 0; i < dataLen / 2; i++) {
              int shift = (3 - ibyte % 4) * 8;
              uint8_t asicCh = buffer[ibyte / 4] >> shift;
              ibyte++;
              shift = (3 - ibyte % 4) * 8;
              uint8_t hitBitSet = buffer[ibyte / 4] >> shift;
              if (m_debug && hitBitSet) std::cout << "ch: " << (unsigned)asicCh << " " <<  std::bitset<8>(hitBitSet) << std::endl;
              // store digit
              digits.appendNew(moduleID, (unsigned)asicCh, hitBitSet);
              ibyte++;
            }
          } else if (febHead.type == 2) {
            unsigned asicCh = 143;
            for (int i = 0; i < dataLen; i++) {
              int shift = (3 - ibyte % 4) * 8;
              uint8_t hitBitSet = buffer[ibyte / 4] >> shift;
              // store digit if hit
              if (hitBitSet & m_bitMask) {
                digits.appendNew(moduleID, asicCh, hitBitSet);
              }
              asicCh--;
              ibyte++;
            }
          } else B2ERROR("ARICHUnpackerModule: Unknown data type" << LogVar("type", febHead.type));
 
        }
 
        if (ceil(ibyte / 4.) != (unsigned)bufferSize)
          B2WARNING("ARICHUnpackerModule: data buffer size mismatch " <<  LogVar("size from copper", bufferSize) << LogVar("size from merger",
                    ceil(
                      ibyte / 4.)));
 
 
        //-- If thscan_mode detected from header: proceed to fill ARICHRawDigit
        //-- If m_rawmode is set to 1: proceed to fill ARICHRawDigit
        if (thscan_mode == 0 && m_rawmode == 0) continue;
        //-- go back to beginning again for second loop in case of thscan
        m_ibyte = begin;
 
        while (m_ibyte < length_all) {
          ARICHRawHeader febHead;
          readFEHeader(buffer, m_ibyte, febHead);
          int type_feb = febHead.type;
          ver = febHead.version;
          boardid = febHead.mergerID;
          febno = febHead.FEBSlot;
 
          vth_thscan = (febHead.vth * 0.0024) - 1.27;
          unsigned int length = febHead.length;
          int evtno = febHead.trigger;
          unsigned int jbyte = 0;
          std::stringstream ss;
          ss << "type=" << type_feb << ", ver=" << ver << " "
             << ", boardid=" << boardid << ", febno=" << febno
             << ", length=" << length << ", evtno=" << evtno << " ";
          bool hasHit = false;
          long long feb_trigno = 0;
          for (int i = 0; i < 10; i++) {
            int val = calbyte(buffer);
            jbyte++;
            if (i > 1 && i < 6) {
              feb_trigno |= (0xff & val) << (5 - i) * 8;
            }
          }
          ARICHRawDigit::FEBDigit feb;
          if (type_feb == 0x02) {//Raw mode
            int ch = 143;
            //B2INFO("raw mode");
            while (jbyte < length) {
              int val = calbyte(buffer);
              if (val != 0) {
                jbyte++;
                ss << "ch# " << ch << "(" << val << ") ";
                hasHit = true;
                if (febno < 0 || febno > 6) {
                  B2ERROR("FEB is bad : " << LogVar("FEB no.", febno) << LogVar("hslb", finesse) << LogVar("type", type_feb) << LogVar("ver",
                          ver) << LogVar("boardid", boardid) << LogVar("febno", febno) << LogVar("length", length) << LogVar("evtno", evtno));
                }
                feb.push_back(ch, val);
              }
              ch--;
              if (ch < 0) break;
            }
          } else if (type_feb == 0x01) { // Suppressed mode
            // The below line is commented since it sometimes causes problem during processing threshold scan data.
            // No harm to comment this line since it is only utilized for threshold scan data.
            //if (length > 144 * 2 + 10) B2FATAL("error " << LogVar("length", length));
            //B2INFO("suppreed mode");
            while (jbyte < length) {
              int ch = calbyte(buffer);
              jbyte++;
              int val = calbyte(buffer);
              jbyte++;
              if (val != 0) {
                ss << "ch# " << ch << "(" << val << ") ";
                hasHit = true;
                if (febno < 0 || febno > 6) {
                  B2ERROR("FEB is bad : " << LogVar("FEB no.", febno) << LogVar("hslb", finesse) << LogVar("type", type_feb) << LogVar("ver",
                          ver) << LogVar("boardid", boardid) << LogVar("febno", febno) << LogVar("length", length) << LogVar("evtno", evtno));
                  return;
                }
                feb.push_back(ch, val);
              }
            }
          }
          rawdigit->addFEB(feb, type, ver, boardid, febno, length, evtno, feb_trigno);
          if (m_debug && hasHit) {
            B2INFO(ss.str());
          }
        }
 
 
 
      }
    } // end of rawData loop
 
//    } // end of regular unpacker
    /*
        // RawUnpacker mode, fill ARICHRawDigit
        if (m_rawmode == 1) {
          for (auto& raw : rawData) {
            for (int finesse = 0; finesse < 4; finesse++) {
              const int* buf = (const int*)raw.GetDetectorBuffer(0, finesse);
              int bufSize = raw.GetDetectorNwords(0, finesse);
              if (bufSize < 1) continue;
              m_ibyte = 0;
              // read merger header
              int type = calbyte(buf);
              int ver = calbyte(buf);
              int boardid = calbyte(buf);
              int febno = calbyte(buf);
              unsigned int length_all = calword(buf);
              unsigned int mrg_evtno = calword(buf);
              ARICHRawDigit* rawdigit = rawdigits.appendNew(type, ver, boardid, febno, length_all, mrg_evtno);
              rawdigit->setCopperId(raw.GetNodeID(0));
              rawdigit->setHslbId(finesse);
              int nfebs = 0;
    //--done
              while (m_ibyte < length_all) {
                int type_feb = calbyte(buf);
                ver = calbyte(buf);
                boardid = calbyte(buf);
                febno = calbyte(buf);
 
                // first line: vth value
                unsigned int vth_int = cal2byte(buf);
                if (vth_int > 0) { vth_thscan = (vth_int * 0.0024) - 1.27; }
                // second line: length
                unsigned int length = cal2byte(buf);
                int evtno = calword(buf);
                unsigned int ibyte = 0;
                std::stringstream ss;
                ss << "type=" << type_feb << ", ver=" << ver << " "
                   << ", boardid=" << boardid << ", febno=" << febno
                   << ", length=" << length << ", evtno=" << evtno << " ";
                bool hasHit = false;
                long long feb_trigno = 0;
                for (int i = 0; i < 10; i++) {
                  int val = calbyte(buf);
                  ibyte++;
                  if (i < 6) {
                    feb_trigno |= (0xff & val) << (5 - i) * 8;
                  }
                }
                ARICHRawDigit::FEBDigit feb;
                nfebs++;
                if (type_feb == 0x02) {//Raw mode
                  int ch = 143;
                  //B2INFO("raw mode");
                  while (ibyte < length) {
                    int val = calbyte(buf);
                    if (val != 0) {
                      ibyte++;
                      ss << "ch# " << ch << "(" << val << ") ";
                      hasHit = true;
                      if (febno < 0 || febno > 6) {
                        B2ERROR("FEB is bad : " << LogVar("FEB no.", febno) << LogVar("hslb", finesse) << LogVar("type", type_feb) << LogVar("ver",
                                ver) << LogVar("boardid", boardid) << LogVar("febno", febno) << LogVar("length", length) << LogVar("evtno", evtno));
                      }
                      feb.push_back(ch, val);
                    }
                    ch--;
                    if (ch < 0) break;
                  }
                } else if (type_feb == 0x01) { // Suppressed mode
                  // The below line is commented since it sometimes causes problem during processing threshold scan data.
                  // No harm to comment this line since it is only utilized for threshold scan data.
                  //if (length > 144 * 2 + 10) B2FATAL("error " << LogVar("length", length));
                  //B2INFO("suppreed mode");
                  while (ibyte < length) {
                    int ch = calbyte(buf);
                    ibyte++;
                    int val = calbyte(buf);
                    ibyte++;
                    if (val != 0) {
                      ss << "ch# " << ch << "(" << val << ") ";
                      hasHit = true;
                      if (febno < 0 || febno > 6) {
                        B2ERROR("FEB is bad : " << LogVar("FEB no.", febno) << LogVar("hslb", finesse) << LogVar("type", type_feb) << LogVar("ver",
                                ver) << LogVar("boardid", boardid) << LogVar("febno", febno) << LogVar("length", length) << LogVar("evtno", evtno));
                        return;
                      }
                      feb.push_back(ch, val);
                    }
                  }
                }
                rawdigit->addFEB(feb, type, ver, boardid, febno, length, evtno, feb_trigno);
                if (m_debug && hasHit) {
                  B2INFO(ss.str());
                }
              }
            }
          }
 
        } // end of raw unpacker
    */
    arichinfo->settrgtype(trgtype);
    arichinfo->setpciedata(m_pciedata);
    if (vth_thscan > -1.27) { arichinfo->setvth_thscan(vth_thscan); }
    arichinfo->setntrack(0);
    arichinfo->setnexthit(0);
    arichinfo->setnhit(0);
    if (thscan_mode > 0 || m_rawmode != 0)
    { arichinfo->setthscan_mode(true); }
    else
    { arichinfo->setthscan_mode(false); }
 
  }

FastTracking()

TVector3 FastTracking ( TVector3  dirf, TVector3  r, double *  refind, double *  z, int  n, int  opt  )

private

Calculates the intersection of the Cherenkov photon emitted from point "r" in "dirf" direction with the detector plane.

Parameters

r	vector of photon emission point
dirf	direction of photon emission
n	number of aerogel layers through which photon passes
refind	array of layers refractive indices
z	array of z coordinates of borders between layers

Definition at line 149 of file ARICHReconstruction.cc.

HitVirtualPosition()

TVector3 HitVirtualPosition ( const TVector3 &  hitpos, int  mirrorID  )

private

Returns the hit virtual position, assuming that it was reflected from mirror.

Parameters

hitpos	vector of hit position
mirrorID	id of mirror from which the photon was reflected

Definition at line 223 of file ARICHReconstruction.cc.

initialize() [1/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

Function is called only once at the beginning of your job at the beginning of the corresponding module. Things that can be done here, should be done here, as it is relatively cheap in terms of CPU time.

Reimplemented from Module.

Definition at line 92 of file arichBtestModule.cc.

initialize() [2/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 73 of file ARICHDigitizerModule.cc.

initialize() [3/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

better use isRequired(), but RawFTSW is not in sim

Reimplemented from HistoModule.

Definition at line 161 of file ARICHDQMModule.cc.

initialize() [4/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 57 of file ARICHFillHitsModule.cc.

  {
 
    StoreArray<ARICHDigit> digits;
    digits.isRequired();
 
    StoreArray<ARICHHit> arichHits;
    arichHits.registerInDataStore();
 
    arichHits.registerRelationTo(digits);
  }

initialize() [5/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 52 of file ARICHMCParticlesModule.cc.

  {
    // Dependecies check
    m_tracks.isRequired();
    m_extHits.isRequired();
 
    m_arichMCPs.registerInDataStore();
    m_tracks.registerRelationTo(m_arichMCPs);
 
  }

initialize() [6/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 69 of file ARICHNtupleModule.cc.

  {
 
    if (m_file) delete m_file;
    m_file = new TFile(m_outputFile.c_str(), "RECREATE");
    if (m_file->IsZombie()) {
      B2FATAL("Couldn't open file '" << m_outputFile << "' for writing!");
      return;
    }
 
    m_tree = new TTree("arich", "ARICH validation ntuple");
 
    m_tree->Branch("evt", &m_arich.evt, "evt/I");
    m_tree->Branch("run", &m_arich.run, "run/I");
    m_tree->Branch("exp", &m_arich.exp, "exp/I");
 
    m_tree->Branch("charge", &m_arich.charge, "charge/S");
    m_tree->Branch("pValue", &m_arich.pValue, "pValue/F");
    m_tree->Branch("d0", &m_arich.z0, "d0/F");
    m_tree->Branch("z0", &m_arich.d0, "z0/F");
 
#ifdef ALIGNMENT_USING_BHABHA
    m_tree->Branch("eop", &m_arich.eop, "eop/F");
    m_tree->Branch("e9e21", &m_arich.e9e21, "e9e21/F");
    m_tree->Branch("etot", &m_arich.etot, "etot/F");
#endif
 
    m_tree->Branch("PDG", &m_arich.PDG, "PDG/I");
    m_tree->Branch("motherPDG", &m_arich.motherPDG, "motherPDG/I");
    m_tree->Branch("primary", &m_arich.primary, "primary/S");
    m_tree->Branch("seen", &m_arich.seen, "seen/S");
    m_tree->Branch("scatter", &m_arich.scatter, "scatter/I");
 
    m_tree->Branch("rhoProd", &m_arich.rhoProd, "rhoProd/F");
    m_tree->Branch("zProd",   &m_arich.zProd,   "zProd/F");
    m_tree->Branch("phiProd", &m_arich.phiProd, "phiProd/F");
    m_tree->Branch("rhoDec", &m_arich.rhoDec, "rhoDec/F");
    m_tree->Branch("zDec",   &m_arich.zDec,   "zDec/F");
    m_tree->Branch("phiDec", &m_arich.phiDec, "phiDec/F");
    m_tree->Branch("status", &m_arich.status, "status/I");
 
    m_tree->Branch("detPhot",  &m_arich.detPhot,  "detPhot/I");
    m_tree->Branch("numBkg",  &m_arich.numBkg,  "e/F:mu:pi:K:p:d");
    m_tree->Branch("expPhot",  &m_arich.expPhot,  "e/F:mu:pi:K:p:d");
    m_tree->Branch("logL",  &m_arich.logL,  "e/F:mu:pi:K:p:d");
 
    m_tree->Branch("recHit",  &m_arich.recHit,  "PDG/I:x/F:y:z:p:theta:phi");
    m_tree->Branch("mcHit",  &m_arich.mcHit,  "PDG/I:x/F:y:z:p:theta:phi");
    m_tree->Branch("winHit",  &m_arich.winHit,  "x/F:y");
    m_tree->Branch("nrec", &m_arich.nRec, "nRec/I");
    m_tree->Branch("nCDC", &m_arich.nCDC, "nCDC/I");
    m_tree->Branch("inAcceptance", &m_arich.inAcc, "inAcc/O");
    m_tree->Branch("photons", "std::vector<Belle2::ARICHPhoton>", &m_arich.photons);
 
    // required input
    m_arichTracks.isRequired();
    m_arichLikelihoods.isRequired();
 
    // optional input
    m_tracks.isOptional();
    m_arichMCPs.isOptional();
    m_arichAeroHits.isOptional();
    m_arichInfo.isOptional();
 
    StoreArray<ARICHHit> arichHits;
    arichHits.isRequired();
 
  }

initialize() [7/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 62 of file ARICHPackerModule.cc.

  {
 
    StoreArray<ARICHDigit> digits(m_inputDigitsName);
    digits.isRequired();
 
    StoreArray<RawARICH> rawData(m_outputRawDataName);
    rawData.registerInDataStore();
  }

initialize() [8/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from HistoModule.

Definition at line 94 of file ARICHRateCalModule.cc.

initialize() [9/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from HistoModule.

Definition at line 54 of file ARICHRawUnpackerModule.cc.

initialize() [10/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 49 of file ARICHRelateModule.cc.

  {
    // Dependecies check
    m_mcParticles.isRequired();
    m_mdstTracks.isRequired();
    m_aeroHits.isRequired();
    m_extHits.isRequired();
    m_extHits.registerRelationTo(m_aeroHits);
  }

initialize() [11/11]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 73 of file ARICHUnpackerModule.cc.

magFieldDistorsion()

void magFieldDistorsion	(	TVector2 &	hit,
		int	copyno
	)

Apply correction to hit position due to non-perpendicular component of magnetic field.

Parameters

hit	local position of simhit
copyno	copy number of hapd

Definition at line 190 of file ARICHDigitizerModule.cc.

mergerClusterHitMap1D()

TH1 * mergerClusterHitMap1D	(	TH1 *	hitMap,
		int	mergerID
	)

Make 1D hit map of specified Merger Board.

Parameters

1D	hit map of all channels
Merfer	board identifier [1:72]

Returns

1D hit map of the merger board

Definition at line 95 of file hitMapMaker.cc.

  {
 
    int m_mergerID = mergerID;
 
    DBObjPtr<ARICHChannelMapping> arichChMap;
    DBObjPtr<ARICHMergerMapping> arichMergerMap;
 
    TH1* m_hitMap = hitMap;
 
    std::vector<int> moduleIDs;
    for (int i = 1; i < 7; i++) {
      moduleIDs.push_back(arichMergerMap->getModuleID(m_mergerID, i));
    }
 
    TH1D* m_mergerHitMap1D = new TH1D("MergerHitMap1D", Form("Hit map in Merger Board %d", m_mergerID), 144 * 6, -0.5, 144 * 6 - 0.5);
    for (int i = 1; i < 7; i++) {
      for (int j = 0; j < 144; j++) {
        int hitsNum = m_hitMap->GetBinContent((moduleIDs[i] - 1) * 144 + i);
        m_mergerHitMap1D->Fill(144 * (i - 1) + j, hitsNum);
      }
    }
    return m_mergerHitMap1D;
  }

mergerClusterHitMap2D()

TCanvas * mergerClusterHitMap2D	(	TH1 *	hitMap,
		int	mergerID
	)

Make display of 6 HAPDs' 2D hit map of the Merger Board.

Parameters

1D	hit map of all channels
Merfer	board identifier [1:72]

Returns

Display of 6 HAPDs' 2D hit map

Definition at line 120 of file hitMapMaker.cc.

moduleDeadMap()

TH2 * moduleDeadMap	(	TH1 *	hitMap,
		int	moduleID
	)

Make chip dead/alive map in HAPD view (2*2 chips)

Parameters

1D	hit map of all channels
Module	identifier [1:420]

Returns

2D dead/alive map on HAPD

Definition at line 62 of file hitMapMaker.cc.

moduleHitMap()

TH2 * moduleHitMap	(	TH1 *	hitMap,
		int	moduleID
	)

Make hit map in HAPD view (12*12 channels)

Parameters

1D	hit map of all channels
Module	identifier [1:420]

Returns

2D hit map on HAPD

Definition at line 38 of file hitMapMaker.cc.

printBits()

void printBits ( const int *  buffer, int  bufferSize  )

private

Unpack raw data given in production format.

Parameters

buffer	raw data buffer
bufferSize	buffer size

Definition at line 529 of file ARICHUnpackerModule.cc.

sectorDeadMap()

TCanvas * sectorDeadMap	(	TH1 *	hitMap,
		int	sector
	)

Make display of 70 HAPDs' 2D dead/alive map of the sector.

Parameters

1D	hit map of all channels
Sector	identifier [1:6]

Returns

Display of 70 HAPDs' 2D dead/alive map

Definition at line 176 of file hitMapMaker.cc.

sectorHitMap()

TCanvas * sectorHitMap	(	TH1 *	hitMap,
		int	sector
	)

Make display of 70 HAPDs' 2D hit map of the sector.

Parameters

1D	hit map of all channels
Sector	identifier [1:6]

Returns

Display of 70 HAPDs' 2D hit map

Definition at line 144 of file hitMapMaker.cc.

terminate() [1/2]

void terminate ( void  )

overridevirtual

Is called at the end of your Module.

Function is called only once at the end of your job at the end of the corresponding module. This function is for cleaning up, closing files, etc.

Reimplemented from Module.

Definition at line 567 of file arichBtestModule.cc.

terminate() [2/2]

void terminate ( void  )

overridevirtual

Termination action.

Clean-up, close files, summarize statistics, etc.

Reimplemented from Module.

Definition at line 336 of file ARICHNtupleModule.cc.