---

GflashHadronShowerModel Class Reference

#include <SimG4Core/GFlash/interface/GflashHadronShowerModel.h>

List of all members.

Public Member Functions
void	DoIt (const G4FastTrack &, G4FastStep &)
	GflashHadronShowerModel (G4String modelName, G4Region *envelope)
G4bool	IsApplicable (const G4ParticleDefinition &)
G4bool	ModelTrigger (const G4FastTrack &)
	~GflashHadronShowerModel ()
Private Member Functions
G4bool	excludeDetectorRegion (const G4FastTrack &fastTrack)
G4bool	isFirstInelasticInteraction (const G4FastTrack &fastTrack)
Private Attributes
G4Navigator *	theGflashNavigator
G4Step *	theGflashStep
G4TouchableHandle	theGflashTouchableHandle
GflashHistogram *	theHisto
GflashHadronShowerProfile *	theProfile

---

Detailed Description

Definition at line 12 of file GflashHadronShowerModel.h.

---

Constructor & Destructor Documentation

GflashHadronShowerModel::GflashHadronShowerModel	(	G4String	modelName,
		G4Region *	envelope
	)

Definition at line 22 of file GflashHadronShowerModel.cc.

References GflashHistogram::instance(), theGflashNavigator, theGflashStep, theGflashTouchableHandle, theHisto, and theProfile.

  : G4VFastSimulationModel(modelName, envelope)
{
  theProfile = new GflashHadronShowerProfile(envelope);
  theHisto = GflashHistogram::instance();

  theGflashStep = new G4Step();
  theGflashNavigator = new G4Navigator();
  theGflashTouchableHandle = new G4TouchableHistory();

  theGflashNavigator->SetWorldVolume(G4TransportationManager::GetTransportationManager()->GetNavigatorForTracking()->GetWorldVolume());

}

GflashHadronShowerModel::~GflashHadronShowerModel

(

)

Definition at line 36 of file GflashHadronShowerModel.cc.

References theGflashStep, and theProfile.

{
  if(theProfile) delete theProfile;
  if(theGflashStep) delete theGflashStep;
}

---

Member Function Documentation

void GflashHadronShowerModel::DoIt	(	const G4FastTrack &	fastTrack,
		G4FastStep &	fastStep
	)

Definition at line 80 of file GflashHadronShowerModel.cc.

References GflashHadronShowerProfile::clearSpotList(), GflashHadronShowerProfile::getEnergySpotList(), GflashHadronShowerProfile::hadronicParameterization(), theGflashNavigator, theGflashStep, theGflashTouchableHandle, and theProfile.

{
  // Kill the parameterised particle:

  fastStep.ProposeTotalEnergyDeposited(fastTrack.GetPrimaryTrack()->GetKineticEnergy());

  // Parameterize shower shape and get resultant energy spots
  theProfile->clearSpotList();
  theProfile->hadronicParameterization(fastTrack);

  std::vector<GflashEnergySpot>& energySpotList = theProfile->getEnergySpotList();

  // Make hits
  G4double timeGlobal = fastTrack.GetPrimaryTrack()->GetStep()->GetPreStepPoint()->GetGlobalTime();
  
  std::vector<GflashEnergySpot>::const_iterator spotIter    = energySpotList.begin();
  std::vector<GflashEnergySpot>::const_iterator spotIterEnd = energySpotList.end();
  
   for( ; spotIter != spotIterEnd; spotIter++){

    // to make a different time for each fake step. (+1.0 is arbitrary)
    timeGlobal += 0.0001*nanosecond;

    // fill equivalent changes to a (fake) step associated with a spot 

    theGflashStep->SetTrack(const_cast<G4Track*>(fastTrack.GetPrimaryTrack()));
    theGflashStep->GetPostStepPoint()->SetGlobalTime(timeGlobal);
    theGflashStep->GetPreStepPoint()->SetPosition(spotIter->getPosition());
    theGflashStep->GetPostStepPoint()->SetPosition(spotIter->getPosition());
    theGflashStep->GetPostStepPoint()->SetProcessDefinedStep(const_cast<G4VProcess*> (fastTrack.GetPrimaryTrack()->GetStep()->GetPostStepPoint()->GetProcessDefinedStep()));

    //put touchable for each energy spot
    theGflashNavigator->LocateGlobalPointAndUpdateTouchable(spotIter->getPosition(),theGflashTouchableHandle(), false);
    theGflashStep->GetPreStepPoint()->SetTouchableHandle(theGflashTouchableHandle);
    theGflashStep->SetTotalEnergyDeposit(spotIter->getEnergy());
    
    // Send G4Step information to Hit/Dig if the volume is sensitive
    // Copied from G4SteppingManager.cc

    G4VPhysicalVolume* aCurrentVolume = theGflashStep->GetPreStepPoint()->GetPhysicalVolume();

    if( aCurrentVolume != 0 ) {
      theGflashStep->GetPreStepPoint()->SetSensitiveDetector(aCurrentVolume->GetLogicalVolume()->GetSensitiveDetector());
      G4VSensitiveDetector* aSensitive = theGflashStep->GetPreStepPoint()->GetSensitiveDetector();
      
      if( aSensitive != 0 ) {
        aSensitive->Hit(theGflashStep);
      }
    }
  }

  fastStep.KillPrimaryTrack();
  fastStep.ProposePrimaryTrackPathLength(0.0);

}

