CgemDigitizerSvc Class Reference

#include <CgemDigitizerSvc.h>

Inheritance diagram for CgemDigitizerSvc:

Public Member Functions
	CgemDigitizerSvc (const std::string &name, ISvcLocator *svcloc)
	~CgemDigitizerSvc ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
StatusCode	setTrack (int layer, int particle, int charge, double p, double trkPosIn[], double trkPosOut[])
StatusCode	setTrack (int layer, std::vector< int >particle, std::vector< int > charge, std::vector< double > p, std::vector< vector< double > > trkPosIn, std::vector< vector< double > > trkPosOut)
int	getNXstrips () const
int	getNVstrips () const
int	getXstripSheet (int n) const
int	getXstripID (int n) const
int	getVstripSheet (int n) const
int	getVstripID (int n) const
double	getXstripQ2 (int n) const
double	getVstripQ2 (int n) const
double	getXstripQ (int n) const
double	getVstripQ (int n) const
double	getXstripT (int n) const
double	getVstripT (int n) const
double	getXfirstT (int n) const
double	getVfirstT (int n) const
double	getQsaturation (int layer, int sheet, int view, int id)
Public Member Functions inherited from ICgemDigitizerSvc
virtual StatusCode	setTrack (int layer, int particle, int charge, double p, double trkPosIn[], double trkPosOut[])=0
virtual StatusCode	setTrack (int layer, std::vector< int >particle, std::vector< int > charge, std::vector< double > p, std::vector< std::vector< double > > trkPosIn, std::vector< std::vector< double > > trkPosOut)=0
virtual int	getNXstrips () const =0
virtual int	getNVstrips () const =0
virtual int	getXstripID (int n) const =0
virtual int	getVstripID (int n) const =0
virtual double	getXstripQ (int n) const =0
virtual double	getVstripQ (int n) const =0
virtual double	getXstripT (int n) const =0
virtual double	getVstripT (int n) const =0

Additional Inherited Members
Static Public Member Functions inherited from ICgemDigitizerSvc
static const InterfaceID &	interfaceID ()

Detailed Description

Definition at line 44 of file CgemDigitizerSvc.h.

Constructor & Destructor Documentation

CgemDigitizerSvc()

CgemDigitizerSvc::CgemDigitizerSvc	(	const std::string &	name,
		ISvcLocator *	svcloc
	)

Definition at line 38 of file CgemDigitizerSvc.cxx.

                                                                           :Service (name, svcloc){
     // declare properties
     declareProperty("IonizationModel", m_ionModel=1);
     declareProperty("DriftAvalancheModel", m_driftAvaModel=1); // 3 mean SamplingGar2, 
     declareProperty("InductionModel", m_inductionModel=1);     // 3 mean InductionGar2,
                                                                // this 2 doesnot extend from base normally and have to be used in conjunction.
     declareProperty("GarfieldDebugging", m_garfDebugging=false);
     declareProperty("SamplingDebugging", m_samplingDebugging=false);
     declareProperty("InductionDebugging", m_debugInduction=false);
     declareProperty("MagConfig", m_magConfig=1);
     declareProperty("LUTFilePath",m_LUTFilePath = "/bes3fs/cgemCosmic/data/CGEM_cosmic_look_up_table_17.root");
     declareProperty("QBranchVsampleDelay",sample_delay = 162.5);
     declareProperty("StoreUnderThrethold",m_storeFlag = false);
     declareProperty("Saturation",m_saturation = true);
     declareProperty("SaveNt",m_saveNt= false);
     declareProperty("SamplingGarGainMultiplier",m_GainMultiplier);
     declareProperty("SamplingGarTransMultiplier",m_TransMultiplier= 1.0);
     declareProperty("SamplingGarDiffuMultiplier",m_DiffuMultiplier= 1.0);
     //cout << "CgemDigitizerSvc::m_ionModel = " << m_ionModel << endl;
     declareProperty("Ngaps_microSector",m_Ngaps_microSector= 40);
     declareProperty("gapShift_microSector",m_gapShift_microSector);
     declareProperty("gap_microSector", m_gap_microSector= 1.1);// in mm
     declareProperty("microSector_width", m_microSector_width);// in rad
     declareProperty("QinGausSigma", m_QinGausSigma);
     declareProperty("ScaleSignalX", m_ScaleSignalX=1.0);
}

