B03Run Class Reference

#include <Doxymodules_biasing.h>

Inheritance diagram for B03Run:

Public Member Functions
	B03Run (const std::vector< G4String > mfdName)
virtual	~B03Run ()
virtual void	RecordEvent (const G4Event *)
G4int	GetNumberOfHitsMap () const
G4THitsMap< G4double > *	GetHitsMap (G4int i)
G4THitsMap< G4double > *	GetHitsMap (const G4String &detName, const G4String &colName)
G4THitsMap< G4double > *	GetHitsMap (const G4String &fullName)
void	DumpAllScorer ()
virtual void	Merge (const G4Run *)

Private Attributes
std::vector< G4String >	fCollName
std::vector< G4int >	fCollID
std::vector< G4THitsMap< G4double > * >	fRunMap

Detailed Description

Definition at line 56 of file Doxymodules_biasing.h.

Constructor & Destructor Documentation

B03Run()

B03Run::B03Run

(

const std::vector< G4String >

mfdName

)

Definition at line 69 of file B03Run.cc.

                                               : G4Run()
{
  G4SDManager* SDman = G4SDManager::GetSDMpointer();
  //=================================================
  //  Initalize RunMaps for accumulation.
  //  Get CollectionIDs for HitCollections.
  //=================================================
  G4int Nmfd = mfdName.size();
  for ( G4int idet = 0; idet < Nmfd ; idet++){  // Loop for all MFD.
    G4String detName = mfdName[idet];
    //--- Seek and Obtain MFD objects from SDmanager.
    G4MultiFunctionalDetector* mfd =
      (G4MultiFunctionalDetector*)(SDman->FindSensitiveDetector(detName));
    //
    if ( mfd ){
        //--- Loop over the registered primitive scorers.
        for (G4int icol = 0; icol < mfd->GetNumberOfPrimitives(); icol++){
            // Get Primitive Scorer object.
            G4VPrimitiveScorer* scorer=mfd->GetPrimitive(icol);
      // collection name and collectionID for HitsCollection,
      // where type of HitsCollection is G4THitsMap in case of primitive scorer.
      // The collection name is given by <MFD name>/<Primitive Scorer name>.
            G4String collectionName = scorer->GetName();
            G4String fullCollectionName = detName+"/"+collectionName;
            G4int    collectionID = SDman->GetCollectionID(fullCollectionName);
            //
            if ( collectionID >= 0 ){
                G4cout << "++ "<<fullCollectionName<< " id " 
                       << collectionID << G4endl;
          // Store obtained HitsCollection information into data members.
          // And, creates new G4THitsMap for accumulating quantities during RUN.
                fCollName.push_back(fullCollectionName);
                fCollID.push_back(collectionID);
                fRunMap.push_back(new G4THitsMap<G4double>(detName
                                                          ,collectionName));
            }else{
                G4cout << "** collection " << fullCollectionName 
                       << " not found. "<<G4endl;
            }
        }
    }
  }
}

~B03Run()

B03Run::~B03Run ( )

virtual

Definition at line 117 of file B03Run.cc.

{
  //--- Clear HitsMap for RUN
  G4int Nmap = fRunMap.size();
  for ( G4int i = 0; i < Nmap; i++){
    if(fRunMap[i] ) fRunMap[i]->clear();
  }
  fCollName.clear();
  fCollID.clear();
  fRunMap.clear();
}

Member Function Documentation

RecordEvent()

void B03Run::RecordEvent ( const G4Event *  aEvent )

virtual

Definition at line 134 of file B03Run.cc.

