missing file...

ECALELFS · Jun 9, 2015 · 17f757d · 17f757d
1 parent 1928742
commit 17f757d
Show file tree

Hide file tree

Showing 2 changed files with 14 additions and 6 deletions.
diff --git a/EcalTiming/plugins/EcalTimingCalibProducer.cc b/EcalTiming/plugins/EcalTimingCalibProducer.cc
@@ -418,15 +418,18 @@ bool EcalTimingCalibProducer::addRecHit(const EcalRecHit& recHit)
 
 void EcalTimingCalibProducer::plotRecHit(const EcalTimingEvent& recHit)
 {
+        std::cout<<"I'm in PlotrecHit"<<std::endl;
 	if(recHit.detid().subdetId() == EcalBarrel) {
 		EBDetId id(recHit.detid());
+                std::cout<<"filling histos for rechit: "<<id.ieta()<<std::endl;
 		// Fill Rechit Energy
 		Event_EneMapEB_->Fill(id.ieta(), id.iphi(), recHit.energy()); // 2D energy map
 		Event_TimeMapEB_->Fill(id.ieta(), id.iphi(), recHit.time()); // 2D time map
 		RechitEnergyTimeEB->Fill(recHit.energy(), recHit.time());
 	} else {
 		// create EEDetId
 		EEDetId id(recHit.detid());
+                std::cout<<"filling histos for rechit: "<<id.ix()<<std::endl;
 		if(id.zside() < 0) {
 			Event_EneMapEEM_->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<EcalRecHitCollection>::const_iterator
+        //std::cout<<"EB recHit handle: "<<ebRecHitHandle->size()<<std::endl;
+        //std::cout<<"EE recHit handle: "<<eeRecHitHandle->size()<<std::endl;
 	for(auto  recHit_itr = ebRecHitHandle->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
@@ -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,