BmnLambdaMisc.h

Go to the documentation of this file.

// @(#)bmnroot/embedding:$Id$
// Author: Pavel Batyuk <pavel.batyuk@jinr.ru> 2019-12-07
 
////////////////////////////////////////////////////////////////////////////////
//                                                                            //
// BmnLambdaMisc                                                              //
//                                                                            //
// Useful instruments to work with mapping files (GEM, SILICON)               //
// (reading, sorting, correspondence e.t.c.)                                  //
//                                                                            //
////////////////////////////////////////////////////////////////////////////////
 
#ifndef BMNLAMBDAMISC_H
#define BMNLAMBDAMISC_H 1
 
#include <BmnCSCDigit.h>
#include <BmnCSCStationSet.h>
#include <BmnGemStripDigit.h>
#include <BmnGemStripStationSet.h>
#include <BmnLambdaEmbeddingQa.h>
#include <BmnSiliconDigit.h>
#include <BmnSiliconStationSet.h>
#include <TClonesArray.h>
#include <TNamed.h>
#include <TSystem.h>
#include <fstream>
#include <iostream>
#include <map>
#include <vector>
 
using namespace std;
 
class MappingInfo : public TObject
{
  public:
    MappingInfo() {};
 
    MappingInfo(Int_t stat, TString map, pair<Int_t, Int_t> ch, pair<Int_t, Int_t> str)
        : channel(-1)
        , strip(-1)
    {
        station = stat;
        mapFile = map;
        channels = ch;
        strips = str;
    }
 
    MappingInfo(Int_t stat, TString map, Int_t ch, Int_t str)
        : channels(make_pair(-1, -1.))
        , strips(make_pair(-1, -1))
    {
        station = stat;
        mapFile = map;
        channel = ch;
        strip = str;
    }
 
    MappingInfo(Int_t stat, TString mapping, map<Int_t, Int_t> strGlobChan)
        : channels(make_pair(-1, -1.))
        , strips(make_pair(-1, -1))
        , channel(-1)
        , strip(-1)
    {
        station = stat;
        mapFile = mapping;
        stripChan = strGlobChan;
    }
 
    virtual ~MappingInfo() {};
 
    Int_t station;
    TString mapFile;
    pair<Int_t, Int_t> channels;
    pair<Int_t, Int_t> strips;
 
    // Additional data members ...
    Int_t channel;
    Int_t strip;
 
    map<Int_t, Int_t> stripChan;
};
 
class BmnLambdaMisc : public TNamed
{
  public:
    BmnLambdaMisc();
 
    virtual ~BmnLambdaMisc();
 
    /* GEM */
    Int_t GemDigiToChannel(BmnStripDigit*, Long_t&);             // GEM-digi ---> channel (serial)
    TString GemDigiToMapping(BmnStripDigit*);                    // GEM-digi ---> corresponding <strip-channel> mapping
    Long_t GemDigiChannelToSerial(pair<BmnStripDigit, Int_t>);   // GEM-digi + channel ---> serial
    Long_t GetGemSerial(Int_t, Int_t, Int_t, Int_t);             // (st, mod, id, channel ---> serial)
 
    /* CSC */
    Int_t CscDigiToChannel(BmnStripDigit*, Long_t&);   // CSC-digi ---> channel (serial)
 
    /* SILICON */
    void SiliconDigiToChannelSampleSerial(BmnStripDigit*,
                                          Int_t&,
                                          Int_t&,
                                          Long_t&);   // SILICON-digi ---> &channel, &sample, &serial
 
    /* Supplementary methods to be used for testing */
    void CheckStripOverlaps();   // To be used for checking whether strips from common serial are overlapped with those
                                 // ones from other serials or not ...
 
  private:
    void GEM();
    void SILICON();
    void CSC();
 
  private:
    /* GEM */
    vector<vector<Int_t>> fGemStatModId;   // vector of (stat, module, id)
    /* SILICON */
    vector<vector<Int_t>> fSiliconStatModLayId;   // vector of (stat, module, layer, id)
    /* CSC */
    vector<vector<Int_t>> fCscStatModLay;   // vector of (stat, module, layer)
 
    /* GEM */
    // Map with commonSerial only ...
    map<vector<Int_t>, vector<Long_t>> corrMapCommonSerial;   // <st, mod, id> ---> <low, high, serial>
 
    // Map with no commonSerial included ...
    map<vector<Int_t>, vector<Long_t>> corrMapNoCommonSerial;   // <st, mod, id> ---> <low, high, serial>
 
    /* SILICON */
    // Map with SILICON serials ...
    map<vector<Int_t>, vector<Long_t>> serialsSilicon;   // <st, mod, layer, id> ---> <low, high, serial>
 
    /* CSC */
    // Map with CSC serials ...
    map<vector<Int_t>, vector<Long_t>> serialsCsc;   // <st, mod, layer> ---> <low, high, serial>
 
    // Geometry sets ...
    BmnGemStripStationSet* fDetectorGEM;
    BmnSiliconStationSet* fDetectorSI;
    BmnCSCStationSet* fDetectorCSC;
 
    UInt_t FindChannelByStrip(TString, Int_t);
 
    void ParseStripChannelMapping(TString, Int_t, map<Int_t, Int_t>&);
 
    ClassDef(BmnLambdaMisc, 1)
};
 
#endif