{
  numberOfEvent++;  // This is an original line.
 
  //=============================
  // HitsCollection of This Event
  //============================
  G4HCofThisEvent* HCE = aEvent->GetHCofThisEvent();
  if (!HCE) return;
  //=======================================================
  // Sum up HitsMap of this Event  into HitsMap of this RUN
  //=======================================================
  G4int Ncol = fCollID.size();
  for ( G4int i = 0; i < Ncol ; i++ ){  // Loop over HitsCollection
    G4THitsMap<G4double>* EvtMap=0;
    if ( fCollID[i] >= 0 ){           // Collection is attached to HCE
      EvtMap = (G4THitsMap<G4double>*)(HCE->GetHC(fCollID[i]));
    }else{
      G4cout <<" Error EvtMap Not Found "<< i << G4endl;
    }
    if ( EvtMap )  {
      //=== Sum up HitsMap of this event to HitsMap of RUN.===
      *fRunMap[i] += *EvtMap;
      //======================================================
    }
  }
 
  
}

GetNumberOfHitsMap()

G4int B03Run::GetNumberOfHitsMap ( ) const

inline

Definition at line 65 of file B03Run.hh.

{return fRunMap.size();}

GetHitsMap() [1/3]

G4THitsMap< G4double > * B03Run::GetHitsMap ( G4int  i )

inline

Definition at line 69 of file B03Run.hh.

{return fRunMap[i];}

GetHitsMap() [2/3]

G4THitsMap< G4double > * B03Run::GetHitsMap	(	const G4String &	detName,
		const G4String &	colName
	)

Definition at line 170 of file B03Run.cc.

                                                                 {
    G4String fullName = detName+"/"+colName;
    G4cout << " getting hits map " << fullName << G4endl;
    return GetHitsMap(fullName);
}

GetHitsMap() [3/3]

G4THitsMap< G4double > * B03Run::GetHitsMap

(

const G4String &

fullName

)

Definition at line 182 of file B03Run.cc.

                                                                {
    G4cout << " getting hits map " << fullName << G4endl;
    G4int Ncol = fCollName.size();
    for ( G4int i = 0; i < Ncol; i++){
        if ( fCollName[i] == fullName ){
            return fRunMap[i];
        }
    }
    return NULL;
}

DumpAllScorer()

void B03Run::DumpAllScorer

(

)

Definition at line 197 of file B03Run.cc.

                          {
 
  // - Number of HitsMap in this RUN.
  G4int n = GetNumberOfHitsMap();
  // - GetHitsMap and dump values.
  for ( G4int i = 0; i < n ; i++ ){
    G4THitsMap<G4double>* RunMap =GetHitsMap(i);
    if ( RunMap ) {
      G4cout << " PrimitiveScorer RUN " 
             << RunMap->GetSDname() <<","<< RunMap->GetName() << G4endl;
      G4cout << " Number of entries " << RunMap->entries() << G4endl;
      std::map<G4int,G4double*>::iterator itr = RunMap->GetMap()->begin();
      for(; itr != RunMap->GetMap()->end(); itr++) {
        G4cout << "  copy no.: " << itr->first
               << "  Run Value : " << *(itr->second) 
               << G4endl;
      }
    }
  }
}

Merge()

void B03Run::Merge ( const G4Run *  aRun )

virtual

Definition at line 220 of file B03Run.cc.

{
  const B03Run * localRun = static_cast<const B03Run *>(aRun);
  //=======================================================
  // Merge HitsMap of working threads
  //=======================================================
  G4int nCol = localRun->fCollID.size();
  for ( G4int i = 0; i < nCol ; i++ ){  // Loop over HitsCollection
    if ( localRun->fCollID[i] >= 0 ){
      *fRunMap[i] += *localRun->fRunMap[i];
    }
  }
  
  G4Run::Merge(aRun);
} 

Member Data Documentation

fCollName

std::vector<G4String> B03Run::fCollName

private

Definition at line 80 of file B03Run.hh.

fCollID

std::vector<G4int> B03Run::fCollID

private

Definition at line 81 of file B03Run.hh.

fRunMap

std::vector<G4THitsMap<G4double>*> B03Run::fRunMap

private

Definition at line 82 of file B03Run.hh.

---

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

• .doxygen/Doxymodules_biasing.h
• extended/biasing/B03/include/B03Run.hh
• extended/biasing/B03/src/B03Run.cc