| 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