~CgemDigitizerSvc()

CgemDigitizerSvc::~CgemDigitizerSvc

(

)

Definition at line 65 of file CgemDigitizerSvc.cxx.

                                   {
}

Member Function Documentation

finalize()

StatusCode CgemDigitizerSvc::finalize ( )

virtual

Definition at line 272 of file CgemDigitizerSvc.cxx.

                                     {
     MsgStream log(messageService(), name());
     log << MSG::INFO << "CgemDigitizerSvc::finalize()" << endreq;
     delete m_pIon;
     delete m_pDriftAndAva;
     delete m_pInduction;
 
     return StatusCode::SUCCESS;
}

getNVstrips()

int CgemDigitizerSvc::getNVstrips ( ) const

inlinevirtual

Implements ICgemDigitizerSvc.

Definition at line 61 of file CgemDigitizerSvc.h.

{return m_nVstrips;}

Referenced by BesCgemDigitizer::Digitize_v3().

getNXstrips()

int CgemDigitizerSvc::getNXstrips ( ) const

inlinevirtual

Implements ICgemDigitizerSvc.

Definition at line 60 of file CgemDigitizerSvc.h.

{return m_nXstrips;}

Referenced by BesCgemDigitizer::Digitize_v3().

getQsaturation()

double CgemDigitizerSvc::getQsaturation ( int  layer, int  sheet, int  view, int  id  )

inline

Definition at line 75 of file CgemDigitizerSvc.h.

{ return Qsaturation[layer][sheet][view][id]; }

Referenced by BesCgemDigitizer::Digitize_v3().

getVfirstT()

double CgemDigitizerSvc::getVfirstT ( int  n ) const

inline

Definition at line 74 of file CgemDigitizerSvc.h.

{return m_vfirstT[n];}

getVstripID()

int CgemDigitizerSvc::getVstripID ( int  n ) const

inlinevirtual

Implements ICgemDigitizerSvc.

Definition at line 65 of file CgemDigitizerSvc.h.

{return m_vstripID[n];}

Referenced by BesCgemDigitizer::Digitize_v3().

getVstripQ()

double CgemDigitizerSvc::getVstripQ ( int  n ) const

inlinevirtual

Implements ICgemDigitizerSvc.

Definition at line 70 of file CgemDigitizerSvc.h.

{return m_vstripQ2[n];}

Referenced by BesCgemDigitizer::Digitize_v3().

getVstripQ2()

double CgemDigitizerSvc::getVstripQ2 ( int  n ) const

inline

Definition at line 68 of file CgemDigitizerSvc.h.

{return m_vstripQ[n];}

getVstripSheet()

int CgemDigitizerSvc::getVstripSheet ( int  n ) const

inline

Definition at line 64 of file CgemDigitizerSvc.h.

{return m_vstripSheet[n];}

Referenced by BesCgemDigitizer::Digitize_v3().

getVstripT()

double CgemDigitizerSvc::getVstripT ( int  n ) const

inlinevirtual

Implements ICgemDigitizerSvc.

Definition at line 72 of file CgemDigitizerSvc.h.

{return m_vstripT[n];}

Referenced by BesCgemDigitizer::Digitize_v3().

getXfirstT()

double CgemDigitizerSvc::getXfirstT ( int  n ) const

inline

Definition at line 73 of file CgemDigitizerSvc.h.

{return m_xfirstT[n];}

getXstripID()

int CgemDigitizerSvc::getXstripID ( int  n ) const

inlinevirtual

Implements ICgemDigitizerSvc.

Definition at line 63 of file CgemDigitizerSvc.h.

{return m_xstripID[n];}

Referenced by BesCgemDigitizer::Digitize_v3().

getXstripQ()

double CgemDigitizerSvc::getXstripQ ( int  n ) const

inlinevirtual

Implements ICgemDigitizerSvc.

Definition at line 69 of file CgemDigitizerSvc.h.

{return m_xstripQ2[n];}//sum of induced charge simulated

Referenced by BesCgemDigitizer::Digitize_v3().

getXstripQ2()

double CgemDigitizerSvc::getXstripQ2 ( int  n ) const

inline

