BesCgemSD Class Reference

#include <BesCgemSD.hh>

Inheritance diagram for BesCgemSD:

Public Member Functions
	BesCgemSD (G4String)
	~BesCgemSD ()
void	Initialize (G4HCofThisEvent *)
void	setGeomPtr (CgemGeomSvc *geomSvc)
G4bool	ProcessHits (G4Step *, G4TouchableHistory *)
void	EndOfEvent (G4HCofThisEvent *)
void	BeginOfTruthEvent (const G4Event *)
void	EndOfTruthEvent (const G4Event *)
Public Member Functions inherited from BesSensitiveDetector
	BesSensitiveDetector (const G4String name)
virtual	~BesSensitiveDetector ()
virtual void	BeginOfTruthEvent (const G4Event *)
virtual void	EndOfTruthEvent (const G4Event *)
virtual void	BeginOfTrack (const G4Track *)
virtual void	EndOfTrack (const G4Track *)

Additional Inherited Members
Protected Member Functions inherited from BesSensitiveDetector
void	GetCurrentTrackIndex (G4int &trackIndex, G4int &g4TrackId) const

Detailed Description

Definition at line 61 of file BesCgemSD.hh.

Constructor & Destructor Documentation

BesCgemSD()

BesCgemSD::BesCgemSD

(

G4String

name

)

Definition at line 80 of file BesCgemSD.cc.

:BesSensitiveDetector(name)
{
    m_F_printTruth  = 0;
    m_F_outputTruth = 0;
    myCgemTruth=NULL;
 
    /* Define hit collection name handled by SD */
    collectionName.insert("BesCgemHitsCollection");
    collectionName.insert("BesCgemTruthCollection");
 
    // get info from JobOptionsSvc.
    ISvcLocator* iface = Gaudi::svcLocator();
    SmartIF<IJobOptionsSvc> jos(iface->service("JobOptionsSvc"));
    const std::vector<const Property*>* ppp = jos->getProperties("BesCgemSD");
 
    string outputFile;
    //outputFile += "/dat/outputMcTruth.txt";
 
    if (ppp) {
        //std::cout << "BesCgemSD size: " << ppp->size() << std::endl;
        for (std::vector<const Property*>::const_iterator it = ppp->begin();
                it != ppp->end(); ++it) {
            if ( (*it)->name() == "outputFile") {
                outputFile = (*it)->toString();
            }
        }
    }
    IMessageSvc* msgSvc;
    Gaudi::svcLocator() -> service("MessageSvc", msgSvc);
    MsgStream log(msgSvc, "BesCgemSD::BesCgemSD():BesSensitiveDetector()");
 
    /* Save McTruth information to compare with  Root RawData Information */
    if (m_F_outputTruth == 1)
    {
        log<< MSG::INFO << "CgemSim::BesCgemSD::Initialize(), McTruth Hit will be output to a txt file!" 
            << endreq
            << "outputFile Path: " << outputFile << endreq;
 
        m_outputTruth.open(outputFile.c_str(), ios::out);
 
        if (!m_outputTruth.good())
        {
            log<< MSG::ERROR << "CgemSim::BesCgemSD::Initialize(), Fail to open McTruth output file: "
                << outputFile << endreq;
            return ;
        }
    } /* End of 'if (m_outputTruth == 1)' */
}

~BesCgemSD()

BesCgemSD::~BesCgemSD

(

)

Definition at line 131 of file BesCgemSD.cc.

{
    m_outputTruth.close();
}

Member Function Documentation

BeginOfTruthEvent()

void BesCgemSD::BeginOfTruthEvent ( const G4Event *  f_event )

virtual

Reimplemented from BesSensitiveDetector.

Definition at line 155 of file BesCgemSD.cc.

{  
    m_collection_truth = new BesCgemHitsCollection(SensitiveDetectorName,collectionName[1]); 
    lvi_ID_hit=0;
 
}

EndOfEvent()

void BesCgemSD::EndOfEvent

(

G4HCofThisEvent *

HCE

)

Definition at line 655 of file BesCgemSD.cc.

