41 #include "CLHEP/Random/RandomEngine.h"
47 #include "G4HCofThisEvent.hh"
48 #include "G4ProcessManager.hh"
49 #include "G4SDManager.hh"
52 #include "G4TransportationManager.hh"
53 #include "G4UserEventAction.hh"
54 #include "G4VProcess.hh"
69 produces<DigiCollectionFP420>().setBranchAlias(alias);
95 std::cout <<
"Creating a DigitizerFP420" << std::endl;
96 std::cout <<
"DigitizerFP420: dn0=" << dn0 <<
" sn0=" << sn0 <<
" pn0=" << pn0 <<
" rn0=" << rn0 << std::endl;
104 std::cout <<
"Destroying a DigitizerFP420" << std::endl;
117 std::cout <<
" ============== DigitizerFP420: start produce= " << std::endl;
128 std::vector<const CrossingFrame<PSimHit> *> cf_simhitvec;
131 cf_simhitvec.push_back(cf_simhit.product());
164 std::cout <<
" ============== DigitizerFP420: MixCollection treatment= " << std::endl;
169 for (isim = allTrackerHits->begin(); isim != allTrackerHits->end(); isim++) {
170 unsigned int unitID = (*isim).detUnitId();
171 int det,
zside, sector, zmodule;
183 double losenergy = (*isim).energyLoss();
185 std::cout <<
" ============== DigitizerFP420: losenergy= " << losenergy << std::endl;
186 std::cout <<
" === for intindex = " << intindex << std::endl;
198 std::cout <<
" ============== DigitizerFP420: put zero to container " << std::endl;
204 for (
int det = 1; det <
dn0; det++) {
205 for (
int sector = 1; sector <
sn0; sector++) {
206 for (
int zmodule = 1; zmodule <
pn0; zmodule++) {
216 inputRange.first = collector.begin();
217 inputRange.second = collector.end();
218 output->putclear(inputRange, detID);
233 std::cout <<
" ============== DigitizerFP420: start loop over det iu " << std::endl;
234 std::cout <<
" ============== DigitizerFP420: SimHitMap.size()= " <<
SimHitMap.size() << std::endl;
263 vector<PSimHit>::const_iterator simHitIter =
SimHitMap[
i].begin();
264 vector<PSimHit>::const_iterator simHitIterEnd =
SimHitMap[
i].end();
265 for (; simHitIter != simHitIterEnd; ++simHitIter) {
266 const PSimHit ihit = *simHitIter;
269 std::cout <<
" ====== DigitizerFP420: unitID= " << unitID <<
"Hit number i= " <<
i << std::endl;
270 int det,
zside, sector, zmodule;
275 std::cout <<
"for Hits iu = " << iu <<
" sector = " << sector <<
" zmodule = " << zmodule
276 <<
" zside = " << zside <<
" det=" << det << std::endl;
284 std::cout <<
" DigitizerFP420: index = " << index <<
" iu = " << iu << std::endl;
289 G4ThreeVector bfield(0., 0., 0.0);
294 std::cout <<
" ============== DigitizerFP420: call run for iu= " << iu << std::endl;
306 std::cout <<
"======= DigitizerFP420: collector size = " <<
collector.size() << std::endl;
321 std::cout <<
" ============= DigitizerFP420:collector "
322 "start!!!!!!!!!!!!!!"
333 unsigned int detID0 = 0;
334 output->put(outputRange, detID0);
339 output->put(outputRange, iu);
377 for (
int det = 1; det <
dn0; det++) {
378 for (
int sector = 1; sector <
sn0; sector++) {
379 for (
int zmodule = 1; zmodule <
pn0; zmodule++) {
384 std::cout <<
"****DigitizerFP420:check2" << std::endl;
389 int newdet, newzside, newsector, newzmodule;
391 std::cout <<
" newdet = " << newdet <<
" newsector = " << newsector <<
" newzmodule = " << newzmodule
392 <<
" newzside = " << newzside << std::endl;
397 outputRange = output->get(iu);
405 for (; sort_begin != sort_end; ++sort_begin) {
406 collector.push_back(*sort_begin);
410 std::cout <<
"====== collector size = " << collector.size() << std::endl;
411 if (!collector.empty()) {
412 std::cout <<
" iu = " << iu <<
" sector = " << sector <<
" zmodule = " << zmodule
413 <<
" zside = " <<
zside <<
" det=" << det <<
" layer = " << layer <<
" orient = " << orient
417 vector<HDigiFP420>::const_iterator simHitIter = collector.begin();
418 vector<HDigiFP420>::const_iterator simHitIterEnd = collector.end();
419 for (; simHitIter != simHitIterEnd; ++simHitIter) {
421 std::cout <<
" strip number=" << istrip.
strip() <<
" adc=" << istrip.
adc() << std::endl;
426 std::cout <<
" ===================================================" << std::endl;
442 std::cout <<
"DigitizerFP420 recoutput" << std::endl;
T getUntrackedParameter(std::string const &, T const &) const
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
FP420DigiMain * stripDigitizer_
std::vector< HDigiFP420 >::const_iterator ContainerIterator
static void unpackFP420Index(const unsigned int &idx, int &det, int &zside, int &station, int &superplane)
static int unpackLayerIndex(int rn0, int zside)
void produce(edm::Event &e, const edm::EventSetup &c) override
std::vector< HDigiFP420 > collector
static void unpackMYIndex(const int &idx, int rn0, int pn0, int sn0, int &det, int &zside, int §or, int &zmodule)
constexpr std::array< uint8_t, layerIndexSize > layer
std::vector< HDigiFP420 > run(const std::vector< PSimHit > &input, const G4ThreeVector &, unsigned int)
static unsigned packMYIndex(int rn0, int pn0, int sn0, int det, int zside, int sector, int zmodule)
~DigitizerFP420() override
static unsigned int packFP420Index(int det, int zside, int station, int superplane)
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
std::pair< ContainerIterator, ContainerIterator > Range
T getParameter(std::string const &) const
static int unpackOrientation(int rn0, int zside)
vstring trackerContainers
unsigned int detUnitId() const