G4bool GflashHadronShowerModel::excludeDetectorRegion

(

const G4FastTrack &

fastTrack

)

[private]

Definition at line 175 of file GflashHadronShowerModel.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), and eta.

Referenced by ModelTrigger().

{
  G4bool isExcluded=false;
  
  //exclude regions where geometry are complicated 
  G4double eta =   fastTrack.GetPrimaryTrack()->GetMomentum().pseudoRapidity() ;
  if(fabs(eta) > 1.30 && fabs(eta) < 1.57) {
    std::cout << "excluding region of eta = " << eta << std::endl;
    return true;  
  }
  //exclude the region where the shower starting point is too close to the end of 
  //the hadronic envelopes (may need to be optimized further!)

  /*
  Gflash::CalorimeterNumber kColor = theProfile->getCalorimeterNumber(fastTrack);

  //@@@ need a proper scale
  const G4double minDistantToOut = 10.;
  

  if(kColor == kHB || kColor == kHE) {
    G4double distOut = fastTrack.GetEnvelopeSolid()->
      DistanceToOut(fastTrack.GetPrimaryTrackLocalPosition(),
                    fastTrack.GetPrimaryTrackLocalDirection());

    if (distOut < minDistantToOut ) {
      std::cout << "excluding region for dsitOut = " << distOut << std::endl;
      isExcluded = true;
    }
  }
  */

  return isExcluded;
}

G4bool GflashHadronShowerModel::IsApplicable

(

const G4ParticleDefinition &

particleType

)

Definition at line 42 of file GflashHadronShowerModel.cc.

{
  return 
    &particleType == G4PionMinus::PionMinusDefinition() ||
    &particleType == G4PionPlus::PionPlusDefinition();
}

G4bool GflashHadronShowerModel::isFirstInelasticInteraction

(

const G4FastTrack &

fastTrack

)

[private]

Definition at line 136 of file GflashHadronShowerModel.cc.

References GflashHistogram::deltaStep, relval_parameters_module::energy, GflashHistogram::energyLoss, GflashHistogram::getStoreFlag(), cmsDriverOptions::isFirst, GflashHistogram::kineticEnergy, GflashHistogram::postStepPosition, GflashHistogram::preStepPosition, and theHisto.

Referenced by ModelTrigger().