Definition at line 67 of file CgemDigitizerSvc.h.

{return m_xstripQ[n];}//measured Q

getXstripSheet()

int CgemDigitizerSvc::getXstripSheet ( int  n ) const

inline

Definition at line 62 of file CgemDigitizerSvc.h.

{return m_xstripSheet[n];}

Referenced by BesCgemDigitizer::Digitize_v3().

getXstripT()

double CgemDigitizerSvc::getXstripT ( int  n ) const

inlinevirtual

Implements ICgemDigitizerSvc.

Definition at line 71 of file CgemDigitizerSvc.h.

{return m_xstripT[n];}

Referenced by BesCgemDigitizer::Digitize_v3().

initialize()

StatusCode CgemDigitizerSvc::initialize ( )

virtual

Definition at line 77 of file CgemDigitizerSvc.cxx.

                                       {
    MsgStream log(messageService(), name());
    log << MSG::INFO << "CgemDigitizerSvc::initialize()" << endreq;
 
    StatusCode sc = Service::initialize();
    if( sc.isFailure() ) return sc;
 
    static IJobOptionsSvc* jobSvc = 0;
    if ( jobSvc == 0 ) {
        sc = service("JobOptionsSvc", jobSvc);
        if ( sc.isFailure() ) {
            std::cout << "Can't get the JobOptionsSvc @ DistBoss::GetPropertyValue()" << std::endl;
            return sc;
        }
    }
 
    // const vector<const Property*>* properties = jobSvc->getProperties(client);
    const vector<const Property*>* properties = jobSvc->getProperties("BesRndmGenSvc");
    // cout << "client of jobSvc: " << jobSvc->getClients() << endl;
    if ( properties == NULL ) {
        std::cout << "In CgemDigitizerSvc::initialize(), can't get client: " << std::endl;
        return StatusCode::FAILURE;
    }
 
    unsigned int randSeed;
    for ( unsigned int i = 0; i < properties->size(); ++i ) {
        if ( properties->at(i)->name() == "RndmSeed" ) {
            cout << "name of property: " << properties->at(i)->name() << endl;
            string strRnd = properties->at(i)->toString();
            sscanf(strRnd.c_str(), "%u", &randSeed);
            cout << "random seed from jobOption: " << randSeed << endl;
            break;
        }
    }
 
    sc = service("CgemGeomSvc", m_geomSvc); // initial pointer of CgemGeomSvc
    if(sc != StatusCode::SUCCESS)
    {
        log << MSG::ERROR << "can not use CgemGeomSvc" << endreq;
        return StatusCode::FAILURE;
    }
 
    gRandom->SetSeed(randSeed); // set random seed
 
    cout << "ionModel = " << m_ionModel << endl;
    if(1 == m_ionModel){
        m_pIon = new IonizationGar();
    } else if(2 == m_ionModel){
        m_pIon = new IonizationGTS();
    }
    // m_pIon = new IonizationGar();
    m_pIon->setDebugging(m_garfDebugging);
    m_pIon->init(randSeed, m_geomSvc, m_magConfig);
 
    if ((3==m_driftAvaModel && 3!=m_inductionModel )||(3!=m_driftAvaModel && 3==m_inductionModel )) {
        cout<<"CgemDigitizerSvc::initialize() Error: DriftAvalancheModel type3 and InductionModel type3 must be used together!"<<endl;
        return StatusCode::FAILURE;
    }
 
    if(1 == m_driftAvaModel){
        m_pDriftAndAva = new SamplingGar();
        SamplingGar *pDriftAndAvaGar=(SamplingGar *) m_pDriftAndAva;
        pDriftAndAvaGar->setGainMultiplier(m_GainMultiplier);
        pDriftAndAvaGar->setTransMultiplier(m_TransMultiplier);
        pDriftAndAvaGar->setDiffuMultiplier(m_DiffuMultiplier);
    } else if(2 == m_driftAvaModel){
        m_pDriftAndAva = new SamplingGTS();
    } else if(3 == m_driftAvaModel){
        m_pDriftAndAva = new SamplingGar2();
        SamplingGar2* pDriftAndAva=(SamplingGar2*)m_pDriftAndAva;
        pDriftAndAva->setMultiElectronMapAddr(&m_hitmap);
        
        // migrate of samGar1;
        pDriftAndAva->setGainMultiplier(m_GainMultiplier);
        pDriftAndAva->setTransMultiplier(m_TransMultiplier);
        pDriftAndAva->setDiffuMultiplier(m_DiffuMultiplier);
        // migrate of indGar1
        pDriftAndAva->setNgapsSect(m_Ngaps_microSector);
        pDriftAndAva->setGapSizeSect(m_gap_microSector);
        pDriftAndAva->setGapShiftSect(m_gapShift_microSector);
        pDriftAndAva->setMicroSectorWidthRad(m_microSector_width);
        
    }
    m_pDriftAndAva->init(m_geomSvc, m_magConfig);
    m_pDriftAndAva->setDebugging(m_samplingDebugging);
 
    if(1 == m_inductionModel){
        InductionGar* indGar = new InductionGar();
        indGar->setNgapsSect(m_Ngaps_microSector);
        indGar->setGapSizeSect(m_gap_microSector);
        indGar->setGapShiftSect(m_gapShift_microSector);
        indGar->setMicroSectorWidthRad(m_microSector_width);
        indGar->setQinGausSigma(m_QinGausSigma);
        indGar->setScaleSignalX(m_ScaleSignalX);
        m_pInduction = indGar;
        //m_pInduction = new InductionGar();
        m_pInduction->setVsampleDelay(sample_delay);
        m_pInduction->setLUTFilePath(m_LUTFilePath);
        m_pInduction->setStoreFlag(m_storeFlag);
        m_pInduction->setSaturation(m_saturation);
    } else if(2 == m_inductionModel){
        m_pInduction = new InductionGTS();
        //m_pInduction = new InductionGar();
    }else if (3 == m_inductionModel){
        m_pInduction = new InductionGar2();
        InductionGar2* indGar=(InductionGar2*)m_pInduction;
        indGar->setQinGausSigma(m_QinGausSigma);
        indGar->setScaleSignalX(m_ScaleSignalX);
 
        m_pInduction->setVsampleDelay(sample_delay);
        m_pInduction->setLUTFilePath(m_LUTFilePath);
        m_pInduction->setStoreFlag(m_storeFlag);
        m_pInduction->setSaturation(m_saturation);
    }
    m_pInduction->init(m_geomSvc, m_magConfig);
    m_pInduction->setDebugOutput(m_debugInduction);
 
    // --- check gaps
    cout<<"CgemDigitizerSvc:"<<endl;
    cout<<"gapShift_microSector="<<m_gapShift_microSector[0]
        <<", "<<m_gapShift_microSector[1]
        <<", "<<m_gapShift_microSector[2]
        <<", "<<m_gapShift_microSector[3]
        <<", "<<m_gapShift_microSector[4]
        <<", "<<m_gapShift_microSector[5]
        <<endl;
    cout<<"microSector_width="<<m_microSector_width[0]
        <<", "<<m_microSector_width[1]
        <<", "<<m_microSector_width[2]
        <<", "<<m_microSector_width[3]
        <<", "<<m_microSector_width[4]
        <<", "<<m_microSector_width[5]
        <<endl;
    cout<<"QinGausSigma="<<m_QinGausSigma[0]
        <<", "<<m_QinGausSigma[1]
        <<endl;
    cout<<"Ngaps_microSector="<<m_Ngaps_microSector<<endl;
    cout<<"gap_microSector="<<m_gap_microSector<<endl;
    
    // --- read saturations
    TFile *LUTFile = TFile::Open(m_LUTFilePath.c_str(),"read");
    TTree *tree = (TTree*)LUTFile->Get("tree");
    Float_t QDC_saturation;
    Int_t layer,sheet,strip_v,strip_x;
    tree->SetBranchAddress("strip_x_boss", &strip_x);
    tree->SetBranchAddress("strip_v_boss", &strip_v);
    tree->SetBranchAddress("layer", &layer);
    tree->SetBranchAddress("sheet", &sheet);
    tree->SetBranchAddress("calib_QDC_saturation", &QDC_saturation);
    for(int i=0;i<tree->GetEntries();i++) {
        tree->GetEntry(i);
        if(strip_x!=-1) {
            Qsaturation[layer][sheet][0][strip_x] = QDC_saturation;
        }
        if(strip_v!=-1) {
            Qsaturation[layer][sheet][1][strip_v] = QDC_saturation;
        }
    }
    for(int i=0;i<832;i++) {
        Qsaturation[2][0][0][i] = 45.;
        Qsaturation[2][1][0][i] = 45.;
    }
    for(int i=0;i<1395;i++) {
        Qsaturation[2][0][1][i] = 45.;
        Qsaturation[2][1][1][i] = 45.;
    }
    LUTFile->Close();
 
    // ntuple
    if(m_saveNt) {
        INTupleSvc* ntupleSvc;
        Gaudi::svcLocator() -> service("NTupleSvc", ntupleSvc);
        NTuplePtr nt(ntupleSvc, "CgemDigitizerSvcNT/digi");
        // NTuplePtr nt(ntupleSvc(), "CgemDigitizerSvcNT/digi");
        if(nt) m_tuple = nt;
        else{
            m_tuple = ntupleSvc->book("CgemDigitizerSvcNT/digi", CLID_ColumnWiseTuple, "noise");
            if(m_tuple){
                m_tuple->addItem("nIonE",          m_ntNIonE);
                m_tuple->addItem("nMultiE",        m_ntNMultiE);
                m_tuple->addItem("nXstrips",       m_ntNxstrips,0,100);
                m_tuple->addItem("nVstrips",       m_ntNvstrips,0,100);
                m_tuple->addIndexedItem("XstripQ",   m_ntNxstrips,m_ntxstripQ);
                m_tuple->addIndexedItem("XstripT",   m_ntNxstrips,m_ntxstripT);
                m_tuple->addIndexedItem("VstripQ",   m_ntNvstrips,m_ntvstripQ);
                m_tuple->addIndexedItem("VstripT",   m_ntNvstrips,m_ntvstripT);
                m_tuple->addIndexedItem("XstripID",  m_ntNxstrips,m_XstripID);
                m_tuple->addIndexedItem("VstripID",  m_ntNvstrips,m_VstripID);
            }
        }
    }
 
    return StatusCode::SUCCESS;
}

queryInterface()

StatusCode CgemDigitizerSvc::queryInterface ( const InterfaceID &  riid, void **  ppvUnknown  )

virtual

Definition at line 68 of file CgemDigitizerSvc.cxx.

                                                                                       {
     if( IID_ICgemDigitizerSvc.versionMatch(riid) ){
      *ppvInterface = static_cast<ICgemDigitizerSvc*> (this);
     } else{
      return Service::queryInterface(riid, ppvInterface);
     }
     return StatusCode::SUCCESS;
}

setTrack() [1/2]

StatusCode CgemDigitizerSvc::setTrack ( int  layer, int  particle, int  charge, double  p, double  trkPosIn[], double  trkPosOut[]  )

virtual

Implements ICgemDigitizerSvc.

Definition at line 282 of file CgemDigitizerSvc.cxx.

                                                                                                                         {
    vector<int> Vparticle             (1,particle);
    vector<int> Vcharge               (1,charge);
    vector<double> Vp                 (1,p);
 
    // awkward way of saying 
    //vector<vector<double> >VtrkPosIn={{trkPosIn[0], trkPosIn[1], trkPosIn[2]}};
    vector<double> posin;
    posin.push_back(trkPosIn[0]);
    posin.push_back(trkPosIn[1]);
    posin.push_back(trkPosIn[2]);
    vector<vector<double> > VtrkPosIn;
    VtrkPosIn.push_back(posin);
 
    vector<double> posout;
    posout.push_back(trkPosOut[0]);
    posout.push_back(trkPosOut[1]);
    posout.push_back(trkPosOut[2]);
    vector<vector<double> > VtrkPosOut;
    VtrkPosOut.push_back(posout);//= {{trkPosOut[0], trkPosOut[1], trkPosOut[2]}};
 
    StatusCode sc=setTrack(layer,Vparticle, Vcharge, Vp, VtrkPosIn, VtrkPosOut);
 
    return sc;
 
    //*******below are original contents of this function.**********
 
    //  //ProcInfo_t info;
    //  //gSystem->GetProcInfo(&info);
    //  //double currentMemory = info.fMemResident/1024;// MB
    //  //cout<<"CgemDigitizerSvc::setTrack() memory 1 "<<currentMemory<<" MB"<<endl;
    //
    //  double time_spent = 0.;
    //  clock_t begin = clock();
    //  //cout<<"setTrack"<<endl;
    //  clear();
    //  m_layer = layer;
    //
    //  //cout<<"before m_pIon->setTrack()"<<endl;
    //  //cout<<"set particle "<<particle<<", p="<<p<<endl;
    //  m_pIon->setTrack(particle, charge, p, trkPosIn, trkPosOut);
    //  //cout<<"after m_pIon->setTrack()"<<endl;
    //  //gSystem->GetProcInfo(&info); currentMemory = info.fMemResident/1024; cout<<"CgemDigitizerSvc::setTrack() memory 2 "<<currentMemory<<" MB"<<endl;
    //
    //  vector<double> ionEx;
    //  vector<double> ionEy;
    //  vector<double> ionEz;
    //  vector<double> ionEt;
    //  int nIonE = m_pIon->getNumberIonE();
    //  for(int i=0; i<nIonE; i++){
    //      ionEx.push_back(m_pIon->getEx(i));
    //      ionEy.push_back(m_pIon->getEy(i));
    //      ionEz.push_back(m_pIon->getEz(i));
    //      ionEt.push_back(m_pIon->getEt(i));
    //  }
    //  //cout<<"after pushing back ionE "<<endl;
    //  if(m_saveNt) m_ntNIonE=nIonE;
    //  clock_t end = clock();
    //  //cout<<" CLOCKS_PER_SEC "<<endl;
    //  //cout<<CLOCKS_PER_SEC<<endl;
    //  //cout<<" calculate time_spent"<<endl;
    //  time_spent += (double)(end - begin) / CLOCKS_PER_SEC;
    //  //cout << "CgemDigitizerSvc::ionization " << time_spent << " seconds" << endl;
    //  // ---------------------------
    //  time_spent = 0.;
    //  begin = clock();
    //
    //  //gSystem->GetProcInfo(&info); currentMemory = info.fMemResident/1024; cout<<"CgemDigitizerSvc::setTrack() memory 3 "<<currentMemory<<" MB"<<endl;
    //  m_pDriftAndAva->setIonElectrons(layer, nIonE, ionEx, ionEy, ionEz, ionEt);
    //
    //  m_nMultiE = m_pDriftAndAva->getNelectrons();
    //  //for(int i=0; i<m_nMultiE; i++){
    //  //  double x = m_pDriftAndAva->getX(i);
    //  //  double y = m_pDriftAndAva->getY(i);
    //  //  double z = m_pDriftAndAva->getZ(i);
    //  //  double t = m_pDriftAndAva->getT(i);
    //  //  (*m_pmultiEx).push_back(x);
    //  //  (*m_pmultiEy).push_back(y);
    //  //  (*m_pmultiEz).push_back(z);
    //  //  (*m_pmultiEt).push_back(t);
    //  //}
    //  m_pmultiEx=& (m_pDriftAndAva->getXContainer());
    //  m_pmultiEy=& (m_pDriftAndAva->getYContainer());
    //  m_pmultiEz=& (m_pDriftAndAva->getZContainer());
    //  m_pmultiEt=& (m_pDriftAndAva->getTContainer());
    //  //std::cout<<m_nMultiE<<"\t";
    //  if(m_saveNt) m_ntNMultiE=m_nMultiE;
    //  end = clock();
    //  time_spent += (double)(end - begin) / CLOCKS_PER_SEC;
    //  //cout << "CgemDigitizerSvc::diffusion " << time_spent << " seconds" << endl;
    //  // --------------------
    //  time_spent = 0.;
    //  begin = clock();
    //
    //  //gSystem->GetProcInfo(&info); currentMemory = info.fMemResident/1024; cout<<"CgemDigitizerSvc::setTrack() memory 4 "<<currentMemory<<" MB"<<endl;
    //  // Induction process
    //  // positionToStrips();
    //  //cout<<" before m_pInduction->setMultiElectrons "<<endl;
    //  m_pInduction->setMultiElectrons(layer, m_nMultiE, *m_pmultiEx, *m_pmultiEy, *m_pmultiEz, *m_pmultiEt);
    //  //cout<<" after  m_pInduction->setMultiElectrons "<<endl;
    //  m_nXstrips = m_pInduction->getNXstrips();
    //  //gSystem->GetProcInfo(&info); currentMemory = info.fMemResident/1024; cout<<"CgemDigitizerSvc::setTrack() memory 5 "<<currentMemory<<" MB"<<endl;
    //  //m_ntNxstrips = m_nXstrips;
    //  //m_ntNvstrips = m_nVstrips;
    //  //  cout<<"read data"<<endl;
    //  for(int i=0; i<m_nXstrips; i++){
    //      //m_ntNxstrips=i;
    //      if(m_saveNt && i<100) {
    //          m_ntxstripQ[i]=m_pInduction->getXstripQ(i);
    //          m_ntxstripT[i]=m_pInduction->getXstripT(i);
    //          m_XstripID[i]=m_pInduction->getXstripID(i);
    //      }
    //      m_xstripSheet.push_back(m_pInduction->getXstripSheet(i));
    //      m_xstripID.push_back(m_pInduction->getXstripID(i));
    //      m_xstripQ.push_back(m_pInduction->getXstripQ(i));
    //      m_xstripQ2.push_back(m_pInduction->getXstripQ_Branch(i));
    //      m_xstripT.push_back(m_pInduction->getXstripT(i));
    //      m_xfirstT.push_back(m_pInduction->getXfirstT(i));
    //
    //  }
    //
    //  m_nVstrips = m_pInduction->getNVstrips();
    //  for(int i=0; i<m_nVstrips; i++){
    //      //m_ntNvstrips=i;
    //      if(m_saveNt && i<100) {
    //          m_ntvstripQ[i]=m_pInduction->getVstripQ(i);
    //          m_ntvstripT[i]=m_pInduction->getVstripT(i);
    //          m_VstripID[i]=m_pInduction->getVstripID(i);
    //      }
    //      m_vstripSheet.push_back(m_pInduction->getVstripSheet(i));
    //      m_vstripID.push_back(m_pInduction->getVstripID(i));
    //      m_vstripQ.push_back(m_pInduction->getVstripQ(i));
    //      m_vstripQ2.push_back(m_pInduction->getVstripQ_Branch(i));
    //      m_vstripT.push_back(m_pInduction->getVstripT(i));
    //      m_vfirstT.push_back(m_pInduction->getVfirstT(i));
    //  }
    //  end = clock();
    //  time_spent += (double)(end - begin) / CLOCKS_PER_SEC;
    //  // cout << "CgemDigitizerSvc::induction " << time_spent << " seconds" << endl;
    //  // cout << "nXstrips = " << m_nXstrips << " ,   m_nVstrips = " << m_nVstrips << endl;
    //  // ---------------
    //  time_spent = 0.;
    //  begin = clock();
    //
    //  if(m_saveNt) {
    //      //m_ntNIonE = nIonE;
    //      //m_ntNMultiE = m_nMultiE;
    //      m_ntNxstrips = m_nXstrips;
    //      m_ntNvstrips = m_nVstrips;
    //      m_tuple->write();
    //  }
    //
    //  end = clock();
    //  time_spent += (double)(end - begin) / CLOCKS_PER_SEC;
    //  //cout << "CgemDigitizerSvc::write " << time_spent << " seconds" << endl;
    //
    //  //std::cout<<"wow test2!\t";
    //  //checkMemorySize();
    //  //gSystem->GetProcInfo(&info);
    //  //double currentMemory2 = info.fMemResident/1024;// MB
    //  //cout<<"CgemDigitizerSvc::setTrack() memory from "<<currentMemory<<" to "<<currentMemory2<<" MB"<<endl;
    //
    //  return StatusCode::SUCCESS;
}

Referenced by BesCgemDigitizer::Digitize_v3(), and setTrack().

setTrack() [2/2]

StatusCode CgemDigitizerSvc::setTrack	(	int	layer,
		std::vector< int >	particle,
		std::vector< int >	charge,
		std::vector< double >	p,
		std::vector< vector< double > >	trkPosIn,
		std::vector< vector< double > >	trkPosOut
	)

---

The documentation for this class was generated from the following files:

• CgemDigitizerSvc.h
• CgemDigitizerSvc.cxx