Ignore:
Timestamp:
Nov 5, 2010, 4:08:39 PM (14 years ago)
Author:
garnier
Message:

update ti head

Location:
trunk/examples/extended/parameterisations/gflash/src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/examples/extended/parameterisations/gflash/src/ExGflashPhysicsList.cc

    r807 r1342  
    3535#include "G4Material.hh"
    3636#include "G4MaterialTable.hh"
     37#include "G4Region.hh"
     38#include "G4RegionStore.hh"
    3739#include "G4ios.hh"
    3840#include <iomanip>   
     
    135137#include "G4PhotoElectricEffect.hh"
    136138
    137 #include "G4MultipleScattering.hh"
     139#include "G4eMultipleScattering.hh"
     140#include "G4MuMultipleScattering.hh"
     141#include "G4hMultipleScattering.hh"
    138142
    139143#include "G4eIonisation.hh"
     
    179183    //electron
    180184      // Construct processes for electron
    181       G4VProcess* theeminusMultipleScattering = new G4MultipleScattering();
     185      G4VProcess* theeminusMultipleScattering = new G4eMultipleScattering();
    182186      G4VProcess* theeminusIonisation = new G4eIonisation();
    183187      G4VProcess* theeminusBremsstrahlung = new G4eBremsstrahlung();
     
    208212    //positron
    209213      // Construct processes for positron
    210       G4VProcess* theeplusMultipleScattering = new G4MultipleScattering();
     214      G4VProcess* theeplusMultipleScattering = new G4eMultipleScattering();
    211215      G4VProcess* theeplusIonisation = new G4eIonisation();
    212216      G4VProcess* theeplusBremsstrahlung = new G4eBremsstrahlung();
     
    236240    //muon 
    237241     // Construct processes for muon+
    238      G4VProcess* aMultipleScattering = new G4MultipleScattering();
     242     G4VProcess* aMultipleScattering = new G4MuMultipleScattering();
    239243     G4VProcess* aBremsstrahlung = new G4MuBremsstrahlung();
    240244     G4VProcess* aPairProduction = new G4MuPairProduction();
     
    258262              (particle->GetParticleName() != "chargedgeantino")) {
    259263     // all others charged particles except geantino
    260      G4VProcess* aMultipleScattering = new G4MultipleScattering();
     264     G4VProcess* aMultipleScattering = new G4hMultipleScattering();
    261265     G4VProcess* anIonisation = new G4hIonisation();
    262266     // add processes
     
    335339 
    336340  DumpCutValuesTable();
    337 }
    338 
    339 
    340 
    341 
    342 
    343 
    344 
    345 
    346 
     341// set cuts for region crystals with default Cuts
     342  G4Region* region = G4RegionStore::GetInstance()->GetRegion("crystals");
     343  region->SetProductionCuts(
     344          G4ProductionCutsTable::GetProductionCutsTable()->GetDefaultProductionCuts());
     345}
     346
     347
     348
     349
     350
     351
     352
     353
     354
  • trunk/examples/extended/parameterisations/gflash/src/ExGflashSensitiveDetector.cc

    r807 r1342  
    5858void ExGflashSensitiveDetector::EndOfEvent(G4HCofThisEvent*HCE)
    5959{
    60         if (HCE);       
     60        if (HCE){;}     
    6161        // No Hadronic
    6262        //    cout<<"IEndOfEvent  " << HCID << endl;
     
    7878        caloHit->SetPos(aStep->GetPreStepPoint()->GetPosition());
    7979        caloHitsCollection->insert(caloHit);
    80         if (ROhist);
     80        if (ROhist){;}
    8181        G4VPhysicalVolume* physVol = theTouchable->GetVolume();
    8282        G4int crystalnum=0;
     
    102102        caloHit->SetPos(aSpot->GetEnergySpot()->GetPosition());
    103103        caloHitsCollection->insert(caloHit);
    104         if (ROhist);
     104        if (ROhist){;}
    105105        //cout <<pCurrentVolume->GetName()   << endl;
    106106        G4int crystalnum=0;
Note: See TracChangeset for help on using the changeset viewer.