{
  G4bool isFirst=false;

  G4StepPoint* preStep = fastTrack.GetPrimaryTrack()->GetStep()->GetPreStepPoint();
  G4StepPoint* postStep = fastTrack.GetPrimaryTrack()->GetStep()->GetPostStepPoint();

  G4String procName = postStep->GetProcessDefinedStep()->GetProcessName();  
  G4ParticleDefinition* particleType = fastTrack.GetPrimaryTrack()->GetDefinition();

  //@@@ this part is still temporary and the cut for the variable ratio should be optimized later

  if((particleType == G4PionPlus::PionPlusDefinition() && procName == "WrappedPionPlusInelastic") || 
     (particleType == G4PionMinus::PionMinusDefinition() && procName == "WrappedPionMinusInelastic")) {

    G4double energy = fastTrack.GetPrimaryTrack()->GetKineticEnergy();

    G4double ratio = 0.0;
    if (energy > 0) {
      ratio = fabs(fastTrack.GetPrimaryTrack()->GetStep()->GetDeltaEnergy()/energy);
    }

    if(ratio > 0.1) isFirst=true;

    //Fill debugging histograms and check information on secondaries -
    //remove after final implimentation

    if(theHisto->getStoreFlag()) {
      theHisto->preStepPosition->Fill(preStep->GetPosition().getRho()/cm);
      theHisto->postStepPosition->Fill(postStep->GetPosition().getRho()/cm);
      theHisto->deltaStep->Fill((postStep->GetPosition() - preStep->GetPosition()).getRho()/cm);
      theHisto->kineticEnergy->Fill(energy);
      theHisto->energyLoss->Fill(fabs(fastTrack.GetPrimaryTrack()->GetStep()->GetDeltaEnergy()/GeV));
    }

 }
  return isFirst;
} 

G4bool GflashHadronShowerModel::ModelTrigger

(

const G4FastTrack &

fastTrack

)

Definition at line 49 of file GflashHadronShowerModel.cc.

References excludeDetectorRegion(), and isFirstInelasticInteraction().

{
  // ModelTrigger returns false for Gflash Hadronic Shower Model if it is not
  // tested from the corresponding wrapper process, GflashHadronWrapperProcess. 
  // Temporarily the track status is set to fPostponeToNextEvent at the wrapper
  // process before ModelTrigger is really tested for the second time through 
  // PostStepGPIL of enviaged parameterization processes.  The better 
  // implmentation may be using via G4VUserTrackInformation of each track, which
  // requires to modify a geant source code of stepping (G4SteppingManager2)

  G4bool trigger = false;

  // mininum energy cutoff to parameterize
  if (fastTrack.GetPrimaryTrack()->GetKineticEnergy() < 1.0*GeV) return trigger;

  // check whether this is called from the normal GPIL or the wrapper process GPIL
  if(fastTrack.GetPrimaryTrack()->GetTrackStatus() == fPostponeToNextEvent ) {

    // Shower pameterization start at the first inelastic interaction point
    G4bool isInelastic  = isFirstInelasticInteraction(fastTrack);
    
    // Other conditions
    if(isInelastic) {
      trigger = (!excludeDetectorRegion(fastTrack));
    }
  }

  return trigger;

}

---

Member Data Documentation

G4Navigator* GflashHadronShowerModel::theGflashNavigator [private]

Definition at line 37 of file GflashHadronShowerModel.h.

Referenced by DoIt(), and GflashHadronShowerModel().

G4Step* GflashHadronShowerModel::theGflashStep [private]

Definition at line 36 of file GflashHadronShowerModel.h.

Referenced by DoIt(), GflashHadronShowerModel(), and ~GflashHadronShowerModel().

G4TouchableHandle GflashHadronShowerModel::theGflashTouchableHandle [private]

Definition at line 38 of file GflashHadronShowerModel.h.

Referenced by DoIt(), and GflashHadronShowerModel().

GflashHistogram* GflashHadronShowerModel::theHisto [private]

Definition at line 41 of file GflashHadronShowerModel.h.

Referenced by GflashHadronShowerModel(), and isFirstInelasticInteraction().

GflashHadronShowerProfile* GflashHadronShowerModel::theProfile [private]

Definition at line 35 of file GflashHadronShowerModel.h.

Referenced by DoIt(), GflashHadronShowerModel(), and ~GflashHadronShowerModel().

---

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

• SimG4Core/GFlash/interface/GflashHadronShowerModel.h
• SimG4Core/GFlash/src/GflashHadronShowerModel.cc

---