]>
Commit | Line | Data |
---|---|---|
fa3d1cc7 | 1 | #include <Riostream.h> |
8fd1d27e | 2 | #include "TVirtualMCApplication.h" |
df14e7fb | 3 | #include "TGeoMaterial.h" |
4 | #include "TGeoManager.h" | |
18e0cabb | 5 | #include <TParticle.h> |
df14e7fb | 6 | #include "TFlukaCerenkov.h" |
a7bb59a2 | 7 | |
fa3d1cc7 | 8 | #include "TFluka.h" |
d566901f | 9 | #include "TFlukaCodes.h" |
a7bb59a2 | 10 | |
7dac99f1 | 11 | #include "Fdimpar.h" //(DIMPAR) fluka include |
12 | #include "Ftrackr.h" //(TRACKR) fluka common | |
d566901f | 13 | #include "Fltclcm.h" //(LTCLCM) fluka common |
14 | #include "Fpaprop.h" //(PAPROP) fluka common | |
18e0cabb | 15 | |
fa3d1cc7 | 16 | #ifndef WIN32 |
17 | # define endraw endraw_ | |
18 | #else | |
19 | # define endraw ENDRAW | |
20 | #endif | |
21 | extern "C" { | |
22 | void endraw(Int_t& icode, Int_t& mreg, Double_t& rull, Double_t& xsco, Double_t& ysco, Double_t& zsco) | |
23 | { | |
4b81a588 | 24 | TFluka* fluka = (TFluka*) gMC; |
d566901f | 25 | // nothing to do if particle in dummy region |
26 | if (mreg == fluka->GetDummyRegion()) return; | |
2bc4c610 | 27 | Int_t verbosityLevel = fluka->GetVerbosityLevel(); |
2bc78182 | 28 | Bool_t debug = (verbosityLevel >= 3)? kTRUE : kFALSE; |
4aba9d66 | 29 | Int_t mlttc = (icode==kKASKADinelarecoil) ? TRACKR.lt2trk : TRACKR.lt1trk; //LTCLCM.mlatm1; |
d566901f | 30 | fluka->SetCaller(kENDRAW); |
4b81a588 | 31 | fluka->SetRull(rull); |
32 | fluka->SetXsco(xsco); | |
33 | fluka->SetYsco(ysco); | |
34 | fluka->SetZsco(zsco); | |
d566901f | 35 | fluka->SetMreg(mreg, mlttc); |
4aba9d66 | 36 | |
022d4275 | 37 | Float_t edep = rull; |
38 | ||
df14e7fb | 39 | if (TRACKR.jtrack == -1) { |
18e0cabb | 40 | // Handle quantum efficiency the G3 way |
2bc4c610 | 41 | if (debug) printf("endraw: Cerenkov photon depositing energy: %d %e\n", mreg, rull); |
df14e7fb | 42 | TGeoMaterial* material = (gGeoManager->GetCurrentVolume())->GetMaterial(); |
df14e7fb | 43 | TFlukaCerenkov* cerenkov = dynamic_cast<TFlukaCerenkov*> (material->GetCerenkovProperties()); |
44 | if (cerenkov) { | |
18e0cabb | 45 | Double_t eff = (cerenkov->GetQuantumEfficiency(rull)); |
46 | if (gRandom->Rndm() > eff) { | |
47 | edep = 0.; | |
48 | } | |
df14e7fb | 49 | } |
50 | } | |
2bc78182 | 51 | |
18e0cabb | 52 | TVirtualMCStack* cppstack = fluka->GetStack(); |
18e0cabb | 53 | |
54 | if (debug) { | |
4aba9d66 | 55 | cout << "ENDRAW For icode=" << icode |
18e0cabb | 56 | << " track=" << TRACKR.ispusr[mkbmx2-1] << " pdg=" << fluka->PDGFromId(TRACKR.jtrack) |
4aba9d66 | 57 | << " edep="<< edep << " mreg=" << mreg << endl; |
18e0cabb | 58 | } |
59 | ||
4aba9d66 | 60 | // check region lattice consistency (debug Ernesto) |
61 | // ***************************************************** | |
62 | Int_t nodeId; | |
63 | Int_t volId = fluka->CurrentVolID(nodeId); | |
64 | Int_t crtlttc = gGeoManager->GetCurrentNodeId()+1; | |
4b1c3976 | 65 | if(debug && mreg != volId && !gGeoManager->IsOutside() ) { |
4aba9d66 | 66 | cout << " endraw: track=" << TRACKR.ispusr[mkbmx2-1] << " pdg=" << fluka->PDGFromId(TRACKR.jtrack) |
67 | << " icode=" << icode << " gNstep=" << fluka->GetNstep() << endl | |
68 | << " fluka mreg=" << mreg << " mlttc=" << mlttc << endl | |
69 | << " TGeo volId=" << volId << " crtlttc=" << crtlttc << endl | |
70 | << " common TRACKR lt1trk=" << TRACKR.lt1trk << " lt2trk=" << TRACKR.lt2trk << endl | |
71 | << " common LTCLCM newlat=" << LTCLCM.newlat << " mlatld=" << LTCLCM.mlatld << endl | |
72 | << " mlatm1=" << LTCLCM.mlatm1 << " mltsen=" << LTCLCM.mltsen << endl | |
73 | << " mltsm1=" << LTCLCM.mltsm1 << " mlattc=" << LTCLCM.mlattc << endl; | |
74 | if( mlttc == crtlttc ) cout << " *************************************************************" << endl; | |
75 | } | |
76 | // ***************************************************** | |
77 | ||
78 | ||
d566901f | 79 | if (icode != kEMFSCOstopping1 && icode != kEMFSCOstopping2) { |
80 | fluka->SetIcode((FlukaProcessCode_t)icode); | |
022d4275 | 81 | fluka->SetRull(edep); |
18e0cabb | 82 | if (icode == kKASKADelarecoil && TRACKR.ispusr[mkbmx2-5]) { |
4aba9d66 | 83 | // Elastic recoil and in stuprf npprmr > 0, |
84 | // the secondary being loaded is actually still the interacting particle | |
85 | cppstack->SetCurrentTrack( TRACKR.ispusr[mkbmx2-4] ); | |
86 | // cout << "endraw elastic recoil track=" << TRACKR.ispusr[mkbmx2-1] << " parent=" << TRACKR.ispusr[mkbmx2-4] | |
87 | // << endl; | |
18e0cabb | 88 | } |
89 | else | |
4aba9d66 | 90 | cppstack->SetCurrentTrack(TRACKR.ispusr[mkbmx2-1] ); |
022d4275 | 91 | (TVirtualMCApplication::Instance())->Stepping(); |
92 | } else { | |
2bc78182 | 93 | // |
d566901f | 94 | // For icode 21,22 the particle has fallen below thresshold. |
2bc78182 | 95 | // This has to be signalled to the StepManager() |
96 | // | |
18e0cabb | 97 | cppstack->SetCurrentTrack( TRACKR.ispusr[mkbmx2-1] ); |
a5b73d69 | 98 | fluka->SetRull(edep); |
d566901f | 99 | fluka->SetIcode((FlukaProcessCode_t) icode); |
022d4275 | 100 | (TVirtualMCApplication::Instance())->Stepping(); |
101 | fluka->SetTrackIsNew(kFALSE); | |
d566901f | 102 | fluka->SetIcode((FlukaProcessCode_t)icode); |
a5b73d69 | 103 | fluka->SetRull(0.); |
022d4275 | 104 | (TVirtualMCApplication::Instance())->Stepping(); |
18e0cabb | 105 | // cppstack->SetCurrentTrack( saveTrackId ); |
106 | ||
022d4275 | 107 | } |
fa3d1cc7 | 108 | } // end of endraw |
109 | } // end of extern "C" | |
110 |