Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // >> 26 // $Id: $ 26 // 27 // 27 /// \file Par02FastSimModelHCal.cc 28 /// \file Par02FastSimModelHCal.cc 28 /// \brief Implementation of the Par02FastSimM 29 /// \brief Implementation of the Par02FastSimModelHCal class 29 30 30 #include "Par02FastSimModelHCal.hh" 31 #include "Par02FastSimModelHCal.hh" 31 << 32 #include "Par02EventInformation.hh" 32 #include "Par02EventInformation.hh" 33 #include "Par02Output.hh" << 34 #include "Par02PrimaryParticleInformation.hh" 33 #include "Par02PrimaryParticleInformation.hh" 35 #include "Par02Smearer.hh" 34 #include "Par02Smearer.hh" >> 35 #include "Par02Output.hh" 36 36 37 #include "G4AnalysisManager.hh" << 37 #include "G4Track.hh" 38 #include "G4Event.hh" 38 #include "G4Event.hh" 39 #include "G4RunManager.hh" 39 #include "G4RunManager.hh" 40 #include "G4SystemOfUnits.hh" << 40 #include "g4root.hh" 41 #include "G4Track.hh" << 41 42 #include "Randomize.hh" 42 #include "Randomize.hh" >> 43 #include "G4SystemOfUnits.hh" 43 44 44 //....oooOO0OOooo........oooOO0OOooo........oo 45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 46 46 Par02FastSimModelHCal::Par02FastSimModelHCal(G << 47 Par02FastSimModelHCal::Par02FastSimModelHCal( G4String aModelName, 47 P << 48 G4Region* aEnvelope, Par02DetectorParametrisation::Parametrisation aType ) : 48 : G4VFastSimulationModel(aModelName, aEnvelo << 49 G4VFastSimulationModel( aModelName, aEnvelope ), fCalculateParametrisation(), 49 fCalculateParametrisation(), << 50 fParametrisation( aType ) {} 50 fParametrisation(aType) << 51 {} << 52 51 53 //....oooOO0OOooo........oooOO0OOooo........oo 52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 54 53 55 Par02FastSimModelHCal::Par02FastSimModelHCal(G << 54 Par02FastSimModelHCal::Par02FastSimModelHCal( G4String aModelName, 56 : G4VFastSimulationModel(aModelName, aEnvelo << 55 G4Region* aEnvelope ) : 57 fCalculateParametrisation(), << 56 G4VFastSimulationModel( aModelName, aEnvelope ), fCalculateParametrisation(), 58 fParametrisation(Par02DetectorParametrisat << 57 fParametrisation( Par02DetectorParametrisation::eCMS ) {} 59 {} << 60 58 61 //....oooOO0OOooo........oooOO0OOooo........oo 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 62 60 63 Par02FastSimModelHCal::Par02FastSimModelHCal(G << 61 Par02FastSimModelHCal::Par02FastSimModelHCal( G4String aModelName ) : 64 : G4VFastSimulationModel(aModelName), << 62 G4VFastSimulationModel( aModelName ), fCalculateParametrisation(), 65 fCalculateParametrisation(), << 63 fParametrisation( Par02DetectorParametrisation::eCMS ) {} 66 fParametrisation(Par02DetectorParametrisat << 67 {} << 68 64 69 //....oooOO0OOooo........oooOO0OOooo........oo 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 70 66 71 Par02FastSimModelHCal::~Par02FastSimModelHCal( << 67 Par02FastSimModelHCal::~Par02FastSimModelHCal() {} 72 68 73 //....oooOO0OOooo........oooOO0OOooo........oo 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 74 70 75 G4bool Par02FastSimModelHCal::IsApplicable(con << 71 G4bool Par02FastSimModelHCal::IsApplicable( const G4ParticleDefinition& aParticleType ) { 76 { << 77 G4bool isOk = false; 72 G4bool isOk = false; 78 // Applicable to all hadrons, i.e. any parti 73 // Applicable to all hadrons, i.e. any particle made of quarks 79 if (aParticleType.GetQuarkContent(1) + aPart << 74 if ( aParticleType.GetQuarkContent(1) + 80 + aParticleType.GetQuarkContent(3) + a << 75 aParticleType.GetQuarkContent(2) + 81 + aParticleType.GetQuarkContent(5) + a << 76 aParticleType.GetQuarkContent(3) + 82 != 0) << 77 aParticleType.GetQuarkContent(4) + 83 { << 78 aParticleType.GetQuarkContent(5) + >> 79 aParticleType.GetQuarkContent(6) != 0 ) { 84 isOk = true; 80 isOk = true; 85 } 81 } 86 return isOk; << 82 return isOk; 87 } 83 } 88 84 89 //....oooOO0OOooo........oooOO0OOooo........oo 85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 90 86 91 G4bool Par02FastSimModelHCal::ModelTrigger(con << 87 G4bool Par02FastSimModelHCal::ModelTrigger( const G4FastTrack& /*aFastTrack*/ ) { 92 { << 93 return true; // No kinematical restrictions 88 return true; // No kinematical restrictions to apply the parametrisation 94 } 89 } 95 90 96 //....oooOO0OOooo........oooOO0OOooo........oo 91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 97 92 98 void Par02FastSimModelHCal::DoIt(const G4FastT << 93 void Par02FastSimModelHCal::DoIt( const G4FastTrack& aFastTrack, 99 { << 94 G4FastStep& aFastStep ) { 100 // G4cout << " ________HCal model triggered << 95 //G4cout << " ________HCal model triggered _________" << G4endl; 101 96 102 // Kill the parameterised particle at the en 97 // Kill the parameterised particle at the entrance of the hadronic calorimeter 103 aFastStep.KillPrimaryTrack(); 98 aFastStep.KillPrimaryTrack(); 104 aFastStep.ProposePrimaryTrackPathLength(0.0) << 99 aFastStep.ProposePrimaryTrackPathLength( 0.0 ); 105 G4double Edep = aFastTrack.GetPrimaryTrack() 100 G4double Edep = aFastTrack.GetPrimaryTrack()->GetKineticEnergy(); 106 101 107 // Consider only primary tracks (do nothing 102 // Consider only primary tracks (do nothing else for secondary hadrons) 108 G4ThreeVector Pos = aFastTrack.GetPrimaryTra 103 G4ThreeVector Pos = aFastTrack.GetPrimaryTrack()->GetPosition(); 109 if (!aFastTrack.GetPrimaryTrack()->GetParent << 104 if ( ! aFastTrack.GetPrimaryTrack()->GetParentID() ) { 110 auto info = (Par02EventInformation*)G4Even << 105 Par02EventInformation* info = (Par02EventInformation*) 111 if (info->GetDoSmearing()) { << 106 G4EventManager::GetEventManager()->GetUserInformation(); >> 107 if ( info->GetDoSmearing() ) { 112 // Smearing according to the hadronic ca 108 // Smearing according to the hadronic calorimeter resolution 113 G4ThreeVector Porg = aFastTrack.GetPrima 109 G4ThreeVector Porg = aFastTrack.GetPrimaryTrack()->GetMomentum(); 114 G4double res = fCalculateParametrisation << 110 G4double res = fCalculateParametrisation-> 115 << 111 GetResolution( Par02DetectorParametrisation::eHCAL, 116 G4double eff = fCalculateParametrisation << 112 fParametrisation, Porg.mag() ); 117 << 113 G4double eff = fCalculateParametrisation-> >> 114 GetEfficiency( Par02DetectorParametrisation::eHCAL, >> 115 fParametrisation, Porg.mag() ); 118 G4double Esm; 116 G4double Esm; 119 Esm = std::abs(Par02Smearer::Instance()- << 117 Esm = std::abs( Par02Smearer::Instance()-> 120 Par02Output::Instance()->FillHistogram(2 << 118 SmearEnergy( aFastTrack.GetPrimaryTrack(), res ) ); >> 119 Par02Output::Instance()->FillHistogram( 2, (Esm/MeV) / (Edep/MeV) ); 121 // Setting the values of Pos, Esm, res a 120 // Setting the values of Pos, Esm, res and eff 122 auto primaryInfo = static_cast<Par02Prim << 121 Par02PrimaryParticleInformation* primaryInfo= 123 (aFastTrack.GetPrimaryTrack()->GetDyna << 122 static_cast<Par02PrimaryParticleInformation*>( 124 ->GetUserInformation()); << 123 ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )-> 125 primaryInfo->SetHCalPosition(Pos); << 124 GetUserInformation() ) ; 126 primaryInfo->SetHCalEnergy(Esm); << 125 primaryInfo->SetHCalPosition( Pos ); 127 primaryInfo->SetHCalResolution(res); << 126 primaryInfo->SetHCalEnergy( Esm ); 128 primaryInfo->SetHCalEfficiency(eff); << 127 primaryInfo->SetHCalResolution( res ); >> 128 primaryInfo->SetHCalEfficiency( eff ); 129 // The (smeared) energy of the particle 129 // The (smeared) energy of the particle is deposited in the step 130 // (which corresponds to the entrance of 130 // (which corresponds to the entrance of the hadronic calorimeter) 131 aFastStep.ProposeTotalEnergyDeposited(Es << 131 aFastStep.ProposeTotalEnergyDeposited( Esm ); 132 } << 132 } else { 133 else { << 134 // No smearing: simply setting the value 133 // No smearing: simply setting the value of Edep 135 ((Par02PrimaryParticleInformation*)(cons << 134 ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* > 136 aF << 135 ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )-> 137 << 136 GetUserInformation() ) )->SetHCalEnergy( Edep ); 138 << 139 -> << 140 ->SetHCalEnergy(Edep); << 141 // The (initial) energy of the particle 137 // The (initial) energy of the particle is deposited in the step 142 // (which corresponds to the entrance of 138 // (which corresponds to the entrance of the hadronic calorimeter) 143 aFastStep.ProposeTotalEnergyDeposited(Ed << 139 aFastStep.ProposeTotalEnergyDeposited( Edep ); 144 } 140 } 145 } 141 } 146 } 142 } 147 143 148 //....oooOO0OOooo........oooOO0OOooo........oo 144 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 145 149 146