{
    IMessageSvc* msgSvc;
    Gaudi::svcLocator() -> service("MessageSvc", msgSvc);
    MsgStream log(msgSvc, "BesCgemSD::EndOfEvent()");
 
    bool printFlag=false;
    m_map_truth.clear();
 
    /* Get the unique ID for this collection */
    /* GetCollectionID() is a heavy operation, shouldn't be called for every event. */
    /* It's available after SD is constructed and registed to G4SDManager. */
    /* It can't be invoked in the constructor. */
    static G4int HCID = -1;
    if (HCID < 0)
    { 
        HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]); 
    }
 
    /* Attach hit collection to G4HCofThisEvent object */
    if(printFlag)
    {
        log<< MSG::INFO << "-------------------------------------------------------------------------------" << endreq;
        log<< MSG::INFO << "CgemSim::BesCgemSD::EndOfEvent(),      Entries of hits  collection is " 
            << m_collection_hit->entries() << endreq;
    }
 
    HCE->AddHitsCollection(HCID, m_collection_hit); 
 
    if (verboseLevel > 0) 
    {
        m_collection_hit->PrintAllHits();
    } 
}

EndOfTruthEvent()

void BesCgemSD::EndOfTruthEvent ( const G4Event *  f_event )

virtual

Reimplemented from BesSensitiveDetector.

Definition at line 164 of file BesCgemSD.cc.

{  
    static G4int HLID=-1;
    if (HLID < 0)
    {
        HLID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[1]); 
    }
 
    G4HCofThisEvent *HCE = f_event->GetHCofThisEvent();
 
    IMessageSvc* msgSvc;
    Gaudi::svcLocator() -> service("MessageSvc", msgSvc);
    MsgStream log(msgSvc, "BesCgemSD::EndOfTruthEvent()");
    //log<< MSG::INFO  << "CgemSim::BesCgemSD::EndOfTruthEvent(), Entries of truth collection is " 
    //  << m_collection_truth->entries() << endreq;
 
    HCE->AddHitsCollection(HLID, m_collection_truth);
 
    /* Print and output McTruth information */
    if (m_F_printTruth==1 or m_F_outputTruth==1)
    {
        /* Print McTruth information */
        if (m_F_printTruth == 1)
        {
            log<< MSG::INFO << "-------------------------------------------------------------------------------" 
                << endreq;
            log<< MSG::INFO << "CgemSim::BesCgemSD::EndOfTruthEvent(), McTruth information: " << endreq;
            log<< MSG::INFO << left << "Track " << left << "Layer " << left << "PDG   " 
                << left << setw(12) << "R_Pre " 
                << left << setw(12) << "R_Post " 
                << left << setw(40) << "pos_Pre " 
                << left << setw(40) << "pos_Post " 
                << left << setw(12) << "P_Pre " 
                << left << setw(12) << "P_Post " 
                << left << setw(16) << "StepL "
                << left << setw(16) <<"Time " 
                << left << "Edep "
                << endreq;
        }
 
 
        for (G4int ii=0; ii<m_collection_truth->entries(); ii++)
        {
            BesCgemHit *lv_truth = (*m_collection_truth)[ii];
 
            /* Print McTruth information */
            if (m_F_printTruth == 1)
            {
                G4cout << left << setw( 6) << lv_truth-> GetTrackID             ()
                    << left << setw( 6) << lv_truth-> GetLayerID             ()
                    << left << setw( 6) << lv_truth-> GetPDGCode             ()
                    << left << setw(12) << lv_truth-> GetPositionOfPrePoint  ().rho()
                    << left << setw(12) << lv_truth-> GetPositionOfPostPoint ().rho()
                    << left <<             lv_truth-> GetPositionOfPrePoint()
                    << left <<             lv_truth-> GetPositionOfPostPoint()
                    << left << setw(12) << lv_truth-> GetMomentumOfPrePoint().mag()
                    << left << setw(12) << lv_truth-> GetMomentumOfPostPoint().mag()
                    << left << setw(16) << lv_truth-> GetStepLength          ()
                    << left << setw(16) << lv_truth-> GetGlobalTime          ()
                    << left << setw(16) << lv_truth-> GetTotalEnergyDeposit  ()
                    << G4endl;
            }
 
            /* Save McTruth information to compare with  Root RawData Information */
            if (m_F_outputTruth == 1)
            {
                m_outputTruth << left << setw( 5) << lv_truth->GetTrackID()                 
                    << left << setw( 3) << lv_truth->GetLayerID()                 
                    << left << setw(12) << lv_truth->GetTotalEnergyDeposit()      
                    << left << setw(12) << lv_truth->GetPositionOfPrePoint().x()  
                    << left << setw(12) << lv_truth->GetPositionOfPrePoint().y()  
                    << left << setw(12) << lv_truth->GetPositionOfPrePoint().z()  
                    << left << setw(12) << lv_truth->GetPositionOfPostPoint().x() 
                    << left << setw(12) << lv_truth->GetPositionOfPostPoint().y() 
                    << left << setw(12) << lv_truth->GetPositionOfPostPoint().z() 
                    << left << setw(12) << lv_truth->GetMomentumOfPrePoint().x()  
                    << left << setw(12) << lv_truth->GetMomentumOfPrePoint().y()  
                    << left << setw(12) << lv_truth->GetMomentumOfPrePoint().z()  
                    << left << setw(12) << lv_truth->GetMomentumOfPostPoint().x() 
                    << left << setw(12) << lv_truth->GetMomentumOfPostPoint().y() 
                    << left << setw(12) << lv_truth->GetMomentumOfPostPoint().z()
                    << G4endl;
            }
        } /* End of 'for (G4int ii=0; ii<m_collection_truth->entries(); ii++)' */
 
        //m_outputTruth.close();
 
    } /* End of 'if (m_F_printTruth==1 or m_F_outputTruth==1)' */
}

