diff --git a/EcalTiming/plugins/EcalTimingCalibProducer.cc b/EcalTiming/plugins/EcalTimingCalibProducer.cc index 0a7c1c6..9a77c3d 100644 --- a/EcalTiming/plugins/EcalTimingCalibProducer.cc +++ b/EcalTiming/plugins/EcalTimingCalibProducer.cc @@ -418,8 +418,10 @@ bool EcalTimingCalibProducer::addRecHit(const EcalRecHit& recHit) void EcalTimingCalibProducer::plotRecHit(const EcalTimingEvent& recHit) { + std::cout<<"I'm in PlotrecHit"<Fill(id.ieta(), id.iphi(), recHit.energy()); // 2D energy map Event_TimeMapEB_->Fill(id.ieta(), id.iphi(), recHit.time()); // 2D time map @@ -427,6 +429,7 @@ void EcalTimingCalibProducer::plotRecHit(const EcalTimingEvent& recHit) } else { // create EEDetId EEDetId id(recHit.detid()); + std::cout<<"filling histos for rechit: "<Fill(id.ix(), id.iy(), recHit.energy()); Event_TimeMapEEM_->Fill(id.ix(), id.iy(), recHit.time()); @@ -488,6 +491,8 @@ EcalTimingCalibProducer::Status EcalTimingCalibProducer::duringLoop(const edm::E // loop over the recHits // recHit_itr is of type: edm::Handle::const_iterator + //std::cout<<"EB recHit handle: "<size()<size()<begin(); recHit_itr != ebRecHitHandle->end(); ++recHit_itr) { addRecHit(*recHit_itr); // add the recHit to the list of recHits used for calibration (with the relative information) } @@ -534,7 +539,8 @@ EcalTimingCalibProducer::Status EcalTimingCalibProducer::duringLoop(const edm::E continue; } - if(_makeEventPlots) plotRecHit(event); + plotRecHit(event); + //if(_makeEventPlots) plotRecHit(event); _timeCalibMap[it.first].add(event); } diff --git a/EcalTiming/test/ecalTime_fromAlcaStream_cfg.py b/EcalTiming/test/ecalTime_fromAlcaStream_cfg.py index 79d408d..6c9917b 100644 --- a/EcalTiming/test/ecalTime_fromAlcaStream_cfg.py +++ b/EcalTiming/test/ecalTime_fromAlcaStream_cfg.py @@ -25,7 +25,7 @@ ### setup any defaults you want options.output="ecalTiming.root" options.secondaryOutput="ntuple.root" -options.files= "" +options.files= "/store/data/Commissioning2015/AlCaPhiSym/RAW/v1/000/244/768/00000/A8219906-44FD-E411-8DA9-02163E0121C5.root" options.maxEvents = -1 # -1 means all events ### get and parse the command line arguments options.parseArguments() @@ -57,10 +57,11 @@ process.load('Configuration/StandardSequences/Reconstruction_cff') process.recoSequence = cms.Sequence(process.calolocalreco )#+ process.egammaCosmics) -process.load('PhiSym.EcalCalibAlgos.ecalPhiSymLocarecoWeights_cff') +#process.load('PhiSym.EcalCalibAlgos.ecalPhiSymLocarecoWeights_cff') +#process.load('RecoLocalCalo.Configuration.ecalLocalRecoSequence_cff') process.load('Configuration.StandardSequences.EndOfProcess_cff') process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff') - +process.load('EcalTiming.EcalTiming.ecalLocalRecoSequenceAlCaStream_cff') ## Raw to Digi process.load('Configuration/StandardSequences/RawToDigi_Data_cff') @@ -169,7 +170,8 @@ process.filter+=process.spashesHltFilter process.reco_step = cms.Sequence(process.caloCosmicOrSplashRECOSequence) else: - process.reco_step = cms.Sequence(process.reconstruction_step_multiFit) + #process.reco_step = cms.Sequence(process.reconstruction_step_multiFit) + process.reco_step = cms.Sequence(process.ecalLocalRecoSequenceAlCaStream) ### Process Full Path if(options.isSplash==0): @@ -194,7 +196,7 @@ recHitFlags = cms.vint32([0]), # only recHits with these flags are accepted for calibration recHitMinimumN = cms.uint32(10), #recHitMinimumN = cms.uint32(2), - minRecHitEnergy = cms.double(1), + minRecHitEnergy = cms.double(0), globalOffset = cms.double(options.offset), produceNewCalib = cms.bool(True), outputDumpFile = process.TFileService.fileName,