Initialize()

void BesCgemSD::Initialize

(

G4HCofThisEvent *

HCE

)

Definition at line 138 of file BesCgemSD.cc.

{
    /* Instantiate hit collection */
    m_collection_hit = new BesCgemHitsCollection(SensitiveDetectorName,collectionName[0]); 
 
    //Add by sunxh (lvi_ID_hit)
    lvi_ID_hit = 0;
    //end add
    //IMessageSvc* msgSvc;
    //Gaudi::svcLocator() -> service("MessageSvc", msgSvc);
    //MsgStream log(msgSvc, "BesCgemSD::Initialize()");
    //log<< MSG::INFO  << "|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||" << endreq;
    //log<< MSG::INFO  << "CgemSim::BesCgemSD::Initialize(), Begin !!!" << endreq;
}

ProcessHits()

G4bool BesCgemSD::ProcessHits	(	G4Step *	f_step,
		G4TouchableHistory *	f_touchable
	)

Definition at line 257 of file BesCgemSD.cc.

{
    bool printFlag=false;
 
    IMessageSvc* msgSvc;
    Gaudi::svcLocator() -> service("MessageSvc", msgSvc);
    MsgStream log(msgSvc, "BesCgemSD::ProcessHits()");
 
    /* Event -> Track -> Step */
    G4Track *lv_track = f_step->GetTrack();
    
    /* Check tracks */
    /* Skip neutral tracks */
    G4double charge = lv_track->GetDefinition()->GetPDGCharge();
    if (charge == 0.) 
    {
        if(printFlag) 
            log<< MSG::WARNING  << "CgemSim::BesCgemSD::ProcessHits(), Track is neutral, skip the track" << endreq;
        return false;
    }
    G4int   lvi_pdg_code    = lv_track->GetDefinition()->GetPDGEncoding();
 
    /* Global Time (time since the current event began ) */
    G4double      lvd_global_time = lv_track->GetGlobalTime();
    if (isnan(lvd_global_time))
    {
        if(printFlag) 
            log<< MSG::WARNING << "CgemSim::BesCgemSD::ProcessHits(), Gloabal time is nan" << endreq;
        return false;
    }
    /* MDC T window is 2 microsecond */
    if (lvd_global_time > 2000)
    {
        if(printFlag) 
            log<< MSG::WARNING << "CgemSim::BesCgemSD::ProcessHits(), Time window is >2000" << endreq;
        return false; 
    }
    
    /* Skip no energy deposit tracks */
    /*  if (lvd_E_deposit == 0.) 
        {
        log<< MSG::WARNING << "CgemSim::BesCgemSD::ProcessHits(), Energy deposition is 0,lose track=" 
        << primTrackIdx << endreq;
        return false;
        }*/
    
    /* Skip no length tracks */
    G4double      lvd_L_step      = f_step->GetStepLength();
    if (lvd_L_step == 0.)
    {
        if(printFlag) 
            log<< MSG::WARNING << "CgemSim::BesCgemSD::ProcessHits(), Step length is 0, skip track" << endreq;
        return false;
    }
    
    /* A G4Step object consists of two points */
    G4StepPoint *lv_step_pre  = f_step->GetPreStepPoint();
    G4StepPoint *lv_step_post = f_step->GetPostStepPoint();
    G4TrackStatus currentTrackStatus = lv_track->GetTrackStatus();
 
    /* Skip delta electron, added by sunxh */
    G4int lvi_RdtElectron = 0;
    G4String process_name="Generator"; 
    if (NULL != lv_track->GetCreatorProcess()){
        process_name = lv_track->GetCreatorProcess()->GetProcessName();
        if (process_name == "hIoni" || process_name == "ionIoni" || process_name == "conv" || process_name == "compt" || process_name == "muIoni" || process_name == "eIoni" || process_name == "phot" || process_name == "muMinusCaptureAtRest" )  lvi_RdtElectron = 1;
        //if(lvi_RdtElectron != 1&&lv_track->GetDefinition()->GetPDGEncoding()!=22){
        //               G4cout << "process_name = " <<  process_name << G4endl;
        //               G4cout << "PDG = " << lv_track->GetDefinition()->GetPDGEncoding()<<"\tTrack ID = "<<lv_track->GetTrackID() <<"\tParent ID = "<<lv_track->GetParentID() << G4endl;
        //}
        //cout<<"particle pdg="<<lvi_pdg_code<<" with CreatorProcess "<<process_name<<endl;
    }
    if(lvi_RdtElectron) return false;
 
    //else {
    //               G4cout << "process_name = null" << G4endl;
    //               G4cout << "PDG = " << lv_track->GetDefinition()->GetPDGEncoding()<<"\tTrack ID = "<<lv_track->GetTrackID() <<"\tParent ID = "<<lv_track->GetParentID() << G4endl;
    //}
    //   G4cout << "PDG = " << lv_track->GetDefinition()->GetPDGEncoding() << G4endl;
    //   else G4cout << "No process_name PDG = " << lv_track->GetDefinition()->GetPDGEncoding() << G4endl;
    //   if (abs(lv_track->GetDefinition()->GetPDGEncoding()) == 11 )  lvi_RdtElectron = 1;
    //   else{
    //      if (abs(lv_track->GetDefinition()->GetPDGEncoding()) == 11) {
    //       G4cout << "No process_name PDG = " << lv_track->GetDefinition()->GetPDGEncoding() << G4endl;
    //      }
    //   }
    //End Add
 
    /* Local variables used to get Hit information */
    /* Get track ID */
    G4int primTrackIdx   = 0;// trk index for the primary track, comments by llwang
    G4int G4TrackId   = 0;// G4 trk id for the current track
    G4int primG4TrackId = 0;// G4 trk id for the primary track
 
    /* In the following we use 'SetTrackID(primTrackIdx)' */
    //primTrackIdx = lv_track->GetTrackID();
    //GetCurrentTrackIndex(G4TrackId, primG4TrackId);
 
    //Add by sunxh
    G4TrackId = lv_track->GetTrackID();
    GetCurrentTrackIndex(primTrackIdx, primG4TrackId);
    //   G4cout<<"PDG = "<<lv_track->GetDefinition()->GetPDGEncoding()<<"\tprocess_name = "<<process_name<<"\tlvi_ID_truth = " << G4TrackId<<"\tprimTrackIdx = "<<primTrackIdx<<"\tlvi_ID_G4Track = "<<primG4TrackId<<G4endl;
    //end
 
    G4int         lvi_ID_layer    = 0;
    //G4int         lvi_pdg_code    = lv_track->GetDefinition()->GetPDGEncoding();
    G4int         lvi_ID_parent   = lv_track->GetParentID();// G4 trk id for the parent
    G4double      lvd_E_deposit   = f_step->GetTotalEnergyDeposit();
    /* Positions in the global coordinate system */
    G4ThreeVector lv3_XYZ_pre     = lv_step_pre ->GetPosition(); 
    G4ThreeVector lv3_XYZ_post    = lv_step_post->GetPosition();
 
    G4ThreeVector lv3_P_pre       = lv_step_pre ->GetMomentum();
    G4ThreeVector lv3_P_post      = lv_step_post->GetMomentum();
 
    if(printFlag) {
        log<< MSG::INFO  << "..............................................................................." << endreq;
        log<< MSG::INFO  << "CgemSim::BesCgemSD::ProcessHits(), Process a hit   at track=" 
            << primTrackIdx << endreq;
    }
 
    
 
    /* Check and Save Hits and MCTruth Hits */
    /* Get pre step point status and its volume */
    /* Volume depth: 0-Gap_D, 1-CgemLayer, 2-Cgem, 3-BES */
    G4StepStatus lv_status_pre = lv_step_pre  ->GetStepStatus();
    G4StepStatus lv_status_post = lv_step_post->GetStepStatus();
 
    /* Retrieve a touchable by creating a handle for it */
    const G4VTouchable *lv_TOUCHABLE = lv_step_pre  ->GetTouchable();
    /* Get the current volume(where the step has just gone through) and its logical name */
    G4String lvs_volume_name = lv_TOUCHABLE ->GetVolume(0)->GetLogicalVolume()->GetName();
 
    /* Save Hits in SD Geometry only! */
    //if ((lv_status_pre == fGeomBoundary) && (lvs_volume_name == "Gap_D_Logic"))
 
 
    //if (lvs_volume_name == "Gap_D_logic")
    if (lvs_volume_name == "Gap_D_logic0" or lvs_volume_name == "Gap_D_logic1" or lvs_volume_name == "Gap_D_logic2")
    {
        G4double trkLen   = lv_track->GetTrackLength();// track length at post point?
 
        lvi_ID_layer = lv_TOUCHABLE->GetVolume(1)->GetCopyNo();
 
        if(printFlag) {
 
            log<< MSG::INFO << "CgemSim::BesCgemSD::ProcessHits(), Track enters the SD volume: " 
                << lvs_volume_name << endreq;
 
            /* CgemLayer copy number */
 
            //    KeyType lv_key(primTrackIdx, lvi_ID_layer);
 
            /* Create Hits */
            log<< MSG::INFO << "CgemSim::BesCgemSD::ProcessHits(), Create a hit   at layer=" 
                << lvi_ID_layer << "  track=" << primTrackIdx << endreq;
        }
 
        BesCgemHit *l_new_hit = new BesCgemHit();    
 
        //Add by sunxh
        l_new_hit->SetRdtEl               ( lvi_RdtElectron );
        l_new_hit->SetHitID               ( lvi_ID_hit      );
        //end 
        l_new_hit->SetTrackID             ( primTrackIdx    );// primary track idx in McParticle Collection
        G4int isSecondary = (G4TrackId==primG4TrackId)? 0:1;
        //l_new_hit->SetCurrentTrackID      ( G4TrackId       );
        l_new_hit->SetIsSecondary         ( isSecondary     );
        l_new_hit->SetLayerID             ( lvi_ID_layer    );
        l_new_hit->SetPDGCode             ( lvi_pdg_code    );
        l_new_hit->SetParentID            ( lvi_ID_parent   );// G4 trk id for the parent 
        l_new_hit->SetGlobalTime          ( lvd_global_time );
        l_new_hit->SetTotalEnergyDeposit  ( lvd_E_deposit   );
        l_new_hit->SetStepLength          ( lvd_L_step      );
 
        // --> modified by wulh for alignment on Apr 5, 2019    begin:
        double x_pre_glo = lv3_XYZ_pre.x();
        double y_pre_glo = lv3_XYZ_pre.y();
        double z_pre_glo = lv3_XYZ_pre.z();
        double x_post_glo = lv3_XYZ_post.x();
        double y_post_glo = lv3_XYZ_post.y();
        double z_post_glo = lv3_XYZ_post.z();
        double r_pre  = sqrt(x_pre_glo*x_pre_glo+y_pre_glo*y_pre_glo);
        double r_post = sqrt(x_post_glo*x_post_glo+y_post_glo*y_post_glo);
 
        CgemGeoAlign* cgemAlign = m_cgemGeomSvc->getAlignPtr();
        HepPoint3D pos_pre(x_pre_glo, y_pre_glo, z_pre_glo);
        HepPoint3D pos_post(x_post_glo, y_post_glo, z_post_glo);
        StraightLine line_glo(pos_pre, pos_post);
        StraightLine line_loc = cgemAlign->StraightLineConversion(lvi_ID_layer, line_glo);
        HepPoint3D crossPoint_pre;
        HepPoint3D crossPoint_post;
 
        double rDrift_inn = m_cgemGeomSvc->getCgemLayer(lvi_ID_layer)->getInnerROfGapD();
        double rDrift_out = m_cgemGeomSvc->getCgemLayer(lvi_ID_layer)->getOuterROfGapD();
 
        double distLineGlo = fabs(line_glo.dr());
        if(fabs(distLineGlo) > rDrift_inn){
             crossPoint_pre = line_loc.xAtR(r_post, -1);
             crossPoint_post = line_loc.xAtR(r_post, 1);
        } else if(r_pre < r_post){
             crossPoint_pre = line_loc.xAtR(r_pre, 1);
             crossPoint_post = line_loc.xAtR(r_post, 1);
        } else{
             crossPoint_pre = line_loc.xAtR(r_pre, -1);
             crossPoint_post = line_loc.xAtR(r_post, -1);
        }
 
        G4ThreeVector lv3_XYZ_pre_align(crossPoint_pre.x(), crossPoint_pre.y(), crossPoint_pre.z());
        G4ThreeVector lv3_XYZ_post_align(crossPoint_post.x(), crossPoint_post.y(), crossPoint_post.z());
 
        l_new_hit->SetPositionOfPrePointAlign  ( lv3_XYZ_pre_align     );
        l_new_hit->SetPositionOfPostPointAlign ( lv3_XYZ_post_align    );
        // <-- modified by wulh for alignment on Apr 5, 2019 --- end.
 
 
        l_new_hit->SetPositionOfPrePoint  ( lv3_XYZ_pre     );
        l_new_hit->SetPositionOfPostPoint ( lv3_XYZ_post    );
        l_new_hit->SetMomentumOfPrePoint  ( lv3_P_pre       );
        l_new_hit->SetMomentumOfPostPoint ( lv3_P_post      );
        l_new_hit->SetCreatorProcess      ( process_name    );
        l_new_hit->SetFlightLengthPostPoint(trkLen);
 
        m_collection_hit->insert(l_new_hit);
 
        //Add by sunxh
        hit_ID_vector.push_back(lvi_ID_hit);
        lvi_ID_hit++;
        //end
        
        
        
        /* Create Truth Hit */
        /* MCTruth Hit only contains the point input and output the volume (not step by step) */
        /* Meaning of 'MCTruth hit' parameters is not the same as 'hit' */
        /* All steps in SD volume will be merged to 1 step  */
        /* Check over that the particle has just entered in the current volume */
        /* It is at the first step in the volume; the preStepPoint is at the boundary */
 
        KeyType lv_key(primTrackIdx, lvi_ID_layer);
 
        if (m_collection_truth)
        {
            int newVer=1;
 
            if(newVer==0) {
                m_map_truth_it = m_map_truth.find(lv_key);
                if (m_map_truth_it == m_map_truth.end())
                {
                    log<< MSG::INFO << "CgemSim::BesCgemSD::ProcessHits(), Create a Truth at layer=" 
                        << lvi_ID_layer << "  track=" << primTrackIdx << endreq;
 
                    BesCgemHit *l_new_truth = new BesCgemHit();
 
                    l_new_truth->SetTrackID             ( primTrackIdx    ); /* SAME */
                    l_new_truth->SetLayerID             ( lvi_ID_layer    ); /* SAME */
                    l_new_truth->SetPDGCode             ( lvi_pdg_code    ); /* SAME */
                    l_new_truth->SetParentID            ( lvi_ID_parent    ); /* SAME */
                    //l_new_truth->SetCurrentTrackID      ( G4TrackId        );
                    l_new_truth->SetGlobalTime          ( lvd_global_time ); /* DIFF */
                    l_new_truth->SetTotalEnergyDeposit  ( lvd_E_deposit   ); /* SAME */
                    l_new_truth->SetStepLength          ( lvd_L_step      ); /* SAME */
                    l_new_truth->SetPositionOfPrePoint  ( lv3_XYZ_pre     ); /* SAME */
                    l_new_truth->SetPositionOfPostPoint ( lv3_XYZ_post    ); /* SAME */
                    l_new_truth->SetPositionOfPrePointAlign(lv3_XYZ_pre_align); /* SAME */
                    l_new_truth->SetPositionOfPostPointAlign(lv3_XYZ_post_align); /* SAME */
                    l_new_truth->SetMomentumOfPrePoint  ( lv3_P_pre       ); /* SAME */
                    l_new_truth->SetMomentumOfPostPoint ( lv3_P_post      ); /* SAME */
 
                    //Add by sunxh
                    int tmp[2000];
                    for(int ii =0; ii < hit_ID_vector.size();ii++){
                        tmp[ii] = hit_ID_vector[ii];
                    }
                    l_new_truth->AddIdentifier(tmp,hit_ID_vector.size());
                    hit_ID_vector.clear();
                    //end
 
 
                    m_collection_truth->insert(l_new_truth);
 
                    m_map_truth[lv_key] = m_collection_truth->entries() - 1;
 
                    if (lv_step_pre->GetStepStatus() == fGeomBoundary)
                    {
                        log<< MSG::INFO << "CgemSim::BesCgemSD::ProcessHits(), Truth starts from boundry!!!"
                            << endreq;
                    }
                    else 
                    {
                        log<< MSG::INFO << "CgemSim::BesCgemSD::ProcessHits(), Truth doesn't start from boundry!!!"
                            << endreq;
                    }
                }//if (m_map_truth_it == m_map_truth.end())
                else
                {
                    if(printFlag) 
                        log<< MSG::INFO << "CgemSim::BesCgemSD::ProcessHits(), Add a step to Truth at layer=" 
                            << lvi_ID_layer << "  track=" << primTrackIdx << endreq;
 
                    G4int    lvi_ID_truth     = (*m_map_truth_it).second;
                    G4int    primTrackIdx_pre = (*m_collection_truth)[lvi_ID_truth]->GetTrackID();
                    G4double lvd_E_pre        = (*m_collection_truth)[lvi_ID_truth]->GetTotalEnergyDeposit();
                    G4double lvd_L_step_pre   = (*m_collection_truth)[lvi_ID_truth]->GetStepLength();
 
                    (*m_collection_truth)[lvi_ID_truth]->SetTotalEnergyDeposit(lvd_E_deposit + lvd_E_pre);
                    (*m_collection_truth)[lvi_ID_truth]->SetStepLength(lvd_L_step + lvd_L_step_pre);
                    (*m_collection_truth)[lvi_ID_truth]->SetPositionOfPostPoint(lv3_XYZ_post);
                    (*m_collection_truth)[lvi_ID_truth]->SetPositionOfPostPointAlign(lv3_XYZ_post_align);
                    (*m_collection_truth)[lvi_ID_truth]->SetMomentumOfPostPoint(lv3_P_post);
                }
 
                /* Check over that the particle is leaving the current volume */
                /* It is at the last step in the volume; the postStepPoint is at the boundary */
                if(printFlag)
                    if (lv_step_post->GetStepStatus() == fGeomBoundary)
                    {
                        log<< MSG::INFO << "CgemSim::BesCgemSD::ProcessHits(), Finish a Truth at layer=" 
                            << lvi_ID_layer << "  track=" << primTrackIdx << endreq;
                    }
            }// old version
            else
            {
                if(myCgemTruth==NULL)
                {
                    myCgemTruth= new BesCgemHit();
                    if(printFlag)
                        cout<<"creat a new myCgemTruth, prePos is Boundary: "<<(lv_status_pre==fGeomBoundary)<<" trkID="<<primTrackIdx<<", LayerID="<<lvi_ID_layer<<endl;
                    //         myCgemTruth->SetTrackID             ( G4TrackId); /* SAME */
                    myCgemTruth->SetTrackID             ( primTrackIdx); /* SAME */
                    myCgemTruth->SetLayerID             ( lvi_ID_layer    ); /* SAME */
                    myCgemTruth->SetPDGCode             ( lvi_pdg_code    ); /* SAME */
                    myCgemTruth->SetParentID            ( lvi_ID_parent    ); /* SAME */
                    myCgemTruth->SetGlobalTime          ( lvd_global_time ); /* DIFF */
                    myCgemTruth->SetTotalEnergyDeposit  ( lvd_E_deposit   ); /* SAME */
                    myCgemTruth->SetStepLength          ( lvd_L_step      ); /* SAME */
                    myCgemTruth->SetPositionOfPrePoint  ( lv3_XYZ_pre     ); /* SAME */
                    myCgemTruth->SetPositionOfPostPoint ( lv3_XYZ_post    ); /* SAME */
                    myCgemTruth->SetPositionOfPrePointAlign(lv3_XYZ_pre_align); /* SAME */
                    myCgemTruth->SetPositionOfPostPointAlign(lv3_XYZ_post_align); /* SAME */
                    myCgemTruth->SetMomentumOfPrePoint  ( lv3_P_pre       ); /* SAME */
                    myCgemTruth->SetMomentumOfPostPoint ( lv3_P_post      ); /* SAME */
                    myCgemTruth->SetCreatorProcess      ( process_name    );
                    myCgemTruth->SetRdtEl               ( lvi_RdtElectron );
                    myCgemTruth->SetIsSecondary         (isSecondary      );
                    myCgemTruth->SetFlightLengthPostPoint(trkLen);
                    //double dphi=lv3_XYZ_pre.phi()-lv3_XYZ_post.phi();
                    //double dX=dphi*87.5;
                    //if(lvi_ID_layer==0) cout<<"dX="<<dX<<" , nXStrips="<<dX/0.65<<endl;
                }
                else
                {
                    G4double lvd_E_pre        = myCgemTruth->GetTotalEnergyDeposit();
                    G4double lvd_L_step_pre   = myCgemTruth->GetStepLength();
                    G4ThreeVector pos_pre_last  = myCgemTruth->GetPositionOfPrePoint();
                    G4ThreeVector pos_post_last = myCgemTruth->GetPositionOfPostPoint();
                    //cout<<"add CGEM truth, last hit: "<<pos_pre_last<<" to "<<pos_post_last<<endl;
                    //cout<<"                new  hit: "<<lv3_XYZ_pre <<" to "<<lv3_XYZ_post<<endl;
                    myCgemTruth->SetTotalEnergyDeposit(lvd_E_deposit + lvd_E_pre);
                    myCgemTruth->SetStepLength(lvd_L_step + lvd_L_step_pre);
                    myCgemTruth->SetPositionOfPostPoint(lv3_XYZ_post);
                    myCgemTruth->SetPositionOfPostPointAlign(lv3_XYZ_post_align);
                    myCgemTruth->SetMomentumOfPostPoint(lv3_P_post);
                    myCgemTruth->SetFlightLengthPostPoint(trkLen);
                    if(printFlag)
                        cout<<"update myCgemTruth: prePos is Boundary: "<<(lv_status_pre==fGeomBoundary)<<", postPos is Boundary:"<<(lv_status_post==fGeomBoundary)<<", trkID="<<primTrackIdx<<", LayerID="<<lvi_ID_layer<<endl;
                }
                if( myCgemTruth && (currentTrackStatus==fStopAndKill || lv_status_post==fGeomBoundary) )
                {
                    //Add by sunxh
                    int tmp[2000];
                    for(int ii =0; ii < hit_ID_vector.size();ii++){
                        tmp[ii] = hit_ID_vector[ii];
                    }
                    myCgemTruth->AddIdentifier(tmp,hit_ID_vector.size());
                    //TArrayI hit_tmp=myCgemTruth->GetIdentifier();
                    //end
 
                    m_collection_truth->insert(myCgemTruth);
                    myCgemTruth=NULL;
 
                    //Add by sunxh
                    hit_ID_vector.clear();
                    //end
 
                    if(printFlag)
                        cout<<"save a myCgemTruth: prePos is Boundary: "<<(lv_status_pre==fGeomBoundary)<<", postPos is Boundary:"<<(lv_status_post==fGeomBoundary)<<", trkID="<<primTrackIdx<<", LayerID="<<lvi_ID_layer<<endl;
                }
            }//if(newVer==0)
        } /* End of 'if (m_collection_truth)' */
    } /* End of 'if (lvs_volume_name == "Gap_D_Logic")' */
 
    return true;
}

setGeomPtr()

void BesCgemSD::setGeomPtr ( CgemGeomSvc *  geomSvc )

inline

Definition at line 68 of file BesCgemSD.hh.

{m_cgemGeomSvc=geomSvc;}

Referenced by BesCgemConstruction::Construct().

---

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

• BesCgemSD.hh
• BesCgemSD